>	with(plots):

>	surface:=implicitplot3d(y^2+z^2=(x^3-6*x^2+9*x+2)^2,x=0..4,y=-6..6,z=-6..6,axes=normal,numpoints=2500,labels=[x,y,z]):

>	curve1:=spacecurve([t,t^3-6*t^2+9*t+2,0,t=0..4],color=blue,thickness=4,axes=normal,labels=[x,y,z]):

>	curve2:=spacecurve([t,0,t^3-6*t^2+9*t+2,t=0..4],color=red,thickness=4,axes=normal,labels=[x,y,z]):

>	display(curve1,curve2);

>	display(surface,curve1,curve2);

>