CLASE: Nro.
8-9
UC 1: USO DE SOFTWARE ESPECIALIZADO PARA CÁLCULO
NUMÉRICO, TEORIA DEL ERROR Y SOLUCIÓN DE ECUACIONES NO
LINEALES
MÉTODOS NUMÉRICOS
Solución de ecuaciones no lineales
Problemas
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
CONTENIDO
Título Programación para la solución de ecuaciones no lineales
Duración 2 horas
Información general Comandos y algoritmos para encontrar la solución de ecuaciones
no lineales
Objetivo Conocer y aplicar los comandos básicos y los algoritmos de
programación en Octave para la solución de ecuaciones no lineales
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
MÉTODOS DE SOLUCIÓN PARA ECUACIONES NO LINEALES
PROBLEMAS
Para resolver cualquier ecuación, recordemos que, sin importar la forma que nos presenten el modelo
matemático, tendremos que transformar a la forma:
𝑓(𝑥) = 0, (1)
Toda vez que, los métodos de solución se basan en esta estructura algebraica.
1) Encuentre todas las soluciones a la siguiente ecuación:
2 −6𝑥+9)
2𝑙𝑜𝑔𝑥 (𝑥 = 32𝑙𝑜𝑔𝑥 √x−1
Es importante aclarar que esta ecuación tiene 2 soluciones exactas, que son: 2,4
Sin embargo, la resolveremos con los métodos estudiados, para analizar y comparar la convergencia
y efectividad de cada método.
Ahora bien, para empezar a usar los métodos estudiados, empezaremos como se mencionó
anteriormente, igualando la ecuación a cero:
2 −6𝑥+9)
2𝑙𝑜𝑔𝑥 (𝑥 − 32𝑙𝑜𝑔𝑥 √x−1 = 0
Adicionalmente, es importante recalcar que el dominio de la ecuación debe conocerse de antemano.
De tal manera que, para esta ecuación, tenemos que:
𝑥 ∈ (0; +∞) − {1}
Resolución con el método de la Bisección
En nuestra función, creada en Octave o Matlab ingresamos los datos:
>> raiz=met_biseccion
Ingrese la función a resolver : 2^(log(x^2-6*x+9)/log(x))-3^(2*log(sqrt(x))/log(x)-1)
f = 2^(log(x^2-6*x+9)/log(x))-3^(2*log(sqrt(x))/log(x)-1)
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
De acuerdo a la gráfica y al dominio de la ecuación, para encontrar la primera raíz, el intervalo a
considerarse podría ser: [1.2; 2.5]
Ingrese a: 1.2
Ingrese b: 2.5
Ingrese el número de iteraciones: 100
Ingrese la tolerancia de error: 0.00001
an cn bn
1.200000 1.850000 2.500000
1.850000 2.175000 2.500000
1.850000 2.012500 2.175000
1.850000 1.931250 2.012500
1.931250 1.971875 2.012500
1.971875 1.992188 2.012500
1.992188 2.002344 2.012500
1.992188 1.997266 2.002344
1.997266 1.999805 2.002344
1.999805 2.001074 2.002344
1.999805 2.000439 2.001074
1.999805 2.000122 2.000439
1.999805 1.999963 2.000122
1.999963 2.000043 2.000122
1.999963 2.000003 2.000043
f(c)=-0.00000610
niter = 15
raiz = 2.000003051757813
Observamos que la solución se alcanza en 15 iteraciones. Para la otra raíz, podemos considerar el
intervalo [2.2;4.5]
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
Resolución con el método de la Secante
Usamos el programa correspondiente y considerando el intervalo anterior, obtenemos:
>> met_secante
Ingrese la ecuación a resolver: 2^(log(x^2-6*x+9)/log(x))-3^(2*log(sqrt(x))/log(x)-1)
Ingrese el primer valor: 1.2
Ingrese el segundo valor: 2.5
Ingrese la tolerancia: 0.00001
n xo x1 x2 error
0 1.2000 2.5000 ---- ----
0 1.2000 2.5000 2.4903 0.6408
1 2.5000 2.4903 1.7817 0.6065
2 2.4903 1.7817 2.1262 0.2196
3 1.7817 2.1262 2.0346 0.0665
4 2.1262 2.0346 1.9949 0.0103
5 2.0346 1.9949 2.0002 0.0004
6 1.9949 2.0002 2.0000 0.0000
f(x2)=-0.0000
niter = 7
ans = 2.000001348149192
Observamos que la solución se alcanza en 7 iteraciones.
Resolución con el método de Newton
Ingresamos los datos en el programa correspondiente:
>> raiz=metod_newton(100,0.00001)
Ingrese la función a resolver: 2^(log(x^2-6*x+9)/log(x))-3^(2*log(sqrt(x))/log(x)-1)
Ingrese el valor en donde empezarán las iteraciones: 1.2
.
.
.
.
warning: passing floating-point values to sym is dangerous, see "help sym"
warning: called from
double_to_sym_heuristic at line 50 column 7
sym at line 379 column 13
subs at line 226 column 9
metod_newton at line 18 column 7
Se terminó las iteraciones
iteraciones = 10
raiz = 1.999999999992994
como vemos se alcanza la convergencia en 10 iteraciones, pero adicionalmente se generan varios
mensajes de advertencia debido a la complejidad de los cálculos que debe realizar Octave.
En conclusión, el método de la secante resultó ser el más eficiente. Sin embargo, en general el método
de Newton es el más eficiente, la desventaja como observamos es la obtención de la derivada dentro
del mismo script.
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
3
2) Encuentre una aproximación a √2 y √3 mediante el método de la bisección.
Recordemos que para encontrar estas aproximaciones debemos trabajar con ecuaciones de la
forma 𝑓(𝑥) = 0
• Entonces para calcular √2, vamos a partir de la ecuación
3
3
𝑥 = √2
Luego, sabemos que
𝑥3 = 2
y finalmente
𝑥3 − 2 = 0
Por lo que, la función a resolver será
𝑓(𝑥) = 𝑥 3 − 2
Así, con el método de Newton, se consigue en 6 iteraciones y con una tolerancia de 10−7 :
3
√2 = 1.259921049894873.
• Para calcular √3, vamos a partir de la ecuación
𝑥 = √3
y siguiendo la misma idea del ejercicio anterior, la función a resolver es
𝑓(𝑥) = 𝑥 2 − 3
Usando también el método de Newton, obtenemos
√3 = 1.732050807568877
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
3) Calcule las raíces reales de 𝑓(𝑥) = −12 − 21𝑥 + 18x 2 − 2.75𝑥 3
4) Emplee el método de la secante con la ecuación de la circunferencia (𝑥 + 1)2 + (𝑦– 2)2 = 16, a
fin de encontrar una raíz real positiva. Haga que el valor inicial sea 𝑥𝑖−1 = 3 y 𝑥𝑖 = 0.5.
Aproxímese a la solución del primer y cuarto cuadrantes. Cuando resuelva para 𝑓(𝑥) en el cuarto
cuadrante, asegúrese de tomar el valor negativo de la raíz cuadrada. ¿Por qué diverge la solución?
En este caso despejamos y de la ecuación dada
𝑦 = ±√16 − (𝑥 + 1)2 + 2
Ahora bien, las gráficas correspondientes a ambas funciones son
Y podemos observar que la única función que tiene raíces es 𝑔 = ±√16 − (𝑥 + 1)2 + 2. Sin
embargo, vemos que diverge la solución si ingresamos 𝑥0 = 3 y 𝑥1 = 0.5, y se debe a que el orden
de las x es incorrecto. La gráfica lo demuestra:
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
5)
Graficando la función, se puede observar que existe un cero, cerca del ángulo de 25°. Entonces,
usando el método de Newton, con valor inicial de 25° (en radianes), obtenemos:
Que es el ángulo buscado en grados sexagesimales.
6) *El producto de 2 números reales positivos es 5. Mientras que al sumar el cubo del primero más
del cuadrado del segundo se obtiene 40. Encuentre estos dos números. Incluir comprobación.
El problema se modela mediante el siguiente sistema de ecuaciones
𝑥𝑦 = 5(1)
{
𝑥 3 + 𝑦 2 = 40(2)
Despejando 𝑦 de la primera ecuación, consecuentemente reemplazando en la segunda e
igualando a cero obtenemos
25
𝑥3 + − 40 = 0
𝑥2
que es la ecuación que debemos resolver.
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
7) Suponga el lector que está diseñando un tanque esférico de almacenamiento de agua para un
poblado pequeño de un país en desarrollo. El volumen del líquido que puede contener se calcula
con
donde V = volumen [pie3], h = profundidad del agua en el tanque [pies], y R = radio del tanque
[pies].
Si R = 3 m, ¿a qué profundidad debe llenarse el tanque de modo que contenga 30 m3? Haga tres
iteraciones del método de Newton-Raphson para determinar la respuesta. Encuentre el error
relativo aproximado después de cada iteración. Observe que el valor inicial de R convergerá
siempre.
8) Un kilogramo mol de CO está contenido en un recipiente a 𝑇 = 215 𝐾 y 𝑃 = 70 𝑏𝑎𝑟. Calcule el
volumen del gas utilizando la ecuación de estado de van der Waals para un gas no ideal, dada
por
𝑎
(𝑃 + ) (𝑣 − 𝑏) = 𝑅𝑇
𝑣2
𝑏𝑎𝑟 𝑚3 𝑏𝑎𝑟 𝑚6 𝑚3
donde 𝑅 = 0.08314 , 𝑎 = 1463 (𝑘𝑔 𝑚𝑜𝑙)2 y 𝑏 = 0,0394 𝑘𝑔 𝑚𝑜𝑙
𝑘𝑔 𝑚𝑜𝑙 𝐾
𝑚3
Determine el volumen específico 𝑣 (𝑘𝑔 𝑚𝑜𝑙) y compare el resultado con el volumen calculado
por la ecuación de los gases ideales, 𝑃𝑣 = 𝑅𝑇. Puede elegir cualquier método de resolución.
Mgtr. Jacobo Vásquez
CLASE: Nro. 8-9
Más problemas….
• La velocidad v de un paracaidista que cae está dada por
donde g = 9.8 m/s2. Para un paracaidista con coeficiente de arrastre de c = 15 kg/s, calcule la
masa m de modo que la velocidad sea v = 35 m/s en t = 9s. Utilice el método de la secante
para determinar m a un nivel de = 0.1%.
REFERENCIAS BIBLIOGRAFICAS
1. Chapra Steven, Canale Raymond, Métodos numéricos para ingenieros, McGraw-Hill, Quinta
edición, 2007.
2. A. Quarteroni, F. Saleri, Cálculo Científico con Matlab y Octave. Springer-Verlag Italia, milano
2006.
3. A. Delgado, J. Nieto, A. Robles, O. Sánchez; Métodos numéricos básicos con Octave, Editorial
Técnica AVICAM, Granada, septiembre 2018
4. [Link]
Mgtr. Jacobo Vásquez
