# File rapid_equilibrium_approximation1.ode
# Figure 2.11, 2.12 Simulation and rapid equilibrium approximation

#assign parameter values
par k1=9, km1=12, k2=2

# equations
# original model
a'=-k1*a+km1*b
b'=k1*a-km1*b-k2*b
c'=-k2*k1/(k1+km1)*c

aux a_tilde=km1/(k1+km1)*c
aux b_tilde=k1/(k1+km1)*c

# initial conditions
init a=0, b=10, c=10

# final time
@ total=3
# set the plotting window size:
@ xlo=0, xhi=3, ylo=0, yhi=10

# run the simulation in the GUI by selecting (I)nitialconds|(G)o
# plot species b, a_tilde and b_tilde by selecting (G)raphic stuff|(A)dd curve 
# and assigning the y-axis

done