Estimado Mauro,
El comando que vamos a usar en el curso es lsode que solo sirve para EDOs de primer orden. Dado el caso de tener un sistema de segundo orden hay que reescribir el problema como de primer orden.
Según mi conocimiento Octave no tiene una rutina que resuelva sistemas de orden superior. La ayuda de Octave solo describe rutina para EDOs y EDAs (https://octave.org/doc/v4.0.0/Differential-Equations.html#Differential-Equations).
Saludos,
Jimena F.