ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC.

VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

METODOS NUMERICOS CAPITULO 1

RESOLUCION DE ECUACIONES
1. Método grafico
Método Gráfico Es un método simple para obtener una aproximación a la raíz de
la ecuación. Consiste en graficar la función y observar dónde cruza el eje x. Este
punto, que representa el valor de x para el cual, ofrece una aproximación inicial
de la raíz.
Un método simple para obtener una aproximación a la raíz de la ecuación f(x) =
0 consiste en graficar la función y observar en donde cruza el eje x. Este punto,
que representa el valor de x para la cual f(x) = 0, proporciona una aproximación
inicial de la raíz. Fig. 1 La aproximación gráfica para determinar las raíces de una
ecuación

Lo esencial en este método es poder construir un modelo gráfico de la ecuación

y luego por inspección estimar una aproximación a la raíz. El mayor
inconveniente de éste método es su poca precisión y exactitud. Sin embargo,
hoy día se cuenta con excelentes herramientas de software para realizar
rápidamente gráficas con un alto grado de realismo Fig. 2 Método gráfico para
hallar raíces.
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

Las interpretaciones gráficas, además de proporcionar aproximaciones iniciales

de la raíz, son herramientas importantes en la comprensión de las propiedades
de las funciones, previendo las fallas de los métodos numéricos.
Las raíces de la ecuación:
𝑓 (𝑥 ) = 0
Pueden determinarse en forma aproximada considerando las abscisas de los
puntos de intersección de la gráfica de la función y = f(x) con el eje x.
Resulta aconsejable a veces sustituir la ecuación dada por una ecuación
equivalente (dos ecuaciones se denominan equivalentes si tienen exactamente
las mismas raíces):
∅(𝑥 ) = 𝜑(𝑥)
Donde las funciones ∅(𝑥 ) 𝑦 𝜑(𝑥) son más sencillas que f(x), posteriormente se
grafica las funciones 𝑦 = ∅(𝑥 ) 𝑒 𝑦 = 𝜑(𝑥) y las raíces deseadas serán entonces
las abscisas de los puntos de intersección de esas gráficas.
EJERCICIO1.
Resuelva gráficamente la siguiente ecuación:
𝑓 (𝑥 ) = 𝑒 −𝑥 − 𝑥
Solución:
Procedemos a realizar una tabla para poder graficar la función:
Para X=0.
𝑦 = 𝑓(0) = 𝑒 0 − 0 →→ 𝑓 (0) = 1
Para X=0.5
𝑦 = 𝑓 (0.5) = 𝑒 0.5 − 0.5 →→ 𝑓(0.5) = 1.1487
Procedemos con los demás valores entonces tenemos:
𝑓 (𝑥 ) = 𝑒 −𝑥 − 𝑥
X
0 1
0.5 1.1487
1 -0.6321
1.5 -1.2769
2 -1.8647
2.5 -2.4179
3 -2.9502
3.5 -3.4698
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

Revisando la gráfica tenemos

X = 0.56

EJERCICIO2
Resuelva gráficamente la siguiente ecuación:
𝑥𝑙𝑜𝑔(𝑥 ) = 1
SOLUCION:
Primeramente, procederemos a rescribir esta ecuación de la siguiente manera
1
log(𝑥) =
𝑥
las raíces pueden hallarse fácilmente, ya que son las abscisas de los puntos de
intersección de:

𝑐𝑢𝑟𝑣𝑎 𝑝𝑜𝑙𝑖𝑛𝑜𝑚𝑖𝑐𝑎 𝑦 = log(𝑥 ) 𝑦 𝑙𝑎 ℎ𝑖𝑝𝑒𝑟𝑏𝑜𝑙𝑎 𝑦 = 1/𝑥

Para graficar las ecuaciones tenemos las siguientes tablas:

x y= log(x) x y = 1/x
0.1 -1 -1 -1
1 0 1 1
2 0.3010 2 0.5
3 0.4771 3 0.3333
4 0.6021 4 0.25
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

5 0.6990 5 0.2
6 0.7782 6 0.1667

Entonces viendo la gráfica el valor aproximado de la única raíz es:

