%File michaelis_menten.ode
%Figure 3.3 Michaelis-Menten kinetics

%declare model parameters
par et=1, k1=30, km1=1, k2=10

%dynamics for full model
s'=-k1*s*(et-c) + km1*c 
c'=-km1*c + k1*s*(et-c)-k2*c 
p'=k2*c

%dynamics for reduced model 
s2'=-k1*k2*et*s2/(km1+k2+k1*s2)

%auxilliary variables for plotting
aux e=et-c
aux p2=5-s2


%set initial condition. State vector is
init s=5, c=0, p=0, s2=5


# set the plotting window size:
@ xlo=0, xhi=1, ylo=0, yhi=5

#set an appropriate stepsize
@ dt=0.01
done