Benemérita Universidad Autónoma de Puebla

Colegio de ingeniería química

ANÁLISIS NUMÉRICO Y PROGRAMACIÓN

Profesor: José Antonio Rivera Márquez

Alumna: Araceli Felipe Martínez

Matricula: 202145264 NRC: 43621

Examen ordinario.
Ejercicio 1
En la ingeniería ambiental una especialidad de la ingeniería civil la ecuación
siguiente se emplea para calcular el nivel de oxigeno c(mg/L) en un rio agua
abajo de la descarga de drenaje.
𝑐 = 10 − 20(𝑒 −0.15𝑥 − 𝑒 −0.5𝑥 )
Donde x es la distancia agua abajo en kilómetros. Calcule la distancia aguas
abajo a la cual el oxigeno se encuentra al mínimo (recomendación: esta dentro
de los 3 y 3.6 km de la descarga) ¿Cuál es la concentración en dicha ubicación?
Respuesta
De la ecuación procedemos a calcular la derivada de C
𝑐 = 10 − 20(𝑒 −0.15𝑥 − 𝑒 −0.5𝑥 )
𝑐 ′ = 0 − 20(−0.15𝑒 −0.15𝑥 + 0.5𝑒 −0.5𝑥 ) = 0
−0.15𝑒 −0.15𝑥 + 0.5𝑒 −0.5𝑥 = 0
0.15𝑒 −0.15𝑥 = 0.5𝑒 −0.5𝑥
𝑒 −0.5𝑥 = 0.3𝑒 −0.15𝑥
𝑒 −0.35𝑥 = 0.3
Aplicamos logaritmo natural (Ln) para eliminar a e
𝑙𝑛−0.35𝑥 = 𝑙𝑛0.3
0.35𝑥(𝑙𝑛𝑒) = 𝑙𝑛0.3
𝑙𝑛0.3
𝑥=
−0.35
𝑥 = 3.44𝐾𝑚
Se calcula el valor del nivel de oxigeno mínimo sustituyendo el valor de x en la
ecuación original.
𝑐 = 10 − 20(𝑒 −0.15𝑥 − 𝑒 −0.5𝑥 )
𝑐 = 10 − 20(𝑒 −0.15(3.44) − 𝑒 −0.5(3.44 )
𝑚𝑔
𝑐 = 1.64
𝐿
Ejercicio 2
Determine las dos raíces positivas de las siguientes ecuaciones simultáneas, por
el método de Newton Raphson. Usar el método gráfico para hacer una propuesta
razonable de los valores iniciales de las dos raíces solicitadas.
𝑦 = −𝑥 2 + 𝑥 + 0.75
𝑦 + 5𝑥𝑦 = 2
Respuesta
𝐹(𝑥, 𝑦) = −𝑥 2 + 𝑥 + 0.75 − 𝑦 = 0
𝐺(𝑥, 𝑦) = 𝑦 + 5𝑥𝑦 − 2 = 0
Resolvemos obteniendo las derivadas parciales para construir una matriz
jacobiana.
𝜕𝐹 𝜕𝐹
= −2𝑥 + 1, = −1
𝜕𝑥 𝜕𝑦
𝜕𝐺 𝜕𝐺
= 5𝑦, = 1 + 5𝑥
𝜕𝑥 𝜕𝑦
Entonces nos queda una matriz de tal manera que:
−2𝑥 + 1 −1
𝐽 =[ ]
5𝑦 1 + 5𝑥
Primer punto de intersección: (𝑥0, 𝑦0) ≈ (0.2, 0.6)
Segundo punto de intersección: (𝑥1, 𝑦1) ≈ (0.7, 0.2)
𝑥𝑛+1 𝑥𝑛 𝐹(𝑥𝑛, 𝑦𝑛 )
[𝑦 ] = [𝑦 ] − 𝐽−1 [ ]
𝑛+1 𝑛 𝐺(𝑥𝑛, 𝑦𝑛 )
Donde nuestros valores iniciales son:
𝑥0 = 0.2, 𝑦0 = 0.6
Calculamos 𝐹(𝑥0, 𝑦0 ) 𝑦 𝐺(𝑥0, 𝑦0 ):
𝐹(0.2, 0.6) = −0.2 2 + 0.2 + 0.75 − 0.6 = 0.59
𝐺(0.2, 0.6) = 0.6 + 5 ∗ 0.2 ∗ 0.6 − 2 = −0.8
Calculamos 𝐽(𝑥0, 𝑦0 ):
−2 ∗ 0.2 + 1 −1 0.6 −1
𝐽(0.2, 0.6) = [ ]=[ ]
5 ∗ 0.6 1 + 5 ∗ 0.2 3 2
Invertimos J
1 2 1 1 2 1
𝐽−1 = [ ]= [ ]
(0.6 ∗ 2) − (−1 ∗ 3) −3 0.6 1.2 + 3 −3 0.6
0.476 0.238
𝐽−1 = [ ]
−0.714 0.143
Y se procede a calcular la corrección donde:
∆𝑥 𝐹(𝑥0, 𝑦0 ) 0.476 0.238 −0.56
[ ] = 𝐽−1 [ ]=[ ][ ]
∆𝑦 𝐺(𝑥0 𝑦0 ) −0.714 0.143 0.8
0.476(−0.56) + 0.238(0.8)
=[ ]
−0.714(−0.56) + 0.143(0.8)
−0.2808 + 0.1904 −0.09044
=[ ]=[ ]
0.4218 + 0.1144 0.5362
Actualizamos los valores
𝑥1 = 𝑥0 +△ 𝑥 = 0.2 − 0.09044 = 0.10956
𝑦1 = 𝑦0 △ 𝑦 = 0.6 + 0.53626 = 1.13626
Repetimos el proceso iterativo hasta alcanzar la convergencia.
𝐒𝐞𝐠𝐮𝐧𝐝𝐚 𝐫𝐚í𝐳:
𝑉𝑎𝑙𝑜𝑟𝑒𝑠 𝑖𝑛𝑖𝑐𝑖𝑎𝑙𝑒𝑠: 𝑥0 = 0.7, 𝑦0 = 0.2
𝐶𝑎𝑙𝑐𝑢𝑙𝑎𝑚𝑜𝑠 𝐽(𝑥0, 𝑦0 ) 𝑦 𝐺(𝑥0, 𝑦0 ):
𝐹(0.7,0.2) = −0.7 2 + 0.7 + 0.75 − 0.2 = 0.56
𝐺(0.7,0.2) = 0.2 + 5 ∗ 0.7 ∗ 0.2 − 2 = −0.8
Calculamos J(𝑥0, 𝑦0 ):
−2 ∗ 0.7 + 1 −1 −0.4 −1
𝐽(0.7,0.2) = [ ]=[ ]
5 ∗ 0.2 1 + 5 ∗ 0.7 1 4.5
Investimos J
1 4.5 1 1 4.5 1
𝐽−1 = [ ]= [ ]
(0.4 ∗ 4.5) − (−1 ∗ 1) −1 −0.4 −1.8 + 1 −1 −0.4
1 4.5 1
= [ ]
−0.8 −1 −0.4
−5.625 −1.25
𝐽−1 = [ ]
1.25 0.5
Calculamos la corrección
∆𝑥 𝐹(𝑥0, 𝑦0 ) −5.625 −1.25 −0.56
[ ] = 𝐽−1 [ ]=[ ][ ]
∆𝑦 𝐺(𝑥0 𝑦0 ) 1.25 0.5 0.8
−5.625(−0.56) + −1.25(0.8)
=[ ]
1.25(−0.56) + 0.5(0.8)
3.15 −1 2.15
=[ ]=[ ]
0.7 −4 −0.3
Actualizamos los valores
𝑥1 = 𝑥0 +△ 𝑥 = 0.7 + 2.15 = 2.85
𝑦1 = 𝑦0 +△ 𝑦 = 0.2 − 0.3 = −0.1
Ejercicio 3
Fueron escogidas tres longitudes de onda a las cuales la absorbancia de tres
especies, Co2+, Ni2+ y Cu2+, difieren significativamente. Las absortividades
de las tres especies a las tres longitudes de onda son mostradas en la tabla.
Una mezcla de los tres iones metálicos dio las siguientes absorbancias a las tres
longitudes de onda: 394 nm, 0.845; 510 nm, 0.388; 808 nm, 1.696, cuando se
midieron con una celda de 1.00 cm de longitud. Proponer un sistema de
ecuaciones lineales y calcular la concentración de los tres iones metálicos en la
mezcla, usando la ley de Beer A= bc, (A=absorbancia, =absortividad molar,
b=longitud de la celda en cm, c= concentración en mol/L.
Nota : poner todos los pasos para plantear el problema correctamente
Respuesta
Para la longitud de onda 394 nm:
0.845 = (0.995𝑐𝐶𝑜 + 6.868𝑐𝑁𝑖 + 0.188𝑐𝐶𝑢) ∗ 1.00
Para la longitud de onda 510 nm:
0.388 = (6.450𝑐𝐶𝑜 + 0.215𝑐𝑁𝑖 + 0.198𝑐𝐶𝑢) ∗ 1.00
Para la longitud de onda 808 nm:
1.696 = (0.469𝑐𝐶𝑜 + 1.179𝑐𝑁𝑖 + 15.052𝑐𝐶𝑢) ∗ 1.00
Y se genera un sistema de ecuaciones que queda de esta manera:
0.995 6.868 0.188 𝑐𝐶𝑜 0.845
[6.450 0.415 0.198 ] [ 𝑐𝑁𝑖 ] = [ 0.388 ]
0.469 1.179 15.052 𝑐𝐶𝑢 1.696
Si aplicamos la eliminación gaussiana se obtiene que:
𝑐𝐶𝑜 0.108 𝑚𝑜𝑙
[ 𝑐𝑁𝑖 ] = [0.073]
𝑙
𝑐𝐶𝑢 0.103
Ejercicio 4
Es frecuente que las reacciones químicas sigan este modelo:
 𝑑𝑐
= −𝑘 ∗ 𝑐 𝑛
𝑑𝑡
Donde c=concentración, t=tiempo, k=tasa de reacción, y n=orden de reacción.
Es posible evaluar valores dados de c y dc/dt, k y n por medio de regresión lineal
del logaritmo de esta ecuación:
𝑑𝑐
log ( ) = 𝑙𝑜𝑔 − 𝑘 ∗ 𝑛𝑙𝑜𝑔𝑐
𝑑𝑡
Use este enfoque y los datos que siguen para estimar los valores de k y n:

Respuesta
𝑑𝑐
ln ( ) = ln(−𝑘) + 𝑛 ∗ ln(𝑐)
𝑑𝑡
𝑑𝑐
𝑦 = 𝑙𝑛 ( ) , 𝑥 = 𝑙𝑛(𝑐) , 𝑎 = 𝑙𝑛(−𝑘) , 𝑏 = 𝑏
𝑑𝑡

Y por ende la ecuación se convierte: 𝑦 = 𝑎 + 𝑏𝑥

𝑃𝑎𝑟𝑎 𝑇 = 10, 𝑐 = 3.52 𝑥 = ln(3.52) ≈ 1.2599
𝑃𝑎𝑟𝑎 𝑇 = 20, 𝑐 = 2.48 𝑥 = ln(2.48) ≈ 0.9091
𝑃𝑎𝑟𝑎 𝑇 = 30, 𝑐 = 1.75 𝑥 = ln(1.75) ≈ 0.5596
𝑃𝑎𝑟𝑎 𝑇 = 40, 𝑐 = 1.23 𝑥 = ln(1.23) ≈ 0.2070
𝑃𝑎𝑟𝑎 𝑇 = 50, 𝑐 = 0.87 𝑥 = ln(0.87) ≈ −0.1393
𝑃𝑎𝑟𝑎 𝑇 = 60, 𝑐 = 0.61 𝑥 = ln(0.61) ≈ −0.4943

2.48 − 3.52
𝑃𝑎𝑟𝑎 𝑇 = 10, 𝑐 = 3.52 𝑦 = ln ( ) ≈ −0.3567
20 − 10
1.75 − 2.48
𝑃𝑎𝑟𝑎 𝑇 = 20, 𝑐 = 2.48 𝑦 = ln ( ) ≈ 0.4700
30 − 20
1.23 − 1.75
𝑃𝑎𝑟𝑎 𝑇 = 30, 𝑐 = 1.75 𝑦 = ln ( ) ≈ −0.4700
40 − 30
 0.87 − 1.23
𝑃𝑎𝑟𝑎 𝑇 = 40, 𝑐 = 1.23 𝑦 = ln ( ) ≈ −0.4700
50 − 40
0.61 − 0.87
𝑃𝑎𝑟𝑎 𝑇 = 50, 𝑐 = 0.87 𝑦 = ln ( ) ≈ −0.4700
60 − 50
0.61 − 0.87
𝑃𝑎𝑟𝑎 𝑇 = 60, 𝑐 = 0.61 𝑦 = 𝑙𝑛 ( ) ≈ −0.4700
70 − 60
𝛴𝑥 = 1.2599 + 0.9091 + 0.5596 + 0.2070 − 0.1393 − 0.4943
= 2.3018
𝛴𝑦 = −2.2618 − 2.6173 − 2.9552 − 3.3242 − 3.6450 − 3.6450
= −18.4485
𝛴𝑥𝑦 = (1.2599 ∗ −2.2618) + (0.9091 ∗ −2.6173)
+ (0.5596 ∗ −2.9552) + (0.2070 ∗ −3.3242)
+ (−0.1393 ∗ −3.6450) + (−0.4943 ∗ −3.6450)
≈ −2.8475 ± 2.3775 ± 1.6525 ± 0.6878 + 0.5079 + 1.8021
= −5.2553
𝛴𝑥 2 = (1.2599)2 + (0.9091)2 + (0.5596)2 + (0.2070)2
+ (−0.1393)2 + (−0.4943)2
= 1.5874 + 0.8264 + 0.3132 + 0.0429 + 0.0194 + 0.2443
= 3.0336
Obtenemos b y a
𝑏 = 0.8437
𝑎 = −3.3985
Obtenemos valores de K y n
𝑘 ≈ 0.0336𝑠 −1
𝑛 ≈ 0.8437
Ejercicio 5
Se sabe que el número de pulgadas de una estructura recién construida que se
hunde en el suelo está dada por:
 𝑦 = 3 − 3𝑒 −𝑎𝑥
Donde x es el número de meses que lleva construida la estructura. Con los
valores

Respuesta
x y y'
2 1.07 0.65752
4 1.88 0.113329
6 2.26 -0.301105
12 2.78 -1.514127
18 2.97 -3.506557
24 2.99 -4.60517

𝑏 ′ = −0.2437
𝑎 = −𝑏 ′ = 0.2437
Ejercicio 6
El área de la superficie horizontal As(z) de un lago a cierta profundidad, se
calcula a partir del volumen por medio de diferenciación:
0
𝑑𝑉(𝑧)
𝐴𝑠 (𝑧) = −
𝑑𝑧
donde V=volumen (m3 ), y z=profundidad (m), se mide a partir de la superficie
en dirección del fondo. Determinar
𝑧
∫ 𝐴𝑠 (𝑧)𝑑𝑧
0

Con base en los siguientes datos

Respuesta

Ejercicio 7
Escriba la siguiente ecuación diferencial como un sistema de ecuaciones
diferenciales ordinarias de primer orden. Resolver el sistema de ecuaciones
usando el método de Euler de t=0 hasta t=1
𝑦 ′′ − 2𝑦 ′ + 2𝑦 = 𝑒 2𝑡 𝑠𝑒𝑛𝑡 0 ≤ 𝑡 ≤ 1
𝑦(0) = −0.4, 𝑦 ′ (0) = −0.6 ℎ = 0.1
Respuesta