Primer Trabajo Metodos

U N A C Facultad de Ingeniera Qumica
Universidad nacional del callao
FACULTAD: Ingeniera Qumica
MTODOS NUMRICOS
Profesor:
Mg. Ronald Portales
Integrantes:
Antezana Vergara, Miguel ngel.
Huanca Ortiz, Elvis.
Nuez Zumaeta, Mercedes.
Pillaca Quispe, Elizabeth.
Ponte Valverde, Erick.
Vidal Villarreal, Edward.
Villafuertes Pelayo, Joselyn.
Ciclo:
2016-V
Bellavista Callao, 18 de enero del 2016
1
Mtodos Numricos
TRABAJO
N 1
GRUPO 5
2
Mtodos Numricos
I. DETERMINAR LAS RACES DE LOS SIGUIENTES EJERCICIOS:

5.1. Determine las races reales de f(x) = -0.5x2 + 2.5x + 4.5:
a) Grficamente
De la tabla y de la grfica, observamos que existen 2

races, una se encuentra entre los valores de x de -2 y
-1 y la otra entre los valores de 6 y 7.
b) Utilizando el mtodo de la biseccin con tres iteraciones para determinar la raz

ms grande. Emplee como valores iniciales xi = 5 y xu = 10. Calcule el error estimado
Ea y el error verdadero Et para cada iteracin.
De acuerdo a la tabla y la grfica de la parte (a), la raz ms grande se encuentra entre
los valores de x de 6 y 7.
Para ste mtodo, usaremos:
1. Primera iteracin:
x a=5
f
3
Mtodos Numricos
x b=10
x a+ x b 15
= =7,5
2
2
x r=
x
(NEGATIVO) Nuevo intervalo: Inferior o izquierdo.
r
2. Segunda iteracin:
x a=5
4
Mtodos Numricos
x b=7.5
x a+ x b
=6.25
2
x r=
x
(POSITIVO) Nuevo intervalo: Superior o derecho
r
xr
nuevo
x r
anterior
100 =20
|6.257.5
6.25 |
relativo =
absoluto =|x ractual x ranterior|=|6.257.5|=1.25
xr
nuevo
100 =
3. Tercera iteracin:
x a=6.25
5
Mtodos Numricos
x b=7.5
x a+ x b
=6.875
2
x r=
x
(NEGATIVO) Nuevo intervalo: Inferior o izquierdo.
r
relativo =
absoluto =|x ractual x ranterior|=|6.8756.25|=0.625
x rnuevo x ranterior
xr
nuevo
100 =9.0909
|6.8756.25
6.875 |
100 =
En Excel:
6
Mtodos Numricos
Raz= 6.40563965
c) Utilizando el mtodo de la secante con tres iteraciones para determinar la raz ms

grande. Emplee como valores iniciales x i = 5 y xu = 10. Calcule el error estimado E a y
el error verdadero Et para cada iteracin.
x 0=5
x 1=10
De acuerdo a la frmula:
Mtodos Numricos
x
x
x
( i1)f (i)
f
f ( i)(x x )
x =x
i1
Primera iteracin:
x
x
x
( 0)f (1 )
20.5(510)
10
=
f
4.5(20.5) = 5.9
f ( 1)(x x )
x 2=x 1
0
relativo =
absoluto =|x actual x anterior|=|5.910|=4.1
x nuevo xanterior
5.910
100 =
100 =69.491525
nuevo
5.9
x
Segunda iteracin:
x
x
x
( 1)f (2)
f
f (2 )(x x )
x 3=x 2
1
1.845(105.9)
= 5.9 20.5(1.845) = 6.2385
8
Mtodos Numricos
relativo =
absoluto =|x actual x anterior|=|6.23855.9|=0.33858
x nuevo xanterior
6.23855.9
100 =
100 =5.42598
nuevo
6.2385
x
Tercera iteracin:
x
x
x
( 2)f (3)
f
f (3 )(x x )
x 4=x 3
2
relativo =
0.6368(5.96.2385)
= 6.2385 1.845(0.6368)
= 6.4169
x nuevo xanterior
6.41696.2385
100 =
100 =2.780158
nuevo
6.4169
x
9
Mtodos Numricos
absoluto =|x actual x anterior|=|6.41696.2385|=0.1784
En Excel:
Raz= 6.405124451
5.2. Determine las races reales de f(x) = 5x3-5x2+6x-2:
a. Grficamente
10
Mtodos Numricos
3500
3000
delta de x
x
1
f(x)
0
1
2
3
4
5
6
7
8
9
2500
2000
1500
-2
4
30
106
262
528
934
1510
2286
3292
1000
500
0
0
-500
10
De la tabla y de la grfica, observamos que las tres

races se encontrarn entre 0 y 1.
b. Utilizando mtodo de la biseccin para localizar la raz ms pequea. Use los

