> restart:with(plots):with(plottools):
> f:=2*arctan((sqrt(5)-1)/2): evalf(180*f/Pi,3):
> s:= array[1..6]:szin:=[yellow,blue,cyan,magenta,green,red]:
> sp:= (r,d,ro,a,b) -> tubeplot(
> [(r+d*sin(5*t))*cos(t),(r+d*sin(5*t))*sin(t),d*cos(5*t)],
> t=0..2*Pi,radius=ro, numpoints=a, tubepoints=b,scaling=constrained):
>   sp1:=sp(10,2,  1   , 80,11): # fájlnév : 6sp.wrl
# sp1:=sp(10,0.5,0.4 , 80, 7): # fájlnév : 6sp2.wrl
# sp1:=sp(10,3  ,0.75,120, 7): # fájlnév : 2sp3.wrl  és  6sp3.wrl
# sp1:=sp(10,7  ,0.3 ,150, 7): # fájlnév : 2sp4.wrl  és  6sp4.wrl
> s[1]:=rotate(sp1,0,0,Pi):
> s[2]:=rotate(sp1,0,f,0):
> s[3]:=rotate(s[2],0,0,2*Pi/5):
> s[4]:=rotate(s[2],0,0,4*Pi/5):
> s[5]:=rotate(s[2],0,0,6*Pi/5):
> s[6]:=rotate(s[2],0,0,8*Pi/5):
> spiral:=seq(plots[display](s[i],color=szin[i]),i=1..6):
> display(spiral,scaling=CONSTRAINED);
> vrml(display(spiral,scaling=CONSTRAINED),
>     `6sp.wrl`, background_color=COLOR(RGB,0,1,1), shininess=0.4, specular_color=COLOR(RGB,0.5,0.4,0.5),
> directional_light=[1.0, COLOR(RGB,1,0,1), [0,0, 15]],
> directional_light=[1.0, COLOR(RGB,1,0,1), [0,0,-15]]);
>