1) Una viga uniforme con un extremo libre de 5m de longitud, y peso constante de 300kg,

esta empotrada en uno de sus extremos, hallar la deflexión *y* en el extremo libre asi
como la curva elástica de la viga.
y W
El = ( L−x )2
1
2 2
2
( 1+ ( y ) )
Si se asume que la viga se flexiona muy poco, la pendiente de la curva elástica es minima
por lo tanto 1+ ( y )2=1 y El= 150000 . De esta manera el sistema a resolver queda:

dy dz W
=z = ( L−x )2
dx dx 2 El

Y(0)=0

Z(0)=0

Y(5)=?
2) La rapidez de flujo de calor (conducción), entre dos puntos sobre un cilindro calentado en
un extremo, está dada por:
dQ dT
=λA
dt dx
Donde λ , es una constante , A es la sección transversal del cilindro, Q el flujo de calor, T la
temperatura , `t` el tiempo y `x` distancia desde el extremo calentado. Como la ecuación
involucra dos variables, la simplificaremos haciendo:
dT 100( L−x)(20−t)
=
dx 100−xt
Donde L es la longitud del cilindro, combinar las dos ecuaciones y calcular el flujo de calor
desde t=0 hasta 25 segundos. La condición inicial es Q(0)=0, y los parámetros λ =0.40
cal/cm-s, A=10 cm^2 , L=20cm y x=2.5cm. Grafique sus resultados.
3) La encuacion diferencial de primer orden que se muestra describe el comportamiento del
flujo gradualmente variado representado por la curva de remanso o eje hidráulico.
Calcular las coordenadas de dicha curva, si se sabe que ; Q es el caudal , T , es el espejo de
agua de la sección hidráulica , A el área de la sección `g` aceleración de la gravedad , SE
pendiente de la línea de energía, So la pendiente de fondo del canal, Utilizar los
siguientes datos ; Q=10m3/s, So= 0.001 , n= 0.016, z= 1, ancho de solera b=3m . Los
valores iniciales son x=0, y= 3.5m y se deben realizar cálculos hasta un tirante final de
2.8m.
Se
dy
=So (
( 1− )
So
) 1
Q= A R 2/ 3 S1 /2
dx 2
QT n
( 1−
g A3 )
 4) El caso de un reactor de mezcla completa con flujo de entrada y salida, es representado
por la siguiente ecuación diferencial:

dc
V =Qcem−Qc
dt
Si se tiene que Cem=50mg/m3 , Q=5m3/min, V=100m3 y Co=10mg/m3 . Si C=Co en t=0 .
Obtener la solución analítica y numérica de la ecuación. Graficar, calcular errores y comparar
resultados
dC
=2.5−0.05 C C=e 0.05 t−e3.688 +50
dt
n h t(min) C(analitico) C(euler) Error(%)(analitico-euler)
0 1 0 10 10 0
11.054471
1 1 1 1 12 8.55336192
11.108370
2 1 2 9 13.9 25.1308598
11.165034
3 1 3 2 15.705 40.6623548
11.224602
4 1 4 7 17.41975 55.1925755
11.287225 19.048762
5 1 5 3 5 68.7639075
11.353058 20.596324
6 1 6 7 4 81.4165235
 11.422267 22.066508
7 1 7 4 2 93.1885097
11.495024 23.463182
8 1 8 5 7 104.115987
24.790023
9 1 9 11.571512 6 114.233228
26.050522
10 1 10 11.651921 4 123.572769

RELACION EULER - ANALITICO

30
25
20
15
10
5
0
0 2 4 6 8 10 12

C(analitico) EULER