valores iniciales xi = 0 y xu = 1. Calcule el error estimado Ea y el error
verdadero Et para cada iteracin.
De acuerdo a la tabla y la grfica de la parte (a), la raz se encuentra entre los
valores de x de 0 y 1.
Entonces, haciendo los respectivos clculos, se tendr los siguientes resultados:
11
Mtodos Numricos
12
Mtodos Numricos
f ( x ) =5 x 3 - 5 x 2 +6 x - 2
x0 =0
x r =0.5
x1 =1
f ( x0 ) =- 2
f ( xr ) =5(0.5)3 - 5(0.5) 2 +6(0.5) - 2 =0.375
f ( xr ) f ( x0 ) =(0.375) ( - 2) =- 0.75
luego :
x0 =0
x r =0.25
error = 0.25 - 0.5 =0.25
x1 =0.5
f ( xr ) =5(0.25)3 - 5(0.25) 2 +6(0.25) - 2 =- 0.734375
f ( xr ) f ( x0 ) =( - 0.734375) (- 2) =1.46875
luego :
x0 =0.25
x r =0.375
error = 0.375 - 0.25 =0.125
x1 =0.5
f ( xr ) =5(0.375)3 - 5(0.375) 2 +6(0.375) - 2 =- 0.189453
f ( x0 ) =- 0.734375
f ( xr ) f ( x0 ) =( - 0.189453) ( - 0.734375) =0.139129
luego :
luego :
x0 =0.375
x r =0.4375
error = 0.4375 - 0.375 =0.0625
x1 =0.5
f ( xr ) =5(0.4375) 3 - 5(0.4375) 2 +6(0.4375) - 2 =0.0866699
f ( x0 ) =- 0.1894531
f ( xr ) f ( x0 ) =(0.0866699) ( - 0.1894531) =- 0.01641
luego :
x0 =0.375
x r =0.40625
error = 0.40625 - 0.4375 =0.015625
x1 =0.4375
f ( xr ) =5(0.40625)3 - 5(0.40625) 2 +6( 0.40625) - 2 =- 0.052459
f ( x0 ) =- 0.1894531
f ( xr ) f ( x0 ) =( - 0.052459) ( - 0.1894531) =0.009938
luego :
x0 =0.40625
x r =0.421875
error = 0.421875 - 0.40605 =0.0078125
x1 =0.4375
f ( xr ) =5(0.421875)3 - 5(0.421875) 2 +6(0.421875) - 2 =0.016780853
f ( x0 ) =- 0.052459
f ( xr ) f ( x0 ) =(0.016780853) ( - 0.052459) =- 0.00088028
Resolviendo en Excel
13
Mtodos Numricos

i
X1
X2
Xr
0
1
0
0
1
0.5
0.25
0.5
0.375
0.5
0.375
0.4375
0.40625
0.4375
0.40625
0.421875
0.414062
5
0.421875
0.417968
75
0.421875
0.417968
75
0.419921
88
10
0.417968
75
0.418945
31
11
0.417968
75
0.418457
03
0.417968
75
0.418090
82
0.418090
82
0.418090
82
0.418090
82
0.418098
45
0.418098
45
0.418100
36
0.418100
36
0.418212
89
0.418212
89
0.418151
86
0.418121
34
0.418106
08
0.418106
08
0.418102
26
0.418102
26
0.418101
31
12
13
14
15
16
17
18
19
20
f(X1)
0.5
0.25
f(Xr)
-2
-2
f(X1)*f(Xr) error ABS

0.375
-0.75
-0.734375
1.46875
0.25
0.189453 0.139129
13
64
0.125
0.086669 0.016419
92
89
0.0625
0.052459 0.009938
72
66
0.03125
0.016780 0.000880
85
32 0.015625
0.017913 0.000939 0.007812
34
73
5
0.000585
0.003906
62 1.049E-05
25
0.375 -0.734375
0.189453
0.4375
13
0.189453
0.40625
13
0.052459
0.421875
72
0.414062 0.052459
5
72
0.417968 0.017913
75
34
0.419921 0.000585 0.008092 -4.7392E88
62
66
06
0.418945 0.000585 0.003752 -2.1974E31
62
3
06
0.418457 0.000585 0.001583 -9.2705E03
62
04
07
0.418212 0.000585 0.000498 -2.9201E89
62
64
07
0.418090 0.000585 -4.3509E82
62
05 2.548E-08
0.418151 -4.3509E- 0.000227 -9.9009E86
05
56
09
0.418121 -4.3509E- 9.2023E- -4.0039E34
05
05
09
0.418106 -4.3509E- 2.4257E- -1.0554E08
05
05
09
0.418098 -4.3509E- -9.6263E- 4.1883E45
05
06
10
0.418102 -9.6263E- 7.3152E- -7.0418E26
06
06
11
0.418100 -9.6263E- -1.1556E- 1.1124E36
06
06
11
0.418101 -1.1556E- 3.0798E- -3.5589E31
06
06
12
0.418100 -1.1556E- 9.6212E- -1.1118E83
06
07
12
14
Mtodos Numricos
0.001953
13
0.000976
56
0.000488
28
0.000244
14
0.000122
07
6.1035E05
3.0518E05
1.5259E05
7.6294E06
3.8147E06
1.9073E06
9.5367E07
4.7684E07
Raz = 0.418100
c. Mediante el mtodo de la secante
De acuerdo a la frmula:
x
x
x
( i1)f (i)
f
f ( i)(x x )
x i +1=x i
i1
Aplicando la frmula anterior, se tiene:
(4) (1 - 0)
4 - (- 2)
x 2 =0.33333
x0 =0 f ( x0 ) =- 2
x 2 =1 -
x1 =1 f ( x1 ) =4
luego :
x1 =1 f ( x1 ) =4
x 3 =0.33333 -
x2 =0.33333 f ( x2 ) =- 0.37037
(- 0.37037) (0.33333 - 1)
- 0.37037 - (4)
x 3 =0.38983
luego :
x2 =0.33333 f ( x2 ) =- 0.37037
x3 =0.38983 f ( x3 ) =- 0.124648
(- 0.124648) (0.38983 - 0.33333)

- 0.124648 - (- 0.37037)
x 4 =0.41848
x 4 =0.38983 -
luego :
x3 =0.38983 f ( x3 ) =- 0.124648
x4 =0.41848 f ( x4 ) =- 0.0017285
(- 0.0017285) (0.41848 - 0.38983)

- 0.0017285 - (- 0.124648)
x 5 =0.41809
x 5 =0.41848 -
luego :
x4 =0.41848 f ( x4 ) =- 0.0017285
x 6 =0.41809 -
x5 =0.41809 f ( x5 ) =- 0.0000125
x 6 =0.41810
(- 0.0000125) (0.41809 - 0.41848)

