EJERCICIOS CAPITULO II

EJERCICIO 1.- En el cálculo hidráulico de tuberías se utiliza la ecuación de

Darcy - Prandtl- Colebrook que proporciona las pérdidas de presión (h) mediante
la expresión:
𝜆. 𝐿. 𝑢2
ℎ=
2. 𝐷. 𝑔
Dónde: D es el diámetro de la tubería (en m), L la longitud de la tubería (en m),
u es la velocidad del fluido que por ella circula (en m. s−1), g es la aceleración
gravitatoria (en m. s−2) y 𝜆 es el coeficiente de fricción que puede estimarse a su
vez mediante la ecuación de Colebrook:
−1⁄ 2,51 𝐾
𝜆 2 = −2. ln⁡( 1−2
+ )
𝑅𝑒. 𝜆 3,71. 𝐷
Donde Re es el número de Reynolds (Re = u·D/µ), µ es la viscosidad cinemática
del fluido (en m2·s−1) y K es la altura de rugosidad (en m). Calcúlese, mediante
el método de aproximaciones sucesivas, el valor del coeficiente de fricción de un
colector tubular recto sin acometidas para el que se sabe que K = 0,25 · 10 −3 m,
D = 0,3 m y por el que se quiere hacer circular un fluido de tal forma que el
número de Reynolds tome el valor Re = 200000.
Solución:
1) Definir la función:
2.51 𝐾
𝑓 (𝜆) = 𝜆−1/2 + 2 ln ( 1 + )=0
3,71. 𝐷
𝑅𝑒. 𝜆 2

DATOS
K 0,00025
D 0,3
Re 200000

Queda la función:
𝑓 (𝜆) = 𝜆−1/2 + 2 ln(0.00001255𝜆−1/2 + 0.00022462) = 0

2) Graficar la función:
Método de la Bisección:
a) Función a trabajar:
1 1
𝑓 (𝜆) = 𝜆−2 + 2 ln (0.00001255𝜆−2 + 0.00022462) = 0

b) Valores límites:
Valor límite inferior: a = 0.001

Valor límite superior: b = 0.005

a b m F(a) F(b) F(m) Error Iteración

 c) Condiciones:
Si f(a)*f(m) > 0 entonces a se cambia por m
Si f(a)*f(m) < 0 entonces b se cambia por m
d) Ecuación para hallar m:
𝑎+𝑏
𝑚=
2
e) Parámetros de control:
𝑚2 −𝑚1
𝐸𝑟𝑟𝑜𝑟 = | | 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 1 × 10−5
𝑚2

Método de Regula Falsi:

a) Función a trabajar:
1 1
𝑓 (𝜆) = 𝜆−2 + 2 ln (0.00001255𝜆−2 + 0.00022462) = 0

b) Valores límites:

Valor límite inferior: a = 0.001

Valor límite superior: b = 0.005
c) Condiciones:
Si f(a)*f(m) > 0 entonces a se cambia por m
Si f(a)*f(m) < 0 entonces b se cambia por m
d) Ecuación para hallar m:
𝑎+𝑏
𝑚=
2
e) Parámetros de control
 𝑚2 −𝑚1
𝐸𝑟𝑟𝑜𝑟 = | | 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 1 × 10−5
𝑚2

a b m F(a) F(b) F(m) Error Iteración

Método de Newton:
a) Función a trabajar:
1 1
𝑓 (𝜆) = 𝜆−2 + 2 ln (0.00001255𝜆−2 + 0.00022462) = 0

Nota: La ecuación es complicada para derivar, así que, se creara

una nueva función

𝑓 (𝜆) = 𝑥 + 2 ln(0.00001255. 𝑥 + 0.00022462) = 0

1
Donde: 𝜆−2 = 𝑥

b) Derivada de la función:
2 × 0.00001255
𝑓 ′ (𝑥 ) = 1 +
0.00001255⁡. 𝑥 + 0.00022462
c) Determinar valor inicial:
 d) Ecuación para hallar m:
𝑓(𝑥)
𝑚=𝑎−
𝑓′(𝑥)
e) Parámetros de control
𝑚2 −𝑚1
𝐸𝑟𝑟𝑜𝑟 = | | 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 1 × 10−5
𝑚2

X=15.55151346
1 1
1
𝑥 = 𝜆−2 𝑥 = 𝜆−2 = 15.551513462

Método de Punto Fijo:

a) Función a trabajar:
1 1
𝑓 (𝜆) = 𝜆−2 + 2 ln (0.00001255𝜆−2 + 0.00022462) = 0

b) Hallar G(x):
1
𝐺 (𝜆 ) = ( 1 )2
−2. ln(0.00001255 . 𝜆−2 + 0.00022462⁡
c) Parámetros de control
𝑚2 −𝑚1
𝐸𝑟𝑟𝑜𝑟 = | | 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 1 × 10−5
𝑚2
Método de la Secante:
a) Función a trabajar:
1 1
𝑓 (𝜆) = 𝜆−2 + 2 ln (0.00001255𝜆−2 + 0.00022462) = 0

