input mp-solid

nb:=48;subh:=12;

figureespace(-10u,-10u,10u,10u);
Initialisation(500,30,20,30);
arcenciel:=true;
incolor:=orange;
Ferme1:=false;
ObjetCylindre1("(3*(cos(u)**3),3*(sin(u)**3),-2+v)",pi,-pi,0,4);
AffichageObjet1;
finespace;
unit:=0.6;
figureespace(-10u,-10u,10u,10u);
Initialisation(500,30,20,50);
drawoptions(withcolor gris);
for k=-4 upto 4:
  draw Projette((-4,k,0))--Projette((4,k,0));
  draw Projette((k,-4,0))--Projette((k,4,0));
endfor;
drawoptions();
Lumiere:=Oeil;
arcenciel:=true;
incolor:=orange;
Ferme1:=false;
ObjetCylindre1("(3*(cos(u)**3),3*(sin(u)**3)+v,2*v)",pi,-pi,0,2);
AffichageObjet1;
drawarrow Projette((0,3,0))--Projette((0,5,4)) withpen pencircle scaled2bp withcolor violet;
label.rt(btex $\overrightarrow v=(0;2;4)$ etex rotated angle(Projette((0,5,4)) -Projette((0,3,0))),1/2[Projette((0,3,0)),Projette((0,5,4))]);
TraceAxes;
finespace;
end