Find y=y(t), set y as an initial condition, and return y.
scipy.integrate.ode
scipy.integrate.ode.set_f_params