b) Valores iniciales:
X0=0.001
X1=0.05
c) Ecuación para hallar x2
𝑋1 − 𝑋0
𝑋2 = 𝑋0 − × 𝑓 (𝑋0 )
𝑓 (𝑋1 ) − 𝑓 (𝑋0 )
d) Parámetros de control
𝑚2 −𝑚1
𝐸𝑟𝑟𝑜𝑟 = | 𝑚2
| 𝑇𝑜𝑙𝑒𝑟𝑎𝑛𝑐𝑖𝑎 = 1 × 10−5

EJERCICIO 2.- Resolver utilizando los métodos numéricos utilizados F(x) = 2x2 - x - 5=0

1. METODO DE BISECCIÓN

Ecuación

𝒇(𝒙) = 𝟐𝒙𝟐 − 𝒙 − 𝟓 = 𝟎

x y
-3 16
-2 5
-1 -2
0 -5
1 -4
2 1
3 10
Iteración 1

a=1

b=2

m1= (1+2)/2=1.5

f(a1)=2*(1)^2-(1)-5=-4

f(b1)=2*(2)^2-(2)-5=1

f(m1)=2*(1.5)^2-(1.5)-5=-2

Error=(2-1)/2=0.5

Iteración 2

a=1.5

b=2

m2=(1.5+2)/2=1.75

f(a2)=2*(1.5)^2-(1.5)-5=-2

f(b2)=2*(2)^2-(2)-5=1

f(m2)=2*(1.75)^(1.75)-5=-0.62

Error= ABS(1.75-1.5)/1.75

=SI (error>=tolerancia,"CONTINUA","FIN")
Método de bisección Tolerancia 0.00001
xi xs xr fxi fxs fxr 𝒇
a b m f(a) f(b) f(m) Error Iteración
1 2 1.5 -4.00000000 1.00000000 -2.00000000 0.50000000 1
1.5 2 1.75 -2.00000000 1.00000000 -0.62500000 0.14285714 2 Continua
1.75 2 1.875 -0.62500000 1.00000000 0.15625000 0.06666667 3 Continua
1.75 1.875 1.8125 -0.62500000 0.15625000 -0.24218750 0.03448276 4 Continua
1.8125 1.875 1.84375 -0.24218750 0.15625000 -0.04492188 0.01694915 5 Continua
1.84375 1.875 1.859375 -0.04492188 0.15625000 0.05517578 0.00840336 6 Continua
1.84375 1.859375 1.8515625 -0.04492188 0.05517578 0.00500488 0.00421941 7 Continua
1.84375 1.8515625 1.84765625 -0.04492188 0.00500488 -0.01998901 0.00211416 8 Continua
1.84765625 1.8515625 1.84960938 -0.01998901 0.00500488 -0.00749969 0.00105597 9 Continua
1.84960938 1.8515625 1.85058594 -0.00749969 0.00500488 -0.00124931 0.00052770 10 Continua
1.85058594 1.8515625 1.85107422 -0.00124931 0.00500488 0.00187731 0.00026378 11 Continua
1.85058594 1.85107422 1.85083008 -0.00124931 0.00187731 0.00031388 0.00013191 12 Continua
1.85058594 1.85083008 1.85070801 -0.00124931 0.00031388 -0.00046775 0.00006596 13 Continua
1.85070801 1.85083008 1.85076904 -0.00046775 0.00031388 -0.00007694 0.00003298 14 Continua
1.85076904 1.85083008 1.85079956 -0.00007694 0.00031388 0.00011847 0.00001649 15 Continua
1.85076904 1.85079956 1.8507843 -0.00007694 0.00011847 0.00002076 0.00000824 16 Fin

Raíz =1.85079956

2. MÉTODO DE LA REGLA FALSA

𝒇(𝒙) = 𝟐𝒙𝟐 − 𝒙 − 𝟓 = 𝟎
 20
x y
15
-3 16
-2 5 10

-1 -2 5
0 -5 y 0
1 -4 -4 -3 -2 -1 0 1 2 3 4
-5
2 1
3 10 -10

Iteración 1

a=1

b=2

f(a1)=2*(1)^2-(1)-5=-4

f(b1)=2*(2)^2-(2)-5=1

m=(1*1-(-4)*2)/(1-(-4))=1.8

f(m1)=2*(1.8)^2-(1.8)-5= -0.32

Error=(2-1)/2=0.5

Iteración 2

a=1

b=1.8

f(a2)=2*(1)^2-(1)-5=-4

f(b2)=2*(1.8)^2-(1.8)-5=-0.32

m=(0.32*-4-(1)*1.8)/(0.32-(-4))=1.869

