Método de Euler y Runge Kutta
Método de Euler y Runge Kutta
Método de Euler y Runge Kutta
ANÁLISIS NUMÉRICO
CURSO: P2
FECHA: 2020-02-04
Método de Euler
Para utilizar el Método de Euler con el fin de generar una solución numérica a un problema de
valor inicial de la forma:
y ' =f ( x , y )
y ( x 0 )= y 0
para así encontrar las coordenadas de los puntos de nuestra solución numérica. Terminamos
este proceso cuando hemos llegado al final del intervalo deseado.
Desde el punto de vista geométrico, tenemos en definitiva que el Método de Euler aproxima a
la función solución por medio de una línea poligonal, la aproximación será tanto peor cuanto
mayor sea en número de pasos, es decir, cuanto más “lejos” nos encontremos del punto inicial
(x0, y0). Por otro lado, el error será evidentemente tanto mayor cuanto más grande sea el
“paso” del método, h.
dx
=f ( x , y ) , y ( 0 )= y 0
dy
La fórmula básicamente calcula el siguiente valor yn+1 usando el yn actual más el promedio
ponderado de cuatro incrementos.
k1 es el incremento basado en la pendiente al principio del intervalo, utilizando y
k2 es el incremento basado en la pendiente en el punto medio del intervalo, usando y +
hk1/2.
k3 es de nuevo el incremento basado en la pendiente en el punto medio, usando y +
hk2/2.
k4 es el incremento basado en la pendiente al final del intervalo, usando y + hk3.
El método es un método de cuarto orden, lo que significa que el error de truncamiento local es
del orden de O(h5), mientras que el error acumulado total es del orden de O(h4).
ECUACIÓN:
dy
=0,1 √ y +0,4 x2
dx
con y ( 2 )=4
GRÁFICA:
Bibliografía
Barker, C. (2017). Numerical Methods for Solving Differential Equations EULER'S METHOD. San
Joaquin DELTA COLLEGE.
Butcher, J. (2005). Runge–Kutta methods for Differential Equations. The University of Auckland.
Ejercicio en clase
Método de Runge- Kutta de cuarto orden
Resolver la siguiente ecuación diferencial:
f ( x , y )=0.1 √ y + 0.4 x2
Con las siguientes condiciones: y(2)= 4 y h= 0.1. Encontrar el valor para y(2.5).
clear all
clc
fprintf('\n \t ECUACÓN DIFERENCIAL POR EL MÉTODO DE RUNGE-KUTTA DE CUARTO
ORDEN\n')
f=input('\n Ingrese la ecuacion diferencial\n','s');
x0=2
x1=2.5
y0=4
n=6;
h=0.1;
xs=x0:h:x1;
disp('-------------------------');
fprintf('\n Nº x0 y(x1)\n');
disp('-------------------------');
for i=1:n
it=i-1;
x0=xs(i);
x=x0;
y=y0;
k1=h*eval(f);
x=x0+h/2;
y=y0+k1/2;
k2=h*eval(f);
x=x0+h/2;
y=y0+k2/2;
k3=h*eval(f);
x=x0+h;
y=y0+k3;
k4=h*eval(f);
y0=y0+(k1+2*k2+2*k3+k4)/6;
fprintf('\n%2.0f%10.6f%10.6f\n',it,x0,y0);
end
disp('-------------------------');
fprintf('\n El punto aproximado y(x1) es = %8.6f\n',y0);
Resultados
El valor aproximado para y(2.5)= 5.40600
Conclusiones
El método de Runge-Kutta de cuarto orden, en análisis numérico es uno de los
métodos más usados para la resolución de ecuaciones diferenciales, mediante
la iteración se aproxima al valor del problema, haciendo uso del promedio de
las cuatro pendientes que posee la ecuación diferencial.