#file Goodwin.ode
#Goodwin oscillator model
#from Goodwin (1965) Adv. in. Enz. Reg. 3 pp. 425-428.
#Figure 7.17A

par a1=360, kappa1=43, k1=1, b1=1
par alpha1=1, beta1=0.6, gamma1=1, delta1=0.8, n=12

#model equations
X' = a1/(kappa1+k1*Z^n) - b1*X
Y' = alpha1*X - beta1*Y
Z' =  gamma1*Y - delta1*Z

# initial condition
init X=0, Y=0, Z=0

# final time
@ total=35
# set the plotting window size:
@ xlo=0, xhi=35, ylo=0, yhi=6

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


done