function dxdt = michaelis_mentin_odefun(t,x)

%For michaelis-mentin equations, p. 10 of EG, eqns [1.10]
%Takes as arguments the current time t (here, irrelevant) as well as the
%state vector x
%x(1)=c
%x(2)=s
%Returns dx/dt

c=x(1);
s=x(2);


%define parameters
k1=1 ;
km1=1  ;  %k_{-1}
k2=1;
e0=.1;

dxdt = [k1.*(e0-c).*s - (km1+k2).*c ; 
        km1.*c - k1.*(e0-c).*s];