

def jac(self, z, t):
    lmbd = self.lmbd
    j = np.zeros((((self.lband + self.uband) + 1), 3), order='F')

    def set_j(ri, ci, val):
        j[(((self.uband + ri) - ci), ci)] = val
    set_j(0, 0, (- lmbd[0]))
    set_j(1, 0, lmbd[0])
    set_j(1, 1, (- lmbd[1]))
    set_j(2, 1, lmbd[1])
    set_j(2, 2, (- lmbd[2]))
    return j

dataset/ETHPy150Open scipy/scipy/scipy/integrate/tests/test_integrate.py CoupledDecay.jac/original