Ecuaciones no lineales

𝑓(𝑥) = 0

Método Gráfico
Usando la gráfica podemos seleccionar un intervalo donde la ecuación presenta una posible solución o
raíz.

Ejemplo 1

𝑥 3 − 3𝑥 2 + 2 = 0
¿Quién es la función? 𝑓(𝑥) = 0

𝑓(𝑥) = 𝑥 3 − 3𝑥 2 + 2

Esta ecuación tiene tres raíces, dos positivas y una negativa

Un intervalo para la raíz negativa [−1,0]

Un intervalo para la raíz negativa [−0.8, −0.6]

Un intervalo para la segunda raíz positiva [2,3]

Ejemplo 2

𝑥 + 2 = tan(𝑥)
¿Quién es la función? 𝑓(𝑥) = 0

𝑥 + 2 − tan(𝑥) = 0
𝑓(𝑥) = 𝑥 + 2 − tan(𝑥)
Veamos la función tangente tiene múltiples intersecciones con el eje x por ende tiene múltiples raíces
vamos a seleccionar un intervalo adecuado para la primera raíz negativa.

Generemos iteraciones mediante el calculo “al ojo” del valor aproximado de la solución

Iteración 1 (Brandon) [-4.4,-4.2] Iteración 2 (Isaac) [-4.4,-4.3] Iteración 3 (Brandon) [-4.32,-4.3]

Solución: 𝑥 ≈ −4.30 𝑥 ≈ −4.31 𝑥 ≈ −4.303

Métodos cerrados
La solución debe de pertenecer al intervalo y generar un cambio de signo

Método muy confiable pero lentos en converger (converger significa aproximar de forma exacta el
valor)

Bisección
 Sea f una función continua, definida en un intervalo [𝑥𝑖 , 𝑥𝑢 ] donde 𝑓(𝑥𝑖 ) ∙ 𝑓(𝑥𝑢 ) < 0,

∃𝑥 ∈ [𝑥𝑖 , 𝑥𝑢 ]: 𝑓(𝑥) = 0
El método de bisección consiste en
𝑥𝑖 + 𝑥𝑢
𝑥𝑟 =
2
Ejemplos

Probar que la solución vive en el intervalo dado 𝑥𝑖 = 0 𝑥𝑢 = 1 𝑓(𝑥𝑖 ) ⋅ 𝑓(𝑥𝑢 ) < 0 , determine quién
es la función 𝑓(𝑥) = 0, para ello debe dejar la ecuación igualada a cero 𝑥 2 − 3𝑥 + 2 − 𝑒 𝑥 = 0,
entonces la función 𝑓(𝑥) = 𝑥 2 − 3𝑥 + 2 − 𝑒 𝑥 , ahora bien determine las imágenes de los valores
iniciales y finales 𝑓(0) = 1, 𝑓(1) = −𝑒 (edite la función en la calculadora y evalué el valor de 𝑥). Vea
que 1 ⋅ −𝑒 < 0 dentro del intervalo hay una solución.

𝑥𝑖 𝑥𝑟 𝑥𝑢 Error
0 𝑥𝑖 + 𝑥𝑢 0 + 1 1 1 -
= =
2 2 2
𝑓(0) = 1 1 𝑓(1) = −𝑒
𝑓 ( ) = −0.8987
2
1 1
𝑓(0) ⋅ 𝑓 ( ) < 0 𝑓 ( ) ⋅ 𝑓(1) > 0
2 2
0 𝑥𝑖 + 𝑥𝑢 0 + 1/2 1 1 1 1
= = −
2 2 4 2 𝐸𝑟 = |2 4| ∗ 100%
1/4
𝑓(0) = 1 1 1
𝑓 ( ) = 0.0284 𝑓 ( ) = −0.8987 𝐸𝑟 = 100%
4 2
1 1 1
𝑓(0) ⋅ 𝑓 ( ) > 0 𝑓( )⋅𝑓( ) < 0
4 4 2
1 3 1 3 1
𝑥𝑟 = −
4 8 2 𝐸𝑟 = |8 4| ∗ 100%
1 𝑓(𝑥) = −0.4393 1 3/8
𝑓 ( ) = 0.0284 𝑓 ( ) = −0.8987
4 2 𝐸𝑟 = 33.3%
3 1 3 1
𝑓( )⋅𝑓( ) < 0 𝑓( )⋅𝑓( ) > 0
8 4 8 2
1 5 3 𝐸𝑟 = 20%
𝑥𝑟 = ≈ 0.3125
4 16 8
 9.8 ⋅ 82 𝑐