f(m2)=2*(1.869)^2-(1.869)-5= 0.1209

Error=(0.1209-1.8)/0.1209=0.0372093

=SI (error>=tolerancia,"CONTINUA","FIN")
Metodo de la regla falsa Tolerancia 0.00001 𝒇
xi xs xr fxi fxs fxr
a b m f(a) f(b) f(m) Error Iteración
1 2 1.8 -4 1 -0.32 0.5 1
1 1.8 1.869565 -4 -0.32 0.12098299 0.0372093 2 Continua
1.86956522 1.8 1.85048 0.12098299 -0.32 -0.00192684 0.0103136 3 Continua
1.85048011 1.8 1.850786 -0.00192684 -0.32 3.1061E-05 0.00016523 4 Continua
1.85078591 1.8 1.850781 3.1061E-05 -0.32 -5.006E-07 2.6632E-06 5 Continua
1.85078098 1.8 1.850781 -5.006E-07 -0.32 8.0682E-09 4.2923E-08 6 Fin
1.85078106 1.8 1.850781 8.0682E-09 -0.32 -1.3004E-10 6.9179E-10 7 Fin

Raíz= 1.850781

3. METODO DE NEWTON RAPHSON

𝒇(𝒙) = 𝟐𝒙𝟐 − 𝒙 − 𝟓 = 𝟎

𝑓 ′ (𝑥) = 4𝑥 − 1

y
20
x y 15
-3 16
10
-2 5
-1 -2 5
0 -5 0
1 -4 -4 -3 -2 -1 0 1 2 3 4
-5
2 1
3 10 -10

Iteración 1
a=1
f(a1)=2*(1) ^-(1)-5=-4
f´(a1)=4*(1)-1=3
m1=1-(-4/3)=2.333
Iteración 2
a=m=2.33
f(a2)=2*(2.33) ^-(2.33)-5=3.555
f´(a2)=4*(1)-1=3=8.333
m2=2.33-(3.555/8.333)=1.90666
Error=ABS((1.90666-2.333)/ 1.90666)

=SI (error>=tolerancia,"CONTINUA","FIN")
 Metodo de Newton Rhapson Tolerancia 0.00001
xi xr fxi
a m f(a) f´(a) Error Iteración
1 2.33333333 -4 3 1
2.33333333 1.90666667 3.55555556 8.33333333 0.22377622 2 Continua
1.90666667 1.85172368 0.36408889 6.62666667 0.02967127 3 Continua
1.85172368 1.85078134 0.00603746 6.4068947 0.00050916 4 Continua
1.85078134 1.85078106 1.776E-06 6.40312535 1.4986E-07 5 Fin
1.85078106 1.85078106 1.5454E-13 6.40312424 1.3077E-14 6 Fin
Raíz = 1.850781

4. METODO DE PUNTO FIJO

𝒇(𝒙) = 𝟐𝒙𝟐 − 𝒙 − 𝟓 = 𝟎

1
𝑥+5 𝑔′(𝑥)=23/2
𝑥=√ √𝑥+5
2

7
6
x x y
5
0 0 1.58113883
4
1 1 1.73205081
3
2 2 1.87082869
3 3 2 2

4 4 2.12132034 1
5 5 2.23606798 0
0 2 4 6 8
6 6 2.34520788

Iteración 1
X1=2
g(x1)=√(2 + 5)/2=1.8708

Iteración 2
X2=g(x)=1.8708
g(x2)=√(1.8708 + 5)/2=1.85348
Error=(1.85348-1.8708)/ 1.85348=0.09356

=SI (error>=tolerancia,"CONTINUA","FIN")
 Tolerancia 0.00001

x(i) g(x) Error Iteración

2 1.87082869 1
1.87082869 1.85348708 0.00935621 2 Continua
1.85348708 1.85114655 0.00126437 3 Continua
1.85114655 1.85083043 0.0001708 4 Continua
1.85083043 1.85078773 2.3071E-05 5 Continua
1.85078773 1.85078196 3.1164E-06 6 Fin
1.85078196 1.85078118 4.2096E-07 7 Fin
1.85078118 1.85078108 5.6863E-08 8 Fin
Raíz = 1.850781

5. METODO DE LA SECANTE
𝒇𝒙𝒊 (𝒙𝒊−𝒙𝒊−𝟏 )
𝒇(𝒙) = 𝟐𝒙𝟐 − 𝒙 − 𝟓 = 𝟎 𝒇𝒊+𝟏 =
𝒇𝒙𝒊 −𝒇(𝒙𝒊−𝟏)

y
20
x y 15
-3 16
10
-2 5
-1 -2 5
0 -5 0
1 -4 -4 -2 0 2 4
-5
2 1
3 10 -10

