implémentation de Newton: inverse de matrice
Bonjour tous,
ma question est à mi chemin entre l'implementation de la methode de Newton et sa compréhension, j'espere que vous pourrez m'aider...
J'ai une systeme de 3 equations non lineaires à 3inconnues et je veux implementer la resolution de Newton pour trouver ces 3inconnues....
1°) Newton pour plusieurs equations à plusieurs inconnues:
de la meme façon que la methode de Newton classique on montre facilement que le resultat est donné par la formule donnée ici dans le dernier paragraphe:
http://fr.wikipedia.org/wiki/M%C3%A9thode_de_Newton
soit
{Xn+1}={Xn}-([F'(Xn)]^-1){F(Xn)}
2°) ma question est comment implementer cela?
je dois calculer un gradient numeriquement par différence fini et ensuite comment puis je inverser la matrice numeriquement?
J'espere que vous pourrez m'expliquer....
