# Selkov glycolysis model
du/dt = nu-f(u,w)
dw/dt = a*f(u,w)-eta*w
f(u,w)=u*w^g/(w^g*u+w^g+1)
param nu=0.0285,eta=0.1,a=1.0,g=2
init u=.5,w=0.5
@ total=800,dt=0.2
@ xplot=u,yplot=w,xhi=1.4,ylo=0
done