function yp = MMspode(t,y,p)
%Michaelis-Menten ODE
ep = p(1); mu = p(2); lam = p(3);
yp1 = -y(1) + (mu + y(1))*y(2);
yp2 = (y(1) - (lam + y(1))*y(2))/ep;
yp = [yp1; yp2];
end

