Méthode de Halley
En analyse numérique, la méthode de Halley est un algorithme de recherche d'un zéro d'une fonction utilisé pour les fonctions d'une variable réelle dérivables deux fois et à dérivée seconde continue (i.e. C2). La méthode, présentée par l'astronome Edmond Halley en 1694[1], est une généralisation de la méthode de Newton, à convergence cubique.
Énoncé
Soit une fonction de classe C² et un zéro de . La méthode de Halley consiste à itérer
à partir d'une valeur proche de .
Au voisinage de , la suite vérifie :
- ,
avec ; ce qui signifie que la convergence est donc (au pire) cubique.
La formulation suivante montre que la méthode de Halley est un raffinement de celle de Newton, et très proche de celle-ci si .
On remarque que dans cette formulation l'expression n'est calculée qu'une fois.
Obtention à partir de la méthode de Newton
La formule se déduit de la méthode de Newton appliquée à la fonction ; en effet, celle-ci s'écrit :
- ,
avec
d'où le résultat. Si , ceci ne s'applique que si peut être prolongée par continuité en .
Application au calcul d'une fonction réciproque
Appliquée à la fonction , cette méthode conduit à une suite récurrente convergeant vers .
Ceci peut être utilisé pour le calcul d'une racine n-ième, avec ; la suite récurrente définie par converge vers .
La méthode de Halley est aussi utilisée pour calculer un logarithme népérien, en inversant la fonction exponentielle.
Notes et références
Voir aussi
Liens internes
Liens externes
- Modèle:MathWorld
- Modèle:En Newton's method and high order iterations, Pascal Sebah et Xavier Gourdon, 2001