Hace años trabajé en este tema. Recuerdo que el mejor procedimiento para la integración era recurrir a los métodos de Runge-Kutta (quizá ahora ya...