Здравствуйте необходимо решить систему из 7 нелинейных уравнения следующего вида
1)

где

2)

где


3)

Где

являются искомыми, все остальные переменные константы
Всего уравнений 7 из них 4 уравнения вида 1), два уравнения вида 2), и одно вида 3).
Пробовал решать систему методом Градиента, а уравнения интегрировать методом Чебышева.
При достаточно высокой сходимости порядка

, невязка составляла все лишь

Такая низкая невязка, как мне кажется, из-за того что Градиентный метод находит лишь локальные минимумы.
Подскажите пожалуйста какой-либо метод подходящей для решения моей системы .