Example CD13-1: Segregation Model with Bimodal Distribution |POLVER05_0
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   #
ca0 =  1  #
k1 =  1  #
k2 =  1  #
tau =  1.26  #
ra =  -k1*ca  #
rc =  k2*cb  #
rb =  k1*ca-k2*cb  #
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.20970*t^4-1224.66963*t^3+2289.84857*t^2-2265.62125*t-925.46463  #
E =  if(t<=1.82)then(E1)else(if(t<=2.8)then(E2)else(if(t<6)then(E3)else(E4)))  #
E3 =  0.00410*t^4-0.07593*t^3+0.52276*t^2-1.59457*t+1.84445  #
E4 =  0  #
t(0)=0
ca(0)=1
cb(0)=0
cc(0)=0
cabar(0)=0
cbbar(0)=0
ccbar(0)=0
t(f)=7