solucion numerica

Solución numérica de las ecuaciones diferenciales

Las leyes que gobiernan los fenómenos de la naturaleza se expresan habitualmente en forma de ecuaciones diferenciales. Lasecuaciones del movimiento de los cuerpos (la segunda ley de Newton) es una ecuación diferencial de segundo orden, como lo es la ecuación que describe los sistemas oscilantes, la propagación delas ondas, la transmisión del calor, la difusión, el movimiento de partículas subatómicas, etc.

Pocas ecuaciones diferenciales tienen una solución analítica sencilla, la mayor parte delas veces es necesario realizar aproximaciones, estudiar el comportamiento del sistema bajo ciertas condiciones. Así, en un sistema tan simple como un péndulo, la amplitud de laoscilación ha de ser pequeña y el rozamiento ha de ser despreciable, para obtener una solución sencilla que describa aproximadamente su movimiento periódico.

Las páginas que vienen acontinuación, se dedicarán al estudio de la resolución de ecuaciones diferenciales por el procedimiento de Runge-Kutta que se aplica de forma directa a una ecuación diferencial de primer orden,pero veremos como se extienden a un sistema de ecuaciones de primer orden, a un ecuación diferencial de segundo orden, y a un sistema de ecuaciones diferenciales de segundo orden.Crearemos una clase base abstracta con una función miembro que describe el procedimiento numérico, y una clase derivada que describirá el sistema físico particular. Para determinar el estadofinal de un sistema físico, cuyo comportamiento está descrito por una ecuación diferencial, se creará un objeto de la clase derivada, y desde éste llamaremos a la función que describeel procedimiento numérico, pasándole en uno de sus parámetros el objeto que describe el estado inicial, devolviendo en el mismo objeto el estado final cuando dicha función retorna.