#file synchronized_oscillators.ode
#Hasty synthetic oscillator model
#from Hasty et al. (2002) Phys. Rev. Lett. 88 148101
#Figure 7.31

par alpha=11
par sigma=2
par gammax=0.2
par gammay=0.012
par ay=0.2
par D=0.015

#model equations
x' = (1+x^2+alpha*sigma*x^4)/((1+x^2+sigma*x^4)*(1+y^4)) - gammax*x + D*(xx-x)
y' = ay*((1+x^2+alpha*sigma*x^4)/((1+x^2+sigma*x^4)*(1+y^4))) - gammay*y
xx' = (1+xx^2+alpha*sigma*xx^4)/((1+xx^2+sigma*xx^4)*(1+yy^4)) - gammax*xx + D*(x-xx)
yy' = ay*((1+xx^2+alpha*sigma*xx^4)/((1+xx^2+sigma*xx^4)*(1+yy^4))) - gammay*yy


# initial condition
init x=0.3963, y=2.3346, xx=0.5587, yy=1.9317

# final time
@ total=500
# set the plotting window size:
@ xlo=0, xhi=500, ylo=0, yhi=3
@ maxstor=100000

# 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