# File chain1.ode
# Model of two-species reaction chain (-> s1 -> s2 ->)
#
# Declaration of model parameters (optional)
par v0=5
par k1=3, k2=2
#
# Declaration of model equations  (the notation s1' can also be used)
ds1/dt = v0 - k1*s1
ds2/dt = k1*s1 - k2*s2
# Set initial conditions (optional)
init s1=1, s2=2
#
# Declaration auxiliary output variables for plotting (optional)
aux R2rate=k1*s1
#
# Set internal computational parameters (optional)
#
# Specify the variables to be plotted and the size of the plotting window
@ xp=t, yp=s1, xlo=0, xhi=5, ylo=0, yhi=3
# Specify the length of time to be simulated
@ total=5
#Allocate memory for storage of each trajectory
@ maxstor=10000
# Specifies the maximum allowable value for the state variables
@ bound=1000
#
# The file must end with the following command
done
#