%y0=[S0 E0 ES0 P0]


function ydot=MMmodel(t,y)
global k1,k2,k3; %k2=k-1, k3=k2
ydot(1)=-k1*y(1)*y(2)+k2*y(3); %substrate
ydot(2)=-k1*y(1)*y(2)+(k2+k3)*y(3);%enzyme
ydot(3)=k1*y(1)*y(2)-(k3+k2)*y(3);%enzyme complex
ydot(4)=k3*y(3);%product

ydot=ydot';
