# File phage_lambda.ode
# model of phage lambda decision switch
# Figure 7.11

#parameters
par K1=1 
par K2=0.1 
par K3=5 
par K4=0.5 
par delta_r=0.02 
par delta_c=0.02
par a=5 
par b=50


#dynamics
r'=(a + 10*a*K1*(r/2)^2)/(1+K1*(r/2)^2+ K1*K2*(r/2)^3+K3*(c/2)+K3*K4*(c/2)^2)- delta_r*r
c'=(b + b*K3*(c/2))/(1+K1*(r/2)^2+ K1*K2*(r/2)^3+K3*(c/2)+K3*K4*(c/2)^2) - delta_c*c


# initial condition
init r=100, c=100

# final time
@ total=1000
@ maxstor=100000
@ dt=0.05
@ bound=500
# set the plotting window size:
@ xlo=0, xhi=120, ylo=0, yhi=250

# 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