> with(plots):
f := (5*cos(y)^2 -1)/3:
cylinderplot([z*theta,theta,cos(z^2)],theta=0..Pi,z=-2..2, color = theta);

>