restart; Digits := 30: with (plots): with(RootFinding):
#Define given parameters
R:=1: #external radius of particles, cm
d:=10^(3): #diffusivity cm^2 per second
alpha:= 1: #fractional void volume
c0:=10^(6): #concentartion of soltion in void volume of solid initially, moles per liter
C0:=0: #concentration of main body of solution initially, moles per liter
k1:=0.5: #constant in adsorption isotherm (ka)
k2:=0.75: #constant in adsorption isotherm (kd)
k:=2.5: #equilbrium constant for adsorption kinetics
n0:=(k1/k2)*c0:#initial amount absrobed on solid, moles per liter
V:=0.1: #volume of external solution, liters
W:=0.1: #weight of absorbant, grams
rho:=2.0: #solid aparrant density, g/cc
delta:=(1/d)*((p+alpha*k2+k1)/(p+alpha*k2));
beta:=W*alpha*d/(rho*V);
