Термодинамическое равновесие и его численное моделирование. Петрусев А.С. - 21 стр.

  Решение
   задачи
 линейного
программи-
  рования
     ↓
   Метод
наименьших
 квадратов
    ↓
 Начальное
приближение
       0
   λi=λi
    ↓
   ...

Примерно такую структуру имеет, например, блок расчёта равновесия
программы «ИВТАНТЕРМО» и ряд зарубежных программ (см. [2]). В
случае расчёта систем с неизвестной температурой применение
процедуры линейного программирования не очень удобно. Кроме того,
данная процедура, как отмечалось, не решает полностью проблему
расчёта малых концентраций. Но отказ от использования процедуры
линейного программирования приводит к необходимости мириться с
плохим начальным приближением и медленной сходимостью итераций.
     Улучшить расчёт малых концентраций можно путём введения в
алгоритм процедуры квазилинейного исключения:

 Начальное
приближение
   λi=λ0
     ↓←←←←←←←←←←←←←←←←←←←←←←
     ↓                            ↑       ↑
Вычисление                    Проверка    ↑
  матрицы                  совместности   ↑
   Якоби                    вырожденных   ↑
  J=DQ/Dλ                    уравнений    ↑
     ↓                            ↑       ↑
Алгоритм Гаусса            Преобразование     ↑
                вырождение
   с проверкой  →→→→→→          системы       ↑
   вырождения    Якобиана      уравнений      ↑
    ↓                                    ↑
  Процедура                                   ↑
демпфирования                                 ↑
    ↓                                    ↑