0% encontró este documento útil (0 votos)
37 vistas52 páginas

Métodos para Hallar Raíces de Ecuaciones

Cargado por

eia
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
37 vistas52 páginas

Métodos para Hallar Raíces de Ecuaciones

Cargado por

eia
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Raices 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.

Salvo en casos muy simples es prácticamente imposible hallar una solución


exacta a estas ecuaciones.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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.

Salvo en casos muy simples es prácticamente imposible hallar una solución


exacta a estas ecuaciones.

En general la teoría matemática nos puede decir si estas soluciones existen


y más o menos donde se encuentran, pero no nos dice como hallarlas.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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.

Salvo en casos muy simples es prácticamente imposible hallar una solución


exacta a estas ecuaciones.

En general la teoría matemática nos puede decir si estas soluciones existen


y más o menos donde se encuentran, pero no nos dice como hallarlas.

Se buscan métodos que permitan hallar soluciones aproximadas y tener


una cota para el error.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

con un error menor que 10−5

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 3 / 30
Método Babilónico

En la tableta YBC 7289 (c. 18001600 AC)

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

paso valor error


0 1 0.5
1 1.5 0.08333333333333348
2 1.4166666666666665 0.002450980392156854
3 1.4142156862745097 2.1238998197947723e-06

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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.

Supongamos por ejemplo que f (a0 ) < 0 y f (b0 ) > 0. Sea L = b0 − a0 .

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

primer aproximación a la raíz. El error es menor que L2 .

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

primer aproximación a la raíz. El error es menor que L2 .


Paso 2 Miramos el signo de f (s0 ). Si es positivo tiene que haber una
raíz entre a0 y s0 . En ese caso llamamos a1 = a0 y b1 = s0 .
Si es negativo que haber una raíz entre s0 y b0 . En ese caso
llamamos a1 = s0 y b1 = b0 . Ahora repetimos el paso 1. Es
decir tomamos s1 = b1 +a 2
1
y ahora el error es menor que L4 .

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

primer aproximación a la raíz. El error es menor que L2 .


Paso 2 Miramos el signo de f (s0 ). Si es positivo tiene que haber una
raíz entre a0 y s0 . En ese caso llamamos a1 = a0 y b1 = s0 .
Si es negativo que haber una raíz entre s0 y b0 . En ese caso
llamamos a1 = s0 y b1 = b0 . Ahora repetimos el paso 1. Es
decir tomamos s1 = b1 +a 2
1
y ahora el error es menor que L4 .
Paso 3 Repetir hasta que el error sea menos que el pedido.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 9 / 30
Convergencia

Siguiendo los pasos del método obtenemos tres sucesiones an , bn y sn tales


que
1 an ≤ an+1 < b0

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia

Siguiendo los pasos del método obtenemos tres sucesiones an , bn y sn tales


que
1 an ≤ an+1 < b0
2 a0 < bn+1 ≤ bn

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia

Siguiendo los pasos del método obtenemos tres sucesiones an , bn y sn tales


que
1 an ≤ an+1 < b0
2 a0 < bn+1 ≤ bn
3 an < sn < bn

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia

Siguiendo los pasos del método obtenemos tres sucesiones an , bn y sn tales


que
1 an ≤ an+1 < b0
2 a0 < bn+1 ≤ bn
3 an < sn < bn
4 |bn − an | → 0

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Convergencia

Siguiendo los pasos del método obtenemos tres sucesiones an , bn y sn tales


que
1 an ≤ an+1 < b0
2 a0 < bn+1 ≤ bn
3 an < sn < bn
4 |bn − an | → 0
5 Entonces existe un único c tal que an → c, bn → c y por lo tanto
sn → c. Como f es continua f (c) = 0

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 10 / 30
Error

Que ocurre con el error


Sea L = b0 − a0 entonces
1 En el paso 0 el error es menor que L/2.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Error

Que ocurre con el error


Sea L = b0 − a0 entonces
1 En el paso 0 el error es menor que L/2.
2 En el paso 1 el error es menor que L/4.

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Error

Que ocurre con el error


Sea L = b0 − a0 entonces
1 En el paso 0 el error es menor que L/2.
2 En el paso 1 el error es menor que L/4.
3 En el paso n el error es menor que L/2n+1 .

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 11 / 30
Como acotar el error

En el algoritmo del método de bisección se ve que si r es el verdadero


valor buscado, entonces
L
|r − sn | ≤ n+1
2

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 12 / 30
Como acotar el error

En el algoritmo del método de bisección se ve que si r es el verdadero


valor buscado, entonces
L
|r − sn | ≤ n+1
2

Si queremos asegurar que el error cometido es menor que una constante


ϵ > 0 hay que tomar n tal que
L

2n+1

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 12 / 30
Como acotar el error

En el algoritmo del método de bisección se ve que si r es el verdadero


valor buscado, entonces
L
|r − sn | ≤ n+1
2

Si queremos asegurar que el error cometido es menor que una constante


ϵ > 0 hay que tomar n tal que
L

2n+1

Lo que es equivalente a pedir que


L
< 2n+1
ϵ
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 12 / 30
Formula para el error

o que  
L ln(x )
ln2 −1<n Recordemos que ln2 (x ) =
ϵ ln(2)

Fórmula final para acotar el error


Si f es una función continua, a < b, f (a)f (b) < 0, y llamamos L = b − a.
Si queremos hallar la solución de f (x ) = 0 con un error menor que un
número ϵ > 0 por el método de bisección deberemos calcular sn con
   
