function dxdt = michaelis_mentin_rescaled_odefun(t,x)

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

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


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

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