X=@(x,t)x;
Y=@(x,t)(x-x.*sin(3*x)).*cos(t);
Z=@(x,t)(x-x.*sin(3*x)).*sin(t);
fmesh(X,Y,Z,[0,1,0,2*pi])
xlabel('x');ylabel('y');zlabel('z');
title({'$x=x$','$y=(x-x\sin{3x})\cos t$','$z=(x-x\sin{3x})\sin t$'},'interp','la')