Кто-нибудь проводил расчет газожидкостного равновесия (ГЖР) для газоконденсатной смеси по уравнению состояния и компонентному составу? Там надо вычислять константы равновесия в итерационном цикле, при этом надо на каждой итерации решать итерационным способом уравнение Рашфорда-Райса или Рэкфорда-Райса (Rachford and Rice) для нахождения молярной доли газовой фазы. В книгах пишут что нужно применять итерационный метод Ньютона или Ньютона-Рафсона. А на деле у меня этот метод иногда расходится. Применял также метод деления пополам (бисекции) - то же самое. Просто на каком-нибудь шаге итерационного цикла получаются такие значения констант равновесия что уравнение Рашфорда-Райса не имеет корней и асимптотически приближается к нулю и имеет точку разрыва. Посоветуйте пожалуйста как правильно построить алгоритм.