- 0.0000125 - (- 0.0017285)
Mtodo de la secante utilizando Excel
15
Mtodos Numricos
2
3
4
5
6
1
0.333333 0.6666666
33
7
0.389830 0.0564971
51
8
0.418489 0.0286592
79
8
0.418097
79 -0.000392
0.418100
2.8243E62
06
4
-2
0.370370
37
4
0.124647
6 -0.37037037
0.001728
54 -0.1246476
-1.2544E05 0.00172854
-1.3997E09 -1.2544E-05
CONTINU
1 A
0.666666 CONTINU
67 A
0.056497
18
0.028659
28
CONTINU
A
CONTINU
A
CONTINU
0.000392 A
2.8243E06 FIN
Raz =0.41810062
5.3. Determine las races reales de f(x) = -25182x-90x2+44x3-8x4+0.7x5.

a) Grficamente
Observando se nota que una

y 1, y otra de ellas entre x=16
races por los diferentes
b) Usando el mtodo de
X=15
x
f(x)
-10
48820
-9
93449,7
-8
117462,4
-7
125799,1
-6
122536,8
-5
110972,5
-4
93707,2
-3
72729,9
-2
49501,6
-1
25039,3
0
0
1
-25235,3
2
-50477,6
3
-75645,9
4
-100683,2
5
-125472,5
6
-149752,8
-122917,5
7
-173035,1
X=16
-148099.5
X=17
87797,9
X=18
257061,6
raz se encuentra entre 0

y x=17. Hallaremos tales
mtodos ya estudiados.
Newton Raphson.
Existe
raz
16
Mtodos Numricos
POR EL METODO DE NEWTON-RAPHSON, NOS ACERCAMOS POR 16 (RAIZ

SUPERIOR)
1RA ITERACION:
X17 = X16 -
f ( X 16)
f ( X 16)
X17 = 16 -
148099.5
104034
X17 = 17.4235
2DA ITERACION:
X18 = X17 -
f ( X 17)
f ( X 17)
X18 = 17.4235 -
153396.9175
165052.7509
16.4941117.4235
16.49411
= 0.0563
X18 = 16.49411
3RA ITERACION:
X19 = X18 -
f ( X 18)
f ( X 18)
20047,7689
123216,2901
X18 = 16.49411 -
16.331416.49411
16.3314
= 0.009963
X18 = 16.3314
17
Mtodos Numricos
4TA ITERACION:
X20 = X19 -
f ( X 19)
f ( X 19)
X20 = 16.3314 -
534.2947
116676.7558
16.3268216.3314
16.32682
= 0.00028 = 0,0001 aprox
X20 = 16.32682
5TA ITERACION:
X21 = X20 -
f ( X 20)
f ( X 20)
X21 = 16.32682 -
2.004996
116495.1052
16.32681
16.3268116.32682
= 0.000000612
X21 = 16.32681
En Excel:
i
F(x)
F(x)
16
-36012,8
104034
16,34616375
2261,19192
117260,957
16,32688034
7,35852544
116498,276
16,32681717
7,8738E05
%Er
Ea
16,9264677
0,34616375
-306,288727 -0,01928342
116495,7
-6,3164E93454,31
83
05
83
raz
c) Usando el mtodo de la secante.

18
Mtodos Numricos
Sabiendo que una de las races es 0, simplificamos x y la ecuacin queda asi:

F(x) = -25182 90x+44x2-8x3+0.7x4
Nos dan como valores iniciales:
X0 = 0.5
X1 = 1
Por iteraciones y hasta tener un error absoluto de 0.1
xn+1 = xn -
X nX n1
F ( X n )F ( X n1)
. F(xn)
x0 = 0.5
x1 = 1
x2 = x1 -
X 1 X 0
F ( X 1 )F( 0)
. F(x1) = 1 -
10.5
F (1 )F( 0.5)
. F( 1
X2 = -686.8446
Error = 1.001
x3 = x2 -
X 2 X 1
F ( X 2 )F( X 1)
686.84461
F (686.8446 )F(1)
. F(x2) = -686.8446-
. F(-686.8446)
x3 = 0.99989
Error =687.9199
X 3X 2
x4 = x3 - F ( X 3 )F ( X 2) . F(x3) = 0.999890.99989686.8446
F ( 0.99989 )F (686.8446)
. F(0.99989)
x4 = 0.999781
19
Mtodos Numricos
xn
1
2
3
xn-1
F( X n1)
1
0.5
-686.8446
1
0.99781
-686.8446
Error =1.09607
-25216.95625
-25235.3
9.22337
F ( X n)
-25235.3
9.22337
-25235.2975
error
1.001456
687.9199
1.0961
En Excel:
Mediante iteraciones se obtiene que el valor de X es
0.99978083423
d) Usando el mtodo de la Biseccin:
Usando el mtodo de la biseccin para localizar la raz ms grande con
= 10
% Utilic como valores iniciales X1 = 0.5 y X2 = 1.0,

SOLUCION:
1ra Iteracin:
Xa
= 0.5
Xr
Xb
= 1.0
0.5+ 1.0
2
= 0.75
Evaluando en la funcin el valor de
Xr
= 0.75
3
f (0.75)=2518290 ( 0.75 ) +44 (0.75) 8( 0.75) + 0.7(0.75)
f(0.75)
= -25216.42764
Multiplicamos nuestro
f (0.75)
f (0.5)
Xa
= 0.5 y nuestro
Xr
= 0.75 evaluados en la funcion
= (-25216.42764) ( -25199.47813 ) = 635440816.7
Como nos sale positivo (+) tomamos para la siguiente iteracin el rango de
<
Xr
= 0.75 ;
Xb
= 1.0 >
20
Mtodos Numricos
2da Iteracin:
Xa
= 0.75
Xr
Xb
0.75+1.0
2
= 1.0
= 0.875
Xr
= 0.875
3
f (0.875)=2518290 ( 0.875 ) +44 (0.875) 8(0.875) + 0.7(0.875)

f(0.875) = -25225.7601
funcion
f (0.875)
f (0.75)
Xa
= 0.75
y nuestro
Xr
= 0.875 evaluados en la
= (-25225.7601) ( -25216.42764) = 636103553.3
<
Xr
= 0.875;
Xb
= 1.0 >
3ra Iteracin:
Xa
Xr
= 0.875
Xb
0.875+ 1.0
2
= 1.0
= 0.9375

