clear;close;
u=-4:0.1:4;
v=-4:0.1:4;
[U,V]=meshgrid(u,v);
X=U.*sin(V);
Y=U.*cos(V);
Z=U;
mesh(X,Y,Z);
shading interp;
colormap(jet);
axis equal
clear;close;
u=-4:0.1:4;
v=-4:0.1:4;
[U,V]=meshgrid(u,v);
X=(3+0.4*cos(U)).*cos(V);
Y=(3+0.4*cos(U)).*sin(V);
Z=0.4*sin(V);
mesh(X,Y,Z);
shading interp;
colormap(jet);
axis equal