SEMANA N° 1

ECUACIONES NO LINEALES y TRASCENDENTES

Estas ecuaciones se caracterizan por que son definidas mediante funciones polinomiales,
logaritmos naturales, funciones exponenciales, funciones trigonométricas y que en
muchas ocasiones estas ecuaciones no se pueden resolver por los métodos tradicionales,
por lo tanto, debemos utilizar métodos que resuelvan numéricamente a estás ecuaciones,
obteniendo soluciones aproximadas reales de las soluciones verdaderas con cierto error.

Estas ecuaciones son generadas por funciones, que deben tener ciertas condiciones, que
se encuentran en Teorema de Bolzano. Veamos:
Teorema de Bolzano:
y Se observa f(a) y f(b) tienen signos opuestos, esto indica que y=f(x)
tiene intersección con el eje x, en r o en (r, 0). Así:
𝑦 = 𝑓 (𝑥)
𝑓 (𝑏) y=0
y = f(x) f(x) = 0
¿Por qué?

función Ecuación Asociada

(𝑟
, 0) que se va a estudiar
𝑥 4 +3 𝑥 3 − 𝑥 2+ 8 𝑥 −5=0
𝑓 (𝑎) a r b x
𝑥𝑙𝑛 ( 𝑥 ) − 𝑒 −𝑥 =0
Curso: Análisis Numérico Prof. Walter C. Cerna Molina
1
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 NOTA IMPORTANTE: Para estudiar los métodos numéricos para este tema, es indispensables en primer lugar conocer
la función y=f(x), para constatar las condiciones del Teorema de Bolzano y el que garantiza la intersección de la función
con el eje x en r o en (r, 0), y luego se elige los otros elementos que requiere cada método.

MÉTODO DE BISECCIÓN Aquí r es raíz de ecuación asociada f(x)=0, que en ocasiones no se

pueden hallar. El método consiste, teniendo el primer intervalo [a, b]
Eje y se obtiene el punto medio de la forma:
y = f(x) 𝑎 +𝑏
𝑆𝑖 𝑓 ( 𝑎 ) 𝑓 ( 𝑏 ) < 0 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑟 1=
2
f(b) Punto Medio o
(b, f(b))
raíz aproxima
Esto hace que el [a, b] se parta en 2 subintervalos iguales,
entonces debemos elegir el subintervalo en donde esta r

a r 𝑟1 Eje x
b a b
r
f(a) Eje x
(a, f(a)) Subintervalo izquierdo
elegido [a,
𝑎 +𝑏
: 𝑆𝑖 𝑓 ( 𝑎 ) 𝑓 ( 𝑟 1 ) < 0 𝑒𝑛𝑡𝑜𝑛𝑐𝑒𝑠 𝑟 2=
2
𝑟2 r 𝑏=𝑟 1 Punto Medio o
a Eje x raíz
aproximada
Subintervalo elegido
Curso: Cálculo Numérico Prof. Walter C. Cerna Molina
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 {𝑟 1 , 𝑟 2 ,𝑟 3 , …,𝑟 𝑛 } Aproximen cada
vez mas a:
r

EN GENERAL: Este método de bisección consiste en subdividir los intervalos, teniendo en cuenta el
punto medio del intervalo correspondiente: Luego se analizan las dos posibilidades que se da, en la
elección de los subintervalos de la izquierda o de la derecha:
1. .
Luego hallamos el valor de su punto medio.
2. 𝐸𝑠 𝑠𝑢𝑏𝑖𝑛𝑡𝑒𝑟𝑣𝑎𝑙𝑜 𝑑𝑒𝑟𝑒𝑐h𝑜 𝑒𝑙𝑒𝑔𝑖𝑑𝑜 :𝑆𝑎𝑡𝑖𝑠𝑓𝑎𝑐𝑒 𝑓 ( 𝑟 𝑛 ) 𝑓 ( 𝑏 ) < 0 , 𝑎𝑞𝑢í 𝑐𝑎𝑚𝑏𝑖𝑎𝑎=𝑟 𝑛 .
Luego hallamos el valor de su punto medio.
Ejemplo 1: Resolver la
• Definir la función: , para graficarla es necesario usar software:

Como se observa en el gráfico de la función, tiene dos

intersecciones con el Eje x. Debemos de elegir uno de
ellos, como por ejemplo el que esta en el intervalo [0,
1]

Curso: Análisis Numérico Prof. Walter C. Cerna Molina

3
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 • Hallar los valores de la sucesión: { 𝑟 1 , 𝑟 2 ,𝑟 3 , … ,𝑟 𝑛 } mediante iteraciones

Primera iteración: El primer intervalo [a, b]=[0, 1]. Entonces se determina el valor del punto medio:

𝑎 +𝑏 0+1
𝑟 1= = = 0.5
2 2
Teniendo el valor del punto medio, hace que se subdivida [0, 1] en
dos subintervalos iguales, como se muestra en el gráfico, esto
permite elegir el subintervalo derecho [0.5, 1] que será el intervalo
para siguiente iteración.

Segunda iteración: El segundo intervalo [a, b]=[0.5, 1]. Entonces se determina el valor del punto medio:

𝑎 +𝑏 𝑟 1 +𝑏 0.5+ 1
𝑟 2= = = =0.75
2 2 2
Teniendo el valor del punto medio, hace que se subdivida [0.5, 1] en
dos subintervalos iguales, como se muestra en el gráfico, esto
permite elegir el subintervalo izquierdo [0.5, 0.75] que será el
intervalo para siguiente iteración.

Curso: Análisis Numérico Prof. Walter C. Cerna Molina

4
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 Tercera iteración: El tercer intervalo [a, b]=[0.5, 0.75]. Entonces se determina el valor del punto medio:

𝑎+𝑏 𝑎+𝑟 2 0.5+ 0.75

𝑟3= = = =0.625
2 2 2
Teniendo el valor del punto medio, hace que se subdivida [0.5, 0.75]
en dos subintervalos iguales, como se muestra en el gráfico, esto
permite elegir el subintervalo izquierdo [0.5, 0.625] que será el
intervalo para siguiente iteración.

Cuarta iteración: El cuarto intervalo [a, b]=[0.5, 0.625]. Entonces se determina el valor del punto medio:

𝑎+ 𝑏 𝑎+𝑟 3 0.5+0.625
𝑟4= = = =0.562500
2 2 2
Teniendo el valor del punto medio, hace que se subdivida [0.5, 0.625]
en dos subintervalos iguales, como se muestra en el gráfico, esto
permite elegir el subintervalo derecho [0.562500, 0.625] que será el
intervalo para siguiente iteración.
Otra forma de hallar el intervalo 𝑓 ( 𝑎 ) 𝑓 ( 𝑟 4 ) 𝑓 ( 𝑏)
para
la siguiente iteración, es evaluando
𝑓 ( 0.5 ) 𝑓 ( 0.5625 ) 𝑓 (0.625)
la función Aquí − − +¿
Por lo tanto: [0.5625, 0.625]

Curso: Análisis Numérico Prof. Walter C. Cerna Molina

5
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 Quinta iteración: El quinto intervalo [a, b]=[0.5625, 0.625]. Entonces se determina el valor del punto
medio:

𝑎+𝑏 𝑟 4 +𝑏 0.5625+ 0.625

𝑟5= = = =0.593750
2 2 2

• Presentamos Tabla de resultados

Estos resultados deben ser comparados con

lo resultados adquirido con software

Curso: Cálculo Numérico Prof. Walter C. Cerna Molina

6
IV-Ciclo Ingeniería de Sistemas Mg. en Ingeniería de Sistemas
 Comparar resultados con uso de software
a) Hacer un programa en MatLab. 𝑒𝑐𝑢𝑎𝑐𝑖 ó 𝑛: 𝑥 4 +3 𝑥 3 − 𝑥 2+8 𝑥 − 5=0
function biseccion(f,a,b,n)
Es necesario definir las función asociada
disp(' MÉTODO DE BISECCIÓN');
disp('================================='); 𝑓 ( 𝑥)= 𝑥 4 +3 𝑥 3 − 𝑥2 +8 𝑥 − 5
disp(' it a b r'); Para conocer los valores de a y b,
disp('================================='); se debe graficar la función:
F=inline(f);
F(a)*F(b)<0; i=1; Software MathCad
while i<=n
r=(a+b)/2;
fprintf('%3.0f %10.6f %10.6f %10.6f\n',i,a,b,r);
i=i+1;
if F(a)*F(r)<0
b=r;
else
a=r;
end
pause(0.25)
end Elijamos el intervalo [a, b]=[0, 1] donde a=0
fprintf('La Raiz = %8.6f\n',r) y b=1.

