t=linspace(0,2*pi,37);
r=linspace(0,1,10);
[t,r]=meshgrid(t,r);
[x,y]=pol2cart(t,r);
mesh(x,y,r);
hold on
r=r*pi/2;
[x,y,z]=sph2cart(t,r,1+0*t);
mesh(x,y,z+1)
axis equal
t=linspace(0,2*pi,37);
r=linspace(0,1,10);
[t,r]=meshgrid(t,r);
[x,y]=pol2cart(t,r);
mesh(x,y,r);
hold on
r=r*pi/2;
[x,y,z]=sph2cart(t,r,1+0*t);
mesh(x,y,z+1)
axis equal