All published worksheets from http://sagenb.org
S=Cylindrical('radius', ['azimuth', 'height'])
theta,z=var('theta, z')
plot3d(3, (theta,pi,2*pi), (z, -4, 1), transformation=S)