5) Un cable cuelga de dos soportes A y B, el cable esta con cargas distribuidas cuyas
magnitudes varian con `x` como:
πx
W =Wo(1+ sen ( 2lo ))
Donde Wo= 1000lb/ft , la pendiente del cable (dy/dx = 0) en x=0, que es punto mas bajo
del cable. Este también es el punto del cable donde la tensión tiene un minimo, To La
ecuación diferencial que modela el cable es:
Resuelva el caso usando un método numérico adecuado. Grafique la forma del cable (x,y)
6) Un cono circular recto de altura Ho y radio R dispuesto verticalmente con orificio circular en el
fondo de diámetro Ø . Determinar el tiempo de vaciado en función de la altura h si Ø = 2pulg y es el
diámetro del tubo de salida , R=1m y Ho=1.2m . Si el sistema generado con los valores iniciales y al
EDO es:
3 /2 dh −4.8 Ø 2 Ho 2
h =
dt 576 R 2
h(t=0)=Ho=1.2
 h(t=?)=0

  h euler analitica
0 1.2 0 0
-
172890.32 406136.53
1 2.4 2 9
- -
490510.12 889611.40
2 3.6 6 1
-
979517.80
3 4.8 3 -1890548.6
- -
1662926.8 3348309.1
4 6 1 9
- -
2561291.2 5317056.7
5 7.2 8 2
- -
3693359.8 7845706.2
6 8.4 9 6
- -
5076482.4 10979212.
7 9.6 7 6
- -
6726883.3 14759402.
8 10.8 8 8
- -
8659855.9 19225551.
9 12 5 8
- -
10889906. 24414800.
10 13.2 6 3
-
30362469.
11 14.4 -13430865 6
- -
16295971. 37102305.
12 15.6 1 2
- -
19497944. 44666671.
13 16.8 7 5
- -
23049042. 53086708.
14 18 1 2
15 19.2 - -
26961103. 62392460.
 5 1
- -
31245592. 72612984.
16 20.4 9 9
- -
35913631. 83776444.
17 21.6 6 9
- -
40976027. 95910184.
18 22.8 3 7
-
46443299. -
19 24 4 109040799
-
52325700. -
20 25.2 5 123194190

Chart Title
0
2 4 6 8 6 2 4 6 8 2 2 4 6 8 8 2 4 6 8 4 2
1. 2. 3. 4. 7. 8. 9. 10. 1 13. 14. 15. 16. 1 19. 20. 21. 22. 2 25.
-20000000

-40000000

-60000000

-80000000

-100000000

-120000000

-140000000

7) En un tanque perfectamente agitado se tiene 400 litros de una salmuera en la cual están
disueltos 25kg. De sal común(Na), en cierto momento se hace llegar al tanque un gasto de
80litros/min de una salmuera que contiene 0.5kg de sal común por litro. Si tiene un gasto
de salida de 80l/min determine:
a) Que cantidad de sal hay en el tanque, transcurridos 10 min .
b) Que cantidad de sal hay en el tanque, transcurridos un tiempo muy grande.

Si llamamos *x* a los kg de sal en el tanque después de *t* minutos, la acumulación de sal
en el tanque estará dada por:
 dx
=masa de sal que entra−masa de sal que sale
dt
dx x
dt
=80 ( 0.5 )−80 ( )
400
=40−0.2 x

Usando la condición inicial de que hay 25kg de sal en el tiempo cero tenemos :

dx
=40−0.2 x
dt
X(t=0)=25

X(t=10)=? h t X K1 K2 K3
1 0 25      
1 1 56.7333333 35 31.5 29.4
1 2 82.7123556 28.6533333 25.788 24.0688
1 3 103.980515 23.4575289 21.111776 19.7043243
1 4 121.392048 19.203897 17.2835073 16.1312735
1 5 135.64629 15.7215903 14.1494313 13.2061359
1 6 147.315763 12.870742 11.5836678 10.8114232
1 7 156.869171 10.5368474 9.48316267 8.85095182
1 8 164.690228 8.62616575 7.76354917 7.24597923
1 9 171.093067 7.06195436 6.35575892 5.93204166
1 10 176.334857 5.78138663 5.20324797 4.85636477
Para un tiempo
muy grande

