Universidad Autnoma de

Quertaro
Facultad de Ingeniera
Divisin de Posgrado
Matemticas Avanzadas


Mtodos Numricos Para Solucin de Ecuaciones
Diferenciales


Catedrtico:
M.C. Patricia Spindola



Alumno:
Ricardo Alonso Garca Salas





Santiago de Quertaro, Quertaro a 09 de diciembre de 2010.


2010



Matemticas Avanzadas Mtodos Numricos



Pgina 2


Contenido
Introduccin ........................................................................................................................................ 3
Mtodos Numricos de Solucin de Ecuaciones Diferenciales de Primer Orden .............................. 3
Mtodo de Euler ............................................................................................................................. 3
Mtodo de Euler mejorado ............................................................................................................. 4
Mtodo de Runge-Kutta .................................................................................................................. 5
Solucin de una ecuacin diferencial de Primer orden ...................................................................... 7
Solucin Analtica ............................................................................................................................ 7
Solucin Numrica por el mtodo de Euler .................................................................................... 8
Solucin Numrica por el mtodo de Euler Mejorado ................................................................... 9
Solucin Numrica por el mtodo de Runge-Kutta ...................................................................... 10
Conclusiones ..................................................................................................................................... 10
Bibliografa ........................................................................................................................................ 11



Matemticas Avanzadas Mtodos Numricos



Pgina 3


Introduccin

Dada cierta dificultada para encontrar soluciones exactas a la resolucin de ecuaciones
diferenciales, podemos deducir aproximaciones usando mtodos numricos como herramienta.


En este trabajo se presentan tres mtodos para encontrar la solucin de una ecuacin diferencial
con un problema de valor inicial y(x
0
)=y
0


Los mtodos numricos al ser una solucin aproximando presentan entonces un error. Cada
mtodo presentar una aproximacin diferente dependiendo de que tan preciso sea el mtodo
para la ecuacin diferencial analizada.

Verificaremos entonces cual es el mtodo que ofrece un menor rango de error



Mtodos Numricos de Solucin de Ecuaciones Diferenciales de Primer Orden

Mtodo de Euler
El mtodo es una de las tcnicas ms simples para aproximar soluciones de una ecuacin
diferencial. Este mtodo se aplica para encontrar la solucin a Ecuaciones Diferenciales
Ordinarias, esto es, cuando la funcin involucra slo una variable independiente.

Se llama mtodo de Euler o Mtodo de las Tangentes y presenta un esquema iterativo que
aproxima la solucin de un problema de valor inicial de la forma.

( , )
dy
f t y
dt
= ( )
0 0
y t y = =

Se observa que la pendiente de la recta tangente a la curva ( ) y f x = est dada por '( ) f x y es
aproximadamente igual a la pendiente de la recta secante.


Matemticas Avanzadas Mtodos Numricos



Pgina 4



1 1 i i i i
n n
y y y y
x h x h
+ +
+
=
+


Siempre y cuando h sea pequeo de aqu obtenemos que:
1 1 n
x i i i i
i i
y y y y
f
x h x h
+ +

=
+


Con lo cual podemos usar el punto
0 0
( , ) x y para construir el siguiente punto
1 1
( , ) x y y as
sucesivamente. De esta forma generamos la sucesin de puntos:
0 0
( , ) x y ,
1 1
( , ) x y ,, ( , )
n n
x y

La solucin aproximada est dada por la iteracin:
1
( , )
i i i i
y y h f t y
+
= +

En esta frmula se predice un nuevo valor de y por medio de la pendiente que es igual a la
primera derivada en el valor original de x, este nuevo valor habr de extrapolarse en forma lineal
sobre el tamao de paso h.

El tamao entonces del error en el mtodo de Euler para un problema de valor inicial sobre un
tiempo fijo est en proporcin al tamao de paso, o, en otras palabras, es inversamente
proporcional al nmero de pasos que usamos sobre un intervalo fijo.


Mtodo de Euler mejorado
En el mtodo de Euler se tom como vlida para todo el intervalo la derivada encontrada en un
extremo de ste. Para obtener una exactitud razonable se utiliza un intervalo muy pequeo, a
cambio de un error de redondeo mayor (ya que se realizarn ms clculos).

En el mtodo de Euler modificado trata de evitar este problema utilizando un valor promedio de la
derivada tomada en los dos extremos del intervalo. En lugar de la derivada tomada en un solo
extremo.

El mtodo de Euler modificado conste de dos pasos bsicos:

Se parte de y se utiliza el mtodo de Euler a fin de calcular el valor de Y correspondiente a. Este
valor de y se denotar aqu como ya que solo es un valor transitorio para y. Esta parte del
proceso se conoce como paso Predictor.