L
n = ln2 −1
ϵ

donde para un número real positivo x el número [x ] denota un numero


natural n tal que n − 1 ≤ x < n
. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
Raices de ecuaciones 13 de enero de 2023 13 / 30
Ejemplo 1

Queremos calcular 3 con un error menor a 0.01. Tomamos
f (x ) = x 2 − 3, a = 1, b = 2.

Cantidad de pasos
Debemos calcular sn con

n = [ln2 (1/0,01) − 1] = [6,643856189774725 − 1]


= [5,643856189774725] = 6

n a b f(a) f(b) c f(c) new error


0 1.0 2.0 -2.0 1.0 1.5 -0.75 a=c 0.5
1 1.5 2.0 -0.75 1.0 1.75 0.0625 b=c 0.25
2 1.5 1.75 -0.75 0.0625 1.625 -0.35938 a=c 0.125
3 1.625 1.75 -0.35938 0.0625 1.6875 -0.15234 a=c 0.0625
4 1.6875 1.75 -0.15234 0.0625 1.7188 -0.045898 a=c 0.03125
5 1.7188 1.75 -0.045898 0.0625 1.7344 0.0080566 b=c 0.015625
6 1.7188 1.7344 -0.045898 0.0080566 1.7266 -0.018982 a=c 0.0078125

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

n a b f(a) f(b) c f(c) new error


0 3.0 4.0 0.047127 -0.03837 3.5 -0.019757 b=c 0.5
1 3.0 3.5 0.047127 -0.019757 3.25 0.0058487 a=c 0.25
2 3.25 3.5 0.0058487 -0.019757 3.375 -0.0086811 b=c 0.125
3 3.25 3.375 0.0058487 -0.0086811 3.3125 -0.001877 b=c 0.0625
4 3.25 3.3125 0.0058487 -0.001877 3.2812 0.001867 a=c 0.03125
5 3.2812 3.3125 0.001867 -0.001877 3.2969 -3.4225e-05 b=c 0.015625
6 3.2812 3.2969 0.001867 -3.4225e-05 3.2891 0.00090903 a=c 0.0078125
7 3.2891 3.2969 0.00090903 -3.4225e-05 3.293 0.00043557 a=c 0.0039062
8 3.293 3.2969 0.00043557 -3.4225e-05 3.2949 0.00020021 a=c 0.0019531
9 3.2949 3.2969 0.00020021 -3.4225e-05 3.2959 8.288e-05 a=c 0.00097656

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

El método de bisección tiene problemas cuando el cero está en un mínimo


o máximo local de f . Por ejemplo si tomamos

f (x ) = x 2 cos(2x )

Tenemos el siguiente gráfico


y
0
−1 −0, 5 0 0, 5 1
x
−0, 25

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

Supongamos que f (p) = 0 y que x es un punto próximo a p. Podemos


hacer un desarrollo de Taylor alrededor de x

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

Es el más usado en general.


Teorema de Newton-Raphson
Sea f una función de clase C 2 en [a, b] si existe a < p < b tal que
f (p) = 0 y f ′ (p) ̸= 0 entonces existe δ > 0 tal que si |x0 − p| < δ,
entonces la sucesión
f (xn )
xn+1 = xn − ′
f (xn )
converge a p.

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 .

Sea m = (f (x0 ) − f (sn ))/(x0 − sn ) entonces m es la pendiente de la recta


que une los puntos (sn , f (sn )) con (x0 , 0). Por el Teorema de Lagrange
existe sn < c < x0 tal que

m = (f (x0 ) − f (sn ))/(x0 − sn ) = −f (sn )/(x0 − sn ) = f ′ (c)

como f ′′ (x ) < 0 entonces f ′ (sn ) > f ′ (c)

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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 .

Sea m = (f (x0 ) − f (sn ))/(x0 − sn ) entonces m es la pendiente de la recta


que une los puntos (sn , f (sn )) con (x0 , 0). Por el Teorema de Lagrange
existe sn < c < x0 tal que

m = (f (x0 ) − f (sn ))/(x0 − sn ) = −f (sn )/(x0 − sn ) = f ′ (c)

como f ′′ (x ) < 0 entonces f ′ (sn ) > f ′ (c)

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

Tenemos entonces una sucesión monótona creciente sn acotada


superiormente por x0 . Por lo tanto existe s = lı́mn→∞ sn . Entonces como f
es continua y f ′ (s) ̸= 0

f (sn ) f (s)
0 = lı́m (sn+1 − sn ) = lı́m ′
= ′
n→∞ n→∞ f (sn ) f (s)

por lo tanto f (s) = 0 y s = x0 .

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

sn2 − 2 2s 2 − sn2 + 2 s2 + 2 (sn + (2/sn ))


sn+1 = sn − = n = n =
2sn 2sn 2sn 2
ąEsta es la sucesión que aparece en el método Babilónico!

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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

Se puede ver que f ′′ (x ) < 0, ∀x ∈ [6, 8]. Por lo tanto f ′ (x ) es decreciente.


Como f ′ (6) ≈ −0,7, entonces f ′ (x ) < 0, ∀x ∈ [6, 8].

. . . . . . . . . . . . . . . . . . . .
. . . . . . . . . . . . . . . . . . . .
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 puede ver que f ′′ (x ) < 0, ∀x ∈ [6, 8]. Por lo tanto f ′ (x ) es decreciente.


Como f ′ (6) ≈ −0,7, entonces f ′ (x ) < 0, ∀x ∈ [6, 8].

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

También podría gustarte