Numerical solution of a modified poisson-boltzmann equation in electrolyte solution theory