2
Xr
= 0.9375
3
f (0.9375)=2518290 ( 0.9375 ) +44 (0.9375) 8(0.9375) + 0.7(0.9375)

f(0.9375) = -25230.51955
21
Mtodos Numricos
Xa
funcion
f (0.9375)
f (0.875)
= 0.875
y nuestro
Xr
= 0.9375 evaluados en la
= (-25230.5195) (-25225.76) = 636459033.4
<
Xr
= 0.9375 ;
Xb
= 1.0 >
Y seguimos iterando hasta obtener un error menor al 1 %
En Excel:
#iter
Xa
Xb
Xr
f(Xa)
f(Xr)
f(Xa)*f(Xb)
interv
Ea
(%)
0.5
0.75
-25199.47813
-25216.427
635440816
Super.
0.75
0.875
-25216.427
-25225.7601
4.4767
Super.
14.28
0.875
0.9375
-25225.7601
-25230.5195
636459033
Super.
6.67
0.9375
0.96875
-25230.5195
-25232.9089
636639403
Super.
3.22
0.96875
0.98437
-25232.9089
-25234.1045
636729862
Super.
1.58
0.98437
0.9921875
-25234.1045
-25234.70
636775115
Super.
0.787
5.4. Calcule las races reales de f(x) = -12 - 21x + 18x 2 2.75x3. Adems, determine
la primera raz de la funcin con los mtodos de biseccin, Newton-Raphson y
secante. Utilice los valores iniciales de x0=-1 y x1=0 y un criterio de detencin de
1%.
a) Grficamente
22
Mtodos Numricos
Grficamente las races de la ecuacin son x1= -0.414689,

x2=2.2198183 y x3=4.7403257.
Hallaremos la primera raz (x1= -0.414689) con los diversos
mtodos.
b) Usando el mtodo de la biseccin.

1. Primera iteracin:
x a=1
x
x b=0
23
Mtodos Numricos
x r=
x a+ x b
=0.5
2
f
(POSITIVO) Nuevo intervalo: superior o derecho.
2. Segunda iteracin:
x a=0.5
24
Mtodos Numricos
x b=0
x r=
x a+ x b
=0.25
2
f
(NEGATIVO) Nuevo intervalo: inferior o izquierdo.
25
Mtodos Numricos
xr
nuevo
x r
anterior
relativo =
absoluto =|x ractual x ranterior|=|0.25(0.5)|=0.25
xr
nuevo
100 =
0.25(0.5)
100 =100
0.25
3. Tercera iteracin:
x a=0.5
x
x b=0.25
x
x r=
x a+ x b
=0.375
2
26
Mtodos Numricos
f
relativo =
xr
nuevo
100 =
0.375(0.25)
100 =33.3
0.375
4. Cuarta iteracin:
x a=0.5
x
x b=0.375
x
27
Mtodos Numricos
x r=
x a+ x b
=0.4375
2
f
(POSITIVO) Nuevo intervalo: superior o derecho.
| |
xr
nuevo
100 =
0.4375(0.375)
100 =14.3
0.4375
relativo =
5. Quinta iteracin:
x a=0.4375
x
f
28
Mtodos Numricos
x b=0.375
x r=
x a+ x b
=0.40625
2
relativo =
absoluto =|x r
xr
actual
nuevo
x r
| |
100 =
0.40625(0.4375)
100 =7.6
0.40625
|=|0.40625(0.4375)|=0.03125
anterior
29
Mtodos Numricos
6. Sexta iteracin:
x a=0.4375
x
x b=0.40625
x
x r=
x a+ x b
=0.42187
2
30
Mtodos Numricos
x
(NEGATIVO) Nuevo intervalo: superior o derecho.
r
relativo =
absoluto =|x r
xr
actual
nuevo
x r
100 =
0.42187(0.40625)
100 =3.7
0.42187
|=|0.42187(0.40625)|=0.01562
anterior
En Excel
Xa
Xb
Xr
F(Xa)
F(Xb)
F(Xr)
F(Xa)*F(Xr)
DECIDIR
-1
-0,5
29,75
-12
3,34375
99,4765625
SUPERIOR
-0,5
-0,25
3,34375
-12
-5,58203125
-18,66491699
INFERIOR
100
-0,5
-0,25
-0,375
3,34375
-1,44873047
-4,844192505
INFERIOR
33,33333
-0,5
-0,375
-0,4375
3,34375
0,86309814
2,885984421
SUPERIOR
14,28571
0,
-0,4375
-0,375
-0,40625
-0,3136673
-0,270725662
INFERIOR
7,692307
0,0
-0,4375
-0,40625
-0,421875
-0,3136673
0,26947117
0,232580066
SUPERIOR
3,70370
0,0
-0,421875
-0,40625
-0,4140625
-0,3136673
-0,02340519
-0,006307025
INFERIOR
1,886792
0,0
-0,421875
-0,4140625
-0,41796875
0,12270571
0,033065652
SUPERIOR
0,934579
0,0
0,41796875
-0,4140625
-0,41601563
0,0495685
0,006082338
SUPERIOR
0,469483
0,0
0,8630981
4
0,8630981
4
0,2694711
7
0,2694711
7
0,1227057
1
5,58203125
1,44873047
1,44873047
0,02340519
0,02340519
31
Mtodos Numricos
Er (%)

0,41601563
0,41503906
0,41503906
0,41479492
-0,4140625
-0,41503906
-0,4140625
-0,41455078
-0,41455078
-0,41479492
-0,41455078
-0,41467285
0,0495685
0,0130612
2
0,0130611
3
0,0039407
2
0,02340519
0,02340519
0,00517714
0,00517714
0,01306122
0,000647425
SUPERIOR
0,235294
0,0
-0,00517709
-6,76192E-05
INFERIOR
0,117785
0,0
0,00394072
5,14702E-05
SUPERIOR
0,05885
0,0
-0,00061853
-2,43745E-06
INFERIOR
0,029437
0,0
Raz:-0,41467285
c) Usando el mtodo de Newton Raphson.

