Calcul du rotationnel d'un champ de vecteurs pour V200 et TI-89

Bonsoir,

Voici le code d'un petit prog qui calcule le rotationnel d'un champ de vecteurs de la forme E = (P(x,y,z), Q(x,y,z), R(x,y,z)). Attention aux caractéres spéciaux de TI qui ne sortent pas bien au copier/coller. Pour ceux qui le souhaitent, je peux leur communiquer le source.
Pour info, ce progr calcule le rotationnele n faisant le produit vectoriel de l'opérateur Nabla et du vecteur E.

()
Prgm
¨ Programme de calcul du rotationnel d'un champ de vecteurs
¨ Dominique Lefebvre Novembre 2008
¨ ce programme calcule le rotationnel d'un champ de vecteurs
¨ E(P(x,y,z), Q(x,y,z), R(x,y,z)) pour produit vectoriel du
¨ vecteur Nabla et de E.
ClrIO
DelVar P,Q,R,rotat
Lbl Debut
ToolBar
Title "Def E", Comp
Title "x,y,z", Point
Title "Rot", Rota
Title "Init", Init
Title "Quitter", Fin
EndTBar
¨ Saisie des composantes du champ E
Lbl Comp
ClrIO
Disp "Saisie des composantes de E"
Input "P(x,y,z) = ", P
Input "Q(x,y,z) = ", Q
Input "R(x,y,z) = ", R
Goto Debut
¨ Détermination du point de calcul
Lbl Point
ClrIO
Disp "Saisie des coordonnées du point de calcul"
Prompt x
Prompt y
Prompt z
Goto Debut
¨ Calcul et affichage du rotationnel
Lbl Rota
ClrIO
newMat(1,3) » rotat
¶(R,y) - ¶(Q,z) » rotat[1,1]
¶(P,z) - ¶(R,x) » rotat[1,2]
¶(Q,x) - ¶(P,y) » rotat[1,3]
Disp "Le rotationnel est : ", rotat
Goto Debut
¨ Réinitialisation des coordonnées du point de calcul
Lbl Init
DelVar x,y,z
Disp "Coordonnées initialisées"
Goto Debut
Lbl Fin
ClrIO
Disp "FIn du programme Rot"
EndPrgm

Les deux programmes que tu as posté m'interessent! Tu peux me les envoyer ?