El método iterativo de referencia: linealizar f en el iterado actual y saltar a la raíz de la tangente. Deducción completa por tres caminos y demostración del orden cuadrático con su ecuación del error.
La fórmula
xk+1=xk−f′(xk)f(xk),k=0,1,2,…
Método de Newton-Raphson.
Geométricamente, cada iteración sustituye la curva por su recta tangente en (xk,f(xk)) y toma como siguiente aproximación el punto donde esa tangente corta el eje. Es un método de punto fijo con ϕ(x)=x−f′(x)f(x), de un punto, sin memoria y con derivadas.
Escribimos f mediante el Teorema Fundamental del Cálculo desde xk y aproximamos la integral con un rectángulo (integrando constante f′(xk), la misma idea que en la deducción de Euler):
f(x)=f(xk)+∫xkxf′(t)dt≈f(xk)+f′(xk)(x−xk)
Evaluamos en x=α, donde f(α)=0, y despejamos α:
0≈f(xk)+f′(xk)(α−xk)⇒α≈xk−f′(xk)f(xk)
Esa aproximación de α es el siguiente iterado. Aproximar la integral con cuadraturas más ricas (trapecio, punto medio, Simpson) produce métodos de orden mayor por esta misma vía.
Demostración del orden 2 (ecuación del error)
Sea α raíz simple (f(α)=0, f′(α)=0), ek=xk−α el error y c2=21f′(α)f′′(α). Desarrollamos f(xk) por Taylor en torno a α; el término constante f(α) desaparece:
Restamos α en la fórmula de Newton y sustituimos: los términos lineales en ek se cancelan y queda la ecuación del error cuadrática. El método de Newton tiene orden p=2.
ek+1=xk+1−α=ek−f′(xk)f(xk)=c2ek2+O(ek3)
Más allá de Newton
Newton es el punto de partida de casi todo lo demás: sustituyendo la derivada por diferencias se obtienen la secante y Steffensen; componiendo y añadiendo funciones peso se construyen los métodos de alto orden (Traub, Ostrowski, Jarratt); y su versión vectorial resuelve sistemas de ecuaciones no lineales. También es la herramienta estándar para las ecuaciones implícitas que aparecen en los métodos implícitos para EDO.