Usando la siguiente frmula:
x
x
f ( i)
f ( i)
f (x) =
x
=x
3
i
-12-21x+18x2-2.75xi+1
f (x) =
-21+36x-8.25x2
1. Primera iteracin: x1 = 0
x
x
12
f (1)=0
=0.57142857
21
f ( 1)
x 2=x 1
x nuevo xanterior
0.571428570
100 =
100 =100
nuevo
0.57142857
x
relativo =
absoluto =|x actual x anterior|=|0.571428570|=0.57142857

32
Mtodos Numricos
2. Segunda iteracin: x2 = -0.57142857

x
x
6.39067
f (2)=0.57142857
=0.427055
44.265
f ( 2)
x 3 =x2
relativo =
0.427055(0.57142857 )
x nuevo xanterior
100 =
100 =33.80
nuevo
0.427055
x
absoluto =|x actual x anterior|=|0.427055(0.57142857)|=0.14437
3. Tercera iteracin: x3 = -0.427055

x
x
0.465105
f (3)=0.427055
=0.414776
37.87858
f (3 )
x 4= x3
relativo =
0.414746(0.427055 )
x nuevo xanterior
100 =
100 =2.96073
nuevo
0.414746
x
absoluto =|x actual x anterior|=|0.414776(0.427055 )|=0.01228
4. Cuarta iteracin: x4 = -0.414776
33
Mtodos Numricos
x
x
0.0032408
f (4 )=0.4147762
=0.414689
37.3513
f ( 4)
x 5=x 4
relativo =
0.414689(0.414746 )
x nuevo xanterior
100 =
100 =0.029
nuevo
0.414689
x
absoluto =|x actual x anterior|=|0.414689(0.4147762)|=0.000087
En Excel:
o 1era raz= -0.414689416
d) Usando el mtodo de la secante.

1) Haciendo uso del Mtodo de la Secante :
( X X k 1 ) f ( X k )
X k 1 X k k
f ( X k ) f ( X k 1 )
x0 1
x1 0
y
34
Mtodos Numricos
2) Iterando
a) K=1
(1) (12)
0.28742514
(12) (29.75)
Error | 0.28742514 0 | 0.28742514
X2 0
b) K=2
(0.28742514 0) (4.411734887)
0.4545310087
( 4.411734887) ( 12)
Error | 0.4545310087 0.2874251497 | 0.167105859
X 3 0.28742514
c) K=3
(0.4545310087 0.2874251497) (1.005683224)
1.036704562
(1.005683224) (4.411734887)
Error |1.036704562 0.4545310087 | 1.491235571
X 4 0.4545310087
d) K=4
(1.036704562 0.4545310087) (17.48924475)

0.3734433398
(17.48924475) (1.005683224)
Error | 0.3734433398 1.036704562 | 1.410147902
X 5 1.036704562
e) K=5
(0.37344333398 1.036704562) (1.504190111)

0.5061379459
(1.504190111) (17.48924476)
Error | 0.5061379459 0.3734433398 | 0.1326946061
X 6 0.37344333398
f) K=6
( 0.5061379459 0.3734433398) (3.596623585)

0.4125739495
(3.596623585) (1.504190111)
Error | 0.4125739495 0.5061379459 | 0.09356400338
X 7 0.5061379459
g) K=7
X 8 0.4125739495
(0.4125739495 0.5061379459) (0.4125739495)

0.4145827112
(0.4125739495) (3.596623585)
Error | 0.4145827112 0.4125739495 | 2.008761721*10 3
35
Mtodos Numricos
Y as sucesivamente hasta que k=12 donde se tiene el valor de la raz ms

prxima a 0 teniendo en cuenta E= 0.001
En Excel:
Iteracin
Xi
Xi - Xi-1
f(xi)
f(xi-1)
Ea
Observacion
-17,75
-12
Continuar
-2,08695652
-3,08695652
135,219199
-17,75
3,08695652
Continuar
0,64180058
2,7287571
-18,7904661
135,219199
2,7287571
Continuar
0,30886941
-0,33293117
-16,8500842
-18,7904661
0,33293117
Continuar
-2,58227181
-2,89114122
209,606032
-16,8500842
2,89114122
Continuar
0,09374618
2,67601799
-13,8127451
209,606032
2,67601799
Continuar
-0,07169719
-0,16544337
-10,4008167
-13,8127451
0,16544337
Continuar
-0,57602961
-0,50433242
6,59481812
-10,4008167
0,50433242
Continuar
-0,38033344
0,19569617
-1,25793856
6,59481812
0,19569617
Continuar
10
-0,41168215
-0,03134871
-0,11212021
-1,25793856
0,03134871
Continuar
11
-0,41474967
-0,00306752
0,00225066
-0,11212021
0,00306752
Continuar
-0,41468931
6,0365E-05
-3,887E-06
0,00225066
6,0365E-05
Fin
12
1era raz= -0.41468931
II. EJERCICIOS DE INGENIERA QUMICA:

1. La concentracin de saturacin de oxgeno disuelto en agua dulce se calcula con
la ecuacin (APHA 1992)
1.575701 x 105 6.642308 x 10 7 1.243800 x 1010 8.621949 x 10 11
ln =139.34411+
Ta
T a2
T a3
T a4
Donde
36
Mtodos Numricos
:concentracion de saturacion de oxigenodisuelto en agua dulce a 1 atm
( mgL )
T a=Temperatura absoluta ( K ) . Recuerde el lector que T a=T + 273.15
Donde T= Temperatura ( ). De acuerdo con esta ecuacin, la saturacin

disminuye con el incremento de la temperatura. Para aguas naturales comunes en
climas templados, la ecuacin se usa para determinar que la concentracin de
oxigeno varia de 14.621 mg/L a 0 a 6.413 mg/L a 40 . Dado un valor
de concentracin de oxgeno, puede emplearse esta frmula y el mtodo de
biseccin para resolver para la temperatura en .
SOLUCIN POR EL MTODO DE LA SECANTE
Nos dan las temperaturas iniciales:
T1
= 0 + 273.15 = 273.15 K
T2
= 40 + 273.15 = 313.15 K
= 14.621 mg/L
= 6.413 mg/ L
Por iteraciones y hasta tener un error absoluto de 0.05

