Este programa calcula numéricamente la solución de un sistema de ecuaciones diferenciales de primer orden. El método es extensible a una ecuación diferencial de orden n, que es reducible a un conjunto de n ecuaciones de orden 1 (ver apuntes de ecuaciones diferenciales). Hay que tener en cuenta ciertos aspectos del programa para un cálculo efectivo: 1. modificamos en el programa la línea "parameter (n=7)" al número de ecuaciones que tengamos. 2. Evaluamos los valores iniciales para las variables, que se nombran como y(i), donde el índice "i" varía desde uno hasta el número de ecuaciones. 3. Seleccionar el paso de cálculo, h, generalmente pequeño para garantizar la convergencia 4. Seleccionar el número de pasos máximo, npasos. El tiempo total será npasos*h. 5. Bajar hasta encontrar la subrutina DERIVS, que define el conjunto de ecuaciones diferenciales de primer orden. A las derivadas se les llama dy(i), y también varían desde uno hasta el número de ecuaciones. Es importante tener en cuenta que si las ecuaciones tienen constantes las debemos definir dentro del cuadro de la subrutina DERIVS, y no fuera, para evitar complicaciones. Tal y como muestra el ejemplo. La subrutina que calcular se basa en un método de RUnge-Kutta de cuarto orden. La solución se encuentra en el archivo "salida.dat"