Iteración 0
a=1
f(a)=2(1)^2-(1)-5=-4
Iteración 1
a 1 =2
f(a1)=2(2) ^2-(2)-5=1
m1=2-(1*(2-1)/(1-(-4))=1.8
Error= (1.8-2)/1.8=0.1111
Iteración 2
a 2 =m1=1.8
f(a2)=2(1.8) ^2-(1.8)-5=-0.32
m2=1.8-(-0.32*(1.8-2)/(-0.32-(1.8))=1.8484
Error
ABS(1.8484-1.8)/1.8484=0.02622951

=SI (error>=tolerancia,"CONTINUA","FIN")
 Tolerancia 0.00001

a m f(a) Error Iteración

1 -4 0
2 1.8 1 0.11111111 1 Continua
1.8 1.84848485 -0.32 0.02622951 2 Continua
1.84848485 1.85081809 -0.01469238 0.00126066 3 Continua
1.85081809 1.85078103 0.00023714 2.0025E-05 4 Continua
1.85078103 1.85078106 -1.702E-07 1.4362E-08 5 Fin
1.85078106 1.85078106 -1.97E-12 1.6628E-13 6 Fin
1.85078106 1.85078106 0 0 7 Fin
Raíz = 1.850781

EJERCICIO 3.- Resolver usando los métodos numéricos utilizados 𝒙𝟑 − 𝒙𝟐 = 𝟓

1. METODO DE BISECCION

Ecuación
x3-x2=5 pero f(x)=x3-x2-5=0
 x y
-5 -155
-4 -85
-3 -41 y
-2 -17 150
-1 -7 100
0 -5
50
1 -5
2 -1 0
-6 -4 -2 0 2 4 6
3 13 -50
4 43
-100
5 95
-150

-200
ITERACION 1

a1 =2

b1 =2.5

m1= (2+2.5)/2= 2.25

f(a1)= (2)^3-(2)^2-5= -1

f(b1)= (2.5)^3-(2.5)^2-5= 4.375

f(m1)= (2.25)^3-(2.25)^2-5=1.328125

Error= (2.5-2/2)= 0.25

ITERACION 2

Si f(a)*f (m) < 0 entonces a se cambia por m

Si f (b)*f (m) < 0 entonces b se cambia por m
a2 =2

b2 =2.25

m2= (2+2.25)/2= 2.125

f(a2)= (2)^3-(2)^2-5= -1

f(b2)= (2.25)^3-(2.25)^2-5= 1.328125

f(m2)=(2.125)^3-(2.125)^2-5= 0.08007813

Error=ABS (2.125-2.25/2.125)= 0.05882353

=SI (error>=tolerancia,"CONTINUA","FIN")

…
RAICES
X1 2.11636353
X2 imaginario
X3 imaginario

2. METODO DE REGLA FALSA

Ecuación
x3-x2=5 pero f(x)=x3-x2-5=0
x y
-5 -155
y
-4 -85 150
-3 -41 100
-2 -17
50
-1 -7
0 -5 0
-6 -4 -2 0 2 4 6
1 -5 -50
2 -1 -100
3 13
-150
4 43
5 95 -200

ITERACION 1

a1 =2

b1 =2.5

f(a1)= (2)^3-(2)^2-5= -1
f(b1)= (2.5)^3-(2.5)^2-5= 4.375

m1= 2-((-1*(2.5-2))/(4.375-(-1)))= 2.09302326

f(m1)= (2.09302326)^3-(2.09302326)^2-5= -2.117424

Error= (2.5-2/2)= 0.25

ITERACION 2

Si f(a)*f (m) < 0 entonces a se cambia por m

Si f (b)*f (m) < 0 entonces b se cambia por m
a2 =2

b2 =2.09302326

f(a2)= (2)^3-(2)^2-5= -1

f(b2)= 2.09302326)^3-(2.09302326)^2-5= -2.117424

m2= 2-((-1*(2.09302326-2))/ (-2.117424-(-1)))= 2.11801123

f(m2)= (2.11801123)^3-(2.11801123)^2-5= 0.01536662

Error=ABS (2.11801123-2.09302326/2.11801123)= 0.01179785

=SI (error>=tolerancia,"CONTINUA","FIN")

RAICES
X1 2.11634361
X2 imaginario
X3 imaginario

3. METODO DE NEWTON RAPHSON

Ecuación
x3-x2=5 pero f(x)=x3-x2-5=0
f'(x)=3x2-2x=0
 y
x y 20
-3 -41 10
-2 -17 0
-1 -7 -4 -3 -2 -1 -10 0 1 2 3 4

0 -5 -20
1 -5 -30
2 -1 -40
3 13 -50

ITERACION 1

a1 =2

f(a1)= (2)^3-(2)^2-5= -1

f ‘(a1)= 3*(2)^2-2*(2)= 8

m1= 2-(-1/8)=2.125

ITERACION 2

a 2= m1=2.125

f(a2)= (2.125)^3-(2.125)^2-5= 0.08007813

f ‘(a2)= 3*(2.125)^2-2*(2.125)= 9.296875

m2= 2.125-(0.08007813/9.296875)=2.11638655

Error=ABS (2.11638655-2.125/2.11638655)= 0.00406988

=SI (error>=tolerancia,"CONTINUA","FIN")

RAICES
X1 2.1163433
X2 imaginario
X3 imaginario

4. METODO DE PUNTO FIJO

Ecuación
x3-x2=5 pero f(x)=x3-x2-5=0
g(x)=(5+x^2)1/3

Título del gráfico

x x g(x)
x g(x)
1 1 1.817121
2 2 2.080084 7
3 3 2.410142 6
5
4 4 2.758924
4
5 5 3.107233
3
6 6 3.448217 2
1
0
ITERACION 1 0 2 4 6 8
X1=2

g (x1)= (5+2^2)^(1/3)=2.080084

ITERACION 2

X2=g(x1) =2.080084

G(x2)= (5+2^2.080084)^(1/3)=2.104958

Error =ABS((2.080084-2)/ 2.080084)=0.0385003

=SI (error>=tolerancia,"CONTINUA","FIN")

RAICES
X1 2.116342
X2 imaginario
X3 imaginario

5. METODO DE LA SECANTE
Ecuación
x3-x2=5 pero f(x)=x3-x2-5=0

X Y
Y
-3 -41 20

-2 -17 10

-1 -7 0
-4 -2 -10 0 2 4
0 -5
1 -5 -20
-30
2 -1
-40
3 13
-50

ITERACION 0

X0=2

f(x0)= (2)^3-(2)^2-5= -1

ITERACION 1

X1=2.5

f(x1)= (2.5)^3-(2.5)^2-5= 4.375

m1= 2-((2.5-2)/( 4.375-(-1)))* -1 = 2.09302326

Error =ABS(2.09302326-2.5)/ 2.09302326= 0.19444444

=SI (error>=tolerancia,"CONTINUA","FIN")

ITERACION 2

X2=m1=2.09302326

f(x2)= (2.09302326)^3-(2.09302326)^2-5= -0.2117424

m2= 2.5-((2.09302326-2.5)/( -0.2117424-(4.375)))* 4.375= 2.11181093

Error =ABS(2.11181093-2.09302326)/ 2.11181093= 0.00889648

=SI (error>=tolerancia,"CONTINUA","FIN")
RAICES
X1 2.1163433
X2 imaginario
X3 imaginario
EJERCICIO4: Resolver usando los métodos numéricos utilizados la ecuación de REDLICH –
KWONG, para el CO2

Ecuación Redlich-Kwong

Donde:

Ecuación Polinomica de Redlich-Kwong

Método de la Secante

Determinar el volumen molar del 𝐶𝑂2

R= 8.314 m^3 Pa/mol k

P= 10000 Pa

T= 300 k

Tc= 553.6 k

Pc= 4073000Pa

Solución:

Hallamos a y b

a= 52.31327049 1era Exp. -0.24942

2da Exp. 0.000277601
b= 9.79066E-05
3era Exp. -2.95708E-08
PRIMERA ITERACIÓN

F(X)1= 0.0005+(-0.24942*0.00005)+(0.0002776*0.00005)+(-2.9508)=4.70

F(x)2=0.00025+(-0.24942*0.000025)+(0.0002776*0.000025)+(-2.9508)=2.43
M1=((0.00025)-2.43*(0.00005-0.00025))/(4.7-2.43)=-0.000016
Nuestro nuevo valor para “x” será igual a m1
SEGUNDA ITERACIÓN
F(x)2=0.000016+(-0.24942*0.0000016)+(0.0002776*0.000016)+(-2.9508)=-3.43

M1=((0.000016)-2.43*(0.000025-0.000016))/(2.43-(-3.43)))=-0.0001394