xn+1 = xn -
X nX n1
F ( X n )F ( X n1)
. F(xn)
x0 = 273.15K
x1 = 313.15K
x2 = x1 -
X 1 X 0
F ( X 1 )F( 0)
. F(x1) = 313.15 -
313.15273.15
F ( 313.15 )F (273.15)
. F(
313.15
X2 = 416.5086
Error = 0.10311446673
37
Mtodos Numricos
x3
x2
X 2 X 1
F ( X 2 )F( X 1)
416.5086313.15
F ( 416.5086 )F (313.15)
F(x2)
416.5086
. F(416.5086)
x3 = 401.4785
Error = 0.03629667967
n
xn
xn-1
2
3
416.5086
401.4786
313.15
416.5086
F( X n1)
2.68244
0.805309
F ( X n)
0.80531
-0.41426
error
0.103114
0.036297
Mediante iteraciones se obtiene que el valor de la temperaturA es

401.4786K.
En grado Celsius se obtiene 128.3285C
SOLUCIN POR EL MTODO DE LA BISECCIN
a) Si los valores iniciales son de 0 y 40
, con el mtodo de la
biseccin, Cuantas iteraciones se requeriran para determinar la

temperatura con un error absoluto de 0.05
Solucin :
Nos dan las temperaturas iniciales:
T1
= 273 K
= 14.621 mg/L
T2
= 313 K
= 6.413 mg/ L
Aplicamos el mtodo de la biseccin :
38
Mtodos Numricos
1ra Iteracin:
Xa
= 273 K
Xr
Xb
273 K +313 K
2
= 313 K
= 293 K

f (293)
Xr
= 293 K
139.34411+
1.575701 x 105 6.642308 x 107 1.243800 x 1010 8.621949 x 1011
293
2932
293 3
293 4
funcion
f (273)
f (293)
Xa
= 273 K y nuestro
Xr
= 293 K evaluados en la
= (14.621) (2.2104) = 32.3182
Xr
Xb
<
= 293 K ;
= 313 K >
2da Iteracin:
Xa
= 293 K
Xb
= 313 K
39
Mtodos Numricos
Xr
293 K +313 K
2
= 303 K

f (303)
Xr
= 303 K
139.34411+
1.575701 x 105 6.642308 x 107 1.243800 x 1010 8.621949 x 1011
303
3032
303 3
303 4
funcion
f (293)
f (303)
Xa
= 293 K y nuestro
Xr
= (2.2104) (2.025305) = 4.4767
Xr
Xb
<
= 303 K ;
= 313 K >
3ra Iteracin:
Xa
Xr
Xb
= 303 K
303 K +313 K
2
= 313 K
= 308 K

f (308)
= 308 K
139.34411+
1.575701 x 105 6.642308 x 107 1.243800 x 1010 8.621949 x 1011
2
3
4
308
308
308
308
funcin:
f (303)
Xr
f (308)
Xa
= 303 K y nuestro
Xr
= (2.025305) (1.941104) = 3.93132
40
Mtodos Numricos
f
f (X )
Xr
Xb
<
= 308 K ;
= 313 K > ( X )
b
Y as consecutivamente hasta encontrar un error absoluto de 0.05

A continuacin mostraremos una tabla para apreciar mejor esto: En Excel
2. Mediante el Mtodo de Newton Raphson (R=0.082), determinar el

volumen molar.
Determinar el volumen molar del oxigeno mediante la ecuacin de Van Der
Waals:
P+
a
( V b ) =RT , P=100 atm ,T =700 K
V2
paraun gas que tiene :a=1.36 y b=0.0318

Por la ecuacin del gas ideal:
RT 0.082700
A 100 atm ,V m =
=
=0.574 dm 3mol1
P
100
La forma desarrollada de la ecuacin de Van Der Waals es la siguiente:

41
Mtodos Numricos
V m3 b+
RT
a
ab
V m2 + V m =0
P
P
P
Reemplazando datos anteriores se tiene:
V m3 0.0318+
0.082700
1.36
1.360.0318
V m 2+
V m
=0
100
100
100
Quedando:
3
V m 0.6058 V m +0.0136 V m 0.00043248=0

Derivando la ecuacin anterior:
'
2
f ( V )=3 V m 1.2116 V m+ 0.0136
m
Por el mtodo de Newton Raphson:

f (V )
V m (i +1)=V m , i '
f (V )
m ,i
m ,i
Luego:
V m , 0=0.574 dm mol
V m ,2 =0.574
V m ,3 =0.5841
V m , 4=0.5837
f ( 0.574)
f
'
=0.5841
( 0.574 )
f (0.5841 )
f ' ( 0.5841)
f (0.5837 )
f ' (0.5837 )
=0.5837
=0.5837
V m para eloxgeno es 0.5837 dm3mol1
Xi
F(Xi)
0
0,574
0,00310342
1
0,58412304 0,00011542
2
0,58377272 1,4067E-07
Mtodos Numricos
3
0,5837723 2,0983E-13
4
0,5837723 7,0473E-19
5
0,5837723 7,0473E-19
F(Xi)
Er
Ea
0,3065696
0,32947571 1,7330324 0,01012304
0,32867275 0,0600092 0,00035032
0,32867177 7,3313E-05 4,2798E-07
0,32867177 1,0935E-10 6,3838E-13
0,32867177
0
0
42
Mediante Excel:
3. La ecuacin de estado Redlich-Kwong es:
P+
a
(V b ) =RT
TV (V + b )
Donde :
a=17.19344 y b=0.0221141 para el oxgeno molecular .
Si T =373 K y P=30 atm
a. Mtodo grfico:
delta x
0.2
v
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
f(v)
-0.00067644
-0.03230547
-0.09749716
-0.14825152
-0.13656855
-0.01444824
0.2661094
0.75310438
1.49453668
2.53840632
43
Mtodos Numricos

