# file predator_prey.ode
# model of predator-prey network
# from Balagadde et al. (2008) Molecular Systems Biology 4
# problem 7.8.20 (b)

par k2=0.4
par cm=1e5
par d1=1 
par d2=0.3
par Ka1=10
par beta=2
par D=0.2 
par gamma1=0.1
par gamma2=0.1 
par del1=0.017
par del2=0.11
par k1=0.8

#model equations
c1' = k1*c1*(1-(c1+c2)/cm)-d1*c1*Ka1/(Ka1+a2^beta) - D*c1
c2' = k2*c2*(1-(c1+c2)/cm)-d2*c2*a1^beta/(Ka1+a1^beta) - D*c2
a1' = gamma1*c1-(del1+D)*a1
a2' = gamma2*c2-(del2+D)*a2

# initial condition
init c1=10000, c2=10000, a1=1, a2=1

# final time
@ total=500
# set the plotting window size:
@ xlo=0, xhi=500, ylo=0, yhi=10000
@ bound=100000
@ maxstor=1000000
# run the simulation in the GUI by selecting (I)nitialconds|(G)o
# plot other species by selecting (G)raphic stuff|(A)dd curve and 
# assigning the y-axis

done