Raiz 1
x y
0.00005 -1.63142E-08

0.00006 -1.38124E-08
0.00007 -1.13605E-08 GRÁFICA 1
0.00008 -8.95847E-09 2E-08
0.00009 -6.60625E-09 1E-08
0.0001 -4.30387E-09
0.00011 -2.05131E-09 0
0 0.00005 0.0001 0.00015 0.0002 0.00025
0.00012 1.51439E-10 -1E-08
0.00013 2.30437E-09
-2E-08
0.00014 4.4075E-09
0.00015 6.46083E-09
0.00016 8.46436E-09
0.00017 1.04181E-08
0.00018 1.23221E-08
0.00019 1.41763E-08
0.0002 1.59807E-08

Tolerancia 0.00001
X m f(X) Error Iteración
0.0005 4.70E-08 0
0.00025 -0.00001663 2.43E-08 16.03257014 1 CONTINUAR
-1.66306E-05 0.00013947 -3.43E-08 1.11924217 2 CONTINUAR
0.000139469 0.00012207 4.30E-09 0.14252605 3 CONTINUAR
0.000122071 0.00011924 6.01E-10 0.02373961 4 CONTINUAR
0.00011924 0.00011931 -1.43E-11 0.00054967 5 CONTINUAR
0.000119305 0.00011931 4.52E-14 0.00000174 6 FIN

