# Goldbeter-Lefever
#

dz/dt = qsig - kap*y
dy/dt = (z-y)*(1+y)^2 - kap*y

par qsig=100, kap=115

done