$ cat forumdxdy.py
#!/usr/bin/python
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
import numpy as np
fig = plt.figure()
ax = Axes3D(fig)
u = np.linspace(0, np.pi, 6)
v = np.linspace(0, 0.2*np.pi, 6)
x = np.outer(np.cos(u), np.sin(v))
y = np.outer(np.sin(u), np.sin(v))
z = np.outer(np.ones(np.size(u)), np.cos(v))
print "x = ", x
print "y = ", y
print "z = ", z
ax.plot_surface(x, y, z, rstride=1, cstride=1, color='b')
plt.show()