'Set callable to be called at every successful integration step.
Parameters
solout : callable
``solout(t, y)`` is called at each internal integrator step,
t is a scalar providing the current independent position
y is the current soloution ``y.shape == (n,)``
solout should return -1 to stop integration
otherwise it should return None or 0'

def set_solout(self, solout):

    if self._integrator.supports_solout:
        self._integrator.set_solout(solout)
        if (self._y is not None):
            self._integrator.reset(len(self._y), (self.jac is not None))
    else:
        raise ValueError('selected    integrator    does    not    support    solout,    choose    another    one')
