MÉTODO DE EULER
En matemática y computación, el método de Euler, llamado así en honor a Leonhard
Euler, es un procedimiento de integración numérico para resolver ecuaciones
diferenciales ordinarias(EDO) a partir de un valor inicial dado. El método de Euler es
el más simple de los métodos numéricos para resolver un problema de valor inicial.
El Método de Euler o de las Tangentes constituye el primer y más sencillo ejemplo
de método numérico para la resolución de un problema de valor inicial:
y´ = f(x, y), y(x0) = y0
donde suponemos además que se verifican las hipótesis del Teorema de Picard 1 ,
y en consecuencia existe solución única para el problema.
Interpretando la e.d.o. y´ = f(x, y) como un campo de direcciones en el plano xy y la
condición inicial y(x0) = y0 como un punto (x0, y0) de dicho plano, podemos
aproximar la función solución y(x) por medio de la recta tangente a la misma que
pasa por ese punto:
y(x) ∼= y´ + f(x0, y0)(x − x0) 1
Consideraremos en general que la función f(x, y) es diferenciable en un entorno del
punto (x0, y0). Si bien es cierto que se trata de una condición más restrictiva de lo
estrictamente necesario, en la práctica trabajaremos siempre con funciones de ese
tipo,donde se ha utilizado que la pendiente de dicha tangente es:
m = y´(x0)
y, en consecuencia: m = f(x0, y0). Calculamos así de manera aproximada el valor
de la solución y en el punto de abscisa x1 como:
y(x1) ∼= y1 = y0 + f(x0, y0)(x1 − x0)
y con este punto (aproximado) ya calculado, podemos repetir el método para
obtener otro punto aproximado (x2, y2) de la forma:
y(x2) ∼= y2 = y1 + f(x1, y1)(x2 − x1)
y así sucesivamente. Es habitual en este método tomar abscisas equiespaciadas,
es decir, calcular la solución aproximada en puntos de la forma:
xn = xn−1 + h = x0 + nh,
siendo h el paso del método. De esta forma se obtienen las fórmulas que nos
determinan la soluci´on aproximada en la forma:
xn = xn−1 + h; yn = yn−1 + f(xn−1, yn−1) h
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
a 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.
Ejemplo:
Resolveremos el problema de valor inicial ( y 0 = x √y y(1) = 4 por el m´etodo de
Euler con h = 0.1 para los puntos x = 1.1, 1.2, 1.3, 1.4 y 1.5.
En este problema tenemos h = 0.1, (x0, y0) = (1, 4) y la funci´on f(x, y) es: f(x, y) = x
√y.
Por tanto: yn = yn−1 + xn−1 √ yn−1 h
Dado que el problema se puede resolver también de forma exacta, presentamos en
la tabla y gráfica siguientes los resultados:
i Xi Yi
Sol. Exacta
0 1 4 4
1 1.1 4.2 4.21276
2 1.2 4.42543 4.45210
3 1.3 4.67787 4.71976
4 1.4 4.95904 5.01760
5 1.5 5.27081 5.34766
APLICACIÓN DEL MÉTODO DE EULER PASO A PASO:
Consiste en multiplicar los intervalos que va de a en subintervalos de
ancho ; osea:
de manera que se obtiene un conjunto discreto de puntos:
del intervalo de interés . Para cualquiera de estos
puntos se cumlple que:
La condición inicial , representa el punto por donde
pasa la curva solución de la ecuación de el planteamiento inicial, la cual se denotará
como .
Ya teniendo el punto se puede evaluar la primera derivada de en ese
punto; por lo tanto:
Con esta información se traza una recta, aquella que pasa por y de pendiente
. Esta recta aproxima en una vecinidad de . Tómese la recta
como reemplazo de y localícese en ella (la recta) el valor de y
correspondiente a . Entonces, podemos deducir según la Gráfica A:
Se resuelve para :
Es evidente que la ordenada calculada de esta manera no es igual a ,
pues existe un pequeño error. Sin embargo, el valor sirve para que se aproxime
en el punto y repetir el procedimiento anterior a fin de generar
la sucesión de aproximaciones siguiente:
Método de Euler Mejorado
Este método se basa en la misma idea del método anterior, pero hace un
refinamiento en la aproximación, tomando un promedio entre ciertas pendientes.
La fórmula es la siguiente:
Donde:
Para entender esta fórmula, analicemos el primer paso de la aproximación, con base
en la siguiente gráfica:
En la gráfica, vemos que la pendiente promedio corresponde a la pendiente de la
recta bisectriz de la recta tangente a la curva en el punto de la condición inicial y la
"recta tangente" a la curva en el punto donde es la aproximación obtenida
con la primera fórmula de Euler. Finalmente, esta recta bisectriz se traslada
paralelamente hasta el punto de la condición inicial, y se considera el valor de esta
recta en el punto como la aproximación de Euler mejorada.
