lunes, 12 de mayo de 2014

Metodo de Euler para resolver ecuaciones diferenciales

http://bmi.bmt.tue.nl/sysbio/Education/Excel_Euler_simulation.pdf

Considere el problema de calcular la pendiente de una curva desconocida que comienza en un punto dado y safisface una cierta ecuación diferencial dada. Se puede pensar en la ecuación diferencial como una fórmula que nos permite calcular la pendiente de la recta tangente a la curva en cualquier punto de la curva, siempre que el punto se conozca.
La idea es que a pesar de que la curva es desconocida en principio, su punto de comienzo(al cual denotamos por A0) es conocido. Entonces, de la ecuación diferencial se puede computar la pendiente de la curva en el punto A0 y por lo tanto la recta tangente a la curva.
Ahora, dando un pequeño paso sobre dicha recta, podemos tomarnos un nuevo punto A1 y suponer que dicho punto pertenece a la curva, entonces seguimos el mismo razonamiento aplicado anteriormente y volvemos a calcular la pendiente de la recta tangente a la curva en el punto A1. Luego de varios pasos tendremos formada una curva poligonal A0A1A2A3... En general esta curva que obtenemos al aplicar el método no diverge lejos de la curva original, además el error entre ambas curvas se puede minimizar si se dan pasos muy pequeños al avanzar sobre la recta tangente a la curva y además el intervalo sobre el que trabajamos es finito(aunque las cosas son más complicadas para ecuaciones inestables, como se discute más abajo).

  PVI = 
   \left \{
      \begin{array}{rcl}
          {dy\over dx} = f(x,y)\\
           y(x_0) = y_0\\ 
           y(x_i)=?
      \end{array}
   \right .

No hay comentarios:

Publicar un comentario