# file population_control.ode
# model of population control network
# from You et al. (2004) Nature 428 pp. 868-871 
# problem 7.8.20 (a)

par k=0.97
par Nm=124000000
par d=0.004
par ke=5
par de=2
par va=0.00000048
par da=0.64

#model equations
N' = k*N*(1-N/Nm)-d*E*N
E' = ke*A-de*E
A' = va*N-da*A

# initial condition
init N=1, E=1, A=1

# final time
@ total=100
# set the plotting window size:
@ xlo=0, xhi=100, ylo=0, yhi=30000000
@ bound=1000000000

# 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