Maple V code: with(plots): with(plottools): f := t > t^3/22*t^2+2*t+1: start := spacecurve([0,t,f(t)],t=0..3,thickness=3): pic := n> cylinderplot(f(z),theta=0..n*2*Pi/30,z=0..3): display(start,seq(rotate(pic(n),Pi/2,Pi/2,0),n=1..30), insequence=true,axes=normal,tickmarks=[0,0,0]); 

Comments:
