SOLUCION DE ECUACIONES DIFERENCIALES

“METODO DE EULER MEJORADO”

 

 

En el método de Euler se tomó como válida para todo el intervalo la derivada encontrada en un extremo de éste Fig. . Para obtener una exactitud ra­zonable se utiliza un intervalo muy pequeño, a cambio de un error de redondeo mayor (ya que se realizarán más cálculos).

El método 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 METODO DE EULER MODIFICADO CONSTA DE DOS PASOS BASICOS:

 

1. Se parte de (xo,Yo) Y se utiliza el método de Euler a fin de calcular el valor de Y correspondiente a Xl' Este valor de Y se denotará aquí como YI' ya que sólo es un valor transitorio para Yl' Esta parte del proceso se conoce como paso predictor.

2. El segundo paso se llama corrector, pues trata de corregir la predicción. En el nuevo punto obtenido (XI, Yl) se evalúa la derivada [(xI' YI) usando la ecua­ción diferencial ordinaria del PVI que se esté resolviendo; se obtiene la media aritmética de esta derivada y la derivada en el punto inicial (xo' Yo)

 

 

1/2 [F(xo ,Yo) + F(Xl,YI)] = derivada promedio

 

Se usa la derivada promedio para calcular un nuevo valor de y1, con la ecuación y1=y0+hf(x0,y0), que deberá ser mas exacto que y1

 

 

y se tomara como valor definitivo de y1. Este procedimiento se repite hasta llegar a yn.

El esquema iterativo para este método  quedara en general así:

Primero, usando el paso de predicción resulta:

 

 .

Una vez obtenida yi+1 se calcula f(xi+1,yi+1), la derivada en el punto (xi+1,yi+1), y se promedia con la derivada previa (xi,yi) para encontrar la derivada promedio

 

Se sustituye f(xi,yi) con este valor promedio en la ecuación de iteración de euler y se obtiene:

 

 

 

 

EJEMPLO:

 

 RESUELVA EL SIGUIENTE:

 

           

           PVI       POR EL METODO DE EULER MEJORADO.

 

SOLUCION:

 

Al utilizar 5 intervalos se tiene:

PRIMERA ITERACION:

 

Primer paso:   y1=y0+hf(x0,y0)= 2+0.2 (0-2)=1.6

 

Segundo paso:

 

                           Y (0.2)= y1=2+0.2 (-1.7)= 1.66

 

SEGUNDA ITERACION:

Primer paso: y2=y1+hf(x1,y1)=1.66+0.2 (0.2-1.66)=1.368

Segundo paso:

 

                            Y (0.4)=y1=1.66+0.2 (-1.214)=1.4172

 

Al continuar con los cálculos se obtiene:

 

                                   Y5=1.08509

                                   Y5’= 1.11222

ALGORITMO PARA EL METODO DE EULER MEJORADO:

 

Para obtener la aproximación YF a la solución de un PVI, proporcionar la función F(X,Y) y los

 

            DATOS:         La condición inicial X0,Y0, el valor XF donde se desea conocer el valor de YF y el numero N de subintervalos por emplear.

 

            RESULTADOS: aproximación a YF:Y0

 

PASO 1.        Hacer H= (HF-X0)/N

PASO 2.        Hacer I=1

PASO 3.        Mientras I<N, repetir los pasos 4 a 7.

                        PASO 4.        Hacer y1= y0+H*(X0,Y0)

                        PASO 5.        Hacer y0=y0+ H/2*(F(X0,Y0)+F(X0+H,Y1))

                        PASO 6.        Hacer x0=X0+H

                        PASO 7.        Hacer I=I+1

PASO 8.        IMPRIMIR Y0 Y FINALIZAR.