Métodos para Hallar Raíces de Ecuaciones
Métodos para Hallar Raíces de Ecuaciones
13 de enero de 2023
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 1 / 30
Los métodos que vamos a ver se usan para hallar raíces de ecuaciones del
tipo
f (x ) = 0
donde f es una función.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 2 / 30
Los métodos que vamos a ver se usan para hallar raíces de ecuaciones del
tipo
f (x ) = 0
donde f es una función.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 2 / 30
Los métodos que vamos a ver se usan para hallar raíces de ecuaciones del
tipo
f (x ) = 0
donde f es una función.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 2 / 30
Los métodos que vamos a ver se usan para hallar raíces de ecuaciones del
tipo
f (x ) = 0
donde f es una función.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 2 / 30
La raíz cuadrada de 2
Problema
Queremos hallar la solución de
x2 − 2 = 0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 3 / 30
Método Babilónico
se encuentra la solución
x = 1,41421296
que es valida hasta el 6 decimal. . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 4 / 30
Idea
Si x 2 = 2 entonces x = 2/x , Si x0 es un número, entonces el promedio
entre x0 y 2/x0 debería estar más cerca de la solución que x0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 5 / 30
Idea
Si x 2 = 2 entonces x = 2/x , Si x0 es un número, entonces el promedio
entre x0 y 2/x0 debería estar más cerca de la solución que x0
Sucesión
Tomamos x0 = 1 y
2
xn + xn xn 1
xn+1 = = +
2 2 xn
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 5 / 30
Idea
Si x 2 = 2 entonces x = 2/x , Si x0 es un número, entonces el promedio
entre x0 y 2/x0 debería estar más cerca de la solución que x0
Sucesión
Tomamos x0 = 1 y
2
xn + xn xn 1
xn+1 = = +
2 2 xn
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 5 / 30
Algorithm 1 Método Babilónico
procedure Babilonic(xin, cota)
x ← xin ▷ Valor inicial
ϵ ← cota ▷ cota del error
while True do
y ← x /2 + 1/x
err ← |y − x |
if err < ϵ then
Break
end if
x ←y ▷ Nuevo valor de x
end while
Return y
end procedure
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 6 / 30
Otros métodos
Método de incrementos
Sea f (x ) = x 2 − 2. Como f (1) = −1 y f (2) = 2 sabemos que hay una
solución entre 1 y 2 (Teorema de Bolzano). Que hacemos?
Calculamos f (1,001), f (1,00002), f (1,00003), . . . hasta que cambie de
signo y tomamos como solución la última de signo negativo.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 7 / 30
Otros métodos
Método de incrementos
Sea f (x ) = x 2 − 2. Como f (1) = −1 y f (2) = 2 sabemos que hay una
solución entre 1 y 2 (Teorema de Bolzano). Que hacemos?
Calculamos f (1,001), f (1,00002), f (1,00003), . . . hasta que cambie de
signo y tomamos como solución la última de signo negativo.
Uso general
Se puede usar en general.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 7 / 30
Otros métodos
Método de incrementos
Sea f (x ) = x 2 − 2. Como f (1) = −1 y f (2) = 2 sabemos que hay una
solución entre 1 y 2 (Teorema de Bolzano). Que hacemos?
Calculamos f (1,001), f (1,00002), f (1,00003), . . . hasta que cambie de
signo y tomamos como solución la última de signo negativo.
Uso general
Se puede usar en general.
Problema
Hay que hacer MUCHOS pasos en general. En el caso de la raíz alrededor
de 50000.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 7 / 30
Método de bisección
Teorema (Bolzano)
Supongamos que f es una función continua y que encontramos números
a0 < b0 tales que f (a0 ) y f (b0 ) tienen distinto signo. Entonces existe
a0 < c < b0 tal que f (c) = 0.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 8 / 30
Método de Bisección algoritmo
Algoritmo
Vamos a definir una sucesión sn que converge a la solución.
Paso 1 Tomamos s0 = b0 +a
2 . El promedio entre a0 y b0 . Es la
0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 9 / 30
Método de Bisección algoritmo
Algoritmo
Vamos a definir una sucesión sn que converge a la solución.
Paso 1 Tomamos s0 = b0 +a
2 . El promedio entre a0 y b0 . Es la
0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 9 / 30
Método de Bisección algoritmo
Algoritmo
Vamos a definir una sucesión sn que converge a la solución.
Paso 1 Tomamos s0 = b0 +a
2 . El promedio entre a0 y b0 . Es la
0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 9 / 30
Convergencia
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Error
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Error
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Error
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Como acotar el error
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 12 / 30
Como acotar el error
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 12 / 30
Como acotar el error
o que
L ln(x )
ln2 −1<n Recordemos que ln2 (x ) =
ϵ ln(2)
Cantidad de pasos
Debemos calcular sn con
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 14 / 30
Ejemplo 2
Hallar la solución de exp(−x ) ∗ (3,2 ∗ sin(x ) − 0,5 ∗ cos(x )) = 0, con un
error menor a 0,001. Tomamos a = 3, b = 4
Pasos a realizar
Tenemos n = [ln2 (1/0,001) − 1]] = [8,965784284662087] = 9
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 15 / 30
El método de bisección
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 16 / 30
El método de bisección
−10 −5 0 5 10
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 17 / 30
Limitaciones del método
f (x ) = x 2 cos(2x )
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 18 / 30
Metodo de Newton Raphson
2.25 4 x
Figura: Idea
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 19 / 30
Deducción del método
f ′′ (x )
f (p) = f (x ) + f ′ (x )(p − x ) + (p − x )2 + O(|p − x |3 )
2
Como (p − x )2 es mucho más cercano a 0 que (p − x ) podemos
despreciarlo en primera aproximación. Entonces se obtiene
0 = f (x ) + f ′ (x )(p − x )
Despejando p tenemos
f (x )
p=x− .
f ′ (x )
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 20 / 30
Método de Newton Rhaphson II
Uso práctico
El método de Newton Raphson tiene algunas desventajas. Hay que hallar
primero una raíz aproximada y en general hay que ver que f ′ (x ) ̸= 0 en el
intervalo en consideración.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 21 / 30
Teorema
Sea f una función de clase C 2 en [a, b] si existe a < x0 < b tal que
f (x0 ) = 0, y los signos de f ′ y f ′′ son constantes entonces la sucesión
generada por el método de Newton Raphson converge a cero único de f en
(a, b) en forma monótona si se elige como valor inicial s0 al extremo de
[a, b] donde sign(f ′′ ) y sign(f ) son iguales en ese extremo.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 22 / 30
Explicación gráfica
Suposiciones
Suponemos que f ′ (x ) > 0 y que f ′′ (x ) < 0 para todo x ∈ [a, b]. Tomamos
s0 = a.
sn sn+1
x0
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 23 / 30
Demostración
Por inducción
Si sn < x0 entonces f (sn ) < 0, por lo tanto −f (sn )/f ′ (sn ) > 0. Entonces
sn+1 = sn + (−f (sn )/f ′ (sn )) es sn mas algo positivo. Por lo tanto
sn+1 > sn .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 24 / 30
Demostración
Por inducción
Si sn < x0 entonces f (sn ) < 0, por lo tanto −f (sn )/f ′ (sn ) > 0. Entonces
sn+1 = sn + (−f (sn )/f ′ (sn )) es sn mas algo positivo. Por lo tanto
sn+1 > sn .
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 24 / 30
Demostración
Por inducción
Si sn < x0 entonces f (sn ) < 0, por lo tanto −f (sn )/f ′ (sn ) > 0. Entonces
sn+1 = sn + (−f (sn )/f ′ (sn )) es sn mas algo positivo. Por lo tanto
sn+1 > sn .
Finalmente
f (sn ) f (sn )
x0 = sn − > sn − ′ = sn+1
m f (sn ) . . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 24 / 30
Continuación
f (sn ) f (s)
0 = lı́m (sn+1 − sn ) = lı́m ′
= ′
n→∞ n→∞ f (sn ) f (s)
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 25 / 30
Convergencia
Teorema
Supongamos que el método de Newton-Raphson genera una sucesión xn
que converge a una raíz p de f . Sea En = |xn − p| entonces
1 Si p es una raíz simple entonces
f ′′ (p) 2
En+1 ≃ E
2f ′ (p) n
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 26 / 30
Convergencia
Teorema
Supongamos que el método de Newton-Raphson genera una sucesión xn
que converge a una raíz p de f . Sea En = |xn − p| entonces
1 Si p es una raíz simple entonces
f ′′ (p) 2
En+1 ≃ E
2f ′ (p) n
2 Si p es una raíz de orden m entonces
m−1
En+1 ≃ En
m
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 26 / 30
Ejemplo
√
Cálculo de 2
Vamos a calcular la raíz cuadrada de 2. Tomamos f (x ) = x 2 − 2 y
[a, b] = [1, 2]. Entonces f es de clase C 2 , f ′ (x ) = 2x > 0 en [1, 2] y como
f (2)f ′′ (2) > 0. Tomamos s0 = 2. La sucesión de Newton-Raphson es
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 27 / 30
Ejemplo 2
Problema
Hallar una solución (aproximada) de la ecuación
f (x ) = 1,8x + 10 ln(1 − 0,1x ) = 0.
Solución
Vale que f (6) ≈ 1,6370 y f (8) ≈ −1,69. Por lo tanto hay una solución
entre 6 y 8.
2
1.5
1
1.8*x+10*log(1-0.1*x)
0.5
-0.5
-1
-1.5
-2
6 6.5 7 7.5 8
x
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 28 / 30
Continuación
Vale que
f ′ (x ) = 1,8 − 1,0/(1 − 0,1 ∗ x )
y
0,1
f ′′ (x ) = −
(1 − 0,1 x )2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 29 / 30
Continuación
Vale que
f ′ (x ) = 1,8 − 1,0/(1 − 0,1 ∗ x )
y
0,1
f ′′ (x ) = −
(1 − 0,1 x )2
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 29 / 30
Continuación
Vale que
f ′ (x ) = 1,8 − 1,0/(1 − 0,1 ∗ x )
y
0,1
f ′′ (x ) = −
(1 − 0,1 x )2
Se cumplen las hipótesis del Teorema. Como f (8) tiene el mismo signo
que f ′′ (8) tomamos como s0 = 8.
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 29 / 30
Conclusión
sn f (sn ) En 2
En /En−1
8 -1.694379124341
7.470506523 -0.298748430165 0.52949347
7.331770623 -0.014512831490 0.13873590 0.4948
7.324319753 -3.891616602e-5 0.00745086 0.3871
7.324299667 -2.817905909e-10 2.008713e-05 0.3618
7.324299666 0.0 -1.454525389e-10 0.3605
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 30 / 30