36 = (1 − 𝑒 −82∗4 )
𝑐
803.6 2𝑐
36 = (1 − 𝑒 −41 )
𝑐
Probar que la solución vive en el intervalo dado 𝑐𝑖 = 3 𝑐𝑢 = 5 𝑓(𝑐𝑖 ) ⋅ 𝑓(𝑐𝑢 ) < 0 , determine quién es
2𝑐
803.6
la función 𝑓(𝑐) = 0, para ello debe dejar la ecuación igualada a cero (1 − 𝑒 −41 ) − 36 = 0,
𝑐
2𝑐
803.6
entonces la función 𝑓(𝑥) = (1 − 𝑒 −41) − 36, ahora bien determine las imágenes de los valores
𝑐
iniciales y finales 𝑓(3) = 0.4666, 𝑓(5) = −1.2144 (edite la función en la calculadora y evalué el valor
de 𝑐). Vea que 0.4666 ⋅ −1.2144 < 0 dentro del intervalo hay una solución.

El coeficiente de resistencia sería 3.5626 aproximadamente

 𝑏−𝑎
< 10−𝑛
2𝑘
5 − −2
< 10−4
2𝑘
7
< 2𝑘
10−4
7
log 2 ( −4 ) < 𝑘
10
16.09 < 𝑘
𝑘 = 17 se deben hacer 17 iteraciones

Método de la Regla Falsa

Método cerrado 𝑓(𝑥𝑖 ) ⋅ 𝑓(𝑥𝑢 ) < 0 la raíz este dentro del intervalo

Gráficamente traza rectas y el corte con el eje x es la nueva iteración

𝑓(𝑥𝑢 )(𝑥𝑖 −𝑥𝑢 )
Analíticamente 𝑥𝑟 = 𝑥𝑢 −
𝑓(𝑥𝑖 )−𝑓(𝑥𝑢 )

Ejemplo 1

𝑓(𝑥) = 0

𝑥2𝑒 𝑥 − 1 = 0
𝑓(𝑥) = 𝑥 2 𝑒 𝑥 − 1
Se trabajará en forma gráfica para determinar una posible aproximación se realiza 4 iteraciones

para la primera iteración se traza

una recta entre los valores 𝑥1 = −1
y 𝑥𝑢 = 1 el corte que genera entre
esa recta y el eje 𝑥 es el nuevo valor
aproximado 𝑥𝑟 , se busca la imagen
correspondiente en la curva sería el
valor de 𝐷.

para ello se traza una primera recta entre los

Para la segunda iteración se traza la
recta entre los valores 𝐷 y 𝑥𝑢 = 1 el
corte que genera entre esa recta y el
eje 𝑥 es el nuevo valor aproximado 𝑥𝑟2 ,
se busca la imagen correspondiente en
la curva sería el valor de 𝐹.

Para la tercera iteración se traza la

recta entre los valores 𝐹 y 𝑥𝑢 = 1
el corte que genera entre esa
recta y el eje 𝑥 es el nuevo valor
aproximado 𝑥𝑟3 , se busca la
imagen correspondiente en la
curva sería el valor de 𝐻.

Para la cuarta iteración se

traza la recta entre los valores
𝐻 y 𝑥𝑢 = 1 el corte que genera
entre esa recta y el eje 𝑥 es el
nuevo valor aproximado 𝑥𝑟4
vea que cada iteración se va
acercando al valor 𝐶, que es la
raíz de la función.
 Ahora bien, resolver en forma analítica (empleando la fórmula)

𝑓(𝑥𝑢 )(𝑥𝑖 − 𝑥𝑢 )
𝑥𝑟 = 𝑥𝑢 −
𝑓(𝑥𝑖 ) − 𝑓(𝑥𝑢 )

𝑓(𝑥) = 𝑥 2 𝑒 𝑥 − 1

# ite 𝑥𝑖 𝑥𝑟 𝑥𝑢 𝑓(𝑥𝑖 ) 𝑓(𝑥𝑟 ) 𝑓(𝑥𝑢 ) Error %

1 -1 -0.935030 2 -0.632120 -0.656779 28.556224 -
2 -0.935030 -0.869043 2 -0.656779 -0.683289 28.556224 7.59%
3 -0.869043 -0.801997 2 -0.683289 -0.711568 28.556224 8.35%
4 -0.801997 -0.733873 2 -0.711568 28.556224 9.28%

