1.Velocidadterminaldeunobjetoencadalibre.

Lavelocidadvdeunparacaidistaquecaeestdadapor:

( ) t
mg
v=
(1e m )
c

Donde:
DATOS
g(m/s^2)
c(Kg/s)
v(m/s)
t(s)

9.81
15
35
9
0.01

a) Mediante la ayuda de un grfica haga un anlisis previo para determinar el intervalo

Intervalo inicial [a0, b0].

y(v)

20

40

0.53
-1.66
60
80

100

120

140

160

-10
-20
-30
-40
X(m)

Figura 1.1: Muestra el punto en donde la grafica atraviesa el 0 en el eje de las x y es donde esta la raz, se
toman los puntos a y b porque son los que se aproximan mas a este

b) Una vez hallado su intervalo inicial y junto con la tolerancia preestablecida resuelva el
problema para determinar la masa mdel paracaidista.
Del inciso (a) podemos obtener los valores de a0 y b0 con los cuales podemos calcular la
primera aproximacin,
x (r )=

ba
2

y posteriormente obtener las funciones f(a), f(b), f[x(r)] y el producto de f(a)*f[x(r)]. Para
obtener a1 y b1 y as sucesivamente hasta n, ocupamos una funcin de restriccin:
a1=si f ( a )f [ x ( r ) ]> 0 el valor sera x ( r ) , sin o sera a 0
b1=si f ( a )f [ x ( r ) ]< 0 el valor sera x ( r ) , sino serab 0
c) Reporte sus resultados parciales en una tabla indicando el nmero de iteraciones n, a, x,
b, f(x) y el error = b- a.
n a

0 56

61 58.5

1 58.5

61 59.75

2 59.75

60.37
61 5

60.37
3 5
60.68
4 75

x(r)

f(a)
31.712953
38
31.193589
88

30.919959
60.687 30.779716
61 5
36
60.843 30.708745
61 75
97

f(b)
30.637212
52
30.637212
52
30.637212
52
30.637212
52
30.637212
52

Resultados
g(m/s^2)
c(Kg/s)
v(m/s)
t(s)
m(kg)

f(x(r))
31.193589
88

f(a)*f(x(r))

error

30.919959
30.779716
36
30.708745
97
30.673049
49

964.504520 0.0209205
2
02

989.240861
7

951.707567 0.0103519
8
67
945.206490 0.0051493
6
31
941.930884 0.0025680
9
53

9.81
15
35
9
60.375

Para calcular el erro relativo porcentual se utilizo la siguiente formula:

x( r)actual x (r )anterior
100
x(r)actual

Se puede observar que la fiabilidad de este metodo es muy alta ya que

2.Equilibrioqumicodeunareaccin.
LaecuacindeequilibrioqumicoparadeterminarlafraccindelCOrestantex,seescribe
como:
k p=

( 1x ) (3+ x)
1
2

x ( x+1) p

1
2

1
2

0< x <1

donde:
kp=Constantedeequilibrio
T=Temperatura

DATOS
x
P=Presin
Kp
x=FraccindeCO
P(atm)
T(K)
a)Haga un anlisis
x0.

?
3.06
5
3000
0.001

previoparadeterminarelpuntoinicial

Paraanalizarcualvaaserelpuntox0 quevamosaconsiderar,utilizocomoreferenciala
restriccindelafraccindeCO (0< x <1) ,graficode0.1a1vslafuncinenx.Deeste
modoseobtienelaraz,queeselpuntoquecruzaelejedesimetra x yeseeselque
consideraremoscomopuntoinicialx0.

Figura2.1:Graficadelaconcentracinxvslafuncinf(x).

b)Unavezhalladoelpuntoinicial,determinelosmolesxconunatoleranciade=0.001.
EnlaFig(2.1)sepuedeapreciarlaraz(x0=0.8).
Conelvalordex0 calculamosf(x0)yf(x0)yutilizandolaformuladeNewtonRaphson
encontramosx(i+1):
x i+1=x i

