Catedra Metodos Numericos 2015 Unsch 15 PDF

CATEDRA 1
Facultad de Ingeniera de Minas, Geologa y Civil

Departamento acadmico de ingeniera de minas y civil
METODOS
NUMERICOS
Ingeniera Civil
Capitulo XV
Solucin de Ecuaciones
Diferenciales Ordinarias
SOLUCIN NUMRICA DE ECUACIONES
DIFRENCIALES ORDINARIAS
3
C0NTENIDO
FORMULACIN DEL PROBLEMA DE VALOR INICIAL

MTODOS NUMRICOS PARA SOLUCIONAR UN PVI
MTODO DE EULER
MTODO DE TAYLOR
MTODO DE EULER MODIFICADO
MTODO DE RUNGE-KUTTA
MTODO DE PREDICCIN Y CORRECCIN
SOLUCIN DE SISTEMAS DE ECUACIONES DIFERENCIALES DE
SEGUNDO ORDEN
24/07/2015 4
PVI
En esta oportunidad formularemos el Problema de

Valor Inicial PVI y analizamos e interpretamos
grficamente su solucin numrica, debemos
destacar que muchas de leyes generales de la
naturaleza se expresan con el lenguaje de las
ecuaciones diferenciales ordinarias que es aplicado
en una diversidad de campos del conocimiento. En
donde una ecuacin diferencial se debe considerar
como la razn de cambio de y con respecto a x.
5
PVI
1.En general una EDO de primer orden esta dado

por:
dy
f ( x, y ) ..(1)
dx
2. Tericamente se dice que la solucin de una EDO
debe contener una constante arbitraria C,
consecuentemente la solucin general de (1) es:
F ( x, y , c ) 0 (2)
6
PVI
Observaciones:
1.La relacin (2) representa una familia de curvas en
el plano xy, en donde cada curva se obtiene para un
valor particular de C.
2. Cada curva representa a una solucin particular de
EDO.
3. Las constantes C son obtenidos analticamente,
exigiendo que la solucin de esa ecuacin pase por
algn punto (x0, y0) esto es:
7
PVI
y ( x0 ) y 0
..(4)
i.e.: que y vale y0 cuando x es x0
Interpretacin Grficamente:
F3 = 0
Y0 F2 = 0, con Y(X0) = Y0
F1 = 0
8
PVI
4. Como se mencion al inicio la gran mayora de las

ecuaciones no pueden resolverse utilizando tcnicas
analticas, lo que obligan a estudiar mtodos numricos.
5. Debemos resaltar que cuando usamos los mtodos
numricos no encontramos soluciones de la forma
F(x,y,c) = 0 pues se trabajan con nmeros y se tiene
resultados numricos. Pero el propsito es determinar
valores de y que correspondan a valores especficos
de x los cual es factible con mtodos numricos.
9
PVI
El problema de valor inicial (P.V.I.) queda

formulado as:
dy
Una ecuacin diferencial de primer orden: dx f ( x, y)
i. Un valor de y en un punto conocido x0
(condicin inicial) y ( x0 ) y 0
ii. El valor xf es donde se quiere conocer el
valor de y(xf )
y (xf ) = yf
10
PVI
Matemticamente.
dy
f ( x, y )
dx
P.V.I. (5)
y ( x0 ) y 0
y( x f ) ?
11
MTODOS NUMRICOS PARA SOLUCIONAR UNA
EDO
MTODO DE EULER
Este mtodo consiste en dividir el intervalo [x0,xf] en n
subintervalos de ancho h esto es: X f X0
h
n
Lo que permite determinar un conjunto de n+1puntos

discretos, i.e. X0, x1 , x2, x3 xi ,xi+1 ... xn-1
x1 x2 x3 ... xi xi+1 ... xn-1 xn
x0 xf
12
EDO
Observando que:
Para cualquier punto se tiene
x1 x0 h x1 x0 h
x 2 x1 h x 2 x1 h x 2 x0 2h
x3 x 2 h x3 x 2 h x3 x0 3h
En general
xi x0 ih i 0,1,2,3,..., n
CONDICIN INICIAL
1. y ( x0 ) y 0 representa el punto P0 ( x0 , y 0 ) , por donde
pasa la curva solucin de la ecuacin PVI. lo que ser
denotado por F(x) = y, en lugar de F(x,y,c1) = 0. 13
EDO
2.Consecuentemente: teniendo el punto P0 podemos

evaluar la primera derivada de F(x) en ese punto P0. Esto
es:
dy
F ' ( x) f ( x0 , y 0 ) ...............................................(6)
dx P0
3. Teniendo esta informacin (6) trazamos una recta la

que pasa por P0 y de pendiente y y0
f ( x0 , y 0 ) : f ( x0 , y 0 ) : .......L3
x x0
que aproxima F(x) en una vecindad de X0.
14
EDO
4.Tomamos la recta L3 en lugar de F(x) y localizamos en

esta recta el valor de y1 que corresponde a x1. Esto es:
x y
x x
1
f (x , y )
0
.............................................(7)
0 0
1 0
x y
.............................................(8)
1 0
f ( x0 , y 0 ) y1 y 0 f ( x0 , y 0 )h
x1 x0
y 1 y 0 hf ( x 0 , y 0 )
y 2 y 1 hf ( x 1 , y 1 )
y1 F ( x1 )
.
.
y i 1 y i hf ( x i , y i )
.
.
y n y n 1 hf ( x n 1 , y n 1 )
15
EDO
Grfico
F(xf)
f(x1)
y1 error f(x0,y0)
y0
P0(x0,y0)
x0 x1 x3 x4 xi xn
x0 x1
xi xi+1
(1) En esencia se trata de aproximar la curva y = F(x) por

medio de una serie de segmentos de lneas rectas.
(2) El mtodo comete un error de truncamiento que es
propio de el.
16
EDO
(3) El error de (2) se puede anular tanto como se quiera,

reduciendo la longitud de h tericamente.
(4) Debido a (3) se comete un error de redondeo ms
alto.
Ejemplos.
Resolver PVI usando Euler
Ejemplo 1
dy
x y
f ( x, y ) x y
dx

y ( x0 ) y 0
y ( 0) 2
y (1) ? y( x f ) ?

17
EDO
Solucin
1. El intervalo de inters [x0,xf] = [0,1]
2. Determinando h: dividimos el intervalo [0,1] en 5
subintervalos 1 0
h 0 .2
5
xi x0 ih
3. Determinar los argumentos: x0 0
x1 x0 1h x1 0 1(0.2) 0.2
x 2 x 0 2 h x 2 0 2 ( 0 .2 ) 0 .4
x3 x0 3h x3 0 3(0.2) 0.6
x 4 x 0 4 h x 4 0 4 ( 0 .2 ) 0 .8
x5 x0 5h x5 0 5(0.2) 1
18
EDO
4. Determinando los valores de yi

y i 1 y1 hf ( xi , y i )
y1 y 0 hf ( x0 , y 0 ) y1 2 0.2 f (0.2) 2 0.2(0 2) 1.6
y 2 y1 hf ( x1 , y1 ) y 2 1.6 0.2 f (0.2,1.6) 1.6 0.2(0.2 1.6) 1.32
y 3 y 2 hf ( x 2 , y 2 ) y 3 1.32 0.2 f (0.4,1.32) 1.32 0.2(0.4 1.32) 1.136
y 4 y 3 hf ( x3 , y 3 ) y 4 1.136 0.2(0.6 1.136) 1.0288
y 5 y 4 hf ( x 4 , y 4 ) y 5 1.0288 0.2(0.8 1.0288) 0.98304
5. Comparando con la solucin analtica

La solucin analtica es: 1.10364
El error absoluto E y y 0.98304 1.10364 0.12060
A
*
5 5
El error relativo E
E
E
0.12060
R
A
0.1092
y5 R
1.10364
19
EDO
Ejemplo 2
Dada la siguiente ecuacin diferencial con la condicin
inicial:
Aproximar
NOTA
Primero observamos que esta ecuacin s puede
resolverse por mtodos tradicionales de ecuaciones
diferenciales. Por ejemplo, podemos aplicar el mtodo
de separacin de variables. Veamos las dos soluciones 20
EDO
Solucin Analtica.
Sustituyendo la condicin inicial:

Por lo tanto, tenemos que la curva solucin real est
dada
Y por lo tanto, el valor real que se pide es:
21
EDO
Solucin Numrica,
Aplicamos el mtodo de Euler y para ello, observamos
que la distancia entre y no es lo suficientemente
pequea. Si dividimos esta distancia entre cinco
obtenemos un valor de y por lo tanto, obtendremos la
aproximacin deseada en cinco pasos.
De esta forma, tenemos los siguientes datos:
22
EDO
Sustituyendo estos datos en la formula de Euler,

tenemos, en un primer paso:
Aplicando nuevamente la formula de Euler, tenemos, en

un segundo paso:
23
EDO
Y as sucesivamente hasta obtener . Resumimos los

resultados en la siguiente tabla
n
0 0 1
1 0.1 1
2 0.2 1.02
3 0.3 1.0608
4 0.4 1.12445
5 0.5 1.2144
24
EDO
Concluimos que el valor aproximado, usando el mtodo

de Euler es:
Puesto que en este caso, conocemos el valor verdadero,
podemos usarlo para calcular el error relativo
porcentual que se cometi al aplicar la formula de Euler.
Tenemos que:
25
EDO
Ejemplo
Aplicar el mtodo de Euler para aproximar , dada la
ecuacin diferencial.
Solucin
Nuevamente vemos que nos conviene dividir en pasos la
aproximacin. As, elegimos nuevamente para obtener
el resultado final en tres pasos. Por lo tanto, aplicamos
el mtodo de Euler con los siguientes datos:
26
EDO
En un primer paso, tenemos que:
Resumimos los resultados en la siguiente tabla:

n
0 1 2
1 1.1 2.3
2 1.2 2.6855
3 1.3 3.1901
27
EDO
MTODO DE TAYLOR
Podemos observar que el mtodo anterior usa los dos
primeros trminos de la serie de Taylor para su primera
iteracin, i.e.
F ( x1 ) y1 F ( x0 ) F ' ( x0 )( x1 x0 ) .(2)
De manera natural se puede pensar que para

determinar y2 se expandi de nuevo F(x) en la serie de
Taylor. As:
F ( x 2 ) y 2 F ( x1 ) F ' ( x1 )( x 2 x1 ) .......................................(2)
28
EDO
Pero se debe resaltar que no disponemos de los valores

exactos de F(x1) y F(x1), los que se usan en la expansin de
Taylor de F(x) alrededor de x1 lo que permite no evaluar la
parte derecha (2) consecuentemente para los otros valores
de x se usa:
y i 1 y i f ( xi , y i )( xi 1 xi ) ......................... (3)
y i 1 F ( xi ) F ' ( xi )( xi 1 xi )
La relacin (3) tiene mucha similitud con la expansin en serie Taylor.
Si aplicamos la informacin acerca de las series de Taylor con la
finalidad de mejorar la exactitud del mtodo de Euler, obtendremos
los llamados Algoritmos de Taylor.
29
EDO
Usemos tres trminos en lugar de dos en la expresin de

F(x1), i.e.
2
(x x )
F ( x ) y F ( x ) F ' ( x )( x x ) F ' ' ( x )
1 1 0 0 1 0 0
1 0
2!
(4)
Pero
dF ' ( x ) df ( x, y ) y h x1 x0
F ' ' ( x)
dx dx
Luego
h 2 df ( x, y )
y1 y 0 hf ( x0 , y 0 ) x0 , y 0
2! dx .(5)
30
EDO
Entonces se sugiere considerar (5) para obtener y2, y3,..., yn

mejorara la exactitud obtenida con (1) consecuentemente
se propone la formula:
h 2 df ( x, y )
y i 1 y i hf ( x, y ) xi , y i (6)
2! dx
La utilidad de la relacin (6) depende de cuan fcil sea la

diferenciacin de f(x,y)
31
EDO
Si f(x,y) es una funcin solo de x, la diferenciacin con

respecto a x es relativamente fcil y la formula propuesta
es muy prctica.
En general f(x,y) es una funcin de x , y, habr que usar

derivadas totales
La derivada total de f(x,y) con respecto a x esta dada por
df ( x, y ) f ( x, y ) f ( x, y ) dy

dx x y dx
32
EDO
Ejemplos. Resolver por el mtodo de Taylor

dy
dx x y

y (0) 2
y (1) ?

1. Clculo de: h = 0.2
2.Clculo de x x ih x
i 0 0 0 , x1 0.2 , x2 0.4 ,x3 0.6 ,x4 0.8 , x5 1
3. Aplicando h 2 df ( x, y )
y i 1 y i hf ( x, y ) ( xi , y i )
2! dx
h 2 df ( x0 , y 0 )
y1 y (0.2) y 0 h( x0 , y 0 )
2! dx
33
EDO
En donde
df ( x , y ) f ( x , y ) f ( x , y )
( x y ) 1 1( x y ) 1 x y ( x 0 , y 0 )
dx x y
h2
y1 y 0 h ( x 0 y 0 ) (1 x 0 y 0 )
2!
(0.2) 2
2 0.2(0 2) (1 0 2) 1.66
2
h2
y 2 y (0.4) y1 h( x1 y1 ) (1 x1 y1 )
2
0.2 2
1.66 0.2(0.2 1.66) (1 0.2 1.66) 1.4172
2
34
EDO
Continuando
0.2 2
y 3 1.4172 0.2(0.4 1.4172) (1 0.4 1.4172) 1.254104
2
0.2 2
y 4 1.254104 0.2(0.6 1.254104 (1 0.6 1.254104) 1.269184
2
0.2 2
y 5 1.269184 0.2(0.8 1.269184) (1 0.8 1.269184) 1.2047308
2
E A 1.010908
E R 0.915976
E % 9.15%
35
EDO
MTODO DE EULER MODIFICADO

En el mtodo de Euler se tom como vlida para todo el
intervalo la derivada encontrada en un extremo.
F(x0,y0)
Y0
X0 h X1
Si queremos obtener una exactitud razonable se toma h

muy pequea, a cambio de un mayor error de redondeo
36
EDO
El mtodo presente trata de evitar tal problema

utilizando un valor promedio de la derivada tomada
en los extremos del intervalo. Constado de 2 pasos:
1 Se inicia de (x0,y0), usar el mtodo de Euler para

determinar y correspondiente a x1, valor que ser
denotado por 1y , puesto que se trata de un valor
transitorio de y1. Este paso se le llama paso predictor.
37
EDO
2 Este paso se llama corrector, pues trata de corregir

la prediccin en el nuevo punto ( x , y ) se evala la
1 1
derivada f ( x , y ) usando la ecuacin diferencial ordinaria

1 1
P.V.I. que se est resolviendo, se obtiene la media

aritmtica de esta derivada y la derivada en el punto
inicial (x0,y0)
Derivada Promedio = 12 f ( x , y ) f ( x , y )
0 0 1 1
Usamos la derivada promedio para calcular el nuevo

valor y1 con la ecuacin de Euler, que ser mas exacto
que y 1
38
EDO
x1 x0
y1 y 0 f ( x0 , y0 f ( x1 , y1 )
2
Que ser el valor definitivo de y1. El proceso se repite
hasta llegar a yn.
Primero: Paso de Prediccin y i 1 y i hf ( xi , y i )
Segundo: Una vez obtenida y i 1 se calcula f ( x , y ) , la i 1 i 1
derivada en el punto ( xi 1 , y i 1 ) y se promedia con la

derivada previa f ( x , x ) para encontrar la derivada
i i
promedio 1
f x , y f x , y
i 1 i 1
Derivada Promedio 2 i i
39
EDO
Solucin
Considerando las mismas condiciones del ejercicio
tenemos:
h=0.2; y0=2; f(x0,y0)=f(0,2)=0-2=-2
Primera iteracin
y 1 y 0 hf ( x0 , y 0 ) 2 0.2(0 2) 1.6
1
2

f ( x 0 , y 0 ) f ( x1 , y1 ) (0 2) (0.2 1.6) 1.7
1
2
derivada promedio
y1 y 0 0.2(1.7) 2 0.2(1.7) 1.66
40
EDO
Segunda interacin
1 y 2 y1 hf ( x1 , y1 ) 1.66 0.2(0.2 1.66) 1.368
2 1 f ( x , y ) f ( x , y ) 1 (0.2 1.66) (0.4 1.368) 1.214
1 1 2 2
2 2
y ( x 2 ) y 2 1.66 0.2(1.214) 1.4172
Tercera interacin
1 y 3 y 2 hf ( x 2 , y 2 ) 1.4172 0.2(0.4 1.4172) 1.21376
2
1
f ( x2 , y 2 ) f ( x3 , y3 ) 1 (0.4 1.4172) (0.6 1.21376)
2 2
24/07/2015 41
EDO
Ejemplo
Aplicar el mtodo de Euler mejorado, para aproximar
si:
Solucin
Vemos que este es el mismo ejemplo 1 del mtodo
anterior. As que definimos y encontraremos la
aproximacin despus de cinco iteraciones.
42
EDO
A diferencia del mtodo de Euler 1, en cada iteracin

requerimos de dos clculos en vez de uno solo: primero el
de y posteriormente el de .
Para aclarar el mtodo veamos con detalle las primeras dos
iteraciones. Primero que nada, aclaramos que tenemos los
siguientes datos iniciales
43
EDO
En nuestra primera iteracin tenemos:
Ntese que el valor de coincide con el (Euler 1), y es

el nico valor que va coincidir, pues para calcular se
usar y no
Esto lo veremos claramente en la siguiente iteracin:
44
EDO
Ntese que ya no coinciden los valores de (Euler 1) y el de El

proceso debe seguirse hasta la quinta iteracin. Resumimos los
resultados en la siguiente tabla: n
0 0 1
1 0.1 1.01
2 0.2 1.040704
3 0.3 1.093988
4 0.4 1.173192
5 0.5 1.28336
24/07/2015 45
EDO
METODO DE RUNGE-KUTTA
METODO DE RUNGE-KUTTA DE SEGUNDO ORDEN
Estos mtodos que se encuentran relacionados a los
nombres de Runge (1885), Kutta (1901), Heun (1900) y
otros, para solucionar P.V.I .Consiste en obtener un
resultado que se obtendr al utilizar un nmero finito de
trminos de una serie de Taylor de la forma:
h2 h3
yi 1 yi h. f ( xi , yi ) f ' ( xi , y i ) f ' ' ( xi , yi ) ...
2! 3! (1)
46
EDO
Con una aproximacin en la cual se calcula yi1 de una

formula del tipo:
0 f (x, y) 1 f (xi u1h, yi b1h) 2 f (xi u2h, yi b2h) ...

yi1 yi h (2)
p f (xi u p h, yi bp h)
En donde:
, u, b son determinados de modo que si se expandiera con , en
serie de Taylor alrededor de ( xi ,yi ); debemos observar que los
coeficientes de h, h2, h3, etc., coincidiran con los coeficientes de la
ecuacin (1). Supongamos p=1 tendremos
47
EDO
yi 1 yi h 0 xi ; y i 1 . f ( xi uhi ; yi bh) ..............(3)

Observaciones:
1. En esta relacin se evala f(x), en x ; y ( x uh ; y bh) donde
i i i i i
es tal que: xi xi uh xi 1 , para mantener la abscisa del

segundo punto dentro del intervalo de inters, con lo que
. 0 u 1 Grficamente
yi+1
(xi+uh , yi+k0)
yi+1+h f( xi , yi )
(xi,yi)
xi+1
xi
48
EDO
2. b puede ser manejado ms libremente y expresarse y

se puede usar como ordenada arriba o debajo de la
ordenada que da el mtodo de Euler simple.
yi bh yi hf ( xi ; yi ) yi k 0(4)
Con k0 = h f(xi,yi)
3.Queda por determinar 0, 1, , tal que la ecuacin (3)
tenga una aproximacin en potencias de h, cuyos primeros
trminos coinciden con los primeros trminos de ecuacin
(1).
4.Para cumplir con (3) expandimos primero en serie de
Taylor 49
EDO
f ( xi uh, yi k 0 )
f f u 2 h 2 2 f 2 f 2 k 02 2 f
f ( xi yi ) uh k 0
x y

2! x 2
uhk 0
xy

2!y 2
0 h3 .(5)
Todas las derivaciones son evaluadas en x , y i i
Sustituyendo en la ecuacin (3)

yi1
f f u2h2 2 f 2 2k02 2 f 3
yi 0hf (xi , yi ) 1h f (xi , yi ) uh k0
x y 2! 2 y
uhk0
xy 2! y2
0 h

Arreglando en potencias de h, tenemos

f f
y i 1 y i h 0 1 f ( x i , y i ) h 2 1 u f ( xi , yi )
x y
2 2 f 2 f 2 f

h3
3 u
x 2
2u f ( xi , yi )
xy
f ( xi , yi )
2 2
y 2
..(6)
0 h4
2
50
EDO
Para que los coeficientes correspondientes de h, h2

coincidan en las ecuaciones (1) y (6) se requiere que:
1 1
1 u
2 ..(7)
0 1 1 , 1
2
5.Observamos que existen 4 incgnitas para solo tres
ecuaciones y, por tanto se tiene un grado de libertad en la
solucin de la ecuacin (7). Podramos pensar en usar este
grado de libertad para hacer coincidir los coeficientes de
h3. Sin embargo, es obvio que esto es imposible para
cualquier forma que tenga la funcin f(x,y). Existe
entonces un nmero de infinito de soluciones de la
ecuacin (7), pero quizs la ms simple sea : 12 ; u 1
0 1
51
EDO
6. La relacin de (5) conduce a la formula
yi 1 yi
h
f ( xi , yi ) f ( xi h, yi hf ( xi , yi ))
2
o bien
y i 1 y
h
k 0 k1 , con : k 0 f ( xi , yi ) ; k1 f ( xi h, yi hk 0 ) .(8)
2
7. La relacin (8) es conocida como algoritmo de

Runge-Kutta de segundo orden.
52
EDO
Lo de segundo orden por coincidir con los tres primeros

trminos de la serie de Taylor que es la formula de Euler
Modificado.
Este mtodo proporciona mayor exactitud que la de
Euler.
Se puede usar un valor de h no tan pequeo como el
primero .El precio de es la evaluacin f(x,y) dos veces en
cada subintervalo contra uno en el mtodo de Euler.
8.Las formulas de Runge-Kutta de cualquier orden se
puede derivar de manera anloga que la de segundo
orden.
53
EDO
METODO DE RUNGE-KUTTA DE CUARTO ORDEN

y y k 2k 2k k
h .(9)
I 1 i 1 2 3 4
6
k1 f ( xi , yi )
h hk1
k 2 f ( xi , y i )
2 2
h hk
k 3 f ( xi , y i 2 )
2 2
k 4 f ( xi h, yi hk 3 )
54
EDO
9. La ecuacin (9) tiene mucha coincidencia con los

5 primeros trminos de la serie de Taylor lo que
significa gran exactitud sin calculo de derivadas,
pero a cambio, se tiene que evaluar la funcin
f(x,y)cuatro veces en cada subintervalo.
55
EDO
Ejemplo 1
dy
dx x y

P.V .I y (0) 2
y (1) ?

Usando Runge-Kutta de cuarto orden.
56
EDO
Solucin:
Primera Iteracin: Calculo de constantes k1, k2, k3, k4
k1 f ( x0 , y0 ) x0 y0 0 2 2
h hk h hk 0 .2
k 2 f ( xi , y i 1 ) f ( x0 , y 0 1 ) f (0 , 2 0 .2 )
2 2 2 2 2
0 .2
2 0 .2 1 .7
2
h hk h hk 0 .2 0.2(1.7)
k 3 f ( xi , y i 2 ) f ( x 0 , y 0 2 ) f (0 ,2 )
2 2 2 2 2 2
0 .2 0.2(1.7) 10 200 17
2 1.73
2 2 100 100 100
57
EDO
k 4 f ( xi h, yi hk 3 ) f ( x0 h, y0 hk 3 ) f (0 0.2,2 0.2(1.73))
173
0.2 2 1.454
1000
Clculo De y1
y1 y 0
h
k1 2k 2 2k 3 k 4 2 0.2 2 3.4 3.46 1.454 1.6562
6 6
58
EDO
Segunda Iteracin: Calculo de constantes k1, k2, k3, k4

k1 f ( x1 , y1 ) f (0.2,1.6562) 0.2 1.6562 1.4562
h hk 0.2 0.2(1.7)
k 2 f ( x1 , y1 1 ) f (0.2 ,1.6562 )
2 2 2 2
0.2 0.2(1.7)
0.2 1.6562 1.21058
2 2
h hk 0 .2 0.2(1.21058)
k3 f ( x1 , y1 2 ) 0.2 1.6562 1.235142
2 2 2 2
k 4 f ( xi h, yi hk3 ) 0.2 0.2 1.6562 0.2(1.235142) 10091716

Clculo De y2:
y2 y1
h
k1 2k 2 2k3 k 4 1.6562 0.2 1.4562 2(1.2128)... 1.4109
6 6
59
EDO
Continuando se tiene
y3 1.246450474
y4 1.148003885
y5 1.103655714
Observacin:
Los mtodos descritos se llaman tambin mtodos de
un solo paso porque se apoyan y usan (xi,yi) para el
clculo de yi+1.
Estos Mtodos adems se apoyan en puntos xi y xi+1
pero nunca en puntos anteriores a xi.
60
EDO
Ejemplo 2
Usar el mtodo de Runge-Kutta para aproximar dada la
siguiente ecuacin diferencial:
Solucin
Primero, identificamos el mismo ejemplo 1 de los dos
mtodos anteriores. Segundo, procedemos con los
mismos datos:
61
EDO
Para poder calcular el valor de y1 debemos calcular

primeros los valores de k1, k2 ,k3, y k4. Tenemos entonces
que:
62
EDO
Con el fin de un mayor entendimiento de las frmulas,

veamos la siguiente iteracin:
63
EDO
El proceso debe repetirse hasta obtener lo requerido.

Resumimos los resultados en la siguiente tabla:
n
0 0 1
1 0.1 1.01005
2 0.2 1.04081
3 0.3 1.09417
4 0.4 1.17351
5 0.5 1.28403
64
EDO
Concluimos que el valor obtenido con el mtodo de

Runge-Kutta es:
Finalmente, calculamos el error relativo verdadero:
Con lo cual vemos que efectivamente se ha reducido

muchsimo el error relativo. De hecho observamos que
tenemos 6 cifras significativas en la aproximacin
65
EDO
MTODOS DE PREDICTOR-CORRECTOR
Recordemos que en el mtodo de Euler modificado se utiliza
la siguiente relacin
y i 1 y i f xi , y i f xi 1 , y i 1 .(1)
h
2
Obsrvese, que el segundo trmino del miembro de la
derecha recuerda el mtodo de integracin trapezoidal
compuesta, en donde h es el ancho del trapezoide h=xi+1 xi,
y podemos decir que,
(2)
66
EDO
Equivalentemente
..(3)
Que es la ecuacin de correccin del mtodo de Euler

modificado, esto sugiere la obtencin de un esquema
iterativo para la solucin del PVI por medio de la regla de
Simpson u otro mtodo de integracin numrica que usan
mayor numero de puntos.
Considerando esta reflexin se deriva un mtodo
corrector basado en el mtodo de Simpson 1/3
24/07/2015 67
EDO
,....................(4)
Considerando la relacin
,...........................(5)
Tenemos
,..(6)
68
EDO
Entonces se llega a la relacin de correccin,

.......(7)
En donde se debe de obtener con un predictor, a partir
de (x0,y0 ) la ultima relacin tomara la forma de,
,..(8)
Para la primera prediccin es calculada con un predictor

que requiere de y1 y f(x1,y1 ) en consecuencia se requiere
de un paso de inicializacin que muy ben puede ser usado
el mtodo de Runge-Kutta por una sola vez en el proceso
iterativo. 69
EDO
Ejemplo:
Resolver el PVI
dy
dx x y

P.V .I y (0) 2
y (1) ?

Usar el mtodo de prediccin y correccin
70
EDO
Solucin
h=(1-0)/5=0.2,
Primera iteracin
Inicializacin. (Usando Euler modificado obtenemos y1 )
1 y y hf ( x , y ) 2 0.2(0 2) 1.6
1 0 0 0
2 12 f ( x , y ) f ( x , y ) 12 (0 2) (0.2 1.6) 1.7

0 0 1 1
derivada promedio. Luego

y1 y 0 0.2(1.7) 2 0.2(1.7) 1.66
71
, EDO
Prediccin (se usa Euler Modificado para tomar el valor y2)

1 y 2 y1 hf ( x1 , y1 ) 1.66 0.2(0.2 1.66) 1.368
2
1
f ( x1 , y1 ) f ( x2 , y 2 ) 1 (0.2 1.66) (0.4 1.368) 1.214
2 2
y ( x 2 ) y 2 1.66 0.2(1.214) 1.4172
Correccin; usamos la relacin 8
72
, EDO
,
Segunda Iteracin
Prediccin
,
Correccin usamos la relacin 7
73
, EDO
,
Tercera Iteracin
Prediccin
74
3.2. MTODOS NUMRICOS PARA
, SOLUCIONAR UNA EDO
,
Cuarta Iteracin
Prediccin
75
EJERCICIOS
I. Utilizar los mtodos de Euler y de Runge Kutta para

dar solucin a las siguientes ecuaciones diferenciales
con valor frontera.
dy dy x
x y
dy
x y
dy x
dx y
dx dx dx y
(9) (10) y(1) 4 (11) y (0) 1 (12)
y (0) 4
y (0) 2
y (1.5) ? y (0.5) ?
y (1) ? y (0.5) ?
dy dy
dy x y x y2 dy y 2 y
y (2 y )
(13) dx (14) dx (15) dx (16) dx x
y (1) 4 y (1) 0
y ( 0) 3 y (1) 1
y (1.5) ? y (1.5) ? y (1.8) ?
y (0.5) ?
76
EJERCICIOS
,
I. Utilizar los mtodos de Euler y de Runge Kutta para

dar solucin a las siguientes ecuaciones diferenciales
con valor frontera.
dy dy dy dy
1 xsenx 1 y 2 y 2 y 1
dx dx dx dx
(1) y (0) 0
(2) y (0) 4
(3) y ( 0) 1
(4)
y (0) 1
y (1.5) ? y (1) ? y (1) ? y (1) ?
dy 1 y dy
y2 dy
1 y x 1 y dy
1 xy
(5) dx x x (6) dx (7) dx (8) dx
y (1) 1 y (0) 1
y (1) 1
y ( 0) 0
y (2) ? y (1) ?
y (1) ? y (2) ?
77
EJERCICIOS
II .- Estructurar un modelo para las problemticas

siguientes y luego solucionarlo Aplicando Euler y Runge
Kuta.
1.- Un tanque cilndrico de fondo plano con dimetro 2
metros contiene un lquido; de densidad 1.8 kg/l a una
altura H de 4 metros. Se desea saber la altura del lquido
dentro del tanque 10 minutos despus que abre
completamente de la vlvula de salida ubicada en la parte
inferior izquierda, la cual da una gasto de 1 m3/s, donde A
es el rea seccional del tubo de salida que tiene un valor
de 80.5 x 10-4m2, considerar g = 9.81m/s2.
78
EJERCICIOS
2.- Se tiene un tanque esfrico de radio de 8 metros calcular el tiempo

necesario para que el nivel del lquido de dicho tanque pase de 6
metros a 7 metros, la velocidad de salida por el orificio del fondo es v
=5.5 m/s el dimetro de dicho orificio es de 12 cm. Donde a es la
altura de lquido.
3.- En un tanque perfectamente agitado se tiene 500 litros de una
salmuera en la cual este disuelto 30 Kg de sal comn en un momento
determinado se hace llegar al tanque un gasto de 90 l/min de una
salmuera que contiene 1.5 Kg de sal comn por litro si se tiene un
gasto de salida de 90 l/min. Determine:
a.- Que cantidad de sal hay en el tanque transcurrido 20min.
b.- Que cantidad de sal transcurrido un tiempo muy grande.
79
EJERCICIOS
4.- Se hace reaccionar isotrmica mente 300gr de acetato de etilo con 200gr de
hidrxido de sodio en solucin acuosa ajustando el volumen total a 10 litros para
dar acetato de sodio y alcohol etlico de acuerdo con lo siguiente ecuacin
estequiometria:
Acetato de etilo + hidrxido de sodio = acetato de sodio + alcohol etlico
Donde la constante de velocidad de reaccin k esta dado por k = 1.44 x 10-2
Determine la cantidad de acetato de sodio y alcohol etlico presente 40min despus
presentada la reaccin.
5.- Se conecta un inductor de 0.5 henries en serie con una resistencia de 10 ohms
un capacitador de 0.025 faradios y un generador de corriente al terna dad por la
funcin 60 sen 5t voltios t 0.
a.- Establezca una ecuacin diferencial para la carga instantnea en el capacitor.
b.- Encuentre la carga en distintos tiempos
80
EJERCICIOS
6.- Se tiene un tanque de forma cnica de 5 metros de dimetro

superior con 10 metros de altura conteniendo un lquido hasta h
metros de altura, si al momento de llegar el nivel del lquido de 2.
5 metros se hace llegar un gasto de alimentacin de 0.50 m3/s el
nivel de lquido aumentara. Determine el tiempo necesario para
que el nivel se recupere nuevamente a 6 metros.
7.- El tiempo que requiere el tanque del ejercicio anterior para

recuperar su nivel de 2.5 a 6 metros con un gasto de alimentacin
de 0.50 m3/s es aproximadamente 500 s calcule el gasto de
alimentacin que se requiere para reducir este tiempo en la mitad.
81
Muchas Gracias

Catedra Metodos Numericos 2015 Unsch 15 PDF

Cargado por

Copyright:

Formatos disponibles

Catedra Metodos Numericos 2015 Unsch 15 PDF

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

Catedra Metodos Numericos 2015 Unsch 15 PDF

Cargado por

Copyright:

Formatos disponibles

CATEDRA 1

Facultad de Ingeniera de Minas, Geologa y Civil

FORMULACIN DEL PROBLEMA DE VALOR INICIAL

En esta oportunidad formularemos el Problema de

1.En general una EDO de primer orden esta dado

4. Como se mencion al inicio la gran mayora de las

El problema de valor inicial (P.V.I.) queda

Lo que permite determinar un conjunto de n+1puntos

2.Consecuentemente: teniendo el punto P0 podemos

3. Teniendo esta informacin (6) trazamos una recta la

que aproxima F(x) en una vecindad de X0.

4.Tomamos la recta L3 en lugar de F(x) y localizamos en

(1) En esencia se trata de aproximar la curva y = F(x) por

(3) El error de (2) se puede anular tanto como se quiera,

4. Determinando los valores de yi

5. Comparando con la solucin analtica

Sustituyendo la condicin inicial:

Y por lo tanto, el valor real que se pide es:

Sustituyendo estos datos en la formula de Euler,

Aplicando nuevamente la formula de Euler, tenemos, en

Y as sucesivamente hasta obtener . Resumimos los

Concluimos que el valor aproximado, usando el mtodo

En un primer paso, tenemos que:

Resumimos los resultados en la siguiente tabla:

De manera natural se puede pensar que para

Pero se debe resaltar que no disponemos de los valores

Usemos tres trminos en lugar de dos en la expresin de

Entonces se sugiere considerar (5) para obtener y2, y3,..., yn

La utilidad de la relacin (6) depende de cuan fcil sea la

Si f(x,y) es una funcin solo de x, la diferenciacin con

En general f(x,y) es una funcin de x , y, habr que usar

Ejemplos. Resolver por el mtodo de Taylor

MTODO DE EULER MODIFICADO

Si queremos obtener una exactitud razonable se toma h

El mtodo presente trata de evitar tal problema

1 Se inicia de (x0,y0), usar el mtodo de Euler para

2 Este paso se llama corrector, pues trata de corregir

derivada f ( x , y ) usando la ecuacin diferencial ordinaria

P.V.I. que se est resolviendo, se obtiene la media

Usamos la derivada promedio para calcular el nuevo

derivada en el punto ( xi 1 , y i 1 ) y se promedia con la

y1 y 0 0.2(1.7) 2 0.2(1.7) 1.66

A diferencia del mtodo de Euler 1, en cada iteracin

En nuestra primera iteracin tenemos:

Ntese que el valor de coincide con el (Euler 1), y es

Ntese que ya no coinciden los valores de (Euler 1) y el de El

Con una aproximacin en la cual se calcula yi1 de una

0 f (x, y) 1 f (xi u1h, yi b1h) 2 f (xi u2h, yi b2h) ...

yi 1 yi h 0 xi ; y i 1 . f ( xi uhi ; yi bh) ..............(3)

es tal que: xi xi uh xi 1 , para mantener la abscisa del

2. b puede ser manejado ms libremente y expresarse y

Sustituyendo en la ecuacin (3)

Arreglando en potencias de h, tenemos

Para que los coeficientes correspondientes de h, h2

6. La relacin de (5) conduce a la formula

7. La relacin (8) es conocida como algoritmo de

Lo de segundo orden por coincidir con los tres primeros

METODO DE RUNGE-KUTTA DE CUARTO ORDEN

9. La ecuacin (9) tiene mucha coincidencia con los