'Set initial conditions y(t) = y.'

def set_initial_value(self, y, t=0.0):

    if isscalar(y):
        y = [y]
    n_prev = len(self._y)
    if (not n_prev):
        self.set_integrator('')
    self._y = asarray(y, self._integrator.scalar)
    self.t = t
    self._integrator.reset(len(self._y), (self.jac is not None))
    return self