f ( xi )
f ( xi )

de esta manera calculamos x hasta que la tolerancia tenga un valor de = 0.001 o

x
0

f(x)
0.8

0.63049314
7
0.44531562
2
4
1

f'(x)
error%
ERROR
0.02924038
3 0.17250266
0.22863565 0.21252443 0.2688480 0.08475342
4
7
5
7
7.23478616 77.7230046 5.8096077 0.26895219
6
5
45
3

3
4
5
6
7
8
9

0.35223139
2
0.28222800
2
0.2314948
36
0.19447270
8
0.16682172
4
0.14561874
9
0.12895290
4

7.87307306 112.467024
2
9
172.075306
8.72992518
8
9.73449421 262.937193
4
5
10.8341682
6 391.818536
11.9961274 565.775684
5
5
13.2004794 792.067807
1
4
14.4348860 1078.09170
7
3

0.1167131
59

0.01163552
9

0.1103354
94

0.00437521
2

0.0987559
1

0.00158541
1

0.0869130
24

0.00057847
1

0.0765331
04

0.00021602
3

0.0678624
45 8.28241E-05
0.0606976
17 3.25505E-05

Paracalcularelerrorrelativoporcentualyelerrorabsolutoseutilizaronlassiguientes
formulas:

a =

x actualx anterior
100
x actual + x anterior

debidoaquecadaiteracinreduceelerroralamitad,laformulageneralquerelacionael
erroryelnumerodeiteraciones,n,es
Ena =

x0
2n

Resulttados
Kp
P(atm
)
T(K)

3.06
5
3000
0.001
-0.231494836

3.Sistemasdeecuacionesnolineales.
AplicarelmtododeNewtonmultivariablepararesolverelsistemalineal.
4 x 19 x 2 +2 x3 =5
2 x 14 x 2+ 6 x3 =3
x 1x 2+3 x 3=4
a)Culeselvectordefuncionesf(x)deesteproblema?
f ( x 1 , x 2 , x 3 ) =4 x 19 x 2 +2 x 35
g ( x1 , x2 , x 3) =2 x1 4 x 2 +6 x 33
h ( x 1 , x 2 , x 3 )=x 1x 2 +3 x3 4
b)CulseralamatrizJacobianadeestesistema?

f ( x 1 , x 2 , x 3 )
x1
g ( x 1 , x 2 , x 3 )
x1
h ( x 1 , x 2 , x 3 )
x1

f ( x 1 , x 2 , x 3 )
x2
g ( x1 , x 2 , x 3 )
x2
h ( x1 , x2 , x3 )
x2

f ( x 1 , x 2 , x 3 )
x3
4 9 2
g ( x 1 , x 2 , x 3 )
= 2 4 6
x3
1 1 3
h ( x 1 , x 2 , x 3 )
x3

c)Estamatrizesconstanteovariable?
Constante.

d)CuleslamatrizinversadelaJacobiana?

f ( x 1 , x 2 , x 3 )
x1
g ( x 1 , x 2 , x 3 )
x1
h ( x 1 , x 2 , x 3 )
x1

f ( x 1 , x 2 , x 3 )
x2
g ( x1 , x 2 , x 3 )
x2
h ( x1 , x2 , x3 )
x2

f ( x 1 , x 2 , x 3 )
x3
g ( x 1 , x 2 , x 3 )
x3
h ( x 1 , x 2 , x 3 )
x3

0.3 1.25 2.3

e)EncuentrelasolucinaestesistemautilizandoelmtodoNewtonmultivariableconel
vectorinicialx0=(0,0,0).
Para resolver este problema se utilizara la ecuacin de Newton Raphson para tres
ecuaciones.

