Scipy odeint lsoda. solve_ivp, which supports five backends (LSODA, DoPri5, Bogacki–Sham...