def f(self, z, t):
        lmbd = self.lmbd
        return np.array([-lmbd[0]*z[0],
                         -lmbd[1]*z[1] + lmbd[0]*z[0],
                         -lmbd[2]*z[2] + lmbd[1]*z[1]])
dataset/ETHPy150Open scipy/scipy/scipy/integrate/tests/test_integrate.py/CoupledDecay.f/original