# 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