H t x K1 k2 k3
1 0 25      
56.733333
1 1 3 35 31.5 29.4
82.712355 28.653333
1 2 6 3 25.788 24.0688
103.98051 23.457528 19.704324
1 3 5 9 21.111776 3
121.39204 17.283507 16.131273
1 4 8 19.203897 3 5
15.721590 14.149431 13.206135
1 5 135.64629 3 3 9
147.31576 11.583667 10.811423
1 6 3 12.870742 8 2
156.86917 10.536847 9.4831626 8.8509518
1 7 1 4 7 2
164.69022 8.6261657 7.7635491 7.2459792
1 8 8 5 7 3
171.09306 7.0619543 6.3557589 5.9320416
1 9 7 6 2 6
176.33485 5.7813866 5.2032479 4.8563647
1 10 7 3 7 7
 180.62613 4.7330285 4.2597256 3.9757439
1 11 7 2 7 6
1 12 184.13926 3.8747726 3.4872954 3.2548090
4 9 2 6
1 13 187.01534 3.1721472 2.8549325 2.6646036
4 4 1 8
1 14 189.36989 2.5969312 2.3372380 2.1814222
5 1 9 1
1 15 191.29748 2.1260210 1.9134189 1.7858576
7 1 1 5
1 16 192.87554 1.7405025 1.5664522 1.4620221
3 4 8 3
1 17 194.16744 1.4248914 1.2824022 1.1969087
4 1 7 8
1 18 195.22508 1.1665111 1.0498599 0.9798693
1 9 2
1 19 196.09093 0.9549837 0.8594853 0.8021863
3 5 8 5
1 20 196.79977 0.7818133 0.7036320 0.6567232
7 7 3 3
1 21 197.38008 0.6400445 0.5760400 0.5376374
4 4 9 2
1 22 197.85516 0.5239831 0.4715848 0.4401458
2 3 2 3
1 23 198.24409 0.4289675 0.3860707 0.3603327
3 2 7 2
1 24 198.56249 0.3511814 0.3160632 0.2949923
7 1 7 9
1 25 198.82316 0.2875005 0.2587504 0.2415004
5 2 7 3
1 26 199.03656 0.2353670 0.2118303 0.1977083
4 9 8 6
1 27 199.21126 0.1926871 0.1734184 0.1618572
7 9 7 4
1 28 199.35429 0.1577465 0.1419719 0.1325071
1 8 2 3
1 29 199.47137 0.1291418 0.1162276 0.1084791
9 7 8 7
1 30 199.56723 0.1057241 0.0951517 0.0888082
6 4 3 8
 8) Calcular el tiempo para que el nivel de un tanque de forma esférica con radio 5m pase de
4m a 3m, la velocidad de salida del orificio de fondo es V=4.895h^1/2 , siendo el diámetro
de dicho orificio de 10 cm. El sistema que resulta luego del análisis hidráulico del problema
y las condiciones iniciales es:
dh −0.122375 √ h
=
dt (10 h−h2 )
h ( 0 ) =4
h(¿?)=3

  Paso h= -0.1
Método de Runge Kutta de 4to Orden
  x y
altura(m
t(s) k1 k2 k3 k4
  )
- - -
-
  4 0 98.05924 98.2564836 98.4395078
98.2564836
4 2 4
- - -
  3.9 9.825411441 98.43950 98.6081419 -98.608142 98.7622062
8 8 3
- - -
-
  3.8 19.68598281 98.76220 98.9015154 99.0258791
98.9015155
6 8 3
- - -
-
  3.7 29.57588526 99.02587 99.1351008 99.2289782
99.1351008
9 4 6
- -
-
  3.6 39.48913961 99.22897 99.3073027 -99.3698591
99.3073028
8 5
  3.5 49.41960708 - - - -
99.36985 99.4164251 99.4164252 99.4467716
 9 8 7
- - -
-
  3.4 59.36097927 99.44677 99.4606616 99.4578502
