No Title  |POLVER05_0 |1
#LEP-WEB17-1g: Segregation Model with Bimodal Distribution (Multiple Reactions)
d(ca)/d(t) =  ra  #
d(cb)/d(t) =  rb  #
d(cc)/d(t) =  rc  #
d(cabar)/d(t) =  ca*E  #
d(cbbar)/d(t) =  cb*E  #
d(ccbar)/d(t) =  cc*E  #
d(cd)/d(t) =  rd  #
d(ce)/d(t) =  re  #
d(cdbar)/d(t) =  cd*E  #
d(cebar)/d(t) =  ce*E  #
k1 =  1  #
k2 =  1  #
k3 =  1  #
E1 =  0.47219*t^4-1.30733*t^3+0.31723*t^2+0.85688*t+0.20909  #
E2 =  3.83999*t^6-58.16185*t^5+366.2097*t^4-1224.66963*t^3+2289.84857*t^2-2265.62125*t+925.46463  #
E3 =  0.00410*t^4-0.07593*t^3+0.52276*t^2-1.59457*t+1.84445  #
rc =  k1*ca*cb  #
re =  k3*cb*cd  #
ra =  -k1*ca*cb-k2*ca  #
rb =  -k1*ca*cb-k3*cb*cd  #
rd =  k2*ca-k3*cb*cd  #
E =  if(t<=1.82)then(E1)else(if(t<=2.8)then(E2)else(E3))  #
t(0)=0
ca(0)=1
cb(0)=1
cc(0)=0
cabar(0)=0
cbbar(0)=0
ccbar(0)=0
cd(0)=0
ce(0)=0
cdbar(0)=0
cebar(0)=0
cao=1
t(f)=6
X=(cao-ca)/cao
d(Xbar) / d(t) = X*E 
Xbar(0) = 0
Scd=ccbar/(cdbar+0.00001)
Sde=cdbar/(cebar+0.00001)