No Title  |POLVER05_0 |1
#Example 17-6b: Maximum Mixedness Model with Bimodal distribution (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 = 6-z #
k2 = 1 #
k1 = 1 #
k3 = 1 #
rc = k1*ca*cb #
re = k3*cb*cd #
E1 = 0.47219*lam^4-1.30733*lam^3+0.31723*lam^2+0.85688*lam+0.20909 #
E2 = 3.83999*lam^6-58.16185*lam^5+366.2097*lam^4-1224.66963*lam^3+2289.84857*lam^2-2265.62125*lam+925.46463 #
E3 = 0.00410*lam^4-0.07593*lam^3+0.52276*lam^2-1.59457*lam+1.84445 #
rb = -k1*ca*cb-k3*cb*cd #
ra = -k1*ca*cb-k2*ca #
rd = k2*ca-k3*cb*cd #
E = if(lam<=1.82)then(E1)else(if(lam<=2.8)then(E2)else(E3)) #
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)=6
X=(cao-ca)/cao
Scd=cc/(cd+0.00001)
Sde=cd/(ce+0.00001)