( )( )(
x1
x2
x3

( i+ 1)

( i+1)

( i+1)

f ( x 1 , x 2 , x 3 )
x1
f ( x1 , x2 , x3 )
x1
g ( x1 , x 2 , x 3 )
= x2 g ( x1 , x2 , x3 )
x1
x3
h ( x1 , x2 , x3 )
h ( x1 , x2 , x3 )
x1

f ( x 1 , x 2 , x 3 )
x2
g ( x 1 , x 2 , x 3 )
x2
h ( x 1 , x 2 , x 3 )
x2

f ( x1 , x2 , x3 )
x3
g ( x 1 , x 2 , x 3 )
x3
h ( x 1 , x 2 , x 3 )
x3

f)Resumaenunatablaelnmerodeiteracionesn,losvaloresdex1,x2,x3yelerror
estimadoatravslanormaeuclidianadelvectordefuncionesf ( X n) ,conunaprecisin
dee =0.001.

x
1
2
3

xi

f(xi)
0
0
0

-5
-3
-4

X(x+i)
0.25
0.15
0.2

1
2
3

0.25
0.15
0.2

-4.95
-1.9
-3.3

-0.0025
-0.055
-0.035

1
2

-0.0025
-0.055

-4.585
-2.995

0.23175
0.20475

-0.035

-4.0525

0.237625

0.23175

-5.4405

0.20475

-1.92975

0.040275
0.108262
5
0.074618
75

1
2
3

0.237625 -3.260125

0.040275 -4.013775
0.108262 2.934112
5
5
0.074618 4.075318
75
75

0.160413
75
0.254968
125
0.278384
688

Sol. Exacta
x1
x2
x3

6.95
2.5
0.15

LanormaEuclidianadelvectordefuncioneses:

A=f ( x n )=

i=1 j=1

aij

