Выбор начального приближения для решения СНАУ

[firefly2038]

Привет всем,

Я решаю задачу по расчету параметров Чемпена-Жуге при моделировании детонации.
Система состоит из трех уравнений. Решить ее нужно методом Ньютона и Бройдена.

Соответственно, эти методы требуют хорошее начальное приближения, а для 3х неизвестных я не знаю как это можно сделать. Подскажите плз как мне найти это приближение.

Система такого вида:
$\frac{p \cdot \eta}{\gamma(T) - 1} - \frac{p_0 \cdot \eta_0}{\gamma_0(T_0) - 1} = 0.5\cdot (p_0 + p) \cdot (\eta_0 \eta) + \operatorname{const} $
$\frac{p - p_0}{\eta_0 - \eta} = \frac{\gamma(T) \cdot p}{\eta}$
$p = \frac{R \cdot T}{\operatorname{const}}$

Функция $\gamma(T)$ - полином 4ой степени.
Соответственно неизвестные это - $T$, $p$ и $rho = \frac{1}{\eta}$

Методы я вроде правильно написал на языке matlab. Для задачи с другими коэффициентами решение находится без проблем.

[firefly2038]

что не существует какого-то правила для отбора начального приближения в 3хмерном случае?