I provide two examples describing the usage of scipy.integrate.odeint which numerically solves systems of first order differential equations.
0:12 Simple first order ODE
3:09 Second order ODE (the axionic field equation)
Further documentation about odeint can be found at
https://docs.scipy.org/doc/scipy/refe...