t=0:pi/20:2*pi;
[x,y,z]=cylinder(2+sin(t),30);
subplot(1,3,1);surf(x,y,z);
[x,y,z]=sphere;
subplot(1,3,1);surf(x,y,z);