> with(plots):
f := (5*cos(y)^2 -1)/3:
cylinderplot(f, x=0..2*Pi,y=-Pi..Pi,style=PATCH, color = f);

>