Python 或 SageMath 中耦合非线性微分方程的解？

Question

How to solve four first order coupled non-linear differential equations using python or SageMath? Please suggest some references if available. Thank you very much for answering.

Answer 1

You can use Python and the scipy module's odeint, import as follows:

from scipy.integrate import odeint

A fun example is that of "Modeling a zombie apocalypse, see: https://scipy-cookbook.readthedocs.io/items/Zombie_Apocalypse_ODEINT.html