Dominio máximo es [0, +∞[

Por gráfica se selecciona un intervalo adecuado para determinar una de las raíces de solución [1,2]

# ite 𝑥𝑖 𝑥𝑟 𝑥𝑢 𝑓(𝑥𝑖 ) 𝑓(𝑥𝑟 ) 𝑓(𝑥𝑢 ) Error <0.001

1 1 1.6941 2 -2.1585 -0.2334 0.9510 -
2 1.6941 1.7543 2 -0.2334 -0.0153 0.9510 0.03
3 1.7543 1.7581 2 -0.0153 -0.0012 0.9510 0.002
4 1.7581 1.7584 2 -0.0012 0.9510 0.0001

𝑓(𝑥𝑢 )(𝑥𝑖 − 𝑥𝑢 )
𝑓(𝑥) = 𝑥 2 |𝑠𝑒𝑛(√𝑥)| − 3 𝑥𝑟 = 𝑥𝑢 −
𝑓(𝑥𝑖 ) − 𝑓(𝑥𝑢 )
Métodos Abiertos
* Son rápidos.

* Tienen probabilidad de no converger.

* No necesariamente el intervalo debe de estar la raíz.

Punto fijo
Fórmula

𝑔(𝑥) = 𝑥
𝑥𝑘 : valor inicial

𝑔(𝑥𝑘 ) = 𝑥𝑘+1
𝑔(𝑥𝑘+1 ) = 𝑥𝑘+2
Ejemplos

# iter 𝑥 Error Función 𝑔(𝑥) = 𝑥 𝑒 −𝑥 = 𝑥

0 1 -
𝑔1 (𝑥) = 𝑒 −𝑥
1 𝑔1 (1) = 𝑒 −1 172%
0.367879 Derivar la función
2 𝑔1 (0.367879) 47%
= 𝑒 −0.367879 𝑔1′ (𝑥) = 𝑒 −𝑥 ⋅ −1 = −𝑒 −𝑥
0.692200
|𝑔1′ (0.5)| = 0.606530
3 𝑔1 (0.6922) 38%
= 𝑒 −0.6922 𝑔2 (𝑥) = −𝑙𝑛𝑥
0.500473
Derivada la función
𝑙𝑛𝑒 −𝑥 = 𝑙𝑛𝑥
4 0.606243 17%
1
5 0.545396 11% 𝑔2′ (𝑥) = − 𝑥 = −𝑙𝑛𝑥
𝑥
6 0.579612 6%
7 0.560115 3% |𝑔2′ (0.5)| = 2
8 0.571143 2%
9 0.564879 1%
10 0.568428 0.62%
Sacar las derivadas
2
1
𝑔1′ (𝑥) = (𝑥 + 1)−3 ⋅ 1 |𝑔1′ (0.2)| = 0.295182 sabemos que tiene un punto fijo y converge
3
1
1 1 −2 1
𝑔2′ (𝑥) = (1 + ) ⋅− |𝑔2′ (0.2)| = 5.1031 no converge no tiene punto fijo
2 𝑥 𝑥2

0.0025 9.35
𝑥 = 1.14 − 2 log ( + 𝑥)
0.1 3𝑥104
0.0025 9.35
𝑔(𝑥) = 1.14 − 2 log ( + 𝑥)
0.1 3𝑥104
|𝑔′ (0.3)| = |0 − 6.23𝑥10−4 | = 0.000623 < 1

# iter 0.0025 9.35

𝑔(𝑥) = 1.14 − 2 log ( + 𝑥)
0.1 3𝑥104

0 4
1 4.301851
2 4.298743
3 4.298775
4 4.298775

1
Recuerde que 𝑥 = se tiene 𝑥 = 4.298775
√𝑓
 1
= 4.298775
√𝑓
1
= √𝑓
4.298775
2
1 2
( ) = (√𝑓)
4.298775
𝑓 = 0.054114
∴El factor de fricción debe de ser de 0.054114.

Nota hacer en Excel este método

Método de Newton Raphson

𝑓(𝑥𝑖 )
𝑥𝑖+1 = 𝑥𝑖 −
𝑓′(𝑥𝑖 )
 Consiste en elegir un punto inicial cualquiera 𝑥0 .

 Calcular 𝑓(𝑥0 )

 Traza una recta tangente en el punto (𝑥0 , 𝑓(𝑥0 )).

 Tomar como segunda aproximación la coordenada en x del punto de intersección de la recta

tangente con el eje de las abscisas.

 Repetir el proceso.