LÓPEZ CASTAÑÓN AMÉRICA RUBÍ

Tarea
Problema 1
Las coordenadas de reacción de dos reacciones que ocurren simultáneamente ε1 y ε2 se calculan
resolviendo el siguiente sistema de ecuaciones:
4 ∗ 𝜀1 3
= 122.97
(1 − 𝜀1 − 2 ∗ 𝜀2 ) ∗ (1 + 2 ∗ 𝜀1 )2
𝜀2 2
= 21.42
(1 − 𝜀1 − 2 ∗ 𝜀2 )

Encuentre los valores de las coordenadas ε1 y ε2 mediante la aplicación del método de

Newton-Raphson. Inicie sus cálculos con los valores ε1(0)= 0.965 y ε2(0)= 0.016. Realice un
máximo de 10 iteraciones o hasta que |fi | ≤ 0.0001
Solución: (0.9651,0.0158)
• Igualando las ecuaciones a cero.
4 ∗ 𝜀1 3
𝑓1 (𝑥, 𝑦) = − 122.97 = 0
(1 − 𝜀1 − 2 ∗ 𝜀2 ) ∗ (1 + 2 ∗ 𝜀1 )2
𝜀2 2
𝑓2 (𝑥, 𝑦) = − 21.42 = 0
(1 − 𝜀1 − 2 ∗ 𝜀2 )2
• Evaluando f1 y f2 en los valores iniciales.

4 ∗. 9653
𝑓1 (𝑥, 𝑦) = − 122.97 = 16.5979
(1 − .965 − 2 ∗ .016) ∗ (1 + 2 ∗ .965)2
. 0162
𝑓2 (𝑥, 𝑦) = − 21.2 = 7.0244
(1 − .965 − 2 ∗ .016)2
• Las derivadas parciales están dadas por:
𝜕𝑓1 −2 ∗ 𝑥 2 ∗ (𝑥 ∗ 𝑦 + 1.5 ∗ (𝑦 − .5))
= 2
𝜕𝑥 (𝑥 + .5)3 ∗ (𝑥 + 2 ∗ (𝑦 − .5))
𝜕𝑓1 . 5 ∗ 𝑥3
=
𝜕𝑦 (𝑦 + .5 ∗ (𝑥 − 1))2 ∗ (𝑥 + .5)2

𝜕𝑓2 −2 ∗ 𝑦 2
= 3
𝜕𝑥 (𝑥 + 2 ∗ (𝑦 − .5))
 LÓPEZ CASTAÑÓN AMÉRICA RUBÍ
𝜕𝑓2 . 25 ∗ 𝑦 ∗ (𝑥 − 1)
=
𝜕𝑦 (𝑦 + .5 ∗ (𝑥 − .1))3

• Evaluando las derivadas parciales

𝜕𝑓1 −2 ∗ 𝑥 2 ∗ (𝑥 ∗ 𝑦 + 1.5 ∗ (𝑦 − .5))
= 2 = 46766
𝜕𝑥 (𝑥 + .5)3 ∗ (𝑥 + 2 ∗ (𝑦 − .5))
𝜕𝑓1 . 5 ∗ 𝑥3
= = 93045.2
𝜕𝑦 (𝑦 + .5 ∗ (𝑥 − 1))2 ∗ (𝑥 + .5)2

𝜕𝑓2 −2 ∗ 𝑦 2
= 3 = 18963
𝜕𝑥 (𝑥 + 2 ∗ (𝑦 − .5))
𝜕𝑓2 . 25 ∗ 𝑦 ∗ (𝑥 − 1)
= = 41481.5
𝜕𝑦 (𝑦 + .5 ∗ (𝑥 − .1))3

• Calculando el jacobiano
𝜕𝑓1 𝑥 ∗ 𝜕𝑓2 𝑦 − 𝜕𝑓1 𝑦 ∗ 𝜕𝑓2 𝑥 = 𝐽
(46766 ∗ 41481.5) − (93045.2 ∗ 18963) = 1.755 ∗ 10^8
• Calculando ∆𝑥 , ∆𝑦
𝑓2 ∗ 𝜕𝑓1 𝑦 − 𝑓1 ∗ 𝜕𝑓2 𝑦
= ∆𝑥 = −.000199
𝐽
𝑓1 ∗ 𝜕𝑓2 𝑥 − 𝑓2 ∗ 𝜕𝑓1 𝑥
= ∆𝑦 = −.000199
𝐽
• Nuevos valores de “x” y “y”
𝑥 + ∆𝑥 = .965 − .000199 = .9648

k x y 𝑓1 𝑓2 𝜕𝑓1 𝜕𝑓1 𝜕𝑓2 𝜕𝑓2 𝐽 ∆𝒙 ∆𝒚

𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦

0 .965 .016 16.5979 7.0244 46766 93045.2 18963 41481.5 1.755*10^8 - -

.000199 .000078

1 . 9648 .01592 1.7595 1.0912 37386.7 74338.2 13416.6 29661 1.1155*10^8 .000259 -
.000154
2 .96506 .01576 .02523 .0285 36340.8 72252.7 12599.8 27920.2 1.0427*10^8 .000013 -
.000007
3 .96507 .01576 0 0
 LÓPEZ CASTAÑÓN AMÉRICA RUBÍ

Problema 2
Para calcular la concentración de un reactivo en dos reactores interconectados se llegó a las
siguientes ecuaciones:

Donde kp1 = 2 y kp2 = 4. Realice 10 iteraciones con el Método de Newton-Raphson

Multivariable (con Jacobiano) para encontrar los valores de C1y C2utilizando como valores
iniciales ( 1 , 0 ).

Solución: ( 0.5 , 0 )

• Evaluando f1 y f2 en los valores iniciales.

𝟏 𝟎
𝒇𝟏 = ∗ 𝒔𝒆𝒏(𝟏 ∗ 𝟎) − =𝟎
𝟒 𝟒∗𝝅
𝟏 𝒆∗𝟎
(𝟏 − ) ∗ (𝒆𝟐∗𝟏 − 𝒆) + = 𝟒. 𝟐𝟗𝟗𝟎
𝟒∗𝝅 𝝅
• Las derivadas parciales están dadas por:
𝜕𝑓1 cos (𝑥 ∗ 𝑦)
= ∗𝑦
𝜕𝑥 2
𝜕𝑓1 cos(𝑥 ∗ 𝑦) 1
= ∗𝑥−
𝜕𝑦 2 4∗𝜋
𝜕𝑓2 (4 ∗ 𝜋 − 1) ∗ 𝑒 2𝑥
=
𝜕𝑥 2∗𝜋
𝜕𝑓2 𝑒
=
𝜕𝑦 𝜋
• Evaluando las derivadas parciales
𝜕𝑓1 cos (1 ∗ 0)
= ∗0=0
𝜕𝑥 2
𝜕𝑓1 cos(1 ∗ 0) 1
= ∗1− = .42042
𝜕𝑦 2 4∗𝜋
𝜕𝑓2 (4 ∗ 𝜋 − 1) ∗ 𝑒 2∗1
= = 13.6021
𝜕𝑥 2∗𝜋
 LÓPEZ CASTAÑÓN AMÉRICA RUBÍ
𝜕𝑓2 𝑒
= = .86526
𝜕𝑦 𝜋
• Calculando el jacobiano
𝜕𝑓1 𝑥 ∗ 𝜕𝑓2 𝑦 − 𝜕𝑓1 𝑦 ∗ 𝜕𝑓2 𝑥 = 𝐽
(0 ∗ .86526) − (. 42042 ∗ 13.6021) = 1.755 ∗ 10^8
• Calculando ∆𝑥 , ∆𝑦
𝑓2 ∗ 𝜕𝑓1 𝑦 − 𝑓1 ∗ 𝜕𝑓2 𝑦
= ∆𝑥 = −.31606
𝐽
𝑓1 ∗ 𝜕𝑓2 𝑥 − 𝑓2 ∗ 𝜕𝑓1 𝑥
= ∆𝑦 = 0
𝐽
• Nuevos valores de “x” y “y”
𝑥 + ∆𝑥 = 1 − .31606 = .68394

k x y 𝑓1 𝑓2 𝜕𝑓1 𝜕𝑓1 𝜕𝑓2 𝜕𝑓2 𝐽 ∆𝒙 ∆𝒚

𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦
0 1 0 0 4.29909 0 .42042 13.6021 .86526 -5.71863 -.31606 0

1 .68394 0 0 1.11254 0 .26239 7.22903 .86526 -1.89684 -.1539 0

2 .53004 0 0 .15493 0 .18544 5.31379 .86526 -.98540 - 0

.029155
3 .50088 0 0 .00443 0 .17087 5.0128 .86526 -.85651 - 0
.000884
4 .5 0 0 0
 LÓPEZ CASTAÑÓN AMÉRICA RUBÍ

Problema 3
Resuelva el siguiente sistema de ecuaciones no lineales utilizando el método de Newton-
Raphson Multivariable:
f(t,c) = 55et + tc2 – 2 = 0

g(t,c) =
Realice 10 iteraciones con x(0) = ( 1 , 1 ) y un criterio de convergencia ε = 0.01
• Evaluando f y g en los valores iniciales.
𝒇 = 𝟓𝟓 ∗ 𝒆𝒕 + 𝒕 ∗ 𝒄𝟐 − 𝟐 = 𝟏𝟒𝟖. 𝟓𝟎𝟔
𝟑
𝒈 = 𝒆𝒕 − 𝟐𝟎 ∗ 𝒄 = −𝟏𝟕. 𝟐𝟖𝟏𝟕
• Las derivadas parciales están dadas por:
𝜕𝑓
= 55 ∗ 𝑒 𝑡 + 𝑐 2
𝜕𝑡
𝜕𝑓
=2∗𝑡∗𝑐
𝜕𝑐
𝜕𝑔 3
= 3 ∗ 𝑡2 ∗ 𝑒𝑡
𝜕𝑡
𝜕𝑔
= −20
𝜕𝑐

• Evaluando las derivadas parciales

𝜕𝑓
= 55 ∗ 𝑒 1 + 12 = 15
𝜕𝑡
𝜕𝑓
= 2 ∗ 1 ∗ 1 = −17
𝜕𝑐
𝜕𝑔 3
= 3 ∗ 12 ∗ 𝑒 1
𝜕𝑡
𝜕𝑔
= −20
𝜕𝑐
• Calculando el jacobiano
𝜕𝑓1 𝑥 ∗ 𝜕𝑔𝑦 − 𝜕𝑓1 𝑦 ∗ 𝜕𝑔𝑥 = 𝐽
(150.506 ∗ −20) − (2 ∗ 8.15485) = −3026.42
• Calculando ∆𝑥 , ∆𝑦
𝑔 ∗ 𝜕𝑓𝑐 − 𝑓 ∗ 𝜕𝑔𝑐
= ∆𝑡 = −.96997
𝐽
 LÓPEZ CASTAÑÓN AMÉRICA RUBÍ
𝑓 ∗ 𝜕𝑔𝑡 − 𝑔 ∗ 𝜕𝑓𝑡
= ∆𝑐 = −1.25959
𝐽
• Nuevos valores de “x” y “y”
𝑡 + ∆𝑡 = 1 − .96997 = .68394
𝑐 + ∆𝑐 = 1 − .96997 = .68394

k t c 𝑓 𝒈 𝜕𝑓 𝜕𝑓 𝜕𝑔 𝜕𝑔 𝐽 ∆𝒙 ∆𝒚
𝜕𝑡 𝜕𝑐 𝜕𝑡 𝜕𝑐

0 1 1 148.506 - 150.506 2 8.15485 - - -.96997 -

17.2817 20 3026.42 1.25959
1 .03003 - 54.6785 6.19173 56.7439 - .002705 - - -.96352 .30946
.25959 .01559 20 1134.88
2 -.93349 .04987 19.6225 -.55409 21.6273 - 1.15895 - - -.90765 -.08030
.09311 20 432.439
3 - - 6.72332 .61054 8.72595 .11205 .019808 - - - .029764
1.84114 .03043 20 174.521 .770879
4 -2.6120 - 2.03625 .01332 4.03625 .00348 3.728*10^- - -80.725 -.50449 .00066
.00066 7 20
5 -3.1165 0 .43714 0 2.43714 0 0 - - -.17936 0
20 48.7428
6 -3.2958 0 .036961 0 2.03696 0 0 - - - 0
20 40.7392 .018145
7 -3.3140 0 .0003 0