A powerful solver based on succesive over-relaxation for the one dimensional modified Poisson–Boltzmann equation, which includes a dielectric response profile, unequal ionic valencies, surface interactions via potentials of mean force and the presence of impurities. The modified PBE is solved in the form of:
You'll need Python3 with the numba package installed.
python3 setup.py install