3
2.5
2
1.5
1
0.5
0
-0.5
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2
Las races se encontrarn entre 1 y 1.2, como se

muestra en la tabla y la grfica.
b. Mtodo de la Secante
Sea la ecuacin:
T =373 K y P=30 atm a=17.19344 y b=0.0221141
Mtodo de la secante por medio de clculos
44
Mtodos Numricos
x0 =1 f ( x0 ) =- 0.01444824
x1 =1.2 f ( x1 ) =0.2661094
(0.2661094) (1.2 - 1)
0.2661094 - (- 0.01444824)
x 2 =1.01029966
x 2 =1.2 -
luego :
x1 =1.2 f ( x1 ) =0.2661094
(- 0.00428043 ) (1.01029966 - 1.2)

- 0.00428043 - (0.2661094)
x 3 =1.01330272
x 3 =1.01029966 -
x2 =1.01029966 f ( x2 ) =- 0.00428043
luego :
x2 =1.01029966 f ( x2 ) =- 0.00428043
x 4 =1.01330272 -
x3 =1.01330272 f ( x3 ) =- 0.00123574
x 4 =1.01452157
(- 0.00123574) (1.01330272 - 1.01029966)

- 0.00123574 - (- 0.00428043)
luego :
(1.0387 10- 5 ) (1.01452157 - 1.01330272)
1.0387 10- 5 - (- 0.00123574)
x3 =1.01330272 f ( x3 ) =- 0.00123574
x 5 =0.41848 -
x4 =1.01452157 f ( x4 ) =1.0387 10- 5
x 5 =1.01451141
En Excel:
1.2
0.2
0.2661094 -0.01444824
1.01029966 -0.18970034 -0.00428043
0.2 continua
0.2661094
0.18970034 continua
1.01330272
0.00300307 -0.00123574 -0.00428043
0.00300307 continua
1.01452157
0.00121885
0.00121885 continua
1.01451141
-1.016E-05 -2.4841E-08
1.0387E-05 -0.00123574
1.0387E-05
1.016E-05 fin
Raz = 1.01451141
45
Mtodos Numricos
4. Resolver:
En un proyecto de ingeniera qumica se requiere que se determine
exactamente el volumen molar y factor de comprensibilidad del amoniaco a
una presin de 120atm y 500K mediante La siguiente ecuacin.
DATOS:
Tc = 405.5K
Pc = 111.3atm
R = 0.082
( P+ Va )( V b) =RT , P=120 atm ,T =500 K

2
27 ( RT )
a=
64
Pc
b=
270 . 08 22405 .52

=4 . 1908
64111 .3
RTc 0 . 082405 .5
=
=0 .00373
8 Pc
8111. 3
a. MTODO DE LA SECANTE
La forma desarrollada de la ecuacin de Van Der Waals es la siguiente:
V m3 b+
RT
a
ab
V m2 + V m =0
P
P
P
Reemplazando los datos tenemos:
V m30 . 379 V m2 +0 . 0349V m0 . 0013026=0
46
Mtodos Numricos
Resolviendo por el mtodo de la secante, tenemos:
GRAFICA DE LA FUNCION
TABLA INTERACION(EXCEL)
I
1
2
3
4
5
6
7
8
9
10
11
x(i-1)
0
1
0.001986
0.003865
0.039736
0.062173
0.1474577
0.0255976
-0.1163496
-0.0120689
0.0403129
F(xi-1)
-0.0013
0.6546
-0.001235
-0.001173
-4.5149166
-3.5745155
-0.00119
-6.4080547
-0.0255976
-5.2045386
-4.4609018
X(i)
1
0.001986
0.003865
0.039736
0.062173
0.1474577
0.0255976
-0.1163496
-0.0120689
0.0403129
0.0808398
F(xi)
0.6545974
-0.001235
-0.0011733
-4.5149166
-3.5745155
-0.00119
-6.4080547
-0.0255976
-5.2045386
-4.4609018
-4.2978955
X
0.001986
0.003864
0.03974
0.062173
0.1474577
0.0255976
-0.1163496
-0.0120689
0.0403129
0.0808398
1.1493886
error
502.5314
0.902734
0.902734
0.360879
0.578365
4.760599
1.220006
4.467217
0.167586
0.501323
0.929667
47
Mtodos Numricos
12
13
14
15
0.0808398
1.1493886
0.0812743
0.0817116
-4.2978955
1.0565657
-4.327573
-4.3579584
1.1493886
0.0812743
0.0817116
0.0189918
1.0565657
-4.327573
-4.3579584
-7.69638E-
0.0812743
0.0817116
0.0189918
0.1635857
13.14210
0.005351
3.302476
0.8839
16
0.0189918
-7.69638E-
0.1635857
4
-0.001358
-0.17015
1.961420
17
18
19
20
21
22
23
24
25
26
0.1635857
-0.17015
0.184393
0.208691
0.9434958
0.2106618
0.2126064
0.3712495
0.231181
0.2447402
4
-0.001358
-0.0231393
-0.0014840
-0.0014366
0.53413074
-0.0014211
-0.0014039
0.0105858
-0.0011345
-8.030288E
-0.17015
0.184393
0.208691
0.9434958
0.2106618
0.2126064
0.3712495
0.231181
0.2447402
0.2775868
-0.0231393
-0.0014840
-0.0014366
0.53413074
-0.0014211
-0.0014039
0.0105858
-0.0011345
-8.030288E
5.7084082E
0.184393
0.208691
0.9434958
0.2106618
0.2126064
0.3712495
0.231181
0.2447402
0.2775868
0.263939
1.922756
0.116428
0.778811
3.478723
0.009146
0.427322
0.605876
0.0553986
0.118328
0.051708
27
28
0.2775868
0.263939
5.7084082E
-1.067119E
0.263939
0.26608853
-1.067119E
-1.05932E
0.26608853
0.26632542
0.008078
8.8947909E
Mediante el mtodo de la secante nos sale que el volumen molar es: 0.26632542
litros/mol.
a. MTODO DE NEWTON RAPSHON

