ECUACIONES DIFERENCIALES
ORDINARIAS
Mtodo de Euler
Ejemplo 1: Aplicar el mtodo de Euler para estimar un valor aproximado de la
solucin en
x = 1,5 de y = x + y2, y(1) = 0, siendo z0 = 0 y el tamao de paso 0,1; 0,05 y
0,025.
Solucin:
Para h = 0,1
Se construye una tabla para evaluar zn+ 1 = zn + hf(xn, zn), siendo xn = x0 + nh.
As, 1,5 = 1 + 0,1n, por lo que n = 5.
n
xn
zn+
f(xn, zn)
1,1
0,1
1,11
1,2
0,211
1,244521
1,3
0,3354521
1,41252811
1,4
0,47670491
1,62724757
1,5
0,63942967
z5 = 0,63942967 y(1,5).
Para h = 0,05, a mano o con una hoja de clculo, se construye una tabla para
evaluar zn+1 = zn + hf(xn, zn), siendo xn = x0 + nh. As, 1,5 = 1 + 0,05n, por lo que
n = 10.
n
xn
zn
f(xn, zn)
1,05
0,05
1,0525
1,1
0,102625
1,11053189
1,15
0,15815159
1,17501193
1,2
0,21690219
1,24704656
1,25
0,27925452
1,32798309
1,3
0,34565367
1,41947646
1,35
0,4166275
1,52357847
1,4
0,49280642
1,64285817
1,45
0,57494933
1,78056673
10
1,5
0,66397766
Ahora z10 = 0,66397766 y(1,5).
Para h = 0,025 se tiene que n = 20 y:
xn
zn
f(xn, zn)
1,025
0,025
1,025625
1,05
0,05064063
1,05256447
1,075
0,07695474
1,08092203
1,1
0,10397779
1,11081138
1,125
0,13174807
1,14235755
1,15
0,16030701
1,17569834
1,175
0,18969947
1,21098589
1,2
0,21997412
1,24838861
1,225
0,25118383
1,28809332
10
1,25
0,28338616
1,33030772
11
1,275
0,31664386
1,37526333
12
1,3
0,35102544
1,42321886
13
1,325
0,38660591
1,47446413
14
1,35
0,42346752
1,52932474
15
1,375
0,46170063
1,58816748
16
1,4
0,50140482
1,65140679
17
1,425
0,54268999
1,71951243
18
1,45
0,5856778
1,79301849
19
1,475
0,63050326
1,87253437
20
1,5
0,67731662
Ahora z20 = 0,67731662 y (1,5)
Ejemplo 2: Aplicar el mtodo de Euler para estimar la solucin en x = 1 de
y = 1 x + 4y, y(0) = 1, siendo z0 = 1 y los tamaos de paso siguientes: 0,1; 0,05,
0,025, 0,0125. Resolver la ecuacin diferencial. Valorar el error cometido y
compararlo con el tamao de paso utilizado.
Solucin:
Resultados de aplicar el mtodo de Euler con distintos tamaos de paso a:
y = 1 x + 4y, y(0) = 1, para aproximar la solucin en x = 1
h = 0,1
h = 0,05
h = 0,025
h = 0,0125
Exacto
34,411490
45,588399
53,807866
60,037126
Error
global
30,48...
19,309...
11,09...
4,86...
64,897803
Ejemplo 3: Aplicar el mtodo de Euler para calcular el valor aproximado en x =
2 de la solucin del problema de valor inicial
y'=
1
2 y 2
2
1+ x
Usando como y(0)= 0 y tamaos de paso h = 0,2. Hallar el error.
Solucin:
Resolviendo con matlab
En este caso n=10
Tenemos como resultado aproximado de y (2)=0.4068
Ejemplo 4: aplicar el mtodo de Euler para aproximar la solucin del problema
de valor inicial
y ' =cos 2 x+ sen 3 x ; 0 x 1 ; y ( 0 )=1, con h=0.25
La solucin exacta de este problema de valor inicial es
1
1
4
y ( x )= sen 2 x co 0 s 3 x+
2
3
3
En cada paso calcule el error de las aproximaciones.
Solucin:
x i=x 0+(i)h
w i+1=wi +h(cos ( 2 x i ) + sen ( 3 xi ) )
N
puntos
i
0
1
2
3
4
(xi,wi) puntos de sol.
aprox.
xi
Wi (aproximada)
0
1
0,25
1,25
0,5
1,63980533
0,75
2,02425465
1
2,23645725
Ejemplo 5:
y '=senxln y
y (0.13)=0.32
y ( 0.14 )=? ? ?
El nmero de nodos de integracin es 4
Solucin:
h=
0.140.13
=0.0025
4
Resolviendo con matlab
(xi,yi) puntos de sol.
exacta
Yi (exacta)
ERROR
1
0
1,329149813
0,07914981
1,730489759
0,09068443
2,041472034
0,01721738
2,117979546 -0,11847708
Mtodo de Heun (Euler mejorado)
Ejemplo1 : aplicar el mtodo de Euler modificado a la EDO
'
y =x+ y
Y(0)=1
Utilizando un incremento de longitud h=0.2
Solucin:
x i=x i1 +h , i=1,2, . n
z j+1 = y j +hf ( x j , y j ) , j=0,1,2,
y j +1= y j +
h
f ( x j , y j ) + f ( x j+1 , z j+1 ) ]
2[
xn
0,0
0,2
0,4
0,6
0,8
1,0
yn
1,00000
1,24000
1,57680
2,03170
2,63067
3,40542
exacto
1,00000
1,24000
1,58000
2,04424
2,65108
3,43656
Error (%)
0,00
0,23
0,43
0,61
0,77
0,91
El valor aproximado obtenido para y (1) es 3.40542
Ejemplo 2:
'
2 y +3 yt + y=0
y0=1.2
h=0.3
y1=1.2
Solucin:
y'=
3 yt y
2
t0=0
t1=t0+h=0.3
y 'i+1= yi +
Utilizando
h
f ( y i ,t i ) + f ( y i+1 ,t i+1 ) ]
2[
Reemplazar
y '1=1,2+
0.3 3 ( 1.2 ) ( 0 ) 1.2 3 ( 1.2 ) ( 0.3 )1.2
+
= 0.939
2
2
2
Y2=??
Y2=0.939
Y0=1.2
T1=0.3
T2=t1+0.3=0.6
y '2=1,2+
Ejemplo 3:
Aproximar el valor de y (1)
y ' =x+
0.3 3 (1.2 ) ( 0.3 ) 1.2 3 ( 0.939 ) ( 0.6 ) 0.939
+
=0.83181
2
2
2
1
y , y ( 0 )=3 , h=0.2 , x [ 0,1 ]
5
Solucin:
X1= 0 + 0.2 = 0.2
1
z 1=3+ 0.2 0+ (3 ) =3.12
5
y 1=3+
0.21
1
1
0+ (3 ) + 0 .2+ (3 .12 ) =3.1024
2
5
5
((
)(
Seguir resolviendo con matlab
El valor aproximado para y (1) es -3.1305.
))
Ejemplo 4: resolver el problema de valor inicial
y'=
t y
en [ 0,3 ] y ( 0 ) =1
2
Con 12 iteraciones
Solucin.
h=
ba 30
=
=0.25
n
12
Utilizando matlab
Ejemplo 5: usar el mtodo de Euler mejorado para obtener el valor aproximado
de y(0.5) para la solucin del problema de valor inicial
y ' =( x + y +1)2
y (0)=2, h=0.05
Solucin:
Mtodo de RK-2
Ejemplo 1 : El circuito de la figura tiene una auto inductancia de L=50 H , una
resistencia de 20 Ohm y una fuente de voltaje de V =10 volt . Si el interruptor
se cierra en el instante
L
t=0 , la corriente
I (t)
d
I ( t ) + RI ( t ) =E , I ( 0 )=0
dt
Determinar la corriente para
con h=0.1
0 t 10
d
R E
I=
+ f (I , t)
dt
L L
Ahora desarrollando para RK-2 se tiene:
R
E
I n+
L
L
k 2=hf ( y n+ k 1 ,t n+1 ) =h
R
E
( I + k )+
L n 1 L
1
I n+1=I n + ( k 1 +k 2 )
2
Los clculos para los dos primeros pasos son:
n=0,(t=0.1)
(5)
segundos, mediante el mtodo de RK-2,
SOLUCIN:
En primer lugar debemos reescribir la ecuacin (5)
k 1=hf ( y n , t n )=h
satisface la ecuacin
k 1=0.1
10
( 0 )+ ( ) =0.1 [ (0.4 ) ( 0 ) + ( 0.2 ) ] =0.02
([ 20
)
50
50 ]
k 2=0.1 [ (0.4 )( 0+ 0.02 )+ 0.2 ] =0.0192
1
I 1 =I 0 + ( 0.02+0.0192 )=0.0196
2
n=1,(t=0.2)
k 1=0.1 [ (0.4 )( 0.0196 )+(0.2) ]=0.019216
k 2=0.1 [ (0.4 )( 0.0196+ 0.019216 ) +0.2 ] =0.018447
1
1
I 2 =I 1 + ( k 1+ k 2 )=0.0196+ ( 0.019216+ 0.018447 )=0.0348431
2
2
Desarrollando con matlab
Donde:
X=tiempo
Y=corriente
Ejemplo 2: utilizar el mtodo de RK-2 para aproximar el valor de la solucin de la
siguiente ecuacin diferencial en los puntos x= 0.2, 0.4, 0.6, 0.8 y 1, usando
h=0.2
y ' =2 x+ y ; y ( 0 )=1
Solucin:
h=
10
=0.2
n
n=5
Ejemplo 3:
Resolver el problema de valor inicial
y'=
y , y ( 0 )=4
2 x +1
Con h=0.5. Hallar y (2)
Solucin:
Ejemplo 4:
Usar el mtodo de RK-2 para obtener el valor aproximado de y(0.5) para
la solucin del problema de valor inicial
y ' =( x + y 1)2
y ( 0 )=2
Tomando h=0.05
Solucin:
Ejemplo 5:
Aproximar la solucin en x=1 por el mtodo de RK-2 con tamao de paso
de 0.25 y 0.125 de
y ' =x 2+ y 2
y(0)=1.
Solucin:
h=
10
=0.25
n
n=4
h=
10
=0.125
n
n=8
MTODO DE RK-4
Ejemplo 1: Resolver el PVI
dy
=( x y ) , y ( 0 )=2 , aproximar
dx
y (1) con el mtodo RK-4
SOLUCIN:
El intervalo es [0,1] al dividirlo en cinco sub intervalos se tiene
h=
10
=0.2 .
5
Entonces:
x 0=0.0, x1 =x0 + h=0.00+ 0.2=0.2
x 2=x 1+h=0.2+0.2=0.4
x 3=x 2+ h=0.4 +0.2=0.6
x 4=x 3 +h=0.6+0.2=0.8
x 5=x 4 +h=0.8+0.2=1.0
Ahora calculamos las constantes
k1 , k2 , k3 , k4
k 1=f ( x 1 , y 1 )=( 0.21.6562 )=1.4562
) [(
)(
) [(
)(
hk
h
0.2
1.14562
k 2=f x 1 + , y 1 + 1 = 0.2+
1.6562+0.2
2
2
2
2
hk
h
0.2
1.21058
k 3 =f x 1 + , y 1 + 2 = 0.2+
1.6562+0.2
2
2
2
2
))]=1.21058
))]=1.235142
k 4=f ( x1 +h , y 1 +h k 3 ) =[ ( 0.2+0.2 ) + ( 1.6562+0.2 (1.235142 ) ) ] =1.0091716
Clculo de
y2
y ( 0.4 )= y 2= y 1+
1.6562+
h
[ k +2 k 2 +2 k 3 +k 4 ]
6 1
( 0.26 ) (1.4562+2 (1.21058 )+ 2(1.235142 )1.0091716)=1.410972813
Continuando el procedimiento, en las siguientes iteraciones se tiene
y ( 0.6 )= y 3 =1.246450474
y ( 0.8 )= y 4 =1.148003522
y (1 ) = y 5=1.103725867
Resolviendo con matlab
Ejemplo 2: Aplicar el metodo de Runge-Kutta de orden cuatro para
calcular el valor aproximado de x(1) e y(1) en el problema de
valores iniciales
x ' ( t )=f ( t , x , y )=4 y +cost ; x ( 0 )=0
'
y ( t )=g ( t , x , y )= x ; y ( 0 )=0
Tomando h=0.1
Solucin:
Se utiliza
tk
xk
yk
tk
xk
yk
0.0
0.00000
0.00000
0.6
0.43314
0.15432
0.1
0.09917
0.00498
0.7
0.44223
0.19829
0.2
0.19339
0.01967
0.8
0.42726
0.24196
0.3
0.27792
0.04333
0.9
0.38813
0.28293
0.4
0.34843
0.07478
1.0
0.32571
0.31881
0.5
0.40117
0.11242
Las aproximaciones de x(0) y y(0) son 0.32571 y 0.31881 respectivamente
Ejemplo 3: Aplicar el metodo de Runge-Kutta de orden cuatro con h =
0.1 para obtener un valor aproximado de y(1) en el siguiente
problema de valor inicial.
y ' =2ty , y ( 0 )=1
Solucin:
yk
ERROR
tk
0.0
0.1
1.00000
0.99004
0.00000
4.158341010
0.6
0.7
0.69767
0.61262
6.11067 108
2.15806107
0.2
0.96078
3.91674109
0.8
0.52729
5.06502107
0.3
0.91393
1.12525108
0.9
0.44485
9.70467107
0.4
0.85214
1.64987108
1.0
0.5
0.77880
2.52770109
tk
El valor aproximado de y (0) es 0.3678
yk
0.36788
-
ERROR
1.62525106
Ejemplo 4: Resolver el problema de valor inicial:
y ' =f ( t , y ) =4 e0.8 t 0.5 y , x [ 0,2 ] , y ( 0 ) =2
Tomando h=0.1. Encontrar la solucin aproximada para el punto 2.
Solucin:
La solucin aproximada para y (2) es 14.8439.
Ejemplo 5: aproxima la solucin para el siguiente problema.
'
3t
y =t e 2 y , t [ 0,1 ] , h=0.2
Y(1)=??
Solucin:
