differentiate dy/dt = k(2x^0)(e^-k1t) - (k3y/y) + M