> restart:

> evalc(Re((x+y*I)^3)):

> x:=r*cos(t):

> y:=r*sin(t):

> z:=x^3-3*x*y^2:

> a:=plot3d([x,y,z],r=0..1,t=0..2*Pi,scaling=constrained,grid=[5,50],color=red):

> with(plots):

> with(plottools):

>

> x:=cos(t):

> y:=sin(t):

> z1:=-1:

> z2:=cos(3*t):

> z:=(1-s)*z1+s*z2:

> b:=plot3d([x,y,z],s=0..1,t=0..2*Pi,scaling=constrained,grid=[5,50],color=blue):

> c:=display(a,b):

> vrml(c,"c:/My documents/87124/11/11-10.wrl",background_color=white):

> display(c);

>