Example CD13-1: Segregation Model with Asymmetric RTD |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  #
x =  1-ca/ca0  #
rb =  k1*ca-k2*cb  #
E1 =  -2.104*t^4+4.167*t^3-1.596*t^2+0.353*t-0.004  #
E2 =  -2.104*t^4+17.037*t^3-50.247*t^2+62.964*t-27.402  #
E =  if(t<=1.26)then(E1)else(E2)  #
t(0)=0
ca(0)=1
cb(0)=0
cc(0)=0
cabar(0)=0
cbbar(0)=0
ccbar(0)=0
t(f)=2.52