syms u v;
d=[-5:0.5:5];
[u v]=meshgrid(d);
x=u.*sin(v),y=u.*cos(v),z=u;
surf(x,y,z)