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