(P
a
)(V b) RT
V2
27 R 2Tc 2
(
)
64
Pc
RTc
Pc
Datos: Tc=405.5k; Pc=111.3atm; R=0.082

Entonces reemplazando datos:
48
Mtodos Numricos
a 4.19081073
b 0.0373438903863
Desarrollando:
a
)(V b) RT
V2
a ab
PV Pb 2 RT
V V
3
PV PbV 2 aV ab RTV 2
(P
PV 3 ( Pb RT )V 2 aV ab 0
120V 3 (120 0.0373438903863 0.082 500)V 2 4.19081073492V (4.19081073492)(0.0373438903863) 0

120V 3 45.4812668464V 2 4.19081073492V 0.156501176715 0
En Excel:
Sea la ecuacin:
f (v ) 120V 3 45.4812668464V 2 4.19081073492V 0.156501176715 0
Graficando y tabulando:
f(v)
-0,156501177
78,55304271
786,3000529
2843,084529
6968,906472
13883,76588
24307,66276
49
Mtodos Numricos
De:
f (v ) 120V 3 45.4812668464V 2 4.19081073492V 0.156501176715 0

f ' (v) 360v 2 90.9625336929v 4.19081073492
Vk
f(v)
f'(v)
0,71250031
8
0,52299700
7
0,40025727
3
0,32446672
3
0,28358185
4
0,26846836
8
0,26627923
2
0,26623483
8
0,26623482
78,553042
7
23,145250
3
6,7613625
3
1,9293631
8
0,5142012
7
0,1110268
273,22827
7
122,13639
2
3
4
5
6
7
8
9
1
0
0,26623482
0,0125160
4
0,0002439
5
9,9282E08
1,6404E14
0
55,086990
7
25,456513
6
12,576811
7,3462072
1
5,7173430
3
5,4950437
5
5,4905713
2
5,4905694
9
5,4905694
9
| Vk-Vk-1 |
0,2874996
8
0,1895033
1
0,1227397
3
0,0757905
5
0,0408848
7
0,0151134
9
0,0021891
4
4,4395E05
1,8082E08
2,9976E15
50
Mtodos Numricos
Por lo tanto el volumen molar es:

Vm= 0.26623482Lmol-1
El factor de compresibilidad es:
PV (120) (0.2662348)
RT
(0.082) (500)
Z 0.779224
Z
51
Mtodos Numricos

Primer Trabajo Metodos

Cargado por

Copyright:

Formatos disponibles

Primer Trabajo Metodos

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Primer Trabajo Metodos

Cargado por

Copyright:

Formatos disponibles

U N A C Facultad de Ingeniera Qumica

Universidad nacional del callao

FACULTAD: Ingeniera Qumica

Antezana Vergara, Miguel ngel.

Huanca Ortiz, Elvis.

Nuez Zumaeta, Mercedes.

Pillaca Quispe, Elizabeth.

Ponte Valverde, Erick.

Vidal Villarreal, Edward.

Villafuertes Pelayo, Joselyn.

U N A C Facultad de Ingeniera Qumica

U N A C Facultad de Ingeniera Qumica

I. DETERMINAR LAS RACES DE LOS SIGUIENTES EJERCICIOS:

De la tabla y de la grfica, observamos que existen 2

b) Utilizando el mtodo de la biseccin con tres iteraciones para determinar la raz

U N A C Facultad de Ingeniera Qumica

U N A C Facultad de Ingeniera Qumica

absoluto =|x ractual x ranterior|=|6.257.5|=1.25

U N A C Facultad de Ingeniera Qumica

absoluto =|x ractual x ranterior|=|6.8756.25|=0.625

U N A C Facultad de Ingeniera Qumica

c) Utilizando el mtodo de la secante con tres iteraciones para determinar la raz ms

U N A C Facultad de Ingeniera Qumica

absoluto =|x actual x anterior|=|5.910|=4.1

= 5.9 20.5(1.845) = 6.2385

U N A C Facultad de Ingeniera Qumica

absoluto =|x actual x anterior|=|6.23855.9|=0.33858

U N A C Facultad de Ingeniera Qumica

absoluto =|x actual x anterior|=|6.41696.2385|=0.1784

5.2. Determine las races reales de f(x) = 5x3-5x2+6x-2:

U N A C Facultad de Ingeniera Qumica

De la tabla y de la grfica, observamos que las tres

b. Utilizando mtodo de la biseccin para localizar la raz ms pequea. Use los

U N A C Facultad de Ingeniera Qumica

U N A C Facultad de Ingeniera Qumica

error = 0.25 - 0.5 =0.25

error = 0.375 - 0.25 =0.125

error = 0.4375 - 0.375 =0.0625

error = 0.40625 - 0.4375 =0.015625

error = 0.421875 - 0.40605 =0.0078125

U N A C Facultad de Ingeniera Qumica

f(X1)*f(Xr) error ABS

U N A C Facultad de Ingeniera Qumica

c. Mediante el mtodo de la secante

Aplicando la frmula anterior, se tiene:

(- 0.124648) (0.38983 - 0.33333)

(- 0.0017285) (0.41848 - 0.38983)

(- 0.0000125) (0.41809 - 0.41848)

Mtodo de la secante utilizando Excel

U N A C Facultad de Ingeniera Qumica

5.3. Determine las races reales de f(x) = -25182x-90x2+44x3-8x4+0.7x5.

Observando se nota que una

raz se encuentra entre 0

U N A C Facultad de Ingeniera Qumica

POR EL METODO DE NEWTON-RAPHSON, NOS ACERCAMOS POR 16 (RAIZ

U N A C Facultad de Ingeniera Qumica

= 0.00028 = 0,0001 aprox

c) Usando el mtodo de la secante.

U N A C Facultad de Ingeniera Qumica