Capítulo 2 (Sección 2.
4): Ecuaciones No Lineales
Nota: Conviene leer el capítulo 3, página 62 del Libro Shoichiro Nakamura para
complementar los conceptos. Está en la plataforma Paideia, en la carpeta Bibliografía.
2.4 Método de Newton-Raphson
En caso de una ecuación no lineal f(x) = 0, haciendo uso de una aproximación inicial de
x0 cercana a un cero de f(x), se puede usar la Fórmula de Taylor, en torno de x0,
obteniéndose la expresión:
f ( x) f ( x0 ) + f '( x0 )( x − x0 ),
con un error, por la fórmula del Resto:
( x − x0 )2
Error = f ( ( x)) , con ( x) que se ubica entre x y x0
2!
Si en la fórmula de Taylor obtenida se reemplaza a x por , se tiene:
0 = f ( ) f ( x0 ) + f '( x0 )( − x0 )
De donde despejando a , se tiene:
f ( x0 )
x0 − = Nueva aproximación para = x1
f '( x0 )
Y haciendo uso de la Fórmula de Taylor nuevamente, pero empleando ahora como
aproximación inicial al valor x1 , se tendrá:
f ( x1 )
x1 − = Nueva aproximación para = x2
f '( x1 )
De manera similar, se pueden generar las siguientes aproximaciones. En resumen, se
puede definir el algoritmo de Newton para una variable:
f ( xn )
xn +1 = xn − , n 0
f '( xn )
En base a lo expuesto anteriormente, se puede analizar la convergencia del método de
Newton a través del siguiente teorema:
Teorema:
Sea f C2[a, b], [a, b], con f() = 0 y f ‘() 0, entonces > 0 tal que el método de
Newton genera una sucesión {xn}, que converge a , cero de f(x), para cualquier
aproximación inicial x0 [ - , + ].
Comentarios:
1. El Teorema se establece que la aproximación inicial x0, está cercana a , pero no
nos da una idea práctica para determinarla.
2. El método de Newton se conoce también como el método de las tangentes, ya que
xn+1 = xn - f(xn)/ f ‘(xn), determina xn+1 como el punto de intersección de la recta
tangente trazada a la curva y = f(x) en el punto (xn, f(xn)), con el eje X:
Página 1 de 4
Elizabeth Doig
Y
y=f(x)
(xn, f(xn))
xn+1 X
xn
Gráfica 1
3. La interpretación geométrica, dada en el comentario 2 y las condiciones del
Teorema nos permiten analizar los casos cuando el método de Newton converge
a través del siguiente
Dificultades
Como el teorema exige que f ‘(x) 0 en cercanía de , las dificultades surgen cuando
f ( x) = 0 para algún x cercano a ó f ‘() = 0, que es el caso de un cero múltiple de
f(x), condiciones que no siempre son fáciles de detectar.
Algoritmo de Newton-Raphson
Para encontrar una solución de f(x) = 0 dada una aproximación inicial x0 proceder de la
siguiente manera:
ENTRADA: aproximación inicial x0; tolerancia TOL; número máximo de iteraciones N0.
SALIDA: Solución aproximada o mensaje de fracaso.
Paso 1 Tomar i = 1.
Paso 2 Mientras que i N0 seguir pasos 3-6.
Paso 3 Tomar x = x0 - f(x0) / f’(x0). (calcular xi).
Paso 4 Si x - x0 < TOL entonces
SALIDA; (procedimiento completado satisfactoriamente).
PARAR.
Paso 5 Tomar i = i + 1
Paso 6 Tomar x0 = x. (Redefinido x0).
Paso 7 SALIDA (‘El método fracasó después de N0 =’, N0);
(Procedimiento completado sin éxito).
PARAR.
Ejemplo 2
Dada la ecuación no lineal del ejemplo 1:
f ( x) = 1 + xe- x - x = 0 ;
que tiene una raíz en el intervalo [1.2, 1.4], vamos a resolver la ecuación empleando el
método de Newton.
Página 2 de 4
Elizabeth Doig
Lo primero es verificar si se cumplen las condiciones del teorema, observar que:
f ( x) = 1 + xe- x - x es continua en 1.2, 1.4
f ( x) = e − x (1 − x) − 1 es continua en 1.2, 1.4
Si x0 = 1.3, entonces: f (1.3) = −1.081759538 0
Por lo que el método de Newton si converge dentro del intervalo, por lo tanto, se define
el algoritmo:
f ( xn )
xn +1 = xn − , n 0
f ( xn )
De donde se obtiene el siguiente conjunto de iteraciones:
n xn
0 1.3
1 1.350187985
2 1.349976489
3 1.349976485
4 1.349976485
Tabla 1
Se observa la rapidez de convergencia del método debido a su convergencia cuadrática,
al compararlo con el método de Punto Fijo que posee convergencia lineal.
Observación:
A continuación, se muestra la ventana del software MATLAB, para el programa
“newton1var.m”:
Página 3 de 4
Elizabeth Doig
*Recordar: Criterios de Tolerancia
Para detener el algoritmo se pueden utilizar cualesquiera de los siguientes criterios:
1. xn +1 − x n
xn+1− xn
2.
xn+1
3. f ( xn )
4. Fijar un número máximo de iteraciones n = N, para interrumpir el algoritmo.
Página 4 de 4
Elizabeth Doig