Prof. Walter C. Cerna Molina

Mg. en Ingeniería de Sistemas
Prof. Walter C. Cerna Molina
8
Mg. en Ingeniería de Sistemas
Ahora nos toca comparar resultados:
Los adquiridos con software MatLab

De forma operativa

OK

¡ooh!

OK

Prof. Walter C. Cerna Molina

9
Mg. en Ingeniería de Sistemas
Resolviendo con software las ecuaciones no lineales en forma rápida:
En MatLab:

En MatCad:

En Geogebra: Vista Algebraica Vista CAS o Cálculo Simbólico

Prof. Walter C. Cerna Molina

10
Mg. en Ingeniería de Sistemas
Resolviendo con software MatLab utilizando el programa hecho:
Ejecutando el programa debemos obtener el resultado con una aproximación de 6 cifras decimales, como se desarrollo
con software en forma rápida, usando instrucciones cuyo resultado fue que raíz = 0.579803. Anteriormente se ejecuto el
programa hasta 5 iteraciones o sea n=5, obteniéndose que raíz = 0.593750.
OK

Nota: Se observa que se ha necesitado veinte iteraciones,

para llegar al resultado obtenido con los tres software
usando instrucciones rápidas. Veamos
Prof. Walter C. Cerna Molina
11
Mg. en Ingeniería de Sistemas
 Practica N° 01
1. Graficar las siguientes funciones utilizando software (MatLab o MathCad o Geogebra) y así obtener las
condiciones que se necesitan para utilizar el Método de Bisección.
2. Resolver en forma operativa las ecuaciones asociadas de las funciones dadas en la pregunta 1,
obteniendo así una tabla con cinco iteraciones con una aproximación de 6 cifras decimales..
3. Luego ejecute el programa hecho en MatLab, para resolver las ecuaciones y señale a partir de que
iteraciones se repiten los resultados.
4. Resolver las ecuaciones usando con instrucciones rápidas con los software MatLab, MathCad y
Geogebra. Luego compare resultados, con lo hecho en la pregunta 3.
5. Las ecuaciones de lanzamiento de un Misil balístico donde t esta dado en minutos y el alcance en
kilómetros son:

a) Determinar el tiempo transcurrido hasta el impacto con una precisión de 6 decimales.

b) Determine el alcance del disparo del misil con una precisión de 6 decimales.
c) Determine la función y=f(x) y grafique su trayectoria del misil.
d) Calcule la altura máxima del misil.

Prof. Walter C. Cerna Molina

12
Mg. en Ingeniería de Sistemas