M(BA)=128.17
M(BAX)=128.17
// BA polymerization

[BA] = 1.5 mol/L
[AIBN]=0.01
f(AIBN)=0.64
[LM]=9.68

kiAIBN=4.31e15*exp(-131.7/RT)
kbb=7.41e+7*exp(-32.7/RT)	
kp=2.21e+7*exp(-17.9/RT)	
kpt=1.20e+6*exp(-28.6/RT)
km=2.9e+5*exp(-32.6/RT)
kmt=2.0e+5*exp(-46.1/RT)
kmp=4.01e3*exp(-29.0/RT)
kms=63.0
kt=1.30e10*exp(-8.4/RT)
ktsm=4.2e9*exp(-6.6/RT)
ktmm=1.8e7*exp(-5.6/RT)


temperature = 60.0 °C

// BA polymerization
AIBN --kiAIBN--> I + I
I + BA --kp--> SPR(1) 

SPR(n) --kbb--> MCR(n)
MCR(n) + BA --kpt--> SPR(n,SCB)

MCRL(n) + BA --kpt--> SPR(n+1)
SPR(n) + BA --kp--> SPR(n+1) 
SPR(n) + BA --km--> D(n) + SPR(1)
MCR(n) + BA --kmt--> D(n) + SPR(1)

SPR(n) + D(m,BA) --kmp--> D(n) + MCRL(m,BAX)
SPR(n) + LM --kms--> D(n) + I

SPR(n) + SPR(m) --kt--> D(n+m)
SPR(n) + MCR(m) --ktsm--> D(n+m)
SPR(n) + MCRL(m) --ktsm--> D(n+m)
MCR(n) + MCR(m) --ktmm--> D(n+m)
MCR(n) + MCRL(m) --ktmm--> D(n+m)
MCRL(n) + MCRL(m) --ktmm--> D(n+m)
Plugin TerminationChainLengthDependent { coefficient = kt, ic = 30, as = 0.85, al = 0.22}
Plugin TerminationChainLengthDependent { coefficient = ktsm, ic = 30, as = 0.85, al = 0.22}