# full cell cycle mode from page 271
m'=mu*m*(1-m/mstar)
cycbt'=k1-(k2p+k2pp*cdh1+k2ppp*cdc20a)*cycbt
cdh1'=(k3p+k3pp*cdc20a)*(1-cdh1)/(j3+1-cdh1)-(k4*m*cycb_+k4p*sk)*cdh1/(j4+cdh1)
cdc20t'=k5p+k5pp*(m*cycb_)^n/(j5^n+(m*cycb_)^n)-k6*cdc20t
cdc20a'=k7*iep*(cdc20t-cdc20a)/(j7+cdc20t-cdc20a)-k8*mad*cdc20a/(j8+cdc20a)-k6*cdc20a
iep'=k9*m*cycb_*(1-iep)-k10*iep
ckit'=k11-(k12p+k12pp*sk+k12ppp*m*cycb_)*ckit
sk'=k13p+k13pp*tf-k14*sk
tf'=(k15p*m+k15pp*sk)*(1-tf)/(j15+1-tf)-(k16p+k16pp*m*cycb_)*tf/(j16+tf)
sigma=cycbt+ckit+1/keq
trimer=2*cycbt*ckit/(sigma+sqrt(sigma^2-4*cycbt*ckit))
cycb_=cycbt-trimer
# this  gives you cycb as a plottable quantity
aux cycb=cycb_
global -1 cycbt-trimer-cycbth {m=m/2}
# cycb
par k1=.04,k2p=.04,k2pp=1,k2ppp=2,cycbth=0.1
# cdh1
par k3p=1,k3pp=10,k4p=2,k4=35,j3=.04,j4=.04
# cdc20t
par k5p=.005,k5pp=.2,k6=.1,j5=.3,n=4
# cdc20a
par k7=1,k8=.5,j7=.001,j8=.001,mad=1
# iep 
par k9=.1,k10=.02
# ckit
par k11=1,k12p=.2,k12pp=50,k12ppp=100,keq=1000
# sk
par k13p=0,k13pp=1,k14=1,k15p=1.5,k15pp=.05,k16p=1,k16pp=3,j15=.01,j16=.01
# m
par mu=.005,mstar=10
init m=.55,cycbt=.04,sk=.04,cdh1=.99,cdc20t=.05,cdc20a=.004,iep=.15,ckit=.41
init tf=.045

@ meth=cvode,total=500,tol=1e-6,atol=1e-6,dt=.01
@ poimap=per,poivar=cycb,poipln=.105,poisgn=0
@ xhi=300,yhi=1,ylo=0,yp=cycb,yp2=m,nplots=2
done