init CycBT=0.23, Cdh1=0.01, Cdc20T=0.06, Cdc20A=0.0001, IE=0.1, CKIT=0.06
init SK=0.09, TF=0.03, M=0.45

CycBT' = k1*M - (k2p+k2pp*Cdh1+k2ppp*Cdc20A)*CycBT
Cdh1'  = (k3p+k3pp*Cdc20A)*(1-Cdh1)/(J3+1-Cdh1) - (k4p*SK+k4*CycB)*Cdh1/(J4+Cdh1)
Cdc20T'= k5p + k5pp*CycB^4/(J5^4+CycB^4) - k6*Cdc20T
Cdc20A'= k7*IE*(Cdc20T-Cdc20A)/(J7+Cdc20T-Cdc20A) - k8*Mad1*Cdc20A/(J8+Cdc20A)- k6*Cdc20A
IE'    = k9*(1-IE)*CycB - k10*IE
CKIT'  = k11 - (k12p+k12pp*SK+k12ppp*CycB)*CKIT
SK'    = k13p + k13pp*TF - k14*SK
TF'    = (k15p*M+k15pp*SK)*(1-TF)/(J15+1-TF)-(k16p+k16pp*CycB)*TF/(J16+TF)

M'     = mu*M*(1-M/Mstar)

BB     = CycBT+CKIT+Kdiss
Trimer = 2*CycBT*CKIT/(BB+(BB^2-4*CycBT*CKIT)^0.5)
CycB   = (CycBT-Trimer)




p k1=0.04,k2p=0.04,k2pp=1,k2ppp=1
p k3p=1,k3pp=10,k4p=2,k4=35
p k5p=0.005,k5pp=0.2,k6=0.1
p J3=0.04,J4=0.04,J5=0.3
p k7=1,k8=0.5,J7=0.001,J8=0.001,Mad1=1
p k9=0.1,k10=0.02
p k11=1,k12p=0.2,k12pp=50,k12ppp=100,Kdiss=.001
p k13p=0,k13pp=1,k14=1
p k15p=1.5,k15pp=0.05,k16p=1,k16pp=3,J15=0.01,J16=0.01
p mu=.005,Mstar=10


global -1 {CycBT-0.1} {M=0.5*M}