99.4606617
2 5 4
- -
-
  3.3 69.30676708 -99.45785 99.4380842 99.4011015
99.4380842
1 1
- - -
-
  3.2 79.25028856 99.40110 99.3466308 99.2743912
99.3466309
2 8 9
- - -
-
  3.1 89.1846555 99.27439 99.1840914 99.0754292
99.1840915
1 5 4
- - -
-
  3 99.1027586 99.07542 98.9480910 98.8017513
98.9480911
9 8 4

Gráfica
99.1
89.18 100
79.25
69.31 80
59.36
Tiempo (s)

49.42 60
39.49
29.58 40
19.69
9.83 20
0
0
4.2 4 3.8 3.6 3.4 3.2 3 2.8
Altura (m)

9) La reacción de NO con H2 se conoce como una reacción de 3er orden y su rapidez esta
dada por :
dx
=k ( PNO−2 x )2 ( PH 2−x )
dt
Donde x es la presión parcial del producto, K es la constante de la velocidad , PNO y PH2
representan las presiones parciales del NO y del H2 respectivamente, en la mezcla de
reacción antes de que ocurra la reacción ; PNO-2x es la presión parcial del NO en el reactor
y PH2-x es la presión parcial del H2 en el reactor a cualquier tiempo después que la
reacción inicio. En este problema PNO=359mmHg PH2=400mmHg K=1.2x10^-7mm^-2S^-1
Y la condición inicial es que x=0 , 0≤t≤3min. Encontrar la relación x(t) y comparar el
resultado con la solución analítica

N h ti xi (Euler)
0 18 0 0
1 18 18 111.353184
2 18 36 122.934866
3 18 54 130.59424
4 18 72 136.161489
5 18 90 140.443028
6 18 108 143.863953
 7 18 126 146.674323
8 18 144 149.032739
9 18 162 151.045525
10 18 180 152.78707

PRESIONES PARCIALES
180
160
140
120
100
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200

xi (Euler)

10) Se usa un reparador estándar de aceite para separar una mezcla de agua y aceite. El nivel
de derrame de aceite esta 5ft arriba del extremo inferior de la pared de retención.
Obtener la relación v(t) usando la ecuación :
dv Pw−Po 3 fPw 2
dt
=
Po (g−( )
4 DPo
)v

Donde
D= diámetro promedio del glóbulo de aceite 1/16 in
V= velocidad vertical del glóbulo de aceite
Po=densidad del aceite = 0.6Pw
Pw= densidad del agua
g= aceleración debida a la gravedad = 32.2ft/s2
f= coeficiente de friccion = 0.06
La condición inicial es que en t=0,v=0 t(0;0.03s)

dv
=21.47−1.2 v 2
dt
n H t k1 k2 k3 k4 v (RK) v(EULER)
0 0.003 0         0 0
1 0.003 0.003 21.47 21.468755 21.468756 21.465022 0.064405 0.06441
2 0.003 0.006 21.465022 21.458802 21.458804 21.372812 0.1287415 0.1288051
3 0.003 0.009 21.450111 21.438927 21.438933 21.270784 0.1929799 0.1931553
 4 0.003 0.012 21.425311 21.409186 21.409199 21.159037 0.2570904 0.257431
5 0.003 0.015 21.390685 21.369652 21.369674 21.037689 0.3210439 0.3216025
6 0.003 0.018 21.346317 21.320415 21.320448 20.906872 0.3848114 0.3856401
7 0.003 0.021 21.292304 21.261583 21.26163 20.76674 0.4483641 0.4495147
8 0.003 0.024 21.228764 21.193281 21.193343 20.617459 0.5116739 0.5131973
9 0.003 0.027 21.155828 21.11565 21.115728 20.459212 0.5747128 0.5766592
10 0.003 0.03 21.073646 21.028847 21.028944 20.292199 0.6374535 0.639872

