No Title  |POLVER05_0 |1
k = 0.1 
d(X) / d(z) = -(ra/Ca0+E/(1-F)*X) 
X(0) = 0
lam=30-z 
Ca = Ca0*(1-X) 
E1=0.01*lam - 0.1 
E2=0.3 - 0.01*lam
F1= 0.005*lam^2 - 0.1*lam + 0.5
F2=0.3*lam - 0.005*lam^2 - 3.5 
E=if(lam<=10) then (0) else if (lam>=10 and lam <=20) then (E1) else (E2)
F=if(lam<=10) then (0) else if (lam>=10 and lam <=20) then (F1) else (F2)
EF=E/(1-F)
Ca0=8
ra=if(lam<=10) then (0) else if (lam>=10 and lam <=14) then (-k*Ca^0.6) else (-k*Ca^2.7)
z(0) = 0.01
z(f) = 30