figure; [X,Y,Z]=cylinder(0:25,20); surf(X,Y,Z);hold on surf(2*X,Y,2*Z);

figure; z=cplxgrid(20) x=real(z); y=imag(z); fz=sqrt((x.^2+y.^2)) cplxmap(z,fz)