Méthode de Halley

Soit f une fonction C² et a un zéro de f. La méthode de Halley consiste à itérer

${\displaystyle x_{n+1}=x_{n}-{\frac {2f(x_{n})f'(x_{n})}{2{[f'(x_{n})]}^{2}-f(x_{n})f''(x_{n})}}}$

à partir d'une valeur x₀ proche de a.

Au voisinage de a, la suite vérifie :

${\displaystyle |x_{n+1}-a|<K|x_{n}-a|^{3}}$,

avec K > 0 ; ce qui signifie que la convergence est donc (au pire) cubique.

La formule se déduit par exemple de la méthode de Newton appliquée à la fonction ${\displaystyle g=f/{\sqrt {f'}}}$ :

${\displaystyle x_{n+1}=x_{n}-{\frac {g(x_{n})}{g'(x_{n})}}}$,

avec

${\displaystyle g'(x)={\frac {2{[f'(x)]}^{2}-f(x)f''(x)}{2f'(x){\sqrt {f'(x)}}}},}$

d'où le résultat. Si f′(c) = 0, cela ne s'applique que si g peut être prolongée en c.