import sympy as sp
import sympy.plotting as spl

x, y, z, u, v = sp.symbols("x y z u v")


# z = f(x,y)
spl.plot3d(x**2 + y**2)
spl.plot3d(x**2 + y**2, (x,-5,5), (y,0,5))

# x(u,v), y(u,v), z(u,v)
# u = theta, v = z, r = 1
spl.plot3d_parametric_surface(sp.cos(u), sp.sin(u), v)

# u = theta, v = phi, rho = 1
spl.plot3d_parametric_surface(sp.cos(u)*sp.sin(v), sp.sin(u)*sp.sin(v), sp.cos(v), (v, 0, sp.pi/2), (u,0,sp.pi/2))