SciPy

scipy.integrate.ode.integrate

ode.integrate(t, step=0, relax=0)[source]

Find y=y(t), set y as an initial condition, and return y.