# File asymmetric_network.ode
#Model of asymmetric network from Figure 4.1. 
#This code can be used to generate phase plane in Figures 4.2, 4.3, 4.4A, 
# and 4.5, and the continuation diagram in Figure 4.18

par k1=20, k2=5, k3=5, k4=5, k5=2, n=4

#model equations
s1' = k1/(1+s2^n) - k3*s1 - k5*s1
s2' = k2 - k4*s2 + k5*s1

# initial condition for Figure 4.2
init s1=0, s2=0

# final time
@ total=4
# set the plotting window size:
@ xlo=0, xhi=1.5, ylo=0, yhi=2

# run the simulation in the GUI by selecting (I)nitialconds|(G)o
# plot species s2 by selecting (G)raphic stuff|(A)dd curve and 
# assigning the y-axis

#multiple trajectories can be generated by (G)raphic stuff|(F)reeze and
#running simulations from different initial conditions

#generate a phase plane by choosing (V)iewaxis|(2)D and assigning the axis
#the nullclines can be generated with (N)ullclies|(N)ew
#a direction field can be generated with (D)ir.field/flow|(D)irect field

#the continuation diagram in Figure 4.18 can be produced by accessing Auto,
#choosing k1 as the main parameter, setting the window size and 
#parameter range, and running the continuation.

done