FPS = 30;
t=[0:0.1:16*pi];
for i = 0:FPS*20
    x = cos(t + pi*i/FPS).*t/(8*pi);
    y = sin(t + pi*i/FPS).*t/(8*pi);
    z = 20+sin(i/(4*pi)).*t/4;
    plot3(x,y,z); 
    axis([-4 4 -4 4 10 30]);
    pause(1/FPS);
end