kill(all);

kw = 1e-14$
eq1 :Ca + kw/h - h + C*ka/(h+ka)$
eq2:diff(eq1,h,1);

eq3:-log(10)*h$
eq : eq2*eq3;

print(" buffer index 'beta' = ",%)$

beta_index :log(10)*(kw/h+h+(C*h*ka)/(ka+h)^2);