Résolution d'un système de 3 équations non-linéaires

Bonjour,
pour le compte de mon projet technique, je dois resoudre grace a fortran un systeme d'equation non lineaire a 3 inconnues (theta1, theta2, v)a l'aide de 3 equations de la forme suivante:

k1=(x1*sin(theta1)*cos(theta2)+y1*sin(theta1)*sin(theta2)+z1*cos(theta1))/v
k2=(x2*sin(theta1)*cos(theta2)+y2*sin(theta1)*sin(theta2)+z2*cos(theta1))/v
k3=(x3*sin(theta1)*cos(theta2)+y3*sin(theta1)*sin(theta2)+z3*cos(theta1))/v

Les valeurs de ki, xi, yi et zi sont connues.

En fait, cette resolution fait partie d'un tres grand programme fortran realise avant que je commence ce projet.Je n'y connais vraiment pas grand chose, j'ai tente de lire a peu pres tout les document mis a disposition sur le net mais rien ne colle avec ce que je veux.

Faut-il creer une subroutine? Dans une publi, j'ai lu que les auteurs avaient utilise un algorithm de Newton-Raphson pour resoudre ce genre de systeme.

Si quelqu'un pouvait me sortir de ce vilain cauchemard je lui en serai fort reconnaissant.

Merci d'avance, je me tiens a votre disposition pour tout eclaircissement de ce probleme.

Salut!
Avant de programmer quoi que ce soit, il vaut la peine de réfléchir un peu. Ne vois-tu pas que tes équations sentent à plein nez les coordonnées sphériques et que, moyennant un changement de variable judicieux, la solution est évidente.
Alors bon travail!
Jean-Marc Blanc