I was reading the section 6.7.2 on numerical integration of differential equations and I feel like it might need some updating. The recommended way of doing it is to use scipy.integrate.odeint, but ...