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

