RAÍCES MÚLTIPLES
Una raíz múltiple corresponde a un punto donde la función es tg al eje x, o en el
punto de intersección se corresponde con un punto de inflexión.
Por ejemplo:
𝑓(𝑥) = (𝑥 − 3)(𝑥 − 1)(𝑥 − 1) =(𝑥 − 3)(𝑥 − 1)2 o 𝑓(𝑥) = 𝑥 3 − 5𝑥 2 + 7𝑥 − 3
𝑓(𝑥) = (𝑥 − 3)(𝑥 − 1)(𝑥 − 1)(𝑥 − 1) =(𝑥 − 3)(𝑥 − 1)3 o
𝑓(𝑥) = 𝑥 4 − 6𝑥 3 + 12𝑥 2 − 10𝑥 + 3
f(x)
Raíz doble
x
f(x)
Raíz cuádruple
Raíz triple
X
En general la multiplicidad impar cruza al eje x y la multiplicidad par no lo cruza.
El hecho que la función no cambia de signo en una raíz de multiplicidad par
impide el uso de los métodos que usan intervalos. Y los métodos abiertos
tienen la dificultad que pueden ser divergentes.
Otro posible problema es el hecho que no sólo la función f(x) se aproxima a
cero sino también sus posibles derivadas.
Se puede demostrar que el método de Newton-Raphson y el de la secante
convergen en forma lineal, en vez de manera cuadrática, cuando hay raíces
múltiples.
MÉTODO DE NEWTON-RAPHSON MODIFICADO
Una alternativa sugerida por Ralston y Rabinowitz (1978) es la de definir una
nueva función u(x) que es el cociente de la función f(x) y su derivada f ’(x), esto
es:
𝒇(𝒙)
𝒖(𝒙) = 𝒇 ′(𝒙)
Esta función tiene raíces en las mismas posiciones que la función original, por
lo tanto esta función se sustituye en la fórmula de Newton-Raphson dando
lugar a una forma alternativa de ella:
𝒖(𝒙 ) 𝒇′ (𝒙)𝒇′ (𝒙)−𝒇(𝒙)𝒇′′(𝒙)
𝒙𝒊+𝟏 = 𝒙𝒊 − 𝒖 ′(𝒙𝒊 ) Teniendo en cuenta que 𝒖 ′(𝒙) =
𝒊 [𝒇′ (𝒙)]𝟐
reemplazando Obtenemos:
𝒇(𝒙𝒊 )𝒇′(𝒙𝒊 )
𝒙𝒊+𝟏 = 𝒙𝒊 −
[𝒇′ (𝒙𝒊 )]𝟐 − 𝒇(𝒙𝒊 )𝒇′′(𝒙𝒊 )
Que es la fórmula de Newton-Raphson modificado para raíces múltiples.
Ejemplo
Úsense los dos métodos, el estándar y el modificado para calcular la raíz
múltiple de : 𝑓(𝑥) = 𝑥 3 − 5𝑥 2 + 7𝑥 − 3 usando el valor inicial 𝑥0 = 0
𝑓 ′ (𝑥) = 3𝑥 2 − 10𝑥 + 7
𝑥𝑖 3 −5𝑥𝑖 2 +7𝑥𝑖 −3
La fórmula estándar queda entonces: 𝑥𝑖+1 = 𝑥𝑖 − 3𝑥𝑖 2 −10𝑥𝑖 +7
𝑥0 = 0
𝑥1 =0,428571
𝑥2 = 0,685714 𝜀𝑎 % = 38%
𝑥3 = 0,832865 𝜖𝑎 % = 17%
𝑥4 = 0,913328 𝜀𝑎 %= 9 %
𝑓 ′′ (𝑥) = 6𝑥 − 10 con lo que la fórmula de Newton-Raphson modificado queda:
(𝑥𝑖 3 − 5𝑥𝑖 2 + 7𝑥𝑖 − 3)(3𝑥𝑖 2 − 10𝑥𝑖 + 7)
𝑥𝑖+1 = 𝑥𝑖 −
(3𝑥𝑖 2 − 10𝑥𝑖 + 7)2 − (𝑥𝑖 3 − 5𝑥𝑖 2 + 7𝑥𝑖 − 3)(6𝑥𝑖 − 10)
𝑥0 = 0
𝑥1 = 1,105263
𝑥2 = 1,003081 𝜀𝑎 % = 10%
𝑥3 = 1,000002 𝜀𝑎 %= 0,3%