entonceslanormadelamatriz

)(

a 11 a12 a13
4 9 2
=
a 21 a22 a23
2 4 6
1 1 3
a 31 a32 a33

A=f ( x n )= a112 +a122+ a132+ a212+ a22 2+ a232+ a312+ a322 +a 332

A=f ( x n )= 2 42 +9 2+22 +22 + 42 +62 +12 +12+ 32

Precisin

0.001
7.71517E05

4.Columnadeabsorcindeseisplatos.
Resuelvaestesistemadeecuacioneslinealesparalascondicionesdadasenelproblema
usandoelmtododeGauss.
Elconjuntocompletodeecuacionesparalacolumnadeseisplatoses,
( L+Ga ) x1 +Ga x 2=L x 0
L x 1( L+Ga ) x 2 +Ga x 3=0
L x 2( L+Ga ) x 3 +Ga x 4 =0
L x 3( L+Ga ) x 4+ Ga x 5=0
L x 4 ( L+Ga ) x 5+ Ga x 6=0
L x 5( L+Ga ) x6 =Ga

y 7b
a

Unconjuntotpicodeparmetrosparaestesistemason:
a=0.72
b=0.0,G=66.7kmol/min,
L=40.8kmol/min
x0=0.0(lquidopuro)y
y7=0.3(kmoldesoluto/kmoldeinertes)
Lasolucinsimultneapuedeexpresarseennotacinmatricialcomo

sustituyendolosvaloresconocidoslamatrizsepuedeexpresardelasiguientemanera.

-88.824
40.8
0
0
0
0

48.024
-88.824
40.8
0
0
0

0
48.024
-88.824
40.8
0
0

0
0
48.024
-88.824
40.8
0

0
0
0
48.024
-88.824
40.8

0
0
0
0
48.024
-88.824

0
0
0
0
0
-20.01

PararesolverelproblemaporelmtododeGaussesnecesarioconvertirlamatrizala
matrizidentidad:
1
0
0
0
0
0

0
1
0
0
0
0

0
0
1
0
0
0

0
0
0
1
0
0

0
0
0
0
1
0

0
0
0
0
0
1

X1
X2
X3
X4
X5
X6

Primeramenteestosehacedividiendoelprimerterminodelaprimerafilaentreelmismo(
88.824/88.824)=1paraquenosdeelprimer1delamatrizidentidadydespuselprimer
terminodelasegundafilamultiplicarlopor(menos)el1delaprimerafilaysumrseloalel
mismo,sea,[(1*40.8)+40.8],paraquenosdeelceroquenecesitamos.

1
0
0
0

0.5406646
85
66.764880
84
40.8
0

48.024
-88.824
40.8

0
48.024
-88.824

0
0
48.024

0
0
0

0
0
0

40.8

-88.824

48.024

40.8

-88.824

0
-20.01

Yassucesivamenteahoraconelvalordearribayabajodel1necesarioparaformarla
matrizidentidadhastalograrladiagonaldepuros1.

0
0
0

0
0
0

0.3889002
8
0.7193003
18
59.476547
04
40.8
0

48.024
-88.824
40.8

0
48.024
-88.824

0
0
48.024

0
0
0

40.8

-88.824

48.024
-88.824

0
48.024

0
0

40.8

-88.824

0
0

0
0

0
0

0.3140153
21
0.5807949
55
0.8074443
19
55.880271
8
40.8

0.2698675
46
0.4991403
23
0.6939247
92
0.8594088
48

-20.01

-20.01

53.760119
01
48.024

40.8

-88.824

0.2410731
09
0.4458828
46
0.6198841
22
0.7677112
94
0.8933015
94
52.377294
97

Unavezlogradolamatrizidentidadseoobtienenlosvaloresdelasxs.
RESULTADOS
x1 0.09209855
x2 0.170343195
x3 0.236817905
x4 0.293293172
x5 0.341273159
x6 0.382035766

0
-20.01

0
0
0
0
0
-20.01

0.092098
55
0.170343
195
0.236817
905
0.293293
172
0.341273
159
0.382035
766

5.Derivacinnumrica.
Evalelaprimeraysegundaderivadasde
f ( x )=sen ( x )
Medianteelesquemadediferenciasfinitashaciadelante.Evaleelerrorcomparandoel
resultadonumricoconelvalorexacto.
FUNCION
(X)
f(xi)

h
x

0.001
/3

DERIVADAS EXACTAS
f(x)
f"(x)

xi

1.0471975
51

0.866025
404

x(i+1)=xi
+h

1.0481975
51

0.866524 0.4991337
971
25

x(i+2)=xi
+2h

1.0491975
51

0.867023 0.4982669
671
5

x(i-3)=xi3h

1.0501975
51

0.867521 0.4973996
504
78

0.5

0.8660254
04
0.8665249
71
0.8670236
71
0.8675215
04

Formulas por diferencias finitas hacia delante; se presentan dos versiones para cada
derivada.LaultimaversinincorporamastrminosdelaseriedeexpansindeTaylory,en
consecuencia,esmasexacta.
Primeraderivada:
f 1 ( x i )=

f ( x i+1 ) f ( xi )
h

f 2 ( x i )=

f ( x i+2 ) + 4 f ( x i+1 ) 3 f ( x i )
2h

Segundaderivada:

f } rsub {1} left ({x} rsub {i} right ) = {f left ({x} rsub {i+2} right ) +2f left ({x} rsub {i+1} right ) +f( {x} r

f } rsub {2} left ({x} rsub {i} right ) = {-f left ({x} rsub {i+3} right ) +4f left ({x} rsub {i+2} right ) -5f left (

DERIVADAS
NUMERICAS
PRIMERAD DERIVADA
F(1)
F(2)
0.4995669 0.5000001
04
66
SEGUNDA DERIVADA
F"(1)
F"(2)
0.8665248 0.8660261
98
98

Elerrorcomparandoelresultadonumricoconelvalorexactonoarroja.
=

f (x) exactaf (x)numerico

f (x)} rsub {exacta}}

f (x)} rsub {exacta} - {f ( x)numerico

=PRIMERA DERIVADA
F(1)
F(2)
0.0008661 3.32899E92
07
=SEGUNDA DERIVADA
F"(1)
F"(2)
0.0005767 9.16926E67
07

Sepuedeapreciarquelasderivadasnumricassonextremadamenteexactasparaestecaso.
6.Integracinnumrica.
Conelmtododeltrapezoideintegrenumricamentelasiguientefuncin.

( 8+ 3 sen ( x ) ) dx
0

desdea=0hastab=,conn=5.Empleandolaregladeltrapezoidal,recordandoelvalor
correctodelaintegral(31.1328),calculadaconWolframAlpha.
Hay(n+1)puntosbaseigualmenteespaciados(x0,x1,x2,x3,x4,xn).Enconsecuencia,
haynsegmentosdeigualanchura.Eltamaodepasoentoncessecalcula.

h=

ba
n

n
xo=a

5
0
3.141592
654
0.628318
531

xi=b
h=(ba)/n

Paraesteproblemasedebeutilizarlareglatrapezoidaldemltipleaplicacin:
n1

f ( x 0 ) +2 f ( x i ) +f ( x n )
I =(ba)

i=1

2n
n1

donde (ba) eslaanchuray

f ( x 0 ) +2 f ( x i ) +f ( x n )
i=1

eslaalturapromedio.De

2n
acuerdoconestaecuacin,lospuntosinterioresdandosveceselpesodelosdospuntos
extremos f ( x 0 ) y f ( x n ) .

f(a)=f(xo)

f(h)=f(xi)
f(b)=f(xn)
8 9.763355757
8

Valor
correcto

Valor
numerico

error
0.641502
31.1328 11.16104559
03

Conestemtodoseobtieneunamalaaproximacinalvalorcorrectodelaintegral,yaque
elerroresdel(64%).
7.Sistemasdeecuacionesdiferencialesordinarias.
UseelmtododeEulerpararesolver:
dy
=2 y
dx

dz y z
=
dx
3

enelrangodex=0a1,conuntamaodepasode0.2,ylascondicionesinicialesy(0)=2,
yz(0)=4.
Paraesteproblemanecesitamosresolverlasintegralesyobtenerlassolucionesexactasy
podercompararlasdespusconelmtododeEuler.
y=e2 x
z=

3
yx

usandolassiguientesexpresiones
y i+1= yi + f ( xi , y i ) h
z i+1=z i +f ( x i , z i ) h
seprediceunnuevovalordeyyzpormediodelapendiente,quehabrdeinterpolarse
enformalinealsobreeltamaodepasoh.
EstaformulaesconocidacomoelmtododeEuler(oEulerCouchyodepuntomedio).
Paraejemplificarutilizaresoloalayporquedeigualmanerasehaceelmismo
procedimientoparaz.
UtilizandolaecuacindeEulerenx=0.2queeslaprimeraiteracineigualaltamaode
pasotenemos,
y ( 0.2 )= y ( 0 ) +f (0,2)(0,2)

dondey(0)=4ylapendienteestimadaenx=0es
f ( 0,2 )=2 ( 2 )=4
porlotanto
y ( 0.2 )=2+ (4 )( 0.2 ) =1.2
enlasiguientetablasepuedeapreciarlosvaloresrealesconlosresultadosobtenidosporel
mtododeEuler,porlovistonoesmuyprecisoparaestecasoyaqueelerroesmuy
grande.
laformaenquesecalculoelerrosemuestraenlasiguienteformula.

valor verdaderovalor aproximado

y(verdader
o)

z(verdader
o)

y(euler)

0.2

0.6703200
46

1.2

0.4

0.4493289
64

1.8202684
14

0.6

0.3011942
12

1.8795223
15

0.8

0.2018965
18

1.9192413
93

0.1353352
83

1.9458658
87

z(euler)

Error

f(x,y)

22.377370 1.86666666
46
7
16.691556 4.34577848
96
2
16.600584 1.53352820
61
5
18.573871 0.64346789
59
8
22.167168 0.43343365
3
9

20.517542
03
21.081958
11
18.202669
24
19.293921
59
22.673006
63

f(x,z)
10.6666666
7
111.886852
3
41.7288924
1
27.6676410
2
23.2173394
9

-4
1.34064009
2
0.89865792
8
0.60238842
4
0.40379303
6
0.27067056
6 22.1671683