T v(EULER) v(RK)
0 0 0 11) Una puerta cuelga de un soporte sin friccion.,
0.003 0.064405 0.06441 inclinado como muestra la fig. La ecuación
diferencial que expresa la oscilación de la puerta
0.006 0.1287415 0.1288051
es :
0.009 0.1929799 0.1931553
0.012 0.2570904 0.257431
0.015 0.3210439 0.3216025
0.018 0.3848114 0.3856401 d2 ∅ 3 g
0.021 0.4483641 0.4495147 + sin∝ sin ∅=0
dt2 2b
0.024 0.5116739 0.5131973
0.027 0.5747128 0.5766592
0.03 0.6374535 0.639872
En la que g=32.2ft/s^2 , sean los
valores iniciales ∅ o=90º , (d∅ /dt) t=0
=4 pi rad/s .Resolver la ecuación en
el intervalo 0≤t≤3s
12) Sea *x* la temperatura al tiempo *t* de un cuerpo inmerso en un medio cuya
temperatura esta descrita por la expresión:
G(t)=-100(t-4)
La temperatura satisface la ecuación diferencial:
dx
+ kx=kG(t )
dt
 Donde k=0.045 y es una constante de proporcionalidad y la condición inicial es que en
t=0 , x=1000ºF . Obtener la relación x-t con el procedimiento de Runge-kutta y comparar
el resultado con la solución analítica .

x=400− ( 100k )( kt−1)+(600− 100k ) exp ⁡(−kt )

2000
1000
Temperatura(°F)

0
0 10 20 30 40 50 60 70
-1000
Analíica
-2000 Runge-Kutta
-3000
-4000
Tiempo(min)

13) En un recipiente que contiene 120lb de agua se vacian 20lb de azúcar. La concentración
porcentual *c* de la solución a cualquier tiempo t esta expresado como:
dc k
( 120−1.212 c ) = (200−14 c)(100−4 c)
dt 3
 Donde k , el ceoficiente de transferencia de masa , es igual a 0.05889. La condición inicial
es que en t=0 , c=0 . Encontrar la relación c(t) por los métodos Runge kutta y de Euler y
comparar el resultado con la solución analítica representado por:

h Ti Ci(Analt.) k1 k2 k3 k4 Y.R.K Euler

10 0 0         0 0

10 10 2.7672667 3.27166667 -0.19654935 3.59847064 3.43020066 22.5095165 32.7166667

10 20 4.79499224 2.41337613 -0.06007394 2.50195746 0.47914541 15.7277476 52.1753754
10 30 6.33814941 1.84605379 0.03050626 1.8055957 -0.21921522 13.6267297 251.580288
10 40 7.54670781 1.45160055 0.08350312 1.35065753 -0.32833779 12.9907895 -2944.70393
10 50 8.51465173 1.1663599 0.11167115 1.04186027 -0.28968551 12.8529365 23240.6524
10 60 9.30388287 0.95352884 0.1243969 0.8243278 -0.22468911 12.8918003 -188095.452
10 70 9.95683947 0.79061561 0.12777576 0.66591473 -0.16554655 12.9912996 1517382.69
10 80 10.503593 0.66324039 0.12562981 0.54718347 -0.11871348 13.1070952 -12245770
10 90 10.9660502 0.56185895 0.12033241 0.45596792 -0.08347926 13.2218936 98822418
10 100 11.3605504 0.47993313 0.1133534 0.38439536 -0.05754921 13.329186 -797494193
10 110 11.6995299 0.41286118 0.10560044 0.3272119 -0.03866523 13.4269181 6435751278
10 120 11.9926228 0.35732743 0.09762996 0.28081173 -0.02499378 13.5148916 -5.1936E+10
10 130 12.2474091 0.31089327 0.0897777 0.2426568 -0.01514282 13.5936552 4.1912E+11
10 140 12.4699341 0.27173234 0.08223999 0.21091945 -0.00808393 13.6640212 -3.3823E+12
10 150 12.6650784 0.23845455 0.07512493 0.18425638 -0.00306535 13.7268538 2.7295E+13
10 160 12.8368252 0.20998644 0.06848485 0.16166185 0.00046095 13.7829797 -2.2027E+14
10 170 12.9884569 0.18548817 0.06233705 0.14237019 0.00289489 13.8331544 1.7776E+15
10 180 13.1227025 0.16429484 0.05667716 0.12578961 0.00452919 13.8780529 -1.4345E+16
1−e−0.0098 t
c= 100
7−4 e−0.0098t

