No Title  |POLVER05_0 |1
#Example LEP-WEB17-1e: Maximum Mixedness Model with Asymmetric RTD
d(ca)/d(z) =  -(-ra+(ca-cao)*EF)  #
d(cb)/d(z) =  -(-rb+(cb-cbo)*EF)  #
d(cc)/d(z) =  -(-rc+(cc-cco)*EF)  #
d(F)/d(z) =  -E  #
cao =  1  #
k1 =   1  #
k2 =   1  #
tau =   1.26  #
ra =   -k1*ca  #
rc =   k2*cb  #
rb =   k1*ca-k2*cb  #
cbo =  0  #
cco =  0  #
lam =  2.52-z  #
E1 =  -2.104*lam^4+4.167*lam^3-1.596*lam^2+0.353*lam-0.004  #
E2 =  -2.104*lam^4+17.037*lam^3-50.247*lam^2+62.964*lam-27.402  #
EF =  E/(1-F)  #
E =  if(lam<=1.26)then(E1)else(E2)  #
z(0)=0
ca(0)=1
cb(0)=0
cc(0)=0
F(0)=.99
z(f)=2.52
X=1-ca/cao