No Title  |POLVER05_0 |1
#LEP-WEB17-1h: Maximum Mixedness Model with Asymmetric RTD (Multiple Reactions) 
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 #
d(cd)/d(z) = -(-rd+(cd-cdo)*EF) #
d(ce)/d(z) = -(-re+(ce-ceo)*EF) #
cbo = 1 #
cao = 1 #
cco = 0 #
cdo = 0 #
ceo = 0 #
lam = 2.52-z #
k2 = 1 #
k1 = 1 #
k3 = 1 #
rc = k1*ca*cb #
re = k3*cb*cd #
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  #
rb = -k1*ca*cb-k3*cb*cd #
ra = -k1*ca*cb-k2*ca #
rd = k2*ca-k3*cb*cd #
E = if(lam<=1.26)then(E1)else(E2)#
EF = E/(1-F) #
z(0)=0
ca(0)=1
cb(0)=1
cc(0)=0
F(0)=0.99
cd(0)=0
ce(0)=0
z(f)=2.52
X=1-ca/cao
Scd=cc/(cd+0.000001)
Sde=cd/(ce+0.000001)