No Title  |POLVER05_0 |1
d(X) / d(W) = -raprime/Fa0 
Ca = Ca0*(1-X) 
k = kc*kr/ (kc+kr) 
raprime=acprime*ra
ra = -k*Ca 
Sc = nu/DAB 
Re = dp*U/nu
nu=mu/rho
phi=0.4
y=1
Reprime=Re/((1-phi)*y)
X(0) = 0
Shprime = (Reprime^(1/2))*(Sc^(1/3)) 
Ca0 = 1 
kr = 0.5 
dp = 10^-3 
U = 1 
mu = 10^-3 
DAB = 10^-8 
rho = 1000 
kc = Shprime*DAB*(1-phi)*y/(dp*phi) 
Fa0=Ca0*U*Ac
Ac=5*10^-2
W(0) = 0
W(f) = 100
WA=-ra
acprime=6*(1-phi)/(dp*rhoc)
rhoc=3000