lundi 8 août 2011

Maple vs PStricks

Il est intéressant de comparer les réalisations obtenues avec Maple et PStricks. Il est évident que les calculs sont bien plus rapides avec Maple et les calculs littéraux étant effectués par le logiciel, les erreurs sont minimisées. L'avantage pour PStricks est que le dessin est plus facilement personnalisable et que les possibilités d'enrichissement du schéma (couleurs, notations etc.) sont incomparables. À vous de décider...
Le code ci-dessous est adapté de celui(excellent) de la page de Philippe Dumas  : code Maple.
J'ai modifié le paramétrage du cône que j'ai trouvé inutilement compliqué.

>    cone_param:=[r*cos(psi),r*sin(psi),r/tan(theta)];


>    plane_param:=[v,h*cos(omega)-u*sin(omega),h*sin(omega)+u*cos(omega)];
>    plane_cart:=y*cos(omega)+z*sin(omega)-h; 


>    equ:=subs(x=cone_param[1],y=cone_param[2],z=cone_param[3],plane_cart);

>    rr:=solve(equ,r);
>    conic:=subs(r=rr,cone_param); 
Theta:=Pi/4:Omega:=Pi/3:H:=4:
m:=20:
picture[cone]:=plot3d(subs(theta=Theta,h=H,cone_param),r=-m..m,psi=-Pi..Pi,shading=xyz):
>    picture[plane]:=plot3d(subs(omega=Omega,h=H,plane_param),u=-m..1.5*m,v=-m..m,grid=[9,9],color=[0.9,0.9,1]):
>    picture[conic]:=plots[spacecurve](subs(theta=Theta,omega=Omega,h=H,conic),psi=-Pi..Pi,color=blue,thickness=3):
>    plots[display]({seq(picture[i],i={cone,plane,conic})},scaling=constrained,axes=boxed);


Aucun commentaire:

Enregistrer un commentaire