Gutarra Rios, Juan Peceros Lecca, Brenda
Gutarra Rios, Juan Peceros Lecca, Brenda
Gutarra Rios, Juan Peceros Lecca, Brenda
PECEROS LECCA,
BRENDA
Las ecuaciones diferenciales ordinarias que
modelan una realidad especca, aumentan su
complejidad en la medida que se aproximen cada
vez ms al comportamiento del objeto o fenmeno
en estudio, razn por la cual en la mayoria de
los casos hallar su solucin por mtodos
analticos es imposible lo que nos lleva a
utilizar los mtodos numricos.
El trabajo de investigacin tiene
como objetivo; la resolucin de
ecuaciones Diferenciales ordinarias
con Matlab, en el que se desarrollan
los mtodos numricos: De Euler, de
Taylor de orden superior y Runge-
Kutta.
El mtodo de Runge-Kutta es un mtodo genrico de
resolucin numrica de ecuaciones diferenciales.
El mtodo de Runge-Kutta no es slo un nico mtodo,
sino una importante familia de mtodos iterativos,
tanto implcitos como explcitos, para aproximar las
soluciones de ecuaciones diferenciales ordinarias;
estas tcnicas fueron desarrolladas alrededor de 1900
por los matemticos alemanes Carl David Tolm Runge y
Martin Wilhelm Kutta
Mtodo de Runge-Kutta
de Cuarto Orden
+ = + + + + ,
= , , , . . .
Dnde:
= ( , )
= + , +
= + , +
= + , +
ALGORITMO
Este algoritmo calcula la solucin del problema de valor inicial en puntos
equidistantes 1 = 0 + , 2 = 0 + 2, 3 = 0 + 3, , = 0 + , aqu
tiene una solucin nica en [0 , ].
Dnde:
Solucin:
= 2.681575756 + 207.6315151
+1 = +
1era iteracin:
1
0.5 = 174 + + 22 + 23 + 4 0.5
6 1
1 = 0,174 = 258.9626664
2 = 0.25,109.2593334 = 85.35566446
3 = 0.25,152.6610839 = 201.7407463
4 = 0.5,73.12962685 = 11.52888069
Reemplazando en la ecuacin
1
0.5 = 174 + + 22 + 23 + 4 (0.5)
6 1
0.5 = 105.5311163
2da iteracin:
1 = 0.5 +
1
1 = 105.5311163 + + 22 + 23 + 4 (0.5)
6 1
1 = 75.35816778
2 = 24.83850894
3 = 58.70658203
4 = 3.354905671
Reemplazando en la ecuacin:
1
1 = 105.5311163 + + 22 + 23 + 4 (0.5)
6 1
1 = 85.60666263
2. Un Balance de masa para un qumico en un reactor completamente mezclado puede
escribirse como:
2
=
Donde:
V=volumen(10 3 ) c=concentracin F=flujo msico(200 )
3 3
Q=flujo volumtrico(1 ) k=flujo de reaccin de 2do orden(0.1 .)
Condicin inicial C(0)=0, hallar C(0.5) para un paso h de 0.5
Solucin:
Reemplazando los valores en la funcin queda:
200 1 (0.1)(10) 2
=
10
Iterando:
0.5 = 0 +
1
0.5 = 0 + 1 + 22 + 23 + 4 (0.5)
6
1 = 20
2 = 17
3 = 17.76875
4 = 11.21835059
Reemplazando:
1
0.5 = 0 + + 22 + 23 + 4 (0.5)
6 1
0.5 = 8.396320883