SEGUNDA RAIZ

PRIMERA ITERACIÓN

F(X)1= 0.2+(-0.24942*0.2)+(0.0002776*0.2)+(-2.9508)=-1.92

F(x)2=0.4+(-0.24942*0.4)+(0.0002776*0.4)+(-2.9508)=2.42

M1=((0.2)-2.43*(0.4-0.2))/(2.42-(-1.92))=0.2147
 Nuestro nuevo valor para “x” será igual a m 1
SEGUNDA ITERACIÓN
F(x)2=0.2147+(-0.24942*0.2147)+(0.0002776*0.2147)+(-2.9508)=-1.54

M1=((0.2147)-2.43*(0.4-0.2147))/((-1.54-2.42))=-0.2257

Raiz 2
x y
0.2 -0.00192131
0.21 -0.00168016 GRÁFICA 3
0.22 -0.00136289
0.01
0.23 -0.0009635
0.24 -0.000476 0.008

0.25 0.00010562 0.006

0.26 0.00078735 0.004
0.27 0.0015752 0.002
0.28 0.00247517
0
0.29 0.00349325 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
-0.002
0.3 0.00463545
0.31 0.00590776 -0.004

0.32 0.00731619
0.33 0.00886674

X m f(X) Error Iteración

0.2 -1.92E-03 0
0.4 0.21470852 2.42E-02 0.86299081 1 CONTINUAR
0.21470852 0.22579687 -1.54E-03 0.04910762 2 CONTINUAR
0.225796867 0.25753765 -1.14E-03 0.12324716 3 CONTINUAR
0.257537653 0.24648632 6.10E-04 0.04483547 4 CONTINUAR
0.246486322 0.24817296 -1.10E-04 0.00679623 5 CONTINUAR
0.248172962 0.24830440 -7.94E-06 0.00052936 6 CONTINUAR
0.248304405 0.24830248 1.18E-07 0.00000775 7 FIN

Método del punto fijo

PRIMERA ITERACIÓN

F(x)1=(0.0003)^3+(-0.2494*(0.0003)^2)+(0.0002776*0.0003)+(-2.95)=3.13

F´(X)1=(3*(0.0003)^2)-(0.4988*0.0003)+0.0002776=0.0001282
M=0.0003-(3.13/0.0001282)=0.00005597
El valor de M será nuestro nuevo valor de “x”
 SEGUNDA ITERACIÓN

F(x)1=(5.5974)^3+(-0.2494*(5.5974)^2)+(0.0002776*5.5974)+(-2.95)=-1.48

F´(X)1=(3*(5.5974)^2)-(0.4988*5.574)+0.0002776=0.0002496
M=5.5974-(-1.48/0.0002496)=0.00005597
El valor de M será nuestro nuevo valor de “x

Raiz 1
x y

0.00005 -1.63142E-08
0.00006 -1.38124E-08
0.00007 -1.13605E-08 GRÁFICA 1
0.00008 -8.95847E-09
2E-08
0.00009 -6.60625E-09
1E-08
0.0001 -4.30387E-09
0.00011 -2.05131E-09 0
0 0.00005 0.0001 0.00015 0.0002 0.00025
0.00012 1.51439E-10 -1E-08
0.00013 2.30437E-09 -2E-08
0.00014 4.4075E-09
0.00015 6.46083E-09
0.00016 8.46436E-09
0.00017 1.04181E-08
0.00018 1.23221E-08
0.00019 1.41763E-08
0.0002 1.59807E-08

Tolerancia 0.00001
X m f(Xo) f´(Xo) Error Iteración
0.0003 0.00005597 3.13E-08 0.00012822 0
5.59744E-05 0.00011530 -1.48E-08 0.00024969 0.51454259 1 CONTINUAR
0.000115302 0.00011929 -8.77E-10 0.00022012 0.03340371 2 CONTINUAR
0.000119287 0.00011931 -3.95E-12 0.00021814 0.00015195 3 CONTINUAR
0.000119305 0.00011931 -8.19E-17 0.00021813 0.00000000 4 FIN
0.000119305 0.00011931 0.00E+00 0.00021813 0.00000000 5 FIN
0.000119305 0.00011931 0.00E+00 0.00021813 0.00000000 6 FIN

PRIMERA ITERACIÓN

F(x)1=(0.0016)^3+(-0.2494*(0.0016)^2)+(0.0002776*0.0016)+(-2.95)=-0.00000022

F´(X)1=(3*(0.0016)^2)-(0.4988*0.0016)+0.0002776=0.0005128
 M=0.0016-(-0.00000022/0.0005128)=0.001171
