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