Hola estoy empezando a programar con Maxima y tengo una ligera idea de lo que tengo que hacer, pero  Me gustaría que alguien me dijera como puedo programar una función genérica(valga para cualquier sistema de ecuaciones) que resuelva numericamente sistemas de ecuaciones diferenciales ordinarias de orden 1 mediante rk y obtenga los resultados en forma de matriz.
Muchas gracias.