4.00E+15
2.00E+15
0
-2.00E+15 0 20 40 60 80 100 120 140 160 180 200
-4.00E+15
Analitica
-6.00E+15 Runge Kutta
-8.00E+15 Euler
-1.00E+16
-1.20E+16
-1.40E+16
-1.60E+16

14) Un tanque hemisférico de radio R esta inicialmente lleno de agua. En el fondo del tanque
hay un agujero de radio r por el cual escapa el agua bajo la influencia de la gravedad. La
ecuación diferencial que expresa la profundidad del agua como función del tiempo que se
ha obtenido es:
 dy r2 √ 2 g
= =0
dt 2 R y 1/ 2− y 3 /2
Donde g=32.2ft/s^2 R=10ft r=1/12ft. La condicion inicial es que en t=0 , y=0. Encontrar la
relación entre y,t

t euler kuta k1 k2 k3 k4
0 0 0
0.0242500 0.0406924 0.0485000
0.5 2 7 4 0.0401144 0.0401144 0.0351972
0.0418486 0.0727452 0.0319010 0.0319010
1 2 5 0.0351972 8 8 0.0295151
0.0566061 0.1005095 0.0262708
1.5 7 7 0.0295151 0.0277 0.0277 5
0.0262708 0.0251172 0.0251172 0.0241687
2 0.0697416 0.125661 5 5 5 2
0.0818259 0.1490597 0.0241687 0.0233779 0.0233779 0.0227118
2.5 6 2 2 4 4 3
0.0931818 0.1712200 0.0227118 0.0216642
3 7 7 3 0.0221465 0.0221465 5
0.0216642 0.0212516 0.0212516 0.0208983
3.5 0.104014 0.1924816 5 5 5 4
0.1144631 0.2130852 0.0208983 0.0205961 0.0205961 0.0203387
4 7 3 4 8 8 3
0.1246325 0.2332119 0.0203387 0.0201207 0.0201207 0.0199381
4.5 3 1 3 9 9 8
0.1346016 0.2530042 0.0199381 0.0197874 0.0197874 0.0196658
5 2 2 8 6 6 3
0.1444345 0.2725793 0.0196658 0.0195709 0.0195709 0.0195010
5.5 3 6 3 9 9 5
0.1541850 0.2920374 0.0195010 0.0194544 0.0194544 0.0194299
6 6 9 5 5 5 3
0.1639000 0.3114673 0.0194299 0.0194264 0.0194264 0.0194432
6.5 2 4 3 7 7 7
0.1736216 0.3309502 0.0194432 0.0194797 0.0194797 0.0195353
7 6 5 7 1 1 4
0.1833893 0.3505632 0.0195353 0.0196098 0.0196098 0.0197031
7.5 3 4 4 7 7 4
0.3703815 0.0197031 0.0198151 0.0198151 0.0199459
8 0.1932409 1 4 3 3 3
0.2032138 0.3904805 0.0199459 0.0200957 0.0200957 0.0202650
8.5 7 3 3 9 9 5
0.2133463 0.4109381 0.0202650 0.0204541 0.0204541 0.0206638
9 9 4 5 9 9 2
0.0206638 0.0208946 0.0208946 0.0211476
9.5 0.2236783 0.4318365 2 7 7 3
0.2342521 0.4532642 0.0211476 0.0214237 0.0214237 0.0217241
10 2 9 3 3 3 8
 Chart Title
euler kuta
0.5
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
0
0 2 4 6 8 10 12