t = 0:pi/10:2*pi;
[X,Y,Z] = cylinder(2+sin(t));
figure(1)
mesh(X,Y,Z)
axis tight square











