Hola!

Para el que tenga ganas de animársele a las soluciones numéricas y Python les dejo un ejemplo en Google Colab (el link es view only, tal vez les convenga hacer una copia) que se puede usar para resolver (casi) cualquier sistemas de ODEs de 1er orden.

https://colab.research.google.com/drive/1N_-Wi4fwaao-XQJrzXetCLD-u6F1ncxP?usp=sharing

Cualquier cosa pregunten!

Saludos,

Patricio