** Binomial Number mechanics
** from inp_polymer_158.txt


*Time, 
** start, end, console out interval, plot out interval, csv interval, unit
0, 100, 20, 20, 20, seconds
*End


*Element
** 0:element name, 1:initial number, 2:plot color, 3:marker
 ii, 1e18, Orange,  
 iR,  0, Red,  s, 
 Cat, 1e20, Gray, 1, 
 M, 1e23, Green, 4,
*End


*ElementPolymer
** 0: element name, 1: first number, 2:last number, 3:interval, 
** 4:initial Number, 5:type
** ラジカル重合のときは、「M=*」を定義しないこと！
M=*, 1, 120, 1, 0, M=*
M=*_R, 1, 100, 1, 0, M=*_R
*End


*Reaction, 1e23
r1_1, 100, 1, ii, 0.15, 2, iR 
r3_2, 100, 1, M, 1, Cat, 1, iR, 100, 1, M=1_R, 1, Cat  
*End


*ReactionPolymer, 1e23
rP, 200, 1, M=*_R, 1, M, 0.1, 1, M=*_R
rP, 300, 1, M=*_R, 1, M=*_R, 1000, 1, M=*
*End

*Plot, log
ii, iR, 
M,
M=10_R, M=20_R, M=30_R, 
M=10, M=20, M=30, 
M=40, M=50, M=60, 
*End


*comments
2021.5.9
2021.12.4
2022.3.9
2022.4.18