El segundo paso se llama Corrector, pues trata de corregir la Prediccin. En el nuevo punto
obtenido sdkj, se evala la derivada y sando la ecuacin diferencial ordinaria del problema de
valor inicial que se est resolviendo; se obtiene la media aritmtica de esta derivada y la derivada
en el punto inicial.
Matemticas Avanzadas Mtodos Numricos



Pgina 5


Se usa la derivada promedio para calcular un nuevo valor de y1, con la ecuacin
y1=y0+hf(x0,y0) que deber ser ms exacto que y1
| |
1 0
1 0 0 1 1
( , ) ( , )
2
o
x x
y y f x y f x y

= + +

Y se tomar como valor definitivo de y1. Este procedimiento se repite hasta llegar a yn
El esquema iterativo para este mtodo quedar en general as:


1. Usando el paso Predictor resulta:
1
( , )
i i i i
y y h f t y
+
= +

2. Una vez obtenida
1 i
y
+
se calcula
1 1
( , )
i i
f t y
+ +
y se promedia con la derivada previa
( , )
i i
f x y para encontrar la derivada promedio:
1 1
( , ) ( , )
2
i i i i
f t y f t y
+ +
+

3. Se sustituye ( , )
i i
f x y con este valor promedio en la ecuacin de iteracin de Euler y se
obtiene
1 1
1
( , ) ( , )
2
i i i i
i i
f t y f t y
y y
+ +
+
+
= +


Mtodo de Runge-Kutta
En el mtodo de Euler mejorado, usamos un promedio de dos pendientes para determinar cada
valor. En otras palabras, para calcular
1 i
y
+
a partir de
i
y , empleamos el promedio de valores del
lado derecho. Por analoga con la integracin numrica, el mtodo de Euler mejorado es similar
a la regla trapezoidal.

Para la integracin numrica existen algoritmos que por lo general son ms eficientes que la
regla trapezoidal. La regla de Simpson aproxima el rea bajo la grfica usando la interpolacin
parablica y permite una mejor estimacin para la integral. De hecho, la regla de Simpson puede
interpretarse como un promedio ponderado de valores donde se da una importancia doble a los
valores de la funcin en los puntos medios de los sub intervalos, en comparacin a la que se
otorga a los puntos extremos. El mtodo de Runge-Kutta es similar a la regla de Simpson en
tanto que considera un promedio ponderado.

Para calcular el valor de
1 i
y
+
a partir de
i
y , empleamos cuatro pendientes dadas por la funcin
( , ) f t y que define la ecuacin diferencial. Esas pendientes se denominan
k
m ,
k
n ,
k
q y
k
p .
Tambin necesitamos diversas variables intermedias. El control de la nomenclatura resulta en
este caso ms delicado que en el mtodo de Euler mejorado, por lo que esas variables adquieren
muchos tipos diferentes de signos de acentuacin.


Matemticas Avanzadas Mtodos Numricos



Pgina 6


Las cuatro pendientes se determinan sucesivamente como sigue:
La primera pendiente
k
m se calcula igual que en el mtodo de Euler; es decir, ( , )
k k k
m f t y =

En el mtodo de Euler mejorado usamos
k
m para producir una segunda pendiente
correspondiente a un punto con
1 k k
t t t t
+
= = + A . El mtodo de Runge- Kutta hace casi lo
mismo excepto que va slo a la mitad del camino a lo largo del eje t hasta
2
t
k
t t
A
= +

. Esto
significa que utilizamos
k
m para producir un punto ,
k
t y

donde

2
t
k k k
y y m
A
= +

Una vez que hemos determinado este punto, usamos la funcin ( , ) f t y para determinar la
segunda pendiente
k
n por medio de
( , )
k k k
n f t y =



Ahora repetimos el paso previo del algoritmo donde empleamos la pendiente
k
n en lugar de
k
m .
En otras palabras, pasamos de ( , )
k k
t y a la lnea t t =

a lo largo de una lnea de pendiente
k
n .
Obtenemos entonces un nuevo nmero
k
y , donde

2
k k k
t
y y n
A
= +

Dado este punto sobre la lnea t t =

, calculamos la tercera pendiente
k
q mediante
( , )
k k
q f t y
=



Finalmente, obtenemos nuestra cuarta pendiente usando
k
q para producir un punto sobre la lnea
1 k
t t
+
= . Obtenemos
k k k
y y q t
+
= A

Una vez que tememos este cuarto punto, calculamos ah la pendiente con
1
( , )
k k k
p f t y
+
=

Ahora que tenemos las cuatro pendientes, tomamos un promedio ponderado y con ste
calculamos el siguiente paso. Ponderamos al doble cada una de las pendientes que provienen de
los puntos con t t =

, respecto a las otras dos pendientes. En otras palabras, nuestro promedio
ponderado es

2 2
6
k k k k
m n q p + + +

Por lo tanto, el paso que realmente tomamos es
Matemticas Avanzadas Mtodos Numricos



Pgina 7


1
2 2
6
k k k k
k k
m n q p
y y t
+
+ + + | |
= + A
|
\ .

Solucin de una ecuacin diferencial de Primer orden
Solucin Analtica

Sea la ecuacin diferencial ordinaria
1
2
2
4 ( 1) ydy x y dx = +

Rescribiendo la ecuacin
1
2
2
( 1) 4 y ydy xdx

+ =

Integrando
( )
1
2
2
1
1 2 4
2
y ydy xdx

+ =
} }


Entonces resulta
( )
1
2
2 2
1
2 1 2
2
y x c + = +

Simplificando
( )
1
2
2 2
1 2 y x c + = +

Despejando y
2
( )
2
2 2
1 2 y x c + = +

( )
2
2 2
2 1 y x c = +

Sustituyendo las condiciones inciales y(0)=1

2
1 1 2 c c = =

( )
2
2 2
2 2 1 y x =

Despejando y
( )
2
2
2 2 1 y x =
Matemticas Avanzadas Mtodos Numricos



Pgina 8


Solucin Numrica por el mtodo de Euler
De la metodologa anterior usamos en Mtodo de Euler para resolver la ecuacin diferencial con
las condiciones de valor inicial. Los resultados se muestran de forma clara en la siguiente tabla:

METODO DE EULER
h= 0.1


x0= 0


y0= 1
i Xi Yi f(Xi,Yi) Y real E. Absoluto E. relativo
0 0.000 1.000 0.000 1.000 0.000 0.000
1 0.100 1.000 0.566 1.028 0.028 2.732
2 0.200 1.057 1.102 1.110 0.054 4.836
3 0.300 1.167 1.580 1.242 0.075 6.029
4 0.400 1.325 2.005 1.417 0.092 6.500
5 0.500 1.525 2.392 1.632 0.107 6.556
6 0.600 1.764 2.759 1.885 0.121 6.420
7 0.700 2.040 3.118 2.175 0.135 6.211
8 0.800 2.352 3.477 2.502 0.150 5.983
9 0.900 2.700 3.839 2.865 0.165 5.756
10 1.000 3.084 4.205 3.264 0.181 5.538
11 1.100 3.504 4.576 3.702 0.197 5.330
12 1.200 3.962 4.951 4.176 0.214 5.134
13 1.300 4.457 5.329 4.689 0.232 4.947
14 1.400 4.990 5.711 5.240 0.250 4.769
15 1.500 5.561 6.096 5.829 0.268 4.601
Se us un valor de paso h=0.1 la funcin diferencial y su solucin se muestran tambin en la
tabla.

1
2
2
4 ( 1) dy x y
dx y
+
=
( )
2
2
2 2 1 y x =
Matemticas Avanzadas Mtodos Numricos



Pgina 9


Solucin Numrica por el mtodo de Euler Mejorado
De la metodologa anterior usamos en Mtodo de Euler Mejorado para resolver la ecuacin
diferencial con las condiciones de valor inicial. Los resultados se muestran de forma clara en la
siguiente tabla:

MTODO DE EULER MEJORADO
x0= 0




y0= 1


h= 0.1
i xi yi f(x0,y0) yn+1 f(xn+1,yn+1) V.Real E. Absoluto E. relativo
0 0.00 1.00 0.00 1.00 0.000 0.000
1 0.10 1.03 0.56 1.00 0.57 1.03 0.000 0.019
2 0.20 1.11 1.08 1.03 1.12 1.11 0.002 0.155
3 0.30 1.25 1.54 1.11 1.61 1.24 0.005 0.394
4 0.40 1.43 1.95 1.25 2.05 1.42 0.009 0.642
5 0.50 1.65 2.34 1.43 2.44 1.63 0.014 0.831
6 0.60 1.90 2.71 1.65 2.81 1.89 0.018 0.944
7 0.70 2.20 3.08 1.90 3.16 2.18 0.022 0.991
8 0.80 2.53 3.44 2.20 3.52 2.50 0.025 0.991
9 0.90 2.89 3.81 2.53 3.87 2.86 0.028 0.961
10 1.00 3.29 4.18 2.89 4.23 3.26 0.030 0.913
11 1.10 3.73 4.56 3.29 4.60 3.70 0.032 0.856
12 1.20 4.21 4.93 3.73 4.97 4.18 0.033 0.797
13 1.30 4.72 5.32 4.21 5.34 4.69 0.035 0.738
14 1.40 5.28 5.70 4.72 5.72 5.24 0.036 0.681
15 1.50 5.87 6.09 5.28 6.11 5.83 0.037 0.628
Se us un valor de paso h=0.1 la funcin diferencial y su solucin se muestran tambin en la
tabla.

( )
2
2
2 2 1 y x =
1
2
2
4 ( 1) dy x y
dx y
+
=
Matemticas Avanzadas Mtodos Numricos


Pgina
10


Solucin Numrica por el mtodo de Runge-Kutta
De la metodologa anterior usamos en Mtodo de Runge-Kutta para resolver la ecuacin
diferencial con las condiciones de valor inicial. Los resultados se muestran de forma clara en la
siguiente tabla:

MTODO DE RUNGE-KUTTA
x0= 0

y0= 1

h= 0.1

i xi k1 k2 k3 k4 yi V.Real E. Absoluto E. relativo
0 0.000 1.000 1.000 0.000 0.000
1 0.100 0.000 0.028 0.028 0.028 1.023 1.028 0.005 0.449
2 0.200 0.056 0.083 0.082 0.082 1.102 1.110 0.009 0.784
3 0.300 0.108 0.132 0.132 0.132 1.229 1.242 0.012 0.975
4 0.400 0.155 0.176 0.176 0.176 1.402 1.417 0.015 1.056
5 0.500 0.197 0.216 0.216 0.216 1.615 1.632 0.018 1.074
6 0.600 0.235 0.254 0.254 0.254 1.865 1.885 0.020 1.061
7 0.700 0.272 0.291 0.290 0.290 2.153 2.175 0.022 1.034
8 0.800 0.309 0.327 0.327 0.327 2.477 2.502 0.025 1.001
9 0.900 0.345 0.363 0.363 0.363 2.837 2.865 0.028 0.967
10 1.000 0.382 0.400 0.400 0.400 3.234 3.264 0.030 0.932
11 1.100 0.419 0.437 0.437 0.437 3.668 3.702 0.033 0.899
12 1.200 0.456 0.475 0.475 0.475 4.140 4.176 0.036 0.866
13 1.300 0.494 0.513 0.513 0.513 4.650 4.689 0.039 0.835
14 1.400 0.532 0.551 0.551 0.551 5.197 5.240 0.042 0.805
15 1.500 0.570 0.590 0.590 0.590 5.784 5.829 0.045 0.777
Se us un valor de paso h=0.1 la funcin diferencial y su solucin se muestran tambin en la
tabla.





Conclusiones

La siguiente tabla muestra un resumen de los resultados obtenidos por los tres mtodos y
podemos observar cmo se comportan las soluciones numricas con respecto a la solucin
analtica en el grfico.

Matemticas Avanzadas Mtodos Numricos


Pgina
11


Podemos concluir en base a los resultados obtenidos que el Mtodo de Euler mejorado presenta
un error relativo menor comparado con los otros dos mtodos estudiados y el mtodo de Euler
es el que en menor medida se asemeja ms a la funcin solucin.

i xi VALOR REAL EULER
EULER
MEJORADO
RUNGE-
KUTTA
0 0.000 1.000 1.000 1.00 1.000
1 0.100 1.028 1.000 1.03 1.023
2 0.200 1.110 1.057 1.11 1.102
3 0.300 1.242 1.167 1.25 1.229
4 0.400 1.417 1.325 1.43 1.402
5 0.500 1.632 1.525 1.65 1.615
6 0.600 1.885 1.764 1.90 1.865
7 0.700 2.175 2.040 2.20 2.153
8 0.800 2.502 2.352 2.53 2.477
9 0.900 2.865 2.700 2.89 2.837
10 1.000 3.264 3.084 3.29 3.234
11 1.100 3.702 3.504 3.73 3.668
12 1.200 4.176 3.962 4.21 4.140
13 1.300 4.689 4.457 4.72 4.650
14 1.400 5.240 4.990 5.28 5.197
15 1.500 5.829 5.561 5.87 5.784



Bibliografa

Dennis G. Zill.Ecuaciones Diferenciales con Aplicaciones de Modelado Sexta Edicin.

Paul Blanchard, Robert Devaney y Glen Hall.Ecuaciones Diferenciales.1999
0
1
2
3
4
5
6
7
0 0.5 1 1.5 2
E
j
e

Y
Eje X
Valores de Y
VALOR REAL
EULER
EULER MEJORADO
RUNGE-KUTTA