Для вычисления распределения электрического потенциала использовано его представление через бинарную функцию распределения, позволившее выделить коррелированную и некоррелированную его составляющие. Некоррелированная составляющая определена через напряженность электрического поля, которая, в свою очередь, была выражена через интеграл от распределения плотности заряда по пространственной переменной. Это позволило представить некоррелированную часть энергии взаимодействия зарядов через интеграл от квадрата градиента соответствующей части потенциала. Для нахождения коррелированной составляющей использована аппроксимация бинарной функции распределения, широко применяемая для расчетов в рамках среднесферического приближения, а также модифицированное выражение для радиуса Дебая. В Фурье-пространстве получено линейное соотношение между вариациями плотности распределения заряда и электрического потенциала с коэффициентом, зависящим от волнового вектора. В длинноволновом приближении это соотношение переходит в уравнение Пуассона. При включении внешнего поля сформулировано выражение для электрохимического потенциала, приводящее к уравнению для распределения плотности заряда по объему системы.