Jaime Villate me ha contestado y me ha ayudado bastante, pero sigo teniendo un problema, es que necesito una función que resuelva numéricamente cualquier sistema (no uno concreto) de ecuaciones diferenciales ordinarias de orden 1 por rk, y luego ajustarla por splines cuadráticos. Mi problema es cómo introducir la función. Gracias.