Python 3 code
####################################
# Whole sphere generation
# by paramteric equations
# (matplotlib module)
####################################

from mpl_toolkits.mplot3d import Axes3D   # import matplotlib module
import matplotlib.pyplot as plt   # import matplotlib module
import numpy as np    # import numpy module

angle = np.linspace(0, 2 * np.pi, 32)   #  define cylindrical co-ordinates array
theta, phi = np.meshgrid(angle, angle)   # set theta, phi grid

R = 1.0     # set radius value
X = R * np.cos(phi) * np.cos(theta)     # evaluate x co-ordinate array values
Y = R * np.cos(phi) * np.sin(theta)      # evaluate y co-ordinate array values
Z = R * np.sin(phi)     # evaluate z co-ordinate array values

plt.style.use('dark_background')     # set black background color

fig = plt.figure()     # define 3D plot and axes
ax = fig.gca(projection = '3d')

ax.set_xlim3d(-1, 1)        # set x axis limits
ax.set_ylim3d(-1, 1)        # set y axis limits
ax.set_zlim3d(-1, 1)        # set z axis limits

ax.plot_surface(X, Y, Z, cmap = "copper", rstride = 1, cstride = 1)   # plot surface as a whole
ax.axis('off')    # do not plot co-ordinate axes
plt.show()     # show surface a a whole