scipy.integrate.solve_ivp

With the new interface for solving ODE in scipy.integrate, the argument ordering(t,X) is different than the one used in pydy(which is obj(X,t)

It would be nice to map the output of solve_ivp since it supports ode45 and other integrators.
Meanwhile, odeint and ode has been marked as old api so I'm assuming that it indicates we embracing the new interface for future compatibility.