# Tyson et al. (Biophys J 77:2411, 1999)
# Model of circadian rhythms


P Vm=1,Km=.1,Vp=.5,Kp2=0.03,Kp1=10,Kp3=.1
P Keq=200,Jp=0.05,Pcrit=0.1

dP/dt = Vp*M - (Kp2*P+Kp1*P*q)/(Jp+P) - Kp3*P
dM/dt = Vm/(1+(P*(1-q)/(2*Pcrit))^2) - Km*M
q = 2/(1+sqrt(1+8*Keq*P))

done