𝒙 = 𝟐. 𝟓

2. Método de la bisección
Paso1. Elegir valores iniciales a (inferior) y b (superior) comprobar que se cumple
lo siguiente:
𝑓 (𝑎 ) ∗ 𝑓 (𝑏 ) < 0
Paso2. Una aproximación de la raíz se determina mediante:
𝑎+𝑏
𝑚=
2
Paso3. Realice las siguientes evaluaciones para determinaren que subintervalo
se encuentra la raíz.
𝑎) 𝑆𝑖 𝑓 (𝑎) ∗ 𝑓 (𝑚) < 0 𝑐𝑎𝑚𝑏𝑖𝑎𝑟 𝑏 = 𝑚
𝑎) 𝑆𝑖 𝑓 (𝑎) ∗ 𝑓(𝑚) > 0 𝑐𝑎𝑚𝑏𝑖𝑎𝑟 𝑎 = 𝑚
Paso4.
𝑆𝑖 𝑓 (𝑎) ∗ 𝑓 (𝑚) = 0
Entonces m es la raíz aproximada buscada
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

EJERCICIO 1.
𝑓 (𝑥 ) = 4𝑥 2 − 5𝑥 [1, 1.6]
Solución:
Primero se debe verificar que realmente existe solo una solución en este
intervalo, esto se lo hace analizando si existe un cambio de signo en los extremos
del intervalo.
𝑓 (𝑥 ) = 4𝑥 2 − 5𝑥
𝑓 (𝑎) = 𝑓(1) = 4 ∗ (1)2 − 5 ∗ (1) = −1
𝑓 (𝑏) = 𝑓(1.6) = 4 ∗ (1.6)2 − 5 ∗ 1.6 = 2.24
𝑎+𝑏 1 + 1.6
𝑚= → 𝑚= → 𝑚 = 1.3
2 2
Valor a
K a b m f(a) f(b) f(m)
cambiar
0 1.0000 1.6000 1.3000 -1.0000 2.2400 0.2600 b
1 1.0000 1.3000 1.1500 -1.0000 0.2600 -0.4600 a
2 1.1500 1.3000 1.2250 -0.4600 0.2600 -0.1225 a
3 1.2250 1.3000 1.2625 -0.1225 0.2600 0.0631 b
4 1.2250 1.2625 1.2438 -0.1225 0.0631 -0.0311 a
5 1.2438 1.2625 1.2531 -0.0311 0.0631 0.0158 b
6 1.2438 1.2531 1.2485 -0.0311 0.0158 -0.0077 a
7 1.2485 1.2531 1.2508 -0.0077 0.0158 0.0040 b

𝑓(a)𝑓(m) = 0.000030819 ≈ 0

La raíz buscada es:

𝑥 = 1.2508

EJERCICIO 2.
Hallar la solución de la ecuación:
ln(𝑥 + 1) + 𝑥 2 − 4 = 0
Conociendo que la solución de la ecuación se halla en [1, 2]
Solución:
Primero se debe verificar que realmente existe solo una solución en este
intervalo, esto se lo hace analizando si existe un cambio de signo en los extremos
del intervalo.
𝑓 (𝑥 ) = ln(𝑥 + 1) + 𝑥 2 − 4 = 0
𝑓 (𝑎) = 𝑓 (1) = ln(1 + 1) + 12 − 4 = −2.306853
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

𝑓 (𝑏) = 𝑓(2) = ln(2 + 1) + 22 − 4 = 1.098612

Desarrollamos el problema para las dos primeras iteraciones
𝑎+𝑏 1+2
𝑚= = = 1.5
2 2
𝑓 (𝑚) = 𝑓 (1.5) = ln(1.5 + 1) + 1.52 − 4 = −0.833709
Valor a
K a b m f(a) f(b) f(m)
cambiar
0 1.000000 2.00000 1.500000 -2.306853 1.098612 -0.833709 a
1 1.500000 2.00000 1.750000 -0.833709 1.098612 0.074101 b
2 1.500000 1.750000 1.625000 -0.833709 0.074101 -0.394294 a
3 1.625000 1.750000 1.687000 -0.394294 0.074101 -0.163732 a
4 1.687000 1.750000 1.718750 -0.163732 0.074101 -0.015726 a
5 1.718750 1.750000 1.734375 -0.045726 0.074101 0.013960 b
6 1.718750 1.734375 1.726562 -0.045726 0.013960 -0.015940 a
7 1.726562 1.734375 1.730469 -0.015940 0.013960 -0.001005 a
8 1.730469 1.734375 1.732422 -0.001005 0.013960 0.006474 b
9 1.730469 1,732422 1.731445 -0.001005 0.006474 0.002734 b
10 1.730469 1,731445 1.730957 -0.001005 0.002734 0.008640 b

Entonces tomamos como solución el ultimo valor medido después de realizar las
iteraciones

𝑥 = 1.730957
EJERCICIO 3.
Hallar la raíz próxima en [1, 2] de la siguiente ecuación por el método de
bisección.
𝑥 3 + 4𝑥 2 − 10 = 0
Primero se debe verificar que realmente existe solo una solución en este
intervalo, esto se lo hace analizando si existe un cambio de signo en los extremos
del intervalo.
𝑓 (𝑥 ) = 𝑥 3 + 4𝑥 2 − 10 = 0
𝑓(𝑎) = 𝑓 (1) = 13 + 4 ∗ 12 − 10 = −5
𝑓 (𝑏) = 𝑓 (2) = 23 + 4 ∗ 22 − 10 = 14
Desarrollamos el problema para las dos primeras iteraciones
𝑎+𝑏 1+2
𝑚= = = 1.5
2 2
𝑓(𝑚) = 𝑓(1.5) = 1.53 + 4 ∗ 1.52 − 10 = 2.375
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

Valor a
K a b m f(a) f(b) f(m)
cambiar
1 1.00000 2.00000 1.500000 -5.00000 14.00000 2.375000 b
2 1.00000 1.50000 1.250000 -5.00000 2.375000 -1.796875 a
3 1.25000 1.50000 1.375000 -1.796875 2.375000 0.162109 b
4 1.25000 1.37500 1.312500 -1.796875 0.162109 -0.848389 a
5 1.31250 1.37500 1.343750 -0.848389 0.162109 -0.350983 a
6 1.34375 1.37500 1.539375 -0.350983 0.162109 -0.096409 a
7 1.35940 1.37500 1.367200 -0.095998 0.162109 0.032563 b
8 1.35940 1.36720 1.363300 -0.095998 0.032563 -0.031841 a
9 1.36330 1.36720 1.365250 -0.031841 0.032563 0.000330 b
10 1.36330 1.36525 1.364275 -0.031841 0.000330 -0.015763 a

Iteracion 9
f(a)f(m)=1.051*10exp(-5)
Iteracion 10
f(a)f(m)=5.019*10exp(-4)
Podemos observar que en la novena iteración se tiene el valor de:
𝑓 (𝑚) = 0.000330 ≅ 0
Esta más cercano que f(m) en la décima iteración, por lo que tomamos como raíz
más próxima en el intervalo dado es el valor de “n” en el noveno intervalo.

𝑥 = 1.365250

EJERCICIO4.
Sean [a0, b0] [a1, b1]………… [an, bn] intervalos generados por el método de
bisección si Cn es el punto medio de cada intervalo demostrar que:
𝑎𝑛 + 𝑏𝑛
𝐶𝑛 =
2
a) lim 𝑎𝑛 = lim 𝑏𝑛 = lim 𝐶𝑛
𝑛→∞ 𝑛→∞ 𝑛→∞

𝑏𝑜−𝑎𝑜
b) 𝑏𝑛 − 𝑎𝑛 =
2𝑛

SOLUCION:
Para el método de la bisección se sabe que la raíz está dentro de un intervalo.
La raíz debe situarse dentro de:
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO
∆𝑥
± →→→ ∆𝑥 = 𝑋𝑏 − 𝑋𝑎
2
Entonces en nuestro caso el error es el siguiente:
𝑏𝑛 − 𝑎𝑛
𝑒𝑛 = | |
2
a) Aplicando limites tenemos lo siguiente:
𝑎𝑛 + 𝑏𝑛
lim 𝐶𝑛 = lim
𝑛→∞ 𝑛→∞ 2
Trabajando con el error:
𝑏𝑛 − 𝑎𝑛 𝑏𝑛 − 𝑎𝑛
𝑒𝑛 = | | →→ ±𝑒𝑛 =
2 2
Entonces despejando tenemos:
𝑏𝑛 = 𝑎𝑛 ± 2𝑒𝑛 𝑎𝑛 = 𝑏𝑛 ∓ 2𝑒𝑛
Reemplazando en los límites tenemos:
𝑎𝑛 + 𝑏𝑛 𝑎𝑛 + 𝑎𝑛 ± 2𝑒𝑛 2𝑎𝑛 ± 2𝑒𝑛
lim 𝐶𝑛 = lim = lim = lim
𝑛→∞ 𝑛→∞ 2 𝑛→∞ 2 𝑛→∞ 2
lim 𝐶𝑛 = lim (𝑎𝑛 ± 𝑒𝑛)
𝑛→∞ 𝑛→∞

𝑎𝑛 + 𝑏𝑛 𝑏𝑛 + 𝑏𝑛 ∓ 2𝑒𝑛 2𝑏𝑛 ∓ 2𝑒𝑛

lim 𝐶𝑛 = lim = lim = lim
𝑛→∞ 𝑛→∞ 2 𝑛→∞ 2 𝑛→∞ 2
lim 𝐶𝑛 = lim (𝑏𝑛 ∓ 𝑒𝑛)
𝑛→∞ 𝑛→∞

Como 𝑛 → ∞, el error "𝑒𝑛 → 0"

Entonces tenemos lo siguiente:

lim 𝐶𝑛 = lim (𝑎𝑛 ± 𝑒𝑛) = lim (𝑏𝑛 ∓ 𝑒𝑛)

𝑛→∞ 𝑛→∞ 𝑛→∞

lim 𝐶𝑛 = lim 𝑎𝑛 = lim 𝑏𝑛

𝑛→∞ 𝑛→∞ 𝑛→∞

b)
𝑏0 − 𝑎0
𝑏𝑛 − 𝑎𝑛 =
2𝑛
Si el error es igual a:
𝑏𝑛 − 𝑎𝑛 𝑏𝑛 − 𝑎𝑛
𝑒𝑛 = | | →→ ±𝑒 = ( )
2 2
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

Entonces:
𝑏0 − 𝑎0
𝑒0 = | |
2
Para n=1 tenemos:
𝑏1 − 𝑎1
𝑒1 = | |
2
Si:
𝑏0 𝑎0
𝑏1 = 𝑦 𝑎1 =
2 2
Entonces en e1:
𝑏0 − 𝑎0
𝑒1 = | |
22
Para n=2 tenemos:
𝑏2 − 𝑎2
𝑒2 = | |
2

Si:
𝑏1 𝑎1
𝑏2 = 𝑦 𝑎2 =
2 2
Entonces en e2:
𝑏1 − 𝑎1 𝑏0 − 𝑎0
𝑒2 = | 2
|=| |
2 23

Para n=3 tenemos:

𝑏3 − 𝑎3
𝑒3 = | |
2

Si:
𝑏2 𝑎2
𝑏3 = 𝑦 𝑎3 =
2 2
Entonces en e3:
𝑏2 − 𝑎2 𝑏0 − 𝑎0
𝑒3 = | 2
|=| |
2 24
Para n=n tenemos:
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO
𝑏𝑛 − 𝑎𝑛 𝑏0 − 𝑎0
𝑒𝑛 = | | = | 𝑛+1 |
2 2
Despejando tenemos:
𝑏𝑛 − 𝑎𝑛 𝑏0 − 𝑎0
=
2 2𝑛+1
𝑏0 − 𝑎0
𝑏𝑛 − 𝑎𝑛 =
2𝑛
METODO3. METODO DE LA INTERPOLACION LINEAL INVERSA
Sea f(x) una función continua es decir derivable en un intervalo [a, b] tales que
f(a) y f(b) son de signos opuestos.
Se aproximará a la raíz de la funcion mediante una recta L(x) o polinomio de
primer grado:
∞
𝑓 ′ (𝑎 )
𝑓 (𝑥 ) = ∑ (𝑥 − 𝑎)𝑖
𝑖!
𝑖=0

Para nuestro objetivo tomamos nota hasta i=1 así:

𝑓 (𝑥 ) ≈ 𝑓 (𝑎) + 𝑓′(𝑎)(𝑥 − 𝑎)
𝑓 (𝑥 ) ≈ 𝑓 (𝑎 ) + 𝑠 (𝑥 − 𝑎 ) ≡ 𝐿 (𝑥 )
Teorema del valor medio
El teorema del valor medio establece que si una función es continua en el
intervalo cerrado [a, b] y diferenciable en el intervalo abierto (a, b), entonces
existe un punto c contenido en el intervalo (a, b) tal que f'(c) es igual a la razón
de cambio promedio de la función en [a, b]. En otras palabras, la tangente de la
gráfica en algún punto es paralela a la recta secante que pasa por (a, f(a)) y (b,
f(b)).
Entonces de acuerdo a este teorema:
𝑓 (𝑏) − 𝑓(𝑎)
𝑠=
𝑏−𝑎
Entonces en la recta
𝑓 (𝑏 ) − 𝑓 (𝑎 )
𝐿 (𝑥 ) = 𝑓 (𝑎 ) + (𝑥 − 𝑎)
𝑏−𝑎
Si analizáramos directamente L(x) nos referiríamos directamente a una
interpolación lineal, pero se resuelve L(x) en lugar de f(x) por esta razón el
método se conoce como interpolación lineal inversa. Donde se debe hacer
L(c)=0 donde “c” es la raíz de la recta.
𝑓 (𝑏 ) − 𝑓 (𝑎 )
𝐿 (𝑐 ) = 𝑓 (𝑎 ) + (𝑐 − 𝑎)
𝑏−𝑎
ETN – 505 PROGRAMACION Y METODOS NUMERICOS DOCENTE: LIC. VIRGINIA ZOTA
AUX: UNIV. BORIS RODRIGO CUEVAS COYO

𝑓 (𝑏 ) − 𝑓 (𝑎 )
0 = 𝑓 (𝑎 ) + (𝑐 − 𝑎) →→→ 0 = 𝑓(𝑎)(𝑏 − 𝑎) + (𝑓 (𝑏) − 𝑓(𝑎))(𝑐 − 𝑎)
𝑏−𝑎
𝑎𝑓 (𝑏) − 𝑏𝑓(𝑎)
𝑐=
𝑓 (𝑏) − 𝑓(𝑎)
EJERCICO1.
𝑠𝑒𝑎 𝑓 (𝑥 ) = 𝑥 2 − 2𝑥 , 𝑒𝑛 [0, 3]
Solución:
Evaluamos la funcion en 0:
𝑓 (0) = 02 − 20 = −1
Evaluamos la funcion en 3:
𝑓 (3) = 32 − 23 = 1

Reemplazando en la ecuación de “c” tenemos:

0𝑓 (3) − 3𝑓(0) −3(−1) 3
𝑐= = =
𝑓 (3) − 𝑓(0) 1 − (−1) 2
Aplicamos nuevamente en la ecuación:
1𝑓 (3) − 3𝑓 (1) 1 − 3(−1) 4
𝑐= = = =2
𝑓 (3) − 𝑓 (1) 1 − (−1) 2
X=2
El ultimo valor obtenido es el valor de la ecuación, sucede siempre que el
intervalo considerado presente su centro próximo a la raíz.
EJERCICO2.
1
𝑠𝑒𝑎 𝑓(𝑥 ) = , 𝑒𝑛 [−1, 1]
𝑥
Solución:
−1𝑓 (1) − 1𝑓 (−1) −1 − 1(−1) 0
𝑐= = = =0
𝑓 (1) − 𝑓 (−1) 1 − (−1) 2
En cualquier intervalo de (-a, a) resultara c=0 ya que es una funcion discontinua
en el intervalo dado.