# Martiel and Goldbeter

W=y*z/(1+y)
WW=alpha*W*W
phi=(eps+WW)/(1+W*W+WW)
f1=(k1+k2*y)/(1+y)
f2=(k1'+k2'*c*y)/(1+c*y)

dx/dt=q*phi-kit*x
dz/dt=f2*(1-z)-f1*z

p q=2400, kit=2.6, kht=0.18, ke=5.4
p k1=0.36, k2=0.0033, k1'=0.036, k2'=0.67
p eps=0.0003, alpha=3, c=10
p y=0

done