Example 17-8 Using Software to Make Maximum Mixedness Model Calculations |POLVER05_0 d(x)/d(z) = -(ra/cao+E/(1-F)*x) # cao = 8 # k = .01 # lam = 200-z # ca = cao*(1-x) # E1 = 4.44658e-10*lam^4-1.1802e-7*lam^3+1.35358e-5*lam^2-.000865652*lam+.028004 # E2 = -2.64e-9*lam^3+1.3618e-6*lam^2-.00024069*lam+.015011 # F1 = 4.44658e-10/5*lam^5-1.1802e-7/4*lam^4+1.35358e-5/3*lam^3-.000865652/2*lam^2+.028004*lam # F2 = -(-9.30769e-8*lam^3+5.02846e-5*lam^2-.00941*lam+.618231-1) # ra = -k*ca^2 # E = if (lam<=70) then (E1) else (E2) # F = if (lam<=70) then (F1) else (F2) # EF = E/(1-F) # z(0)=0 x(0)=0 z(f)=200