> r:=rho*sin(phi);

> x:=r*cos(theta);

> y:=r*sin(theta);

> z:=rho*cos(phi);

> rho:=2*sin(phi)^2*cos(phi);

> a:=plot3d([x,y,z],phi=0..Pi/2,theta=0..2*Pi,scaling=constrained):

> a;

> with(plottools):

> vrml(a,"c:/042/12/12-9.wrl",background_color=white);

> b:=plot3d([x,y,z],phi=0..Pi/2,theta=0..Pi,scaling=constrained):

> b;

> vrml(a,"c:/042/12/12-9'.wrl",background_color=white);

>