scipy.integrate.ode.integrate

ode.integrate(t, step=0, relax=0)

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