método de euler

Enfoque para encontrar soluciones numéricas de ecuaciones diferenciales ordinarias.

(Figura 1) Ilustración del método de Euler. La curva desconocida está en azul y su aproximación poligonal está en rojo.

En matemáticas y ciencia computacional , el método de Euler (también llamado método de Euler directo ) es un procedimiento numérico de primer orden para resolver ecuaciones diferenciales ordinarias (EDO) con un valor inicial dado . Es el método explícito más básico para la integración numérica de ecuaciones diferenciales ordinarias y es el método de Runge-Kutta más simple . El método de Euler lleva el nombre de Leonhard Euler , quien lo propuso por primera vez en su libro Institutionum calculi integralis (publicado entre 1768 y 1770). ^[1]

El método de Euler es un método de primer orden, lo que significa que el error local (error por paso) es proporcional al cuadrado del tamaño del paso y el error global (error en un momento dado) es proporcional al tamaño del paso. El método de Euler suele servir como base para construir métodos más complejos, por ejemplo, el método predictor-corrector .

Descripción geométrica

Propósito y por qué funciona.

Considere el problema de calcular la forma de una curva desconocida que comienza en un punto dado y satisface una ecuación diferencial dada. Aquí, una ecuación diferencial puede considerarse como una fórmula mediante la cual se puede calcular la pendiente de la recta tangente a la curva en cualquier punto de la curva, una vez calculada la posición de ese punto.

La idea es que si bien la curva se desconoce inicialmente, se conoce su punto de partida, que denotamos por (ver Figura 1). Luego, a partir de la ecuación diferencial, se puede calcular la pendiente de la curva en y, por tanto, la recta tangente. ${\displaystyle A_{0},}$ ${\displaystyle A_{0}}$

Da un pequeño paso a lo largo de esa línea tangente hasta un punto . A lo largo de este pequeño paso, la pendiente no cambia demasiado, por lo que estará cerca de la curva. Si pretendemos que todavía está en la curva, se puede utilizar el mismo razonamiento que para el punto anterior. Después de varios pasos, se calcula una curva poligonal ( ). En general, esta curva no diverge demasiado de la curva desconocida original, y el error entre las dos curvas puede hacerse pequeño si el tamaño del paso es lo suficientemente pequeño y el intervalo de cálculo es finito. ^[2] ${\displaystyle A_{1}.}$ ${\displaystyle A_{1}}$ ${\displaystyle A_{1}}$ ${\displaystyle A_{0}}$ ${\displaystyle A_{0}A_{1}A_{2}A_{3}\dots }$

Proceso de primer orden

