# goldbeter.ode
#
# Dimensional form of the Goldbeter-Lefever model for glycolytic
# oscillations via ADP feedback onto phosphofructokinase 1.
# See p. 21-24 of Keener-Sneyd book.

ATP(0)=0.5
ADP(0)=0.5

# substrate (ATP) influx rate
par nu=1

# product (ADP) removal rate
par eta=120

# ratio k_2/k_(-2)
par k=20

# time scale parameter (in sec)
par tauc=500

# Functions

fun = ATP*(1+k*ADP)^2

# Equations

ATP' = (nu - fun)/tauc 
ADP' = (fun - eta*ADP)/tauc

@ meth=cvode, toler=1.0e-10, atoler=1.0e-10, dt=0.2, total=500, maxstor=20000
@ bounds=10000000, xp=t, yp=ATP, bell=off
@ xlo=0, xhi=500, ylo=0, yhi=4

done