# Goldbeter glycolysis model
du/dt = nu-f(u,w)
dw/dt = a*f(u,w)-eta*w
f(u,w)=u*(1+w)^2
param nu=200,eta=120,a=1.0
init u=.5,w=0.5
@ total=8,dt=0.002
@ xplot=u,yplot=w,xhi=40,yhi=20,xlo=0,ylo=0
done