function [pH] = pHfunc_SecondB(H_plus)

pH = -log10(abs(H_plus));

end