Möbius szalag piros-kék színnel
Adatok
> restart: with(plots):with(plottools):
>  a:=-1: b:=1: # szélesség
>  c:=0: d:=Pi: # hosszúság
>  r:=5:        # kerület
>  emel:=0:     # kiemelkedés az xy síkból
>  s1:= MAGENTA: s2:= MAGENTA: # perem színei 
>  k:=1:  # Csavarás mértéke (Möbius szalag)
>  p:=1:  # Hullámzás mértéke (Henger palást)
>  nev:=`dupla.wrl`:
Számolás
> racs:=round(36*(d-c)/Pi)+1: vonal:=(b-a)/15:
> u:=r*cos(2*t):
> v:=r*sin(2*p*t)/p:
> w:=emel:
> df:=.03:
> x:= u + f*sin(k*(t-df))*cos(2*t): 
> y:= v + f*sin(k*(t-df))*sin(2*t): 
> z:= w + f*cos(k*(t-df)): 
> 
> xb:= u + a*sin(k*t)*cos(2*t): 
> yb:= v + a*sin(k*t)*sin(2*t):
> zb:= w + a*cos(k*t): 
> xj:= u + b*sin(k*t)*cos(2*t): 
> yj:= v + b*sin(k*t)*sin(2*t):
> zj:= w + b*cos(k*t): 
> szalag1:=plot3d([x,y,z],f=a..b,t=c..d,
> grid=[5,racs],scaling=CONSTRAINED,color=RED): 
> 
> x:= u + f*sin(k*(t+df))*cos(2*t): 
> y:= v + f*sin(k*(t+df))*sin(2*t): 
> z:= w + f*cos(k*(t+df)): 
> 
> szalag2:=plot3d([x,y,z],f=a..b,t=c+1.5..d+1.5,
> grid=[5,racs],scaling=CONSTRAINED,color=BLUE): 
> 
> vonal1:=tubeplot({[xb,yb,zb]},t=c..d,radius=vonal,grid=[racs,5],scaling=CONSTRAINED, color=s1,  style=PATCHNOGRID):
> vonal2:=tubeplot({[xj,yj,zj]},t=c..d,radius=vonal,grid=[racs,5],scaling=CONSTRAINED, color=s2,  style=PATCHNOGRID):
Megjelenítés
> display3d(szalag1,szalag2,vonal1,vonal2,style=PATCH, orientation=[-120,60]);

VRML fájl készítés
> vrml(display3d(szalag1,szalag2,color=COLOR(RGB,0.6,0.5,0),style=PATCHNOGRID),nev, background_color=COLOR(RGB,0,1,1), shininess=0.3, specular_color=COLOR(RGB,1,1,1)):
>