# File rapid_equilibrium_approximation2.ode
# Figure 2.10 Rapid equilibrium approximation

#assign parameter values
par k0=5, k1=20, km1=12, k2=2

# equations
# original model
a'=k0 - k1*a + km1*b
b'=k1*a - km1*b - k2*b
# approximate model
c'=k0 - (k2*k1)/(km1+k1)*c
# approximate concentrations
aux a_tilde=km1/(k1+km1)*c
aux b_tilde=k1/(k1+km1)*c

# initial conditions
init a=8, b=4, c=12

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

# 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