Example CD13-1: Maximum Mixedness Model with Bimodal Distribution |POLVER05_0
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 =  6-z  #
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.4646  #
EF =  E/(1-F)  #
E =  if(lam<=1.26)then(E1)else(if(lam<=2.8)then(E2)else(E3))  #
E3 =  0.00410*lam^4-0.07593*lam^3+0.52276*lam^2-1.59457*lam+1.84445  #
z(0)=0
ca(0)=1
cb(0)=0
cc(0)=0
F(0)=.99
z(f)=6