El valor de M será nuestro nuevo valor de “x”

SEGUNDA ITERACIÓN

F(x)1=(0.001171)^3+(-0.2494*(0.001171)^2)+(0.0002776*5.5974)+(-2.95)=-0.00000005

F´(X)1=(3*(0.001171)^2)-(0.4988*0.001171)+0.0002776=-0.0003026
M=-(0.0011713)-(-0.00000005/-0.0003026)=0.001022
El valor de M será nuestro nuevo valor de “x

Raiz 2
x y
0.0002 1.59807E-08
0.0003 3.12888E-08 GRÁFICA 2
0.0004 4.16265E-08 0.0000001
0.0005 4.69999E-08 5E-08
0.0006 4.74148E-08
0
0.0007 4.28773E-08 0 0.0002 0.0004 0.0006 0.0008 0.001 0.0012 0.0014 0.0016 0.0018
-5E-08
0.0008 3.33935E-08
-1E-07
0.0009 1.89692E-08
0.001 -3.89441E-10 -1.5E-07

0.0011 -2.46765E-08 -2E-07

0.0012 -5.3886E-08 -2.5E-07
0.0013 -8.80118E-08
0.0014 -1.27048E-07
0.0015 -1.70989E-07
0.0016 -2.19828E-07

X m f(Xo) f´(Xo) Error Iteración

0.0016 0.00117137 -0.00000022 -0.00051286 0
0.001171371 0.00102259 -0.00000005 -0.00030261 0.14548898 1 CONTINUAR
0.001022595 0.00099885 -0.00000001 -0.00022937 0.02377151 2 CONTINUAR
0.000998851 0.00099821 0.00000000 -0.00021767 0.00063927 3 CONTINUAR
0.000998212 0.00099821 0.00000000 -0.00021736 0.00000046 4 FIN
0.000998212 0.00099821 0.00000000 -0.00021736 0.00000000 5 FIN
 EJERCICIOS CAPITULO III
Ejercicio 4.1: La capacidad calorífica a presión constante se define como:

Donde C, es la capacidad calorífica a presión constante, H Es la entalpia molar, y T es la

temperatura. La siguiente tabla muestra los datos de capacidad de calor frente a
temperatura para el dióxido de carbono
CALCULAR
A. Use las diferencias finitas para calcular la capacidad calorífica del dióxido de
carbono en cada una de las temperaturas dadas. Asegúrese de usar las fórmulas de
"punto final" para la primera y la última entrada.

1. Convertimos la temperatura de °C a °K
Para la entalpía convertiremos de KJ/mol a J/mol
Temperatura Entalpía Temperatura Entalpía
°C (Kj/mol) °K (J/mol)
100 25.186 373 25186
150 27.254 423 27254
200 29.408 473 29408
250 31.64 523 31640
300 33.942 573 33942
350 36.307 623 36307
400 38.732 673 38732
450 41.209 723 41209
500 43.736 773 43736
550 46.307 823 46307
600 48.919 873 48919
650 51.568 923 51568
700 54.252 973 54252
750 58.966 1023 58966
800 59.71 1073 59710
 d(Temperatura) d(entalpía)

1 40.5
1 42.22
1 43.86
1 45.34
1 46.67
1 47.9
1 49.02
1 50.04
1 50.98
1 51.83
1 52.61
1 53.33
1 73.98
1 54.58
1 14.88
1. DIFERENCIAL HACIA ADELANTE

⁡
−𝑌(𝑖 + 2) ⁡+ ⁡4⁡𝑌(𝑖 + 1) − 3𝑌𝑖
Yi =
2∆𝑋
Aplicamos la diferencial hacia adelante en la primera fila acá reemplazamos los valores
de la temperatura y la entalpía en la ecuación ajustada a la curva

2. DIFERENCIALES CENTRALES

𝑌(𝑖 + 1) − ⁡⁡𝑌(𝑖 − 1)
Yi =
2∆𝑋
Estas diferenciales la usamos en el medio de nuestros datos
3. DIFERENCIALES HACIA ATRÁS

𝑌(𝑖 ) − ⁡⁡𝑌(𝑖 − 1)
Yi =
∆𝑋
Estas diferenciales la usamos en las filas anteriores a la final y reemplazamos los datos.

B. Grafique los datos de entalpía y ajústelos a la curva con un polinomio adecuado.

Luego, calcule la capacidad de calor a cada temperatura usando analítica diferenciación.
Compara estos resultados con los de la parte a.
 ECUACIÓN AJUSTADA
y = 0.0135x2 + 30.956x + 11766
R² = 0.9958

Entalpía (J/mol) Valores de la Entalpía obtenidos después de reemplazar en la