Cuando se dan los valores de y , y la derivada de es una función dada de y se denota como . Comience el proceso configurando . A continuación, elija un valor para el tamaño de cada paso a lo largo del eje t y configúrelo (o equivalentemente ). Ahora, el método de Euler se utiliza para encontrar desde y : ^[3] ${\displaystyle t_{0}}$ ${\displaystyle y(t_{0})}$ ${\displaystyle y}$ ${\displaystyle t}$ ${\displaystyle y}$ ${\displaystyle y'(t)=f{\bigl (}t,y(t){\bigr )}}$ ${\displaystyle y_{0}=y(t_{0})}$ ${\displaystyle h}$ ${\displaystyle t_{n}=t_{0}+nh}$ ${\displaystyle t_{n+1}=t_{n}+h}$ ${\displaystyle y_{n+1}}$ ${\displaystyle y_{n}}$ ${\displaystyle t_{n}}$

${\displaystyle y_{n+1}=y_{n}+hf(t_{n},y_{n}).}$

El valor de es una aproximación de la solución en el tiempo , es decir, . El método de Euler es explícito , es decir, la solución es una función explícita de for . ${\displaystyle y_{n}}$ ${\displaystyle t_{n}}$ ${\displaystyle y_{n}\approx y(t_{n})}$ ${\displaystyle y_{n+1}}$ ${\displaystyle y_{i}}$ ${\displaystyle i\leq n}$

Proceso de orden superior

Si bien el método de Euler integra una EDO de primer orden, cualquier EDO de orden puede representarse como un sistema de EDO de primer orden. Cuando se le da la EDO de orden definida como ${\displaystyle N}$ ${\displaystyle N}$

${\displaystyle y^{(N+1)}(t)=f\left(t,y(t),y'(t),\ldots ,y^{(N)}(t)\right),}$

así como , , y , implementamos la siguiente fórmula hasta llegar a la aproximación de la solución a la EDO en el tiempo deseado: ${\displaystyle h}$ ${\displaystyle t_{0}}$ ${\displaystyle y_{0},y'_{0},\dots ,y_{0}^{(N)}}$

${\displaystyle {\vec {y}}_{i+1}={\begin{pmatrix}y_{i+1}\\y'_{i+1}\\\vdots \\y_{i+1}^{(N-1)}\\y_{i+1}^{(N)}\end{pmatrix}}={\begin{pmatrix}y_{i}+h\cdot y'_{i}\\y'_{i}+h\cdot y''_{i}\\\vdots \\y_{i}^{(N-1)}+h\cdot y_{i}^{(N)}\\y_{i}^{(N)}+h\cdot f\left(t_{i},y_{i},y'_{i},\ldots ,y_{i}^{(N)}\right)\end{pmatrix}}}$

Estos sistemas de primer orden pueden manejarse mediante el método de Euler o, de hecho, mediante cualquier otro esquema para sistemas de primer orden. ^[4]

Ejemplo de primer orden

Dado el problema del valor inicial

${\displaystyle y'=y,\quad y(0)=1,}$

Nos gustaría utilizar el método de Euler para aproximarnos . ^[5] ${\displaystyle y(4)}$

Usando un tamaño de paso igual a 1 ( h = 1 )

(Figura 2) Ilustración de la integración numérica de la ecuación. Azul es el método de Euler; verde, el método del punto medio ; rojo, la solución exacta, el tamaño del paso es ${\displaystyle y'=y,y(0)=1.}$ ${\displaystyle y=e^{t}.}$ ${\displaystyle h=1.0.}$

El método de Euler es

${\displaystyle y_{n+1}=y_{n}+hf(t_{n},y_{n}).}$

entonces primero debemos calcular . En esta ecuación diferencial simple, la función está definida por . Tenemos ${\displaystyle f(t_{0},y_{0})}$ ${\displaystyle f}$ ${\displaystyle f(t,y)=y}$

${\displaystyle f(t_{0},y_{0})=f(0,1)=1.}$

Al realizar el paso anterior, hemos encontrado la pendiente de la línea que es tangente a la curva solución en el punto . Recuerde que la pendiente se define como el cambio de dividido por el cambio de , o . ${\displaystyle (0,1)}$ ${\displaystyle y}$ ${\displaystyle t}$ ${\textstyle {\frac {\Delta y}{\Delta t}}}$

El siguiente paso es multiplicar el valor anterior por el tamaño del paso , que aquí tomamos igual a uno: ${\displaystyle h}$

${\displaystyle h\cdot f(y_{0})=1\cdot 1=1.}$

Dado que el tamaño del paso es el cambio en , cuando multiplicamos el tamaño del paso y la pendiente de la tangente, obtenemos un cambio en el valor. Luego, este valor se suma al valor inicial para obtener el siguiente valor que se utilizará para los cálculos. ${\displaystyle t}$ ${\displaystyle y}$ ${\displaystyle y}$

${\displaystyle y_{0}+hf(y_{0})=y_{1}=1+1\cdot 1=2.}$

Los pasos anteriores deben repetirse para encontrar y . ${\displaystyle y_{2}}$ ${\displaystyle y_{3}}$ ${\displaystyle y_{4}}$

${\displaystyle {\begin{aligned}y_{2}&=y_{1}+hf(y_{1})=2+1\cdot 2=4,\\y_{3}&=y_{2}+hf(y_{2})=4+1\cdot 4=8,\\y_{4}&=y_{3}+hf(y_{3})=8+1\cdot 8=16.\end{aligned}}}$

Debido a la naturaleza repetitiva de este algoritmo, puede resultar útil organizar los cálculos en forma de gráfico, como se ve a continuación, para evitar cometer errores.

La conclusión de este cálculo es que . La solución exacta de la ecuación diferencial es , entonces . Aunque la aproximación del método de Euler no fue muy precisa en este caso específico, particularmente debido a un tamaño de paso de valor grande , su comportamiento es cualitativamente correcto como muestra la figura. ${\displaystyle y_{4}=16}$ ${\displaystyle y(t)=e^{t}}$ ${\displaystyle y(4)=e^{4}\approx 54.598}$ ${\displaystyle h}$

Usar otros tamaños de paso

(Figura 3) La misma ilustración para ${\displaystyle h=0.25.}$

Como se sugirió en la introducción, el método de Euler es más preciso si el tamaño del paso es menor. La siguiente tabla muestra el resultado con diferentes tamaños de paso. La fila superior corresponde al ejemplo de la sección anterior y la segunda fila se ilustra en la figura. ${\displaystyle h}$

El error registrado en la última columna de la tabla es la diferencia entre la solución exacta en y la aproximación de Euler. En la parte inferior de la tabla, el tamaño del paso es la mitad del tamaño del paso en la fila anterior y el error también es aproximadamente la mitad del error en la fila anterior. Esto sugiere que el error es aproximadamente proporcional al tamaño del paso, al menos para valores bastante pequeños del tamaño del paso. Esto es cierto en general, también para otras ecuaciones; consulte la sección Error de truncamiento global para obtener más detalles. ${\displaystyle t=4}$

Otros métodos, como el método del punto medio también ilustrado en las figuras, se comportan más favorablemente: el error global del método del punto medio es aproximadamente proporcional al cuadrado del tamaño del paso. Por esta razón, se dice que el método de Euler es un método de primer orden, mientras que el método del punto medio es de segundo orden.

Podemos extrapolar de la tabla anterior que el tamaño de paso necesario para obtener una respuesta correcta con tres decimales es aproximadamente 0,00001, lo que significa que necesitamos 400.000 pasos. Este gran número de pasos conlleva un alto coste computacional. Por esta razón, se emplean métodos de orden superior, como los métodos de Runge-Kutta o los métodos lineales de varios pasos , especialmente si se desea una alta precisión. ^[6]

Ejemplo de orden superior

Para este ejemplo de tercer orden, supongamos que se proporciona la siguiente información:

${\displaystyle {\begin{aligned}&y'''+4ty''-t^{2}y'-(\cos {t})y=\sin {t}\\&t_{0}=0\\&y_{0}=y(t_{0})=2\\&y'_{0}=y'(t_{0})=-1\\&y''_{0}=y''(t_{0})=3\\&h=0.5\end{aligned}}}$

De esto podemos aislar y'''' para obtener la ecuación:

${\displaystyle f\left(t,y,y',y''\right)=y'''=\sin {t}+(\cos {t})y+t^{2}y'-4ty''}$

Usando eso podemos obtener la solución para : ${\displaystyle {\vec {y}}_{1}}$

$${\displaystyle {\vec {y}}_{1}={\begin{pmatrix}y_{1}\\y_{1}'\\y_{1}''\end{pmatrix}}={\begin{pmatrix}y_{0}+h\cdot y'_{0}\\y'_{0}+h\cdot y''_{0}\\y''_{0}+h\cdot f\left(t_{0},y_{0},y'_{0},y''_{0}\right)\end{pmatrix}}={\begin{pmatrix}2+0.5\cdot -1\\-1+0.5\cdot 3\\3+0.5\cdot \left(\sin {0}+(\cos {0})\cdot 2+0^{2}\cdot (-1)-4\cdot 0\cdot 3\right)\end{pmatrix}}={\begin{pmatrix}1.5\\0.5\\4\end{pmatrix}}}$$

${\displaystyle {\vec {y}}_{1}}$ ${\displaystyle {\vec {y}}_{2}}$

$${\displaystyle {\vec {y}}_{2}={\begin{pmatrix}y_{2}\\y_{2}'\\y_{2}''\end{pmatrix}}={\begin{pmatrix}y_{1}+h\cdot y'_{1}\\y'_{1}+h\cdot y''_{1}\\y''_{1}+h\cdot f\left(t_{1},y_{1},y'_{1},y''_{1}\right)\end{pmatrix}}={\begin{pmatrix}1.5+0.5\cdot 0.5\\0.5+0.5\cdot 4\\2+0.5\cdot \left(\sin {0.5}+(\cos {0.5})\cdot 1.5+0.5^{2}\cdot 0.5-4\cdot 0.5\cdot 2\right)\end{pmatrix}}={\begin{pmatrix}1.75\\2.5\\0.9604\end{pmatrix}}}$$

${\displaystyle {\vec {y}}_{i}}$

Derivación

El método de Euler se puede derivar de varias formas.

(1) En primer lugar, está la descripción geométrica anterior.

(2) Otra posibilidad es considerar la expansión de Taylor de la función alrededor de : ${\displaystyle y}$ ${\displaystyle t_{0}}$

${\displaystyle y(t_{0}+h)=y(t_{0})+hy'(t_{0})+{\tfrac {1}{2}}h^{2}y''(t_{0})+O\left(h^{3}\right).}$

La ecuación diferencial establece que . Si esto se sustituye en la expansión de Taylor y se ignoran los términos cuadráticos y de orden superior, surge el método de Euler. ^[7] ${\displaystyle y'=f(t,y)}$

La expansión de Taylor se utiliza a continuación para analizar el error cometido por el método de Euler y se puede ampliar para producir métodos de Runge-Kutta .

(3) Una derivación estrechamente relacionada es sustituir la fórmula de diferencias finitas directas por la derivada,

${\displaystyle y'(t_{0})\approx {\frac {y(t_{0}+h)-y(t_{0})}{h}}}$

en la ecuación diferencial . De nuevo, esto produce el método de Euler. ^[8] ${\displaystyle y'=f(t,y)}$

Un cálculo similar conduce al método del punto medio y al método de Euler hacia atrás .

(4) Finalmente, se puede integrar la ecuación diferencial de a y aplicar el teorema fundamental del cálculo para obtener: ${\displaystyle t_{0}}$ ${\displaystyle t_{0}+h}$

${\displaystyle y(t_{0}+h)-y(t_{0})=\int _{t_{0}}^{t_{0}+h}f{\bigl (}t,y(t){\bigr )}\,\mathrm {d} t.}$

Ahora aproxime la integral mediante el método del rectángulo de la izquierda (con un solo rectángulo):

${\displaystyle \int _{t_{0}}^{t_{0}+h}f{\bigl (}t,y(t){\bigr )}\,\mathrm {d} t\approx hf{\bigl (}t_{0},y(t_{0}){\bigr )}.}$

Combinando ambas ecuaciones, se encuentra nuevamente el método de Euler. ^[9]

Esta línea de pensamiento puede continuarse para llegar a varios métodos lineales de varios pasos .

Error de truncamiento local

El error de truncamiento local del método de Euler es el error cometido en un solo paso. Es la diferencia entre la solución numérica después de un paso, y la solución exacta en el momento . La solución numérica viene dada por ${\displaystyle y_{1}}$ ${\displaystyle t_{1}=t_{0}+h}$

${\displaystyle y_{1}=y_{0}+hf(t_{0},y_{0}).}$

Para la solución exacta, utilizamos la expansión de Taylor mencionada en la sección Derivación anterior:

${\displaystyle y(t_{0}+h)=y(t_{0})+hy'(t_{0})+{\tfrac {1}{2}}h^{2}y''(t_{0})+O\left(h^{3}\right).}$

El error de truncamiento local (LTE) introducido por el método de Euler viene dado por la diferencia entre estas ecuaciones:

${\displaystyle \mathrm {LTE} =y(t_{0}+h)-y_{1}={\tfrac {1}{2}}h^{2}y''(t_{0})+O\left(h^{3}\right).}$

Este resultado es válido si tiene una tercera derivada acotada. ^[10] ${\displaystyle y}$

Esto muestra que para valores pequeños , el error de truncamiento local es aproximadamente proporcional a . Esto hace que el método de Euler sea menos preciso que las técnicas de orden superior, como los métodos de Runge-Kutta y los métodos lineales de varios pasos , para los cuales el error de truncamiento local es proporcional a una potencia mayor del tamaño del paso. ${\displaystyle h}$ ${\displaystyle h^{2}}$

Se puede obtener una formulación ligeramente diferente para el error de truncamiento local utilizando la forma de Lagrange para el término restante del teorema de Taylor . Si tiene una segunda derivada continua, entonces existe una tal que ${\displaystyle y}$ ${\displaystyle \xi \in [t_{0},t_{0}+h]}$

${\displaystyle \mathrm {LTE} =y(t_{0}+h)-y_{1}={\tfrac {1}{2}}h^{2}y''(\xi ).}$^[11]

En las expresiones anteriores para el error, la segunda derivada de la solución exacta desconocida se puede reemplazar por una expresión que involucre el lado derecho de la ecuación diferencial. De hecho, de la ecuación se deduce que ^[12] ${\displaystyle y}$ ${\displaystyle y'=f(t,y)}$

${\displaystyle y''(t_{0})={\frac {\partial f}{\partial t}}{\bigl (}t_{0},y(t_{0}){\bigr )}+{\frac {\partial f}{\partial y}}{\bigl (}t_{0},y(t_{0}){\bigr )}\,f{\bigl (}t_{0},y(t_{0}){\bigr )}.}$

Error de truncamiento global

El error de truncamiento global es el error en un momento fijo , después de tantos pasos que el método deba seguir para llegar a ese momento desde el momento inicial. El error de truncamiento global es el efecto acumulativo de los errores de truncamiento local cometidos en cada paso. ^[13] Se determina fácilmente que el número de pasos es , que es proporcional a , y el error cometido en cada paso es proporcional a (ver la sección anterior). Por tanto, es de esperar que el error de truncamiento global sea proporcional a . ^[14] ${\displaystyle t_{i}}$ ${\textstyle {\frac {t_{i}-t_{0}}{h}}}$ ${\textstyle {\frac {1}{h}}}$ ${\displaystyle h^{2}}$ ${\displaystyle h}$

Este razonamiento intuitivo puede precisarse. Si la solución tiene una segunda derivada acotada y es continua de Lipschitz en su segundo argumento, entonces el error de truncamiento global (denotado como ) está acotado por ${\displaystyle y}$ ${\displaystyle f}$ ${\displaystyle |y(t_{i})-y_{i}|}$

${\displaystyle |y(t_{i})-y_{i}|\leq {\frac {hM}{2L}}\left(e^{L(t_{i}-t_{0})}-1\right)}$

donde es un límite superior de la segunda derivada de en el intervalo dado y es la constante de Lipschitz de . ^[15] O más simplemente, cuando , el valor (tal que se trata como una constante). Por el contrario, donde función es la solución exacta que solo contiene la variable. ${\displaystyle M}$ ${\displaystyle y}$ ${\displaystyle L}$ ${\displaystyle f}$ ${\displaystyle y'(t)=f(t,y)}$ ${\textstyle L={\text{max}}{\bigl (}|{\frac {d}{dy}}{\bigl [}f(t,y){\bigr ]}|{\bigr )}}$ ${\displaystyle t}$ ${\textstyle M={\text{max}}{\bigl (}|{\frac {d^{2}}{dt^{2}}}{\bigl [}y(t){\bigr ]}|{\bigr )}}$ ${\displaystyle y(t)}$ ${\displaystyle t}$

La forma precisa de este límite tiene poca importancia práctica, ya que en la mayoría de los casos el límite sobreestima enormemente el error real cometido por el método de Euler. ^[16] Lo importante es que muestra que el error de truncamiento global es (aproximadamente) proporcional a . Por este motivo se dice que el método de Euler es de primer orden. ^[17] ${\displaystyle h}$

Ejemplo

Si tenemos la ecuación diferencial y la solución exacta , sabemos qué encontrar y para cuándo . ${\displaystyle y'=1+(t-y)^{2}}$ ${\displaystyle y=t+{\frac {1}{t-1}}}$ ${\displaystyle M}$ ${\displaystyle L}$ ${\displaystyle 2\leq t\leq 3}$

$${\displaystyle L={\text{max}}{\bigl (}|{\frac {d}{dy}}{\bigl [}f(t,y){\bigr ]}|{\bigr )}=\max _{2\leq t\leq 3}{\bigl (}|{\frac {d}{dy}}{\bigl [}1+(t-y)^{2}{\bigr ]}|{\bigr )}=\max _{2\leq t\leq 3}{\bigl (}|2(t-y)|{\bigr )}=\max _{2\leq t\leq 3}{\bigl (}|2(t-[t+{\frac {1}{t-1}}])|{\bigr )}=\max _{2\leq t\leq 3}{\bigl (}|-{\frac {2}{t-1}}|{\bigr )}=2}$$

$${\displaystyle M={\text{max}}{\bigl (}|{\frac {d^{2}}{dt^{2}}}{\bigl [}y(t){\bigr ]}|{\bigr )}=\max _{2\leq t\leq 3}\left(|{\frac {d^{2}}{dt^{2}}}{\bigl [}t+{\frac {1}{1-t}}{\bigr ]}|\right)=\max _{2\leq t\leq 3}\left(|{\frac {2}{(-t+1)^{3}}}|\right)=2}$$

$${\displaystyle {\text{error bound}}={\frac {hM}{2L}}\left(e^{L(t_{i}-t_{0})}-1\right)={\frac {0.5\cdot 2}{2\cdot 2}}\left(e^{2(2.5-2)}-1\right)=0.42957}$$

₀ ${\displaystyle 2\leq t\leq 3}$

Estabilidad numérica

(Figura 4) Solución calculada con el método de Euler con tamaño de paso (cuadrados azules) y (círculos rojos). La curva negra muestra la solución exacta. ${\displaystyle y'=-2.3y}$ ${\displaystyle h=1}$ ${\displaystyle h=0.7}$

El método de Euler también puede ser numéricamente inestable , especialmente para ecuaciones rígidas , lo que significa que la solución numérica crece mucho para ecuaciones donde la solución exacta no lo hace. Esto se puede ilustrar usando la ecuación lineal

${\displaystyle y'=-2.3y,\qquad y(0)=1.}$

La solución exacta es , que decae a cero cuando . Sin embargo, si se aplica el método de Euler a esta ecuación con un tamaño de paso , entonces la solución numérica es cualitativamente incorrecta: oscila y crece (ver figura). Esto es lo que significa ser inestable. Si se utiliza un tamaño de paso más pequeño, por ejemplo , entonces la solución numérica decae a cero. ${\displaystyle y(t)=e^{-2.3t}}$ ${\displaystyle t\to \infty }$ ${\displaystyle h=1}$ ${\displaystyle h=0.7}$

(Figura 5) El disco rosa muestra la región de estabilidad del método de Euler.

Si se aplica el método de Euler a la ecuación lineal , entonces la solución numérica es inestable si el producto está fuera de la región ${\displaystyle y'=ky}$ ${\displaystyle hk}$

${\displaystyle {\bigl \{}z\in \mathbf {C} \,{\big |}\,|z+1|\leq 1{\bigr \}},}$

ilustrado a la derecha. Esta región se llama región de estabilidad (lineal) . ^[18] En el ejemplo, entonces, si entonces cuál está fuera de la región de estabilidad y, por lo tanto, la solución numérica es inestable. ${\displaystyle k=-2.3}$ ${\displaystyle h=1}$ ${\displaystyle hk=-2.3}$

Esta limitación, junto con su lenta convergencia del error , significa que el método de Euler no se usa con frecuencia, excepto como un simple ejemplo de integración numérica ^{[ cita requerida ]} . ${\displaystyle h}$

Errores de redondeo

En el paso del método de Euler, el error de redondeo es aproximadamente de la magnitud donde está la máquina épsilon . Suponiendo que los errores de redondeo son variables aleatorias independientes, el error de redondeo total esperado es proporcional a . ^[19] Por lo tanto, para valores extremadamente pequeños del tamaño del paso, el error de truncamiento será pequeño, pero el efecto del error de redondeo puede ser grande. La mayor parte del efecto del error de redondeo se puede evitar fácilmente si se utiliza la suma compensada en la fórmula del método de Euler. ^[20] ${\displaystyle n}$ ${\displaystyle \varepsilon y_{n}}$ ${\displaystyle \varepsilon }$ ${\textstyle {\frac {\varepsilon }{\sqrt {h}}}}$

Modificaciones y ampliaciones

Una modificación simple del método de Euler que elimina los problemas de estabilidad señalados anteriormente es el método de Euler hacia atrás :

${\displaystyle y_{n+1}=y_{n}+hf(t_{n+1},y_{n+1}).}$

Esto difiere del método de Euler (estándar o directo) en que la función se evalúa en el punto final del paso, en lugar del punto inicial. El método de Euler hacia atrás es un método implícito , lo que significa que la fórmula del método de Euler hacia atrás tiene ambos lados, por lo que al aplicar el método de Euler hacia atrás tenemos que resolver una ecuación. Esto hace que la implementación sea más costosa. ${\displaystyle f}$ ${\displaystyle y_{n+1}}$

Otras modificaciones del método de Euler que ayudan con la estabilidad producen el método exponencial de Euler o el método semiimplícito de Euler .

Los métodos más complicados pueden lograr un orden superior (y más precisión). Una posibilidad es utilizar más evaluaciones de funciones. Esto se ilustra con el método del punto medio que ya se menciona en este artículo:

${\displaystyle y_{n+1}=y_{n}+hf\left(t_{n}+{\tfrac {1}{2}}h,y_{n}+{\tfrac {1}{2}}hf(t_{n},y_{n})\right)}$.

Esto conduce a la familia de métodos de Runge-Kutta .

La otra posibilidad es utilizar más valores pasados, como lo ilustra el método de dos pasos Adams-Bashforth:

${\displaystyle y_{n+1}=y_{n}+{\tfrac {3}{2}}hf(t_{n},y_{n})-{\tfrac {1}{2}}hf(t_{n-1},y_{n-1}).}$

Esto conduce a la familia de métodos lineales de varios pasos . Existen otras modificaciones que utilizan técnicas de detección de compresión para minimizar el uso de memoria ^[21]

En la cultura popular

En la película Figuras ocultas , Katherine Goble recurre al método de Euler para calcular el reingreso del astronauta John Glenn desde la órbita terrestre. ^[22]

Ver también

• Método Crank-Nicolson
• El descenso de gradiente utiliza de manera similar pasos finitos, aquí para encontrar mínimos de funciones
• Lista de métodos de Runge-Kutta
• Método lineal de varios pasos
• Integración numérica (para calcular integrales definidas)
• Métodos numéricos para ecuaciones diferenciales ordinarias.

Notas

1. Carnicero 2003, pag. 45; Hairer, Nørsett y Wanner 1993, pág. 35
2. Atkinson 1989, pág. 342; Carnicero 2003, pág. 60
3. Carnicero 2003, pag. 45; Hairer, Nørsett y Wanner 1993, pág. 36
4. Carnicero 2003, pag. 3; Hairer, Nørsett y Wanner 1993, pág. 2
5. Véase también Atkinson 1989, p. 344
6. Hairer, Nørsett y Wanner 1993, pág. 40
7. Atkinson 1989, pág. 342; Hairer, Nørsett y Wanner 1993, pág. 36
8. Atkinson 1989, pág. 342
9. Atkinson 1989, pág. 343
10. Carnicero 2003, pag. 60
11. Atkinson 1989, pág. 342
12. Stoer y Bulirsch 2002, pág. 474
13. Atkinson 1989, pág. 344
14. Carnicero 2003, pag. 49
15. Atkinson 1989, pág. 346; Lakoba 2012, ecuación (1.16)
16. Iserles 1996, pág. 7
17. Carnicero 2003, pag. 63
18. Carnicero 2003, pag. 70; Iserles 1996, pág. 57
19. Carnicero 2003, págs. 74–75
20. Carnicero 2003, págs. 75–78
21. Unni, diputado; Chandra, MG; Kumar, AA (marzo de 2017). "Reducción de memoria para solución numérica de ecuaciones diferenciales mediante detección de compresión". 2017 IEEE 13º Coloquio Internacional sobre Procesamiento de Señales y sus Aplicaciones (CSPA) . págs. 79–84. doi :10.1109/CSPA.2017.8064928. ISBN 978-1-5090-1184-1. S2CID  13082456.
22. Khan, Amina (9 de enero de 2017). "Conozca al matemático de 'Figuras ocultas' que ayudó a enviar a los estadounidenses al espacio". Los Ángeles Times . Consultado el 12 de febrero de 2017 .

Referencias

• Atkinson, Kendall A. (1989). Introducción al análisis numérico (2ª ed.). Nueva York: John Wiley & Sons . ISBN 978-0-471-50023-0.
• Ascher, Uri M.; Petzold, Linda R. (1998). Métodos informáticos para ecuaciones diferenciales ordinarias y ecuaciones algebraicas diferenciales. Filadelfia: Sociedad de Matemáticas Industriales y Aplicadas . ISBN 978-0-89871-412-8.
• Carnicero, John C. (2003). Métodos numéricos para ecuaciones diferenciales ordinarias . Nueva York: John Wiley & Sons . ISBN 978-0-471-96758-3.
• Hairer, Ernst; Norsett, Syvert Paul; Wanner, Gerhard (1993). Resolución de ecuaciones diferenciales ordinarias I: Problemas no rígidos . Berlín, Nueva York: Springer-Verlag . ISBN 978-3-540-56670-0.
• Iserles, Arieh (1996). Un primer curso de análisis numérico de ecuaciones diferenciales. Prensa de la Universidad de Cambridge . ISBN 978-0-521-55655-2.
• Stoer, Josef; Bulirsch, Roland (2002). Introducción al análisis numérico (3ª ed.). Berlín, Nueva York: Springer-Verlag . ISBN 978-0-387-95452-3.
• Lakoba, Taras I. (2012), Método simple de Euler y sus modificaciones (PDF) (Apuntes de la conferencia de MATH334), Universidad de Vermont , consultado el 29 de febrero de 2012
• Unni, MP (2017). "Reducción de memoria para solución numérica de ecuaciones diferenciales mediante detección de compresión". 2017 IEEE 13º Coloquio Internacional sobre Procesamiento de Señales y sus Aplicaciones (CSPA) . IEEE CSPA. págs. 79–84. doi :10.1109/CSPA.2017.8064928. ISBN 978-1-5090-1184-1. S2CID  13082456.

enlaces externos

• Medios relacionados con el método de Euler en Wikimedia Commons
• Implementaciones del método Euler en diferentes lenguajes mediante Rosetta Code
• "Método de Euler", Enciclopedia de Matemáticas , EMS Press , 2001 [1994]