'Integrate from t=t0 to t=t1 using y0 as an initial condition.
Return 2-tuple (y1,t1) where y1 is the result and t=t1
defines the stoppage coordinate of the result.'

def run(self, f, jac, y0, t0, t1, f_params, jac_params):

    raise NotImplementedError('all    integrators    must    define    run(f,    jac,    t0,    t1,    y0,    f_params,    jac_params)')