25190.8295 ecuación ajustada de la gráfica.
27275.9295
29428.5295
31648.6295
33936.2295
36291.3295
38713.9295
41204.0295
43761.6295
46386.7295
49079.3295
51839.4295
52667.0295
57562.1295
60524.7295
 Cálculo de la capacidad de calor a cada temperatura
d(Temperatura) d(entalpía)
Cp(J/mol K)
1 40.5 40.5
1 42.22 42.22
1 43.86 43.86
1 45.34 45.34
1 46.67 46.67
1 47.9 47.9
1 49.02 49.02
1 50.04 50.04
1 50.98 50.98
1 51.83 51.83
1 52.61 52.61
1 53.33 53.33
1 73.98 73.98
1 54.58 54.58
1 14.88 14.88

2. Compare los resultados de Cp con los de la base de datos nist (estos valores
aparecen en la tabla que se encuentra en la parte inferior). Asegúrese de usar
unidades consistentes para la comparación
Base de Datos

Cp(J/mol K) Cp(J/molK)
40.5 40.461
42.22 42.256
43.86 43.881
45.34 45.355
46.67 46.695
47.9 47.917
49.02 49.034
50.04 50.055
50.98 50.989
51.83 51.843
52.61 52.624
53.33 53.339
73.98 53.994
54.58 54.593
14.88 55.144
Como observamos los valores obtenidos analíticamente a comparación de los
datos proporcionados por la tabla de base de datos son muy parecidos pues no
varían mucho.
Ejercicio 4.2: Evaluar la siguiente integral usando

a. La regla trapezoidal. Experimente con el incremento de Δx para producir buenos resultados.

n es número de divisiones entre mayor sea n más preciso será el resultado

Utilizamos n=2
𝑏−𝑎 𝑏
Δx = 𝑓(𝑏) + 𝑓(𝑎)
𝑛 ∫ 𝑓(𝑥) = (𝑏 − 𝑎)
𝑎 2

1−0
Δx= 2
f(xi)
Δx=0.5
1.2
Xi=a+i Δx 1
X0 0+0(0.5)=0 0.8
X1 0+1(0.5)=0.5 0.6
X2 0+2(0.5)=1 0.4
𝑏
𝑏−𝑎 0.2
∫ 𝑓(𝑥) ≈ [𝑓(𝑥0) + 2𝑓(𝑥1) + ⋯
𝑎 2𝑛 0
+ 2𝑓(𝑥𝑛) + 𝑓(𝑛)] 0 0.2 0.4 0.6 0.8 1 1.2
1
1−0
∫ √𝑥 ≈ [𝑓(0) + 2𝑓(0.5) + 𝑓(1)]
0 2(2)
1
1−0
∫ √𝑥 ≈ [0 + 1.414212 + 1]
0 2(2)
1 1
∫0 √𝑥 ≈ 4 [2.414212]
1
∫ √𝑥 ≈ 0.60553
0

Utilizamos n=3
𝑏−𝑎
Δx = 𝑛
1−0
Δx= 3

Δx=0.33
Xi=a+i Δx

𝑏
𝑏−𝑎
∫ 𝑓(𝑥) ≈ [𝑓(𝑥0) + 2𝑓(𝑥1) + ⋯ + 2𝑓(𝑥𝑛) + 𝑓(𝑛)]
𝑎 2𝑛
 1
1−0
∫ √𝑥 ≈ [𝑓(0) + 2𝑓(0.33) + 2𝑓(0.66) X0 0+0(0.33)=0
0 2(3)
+ 𝑓(0.99)] X1 0+1(0.33)=0.33
1
1−0 X2 0+2(0.33)=0.66
∫ √𝑥 ≈ [0 + 1.148912 + 1.624807
0 2(3)
+ 0.99498] X3 0+3(0.33)=0.99
1 1
∫0 √𝑥 ≈ 6 [3.76871]
1
∫ √𝑥 ≈ 0.628117
0

Utilizamos n=4
𝑏−𝑎
Δx = 𝑛

1−0
Δx= 4

Δx=0.25
Xi=a+i Δx

X0 0+0(0.25)=0
𝑏
𝑏−𝑎 X1 0+1(0.25)=0.25
∫ 𝑓(𝑥) ≈ [𝑓(𝑥0) + 2𝑓(𝑥1)
𝑎 2𝑛
+ ⋯ + 2𝑓(𝑥𝑛) X2 0+2(0.25)=0.5
+ 𝑓(𝑛)]
X3 0+3(0.25)=0.75
1
1−0
∫ √𝑥 ≈ [𝑓(0) + 2𝑓(0. .25) X4 0+4(0.25)=1
0 2(4)
+ 2𝑓(0.5)
+ 2𝑓(0.75) + 𝑓(1)]
1
1−0
∫ √𝑥 ≈ [0 + 1 + 1.414213 + 1.73205 + 1]
0 2(4)
1 1
∫0 √𝑥 ≈ 8 [5.146263]
1
∫ √𝑥 ≈ 0.643283
0