Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
RUNGE KUTTA
(PARA ORDEN: SEGUNDO, TERCERO Y CUARTO)
Jhon Jaimes Figueroa, Brandon Martínez Costa, Angélica Popayán Páez, Sara Sánchez Carrillo,
Departamento de Química, Facultad de Ingenierías y Arquitectura, Universidad de Pamplona, Pamplona,
03 de Junio de 2017
____________________________________________________________________________________
SEGUNDO ORDEN:
La versión de segundo orden es:
Dónde:
los valores de a1, a2, p1 y q11 se evalúan al igualar la ecuación de segundo orden con la expansión de la
serie de Taylor hasta el término de segundo orden, Al hacerlo, desarrollamos tres ecuaciones para
evaluar las cuatro constantes desconocidas. Las tres ecuaciones son:
Entonces se resuelven de manera simultánea las ecuaciones, obteniendo:
Debido a que podemos elegir un número infinito de valores para a2, hay un número infinito de métodos
Runge Kutta de segundo orden. Sin embargo, se obtienen diferentes resultados cuando la solución es
más complicada.
MÉTODO DE HEUN CON UN SOLO CORRECTOR (a 2 = 1/2)
Si suponemos que a2 es 1/2, se puede obtener que a1 = 1/2 y p1 = q11 = 1. Estos parámetros, al
sustituirse en la ecuación original dan:
Dónde:
Formando líderes para la construcción de un
1
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
Ejemplo:
Resolver problema de valor inicial
dy
=( y +1 ) ( x +1 ) cos ( x +2 x ) y ( 0 ) =4
2
dx
En el intervalo 0 ≤ x ≤ 2, con h = 0.5.
Tenemos que x 0=0 , y 0=4 ,Para x 1=0. 5 la ordenada correspondiente será:
y 1=4+ ( 12 k + 12 k )∗0.5
1 2
Con:
k 1=f (x 0 , y 0 )=( 4 +1 ) ( 0+1 ) cos ( 02 +2 ( 0 ) ) =5
k 2=f ( x 0 +h , y 0 +k 1∗h )
k 2=( ( 4+ ( 5 )∗( 0.5 ) ) +1 )( ( 0+ ( 0.5 ) ) +1 ) cos ( ( 0+ ( 0.5 ) ) + 2 ( 0+ ( 0.5 ) ) )
2
k 2=3.547
Y así:
y 1=4+ ( 12 ( 5) + 12 ( 3.547) )∗0.5
y 1=6.137
De manera análoga se determinan los puntos:
(x 2 ; y 2 )=(1; 2.613);(x 3 ; y 3)=(1.5 ; 0.836);(x 4 ; y 4 )=(2 ;1.095)
EL MÉTODO DEL PUNTO MEDIO (a2 = 1)
Si suponemos que a2 es 1, entonces a1 = 0, p1 = q11 = 1/2, y la ecuación se convierte en:
Dónde:
Formando líderes para la construcción de un
2
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
Ejemplo:
Resolver problema de valor inicial
dy
=( y +1 ) ( x +1 ) cos ( x +2 x ) y ( 0 ) =4
2
dx
En el intervalo 0 ≤ x ≤ 2, con h = 0.5.
Tenemos que x 0=0 , y 0=4 ,Para x 1=0. 5 la ordenada correspondiente será:
y 1=4+ k 2∗0.5
Con:
k 1=f (x 0 , y 0 )=( 4 +1 ) ( 0+1 ) cos ( 02 +2 ( 0 ) ) =5
( 1
2
1
k 2=f x 0+ h , y 0 + ∗k 1∗h
2 )
(( ) )(( ) ) (( ) ( ))
2
1 1 1 1
k 2= 4 + ∗( 5 )∗( 0.5 ) +1 0+ ∗( 0.5 ) +1 cos 0+ ∗( 0.5 ) +2 0+ ∗( 0.5 )
2 2 2 2
k 2=6.609
Y así:
y 1=14 + ( 6.609 )∗0.5
y 1=7.304
De manera análoga se determinan los puntos:
(x 2 ; y 2 )=(1; 3.468);(x 3 ; y 3)=(1.5 ; 1.932);(x 4 ; y 4)=(2 ; 7.048)
MÉTODO DE RALSTON (a2 = 2/3)
Al seleccionar a2 = 2/3 se obtiene un mínimo en el error de truncamiento para los algoritmos Runge
Kutta de segundo orden. Con esta versión, a1 = 1/3 y p1 = q11 = 3/4 y da:
Dónde:
Formando líderes para la construcción de un
3
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
Ejemplo:
Resolver problema de valor inicial
dy
=( y +1 ) ( x +1 ) cos ( x +2 x ) y ( 0 ) =4
2
dx
En el intervalo 0 ≤ x ≤ 2, con h = 0.5.
Tenemos que x 0=0 , y 0=4 ,Para x 1=0. 5 la ordenada correspondiente será:
y 1=4+ ( 13 k + 23 k )∗0.5
1 2
Con:
k 1=f (x 0 , y 0 )=( 4 +1 ) ( 0+1 ) cos ( 02 +2 ( 0 ) ) =5
( 3
4
3
k 2=f x 0+ h , y 0 + ∗k 1∗h
4 )
(( ) )(( ) ) (( ) ( ))
2
3 3 3 3
k 2= 4 + ∗( 5 )∗( 0.5 ) +1 0+ ∗( 0.5 ) +1 cos 0+ ∗( 0.5 ) +2 0+ ∗( 0.5 )
4 4 4 4
k 2=5.945
Y así:
y 1=4+ ( 13 ( 5) + 23 ( 5.945) )∗0.5
y 1=6.815
De manera análoga se determinan los puntos:
(x 2 ; y 2 )=(1; 2.771);(x 3 ; y3 )=(1.5 ; 1.471);(x 4 ; y 4 )=(2 ; 3.944)
TERCER ORDEN:
Para n = 3, es posible efectuar un desarrollo similar al del método de segundo orden. El resultado de tal
desarrollo genera seis ecuaciones con ocho incógnitas. Por lo tanto, se deben dar a priori los valores de
dos de las incógnitas con la finalidad de establecer los parámetros restantes. Una versión común que se
obtiene es:
Formando líderes para la construcción de un
4
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
Dónde:
Ejemplo:
Resolver el problema de valor inicial
dx
=
(
dy 2 yx +1
y
2
y ( 0 )=1
)
En el intervalo 0 ≤ x ≤ 2, con h = 0.5.
Tenemos que x 0=0 , y 0=1 ,Para x 1=0. 5 la ordenada correspondiente será:
0.5
y 1=1+ ∗(k 1 + 4 k 2 +k 3)
6
Con:
k 1=f (x 0 , y 0 )=
( 2 ( 1 )( 0 )+ 1
1
2 )
=1
( 1
2
1
k 2=f x 0+ h , y 0 + ∗k 1∗h
2 )
( ( )( )
)
1 1
2 1+ ∗( 5 )∗( 0.5 ) 0+ ∗( 0.5 ) +1
2 2
k 2=
( )
2
1
1+ ∗( 5 )∗( 0.5 )
2
k 2=0.65
k 3=
( 2 ( 1−( 1 )∗( 0.5 ) +2∗( 0.65 )∗( 0.5 ) ) ( 0+ ( 0.5 ) ) +1
( 1−( 1 )∗( 0.5 ) +2∗( 0.65 )∗( 0.5 ) )
2
)
k 3=1.626
Y así:
0.5
y 1=1+ ∗(1+4 ( 0.65 ) +1.626)
6
y 1=1.435
Formando líderes para la construcción de un
5
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
De manera análoga se determinan los puntos:
(x 2 ; y 2 )=(1; 1.984 );( x3 ; y 3 )=(1.5 ; 2.663);(x 4 ; y 4 )=(2 ; 3.481)
CUARTO ORDEN:
El más popular de los métodos Runge Kutta es el de cuarto orden. Como en el caso de los
procedimientos de segundo orden, hay un número infinito de versiones. La siguiente, es la forma
comúnmente usada y, por lo tanto, le llamamos método clásico Runge Kutta de cuarto orden:
Dónde:
Ejemplo:
Resolver por un método de Runge-Kutta de cuarto orden el problema de valor inicial:
2
y '=x −3 y ; y ( 0 )=1
En el intervalo 0 ≤ x ≤ 0.4, con h = 0.1.
Tenemos que x 0=0 , y 0=1 , y f (x , y)=x 2−3 y .Para x 1=0.1 la ordenada correspondiente será:
0.1
y 1=1+ ∗(k 1 +2 k 2 +2 k 3 +k 4)
6
Con:
2
k 1=f (x 0 , y 0 )=0 −3 ( 1 ) =−3
( )( ) ( )
2
h h 0.1 0.1
k 2=f x 0+ , y 0+ ∗k 1 = 0+ −3 1+ (−3) =−2.5475
2 2 2 2
( )( ) ( )
2
h h 0.1 0.1
k 3=f x 0+ , y 0+ ∗k 2 = 0+ −3 1+ (−2.5475) =−2.61538
2 2 2 2
k 4=f ( x 0+ h , y 0+ h∗k 3 )=( 0+ 0.1 )2−3 ( 1+0.1 (−2.61538 ) ) =−2.20539
Y así:
Formando líderes para la construcción de un
6
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
0.1
y 1=1+ ∗( (−3 )+ 2 (−2.5475 )+ 2 (−2.61538 )+ (−2.20539 ) )
6
y 1=0. 741148
De manera análoga se determinan los puntos:
(x 2 ; y 2 )=(0 .2 ; 0 . 551151); (x 3 ; y 3 )=(0 .3 ; 0 . 413894 ); (x 4 ; y 4 )=(0 . 4 ; 0 . 317435)
La solución exacta es:
y=
25 −3 x 1 2 2
27
e + x − x+
3 3 (
2
9 )
De manera que:
y ( 0.1 )=0.741127 ; y ( 0.2 )=0.551121 ; y ( 0.3 ) =0.413860 ; y ( 0.4 )=0.317402;
Formando líderes para la construcción de un
7
nuevo país en paz
Universidad de Pamplona
Pamplona - Norte de Santander - Colombia
Tels: (7) 5685303 - 5685304 - 5685305 - Fax: 5682750 - [Link]
Formando líderes para la construcción de un
8
nuevo país en paz