Решение систем алгебраических и трансцендентных уравнений, страница 20

12. Формирование накопителя корней сеансов решения NKSRnp1, присоединяющего корни текущего сеанса решения к корням, вычисленным в предыдущих сеансах (рис. 6 .8, Е).

13. Ввод данных сеанса решения в протокол решения системы ПРСУ (табл. 6.1) и обновление эталонного вектора Rxet0 по выходной таблице Rxet этапа 11.

На рис. 6. 7, Ж, З, И,К показано также вычисление невязки для найденных в сеансе решения корней системы уравнений. Это вычисление может понадобиться для оценки точности вычисления корней и оценки представления корней с определенным количеством знаков в окончательных расчетах (величины компонент невязки должны быть меньше TOL) или для выяснения того, что найденные величины действительно являются корнями системы (при отладке решателя или при испытании точек из графического решения системы). Невязка может вычисляться для одного решения (используются по одному корню для каждой неизвестной) и для группы решений (для каждой неизвестной используется вектор корней, формируемый по выходной таблице решателя системы уравнений).

Для одного решения неважно представление выражений невязок: в обычном или векторизованном виде. Для группы решений выражения невязки следует векторизовать: отсутствие векторизации может привести по результатам расчета к ложному выводу, что исследуемые решения не содержат корней системы уравнений (см. вычисление невязок на рис. 6.7, И и К) или полученные неправильно данные – корни системы уравнений, хотя в действительности они таковыми не являются (см. рис. 6.7, З и Ж).

Вычисление невязки для группы решений системы нелинейных уравнений  в виде выходных таблиц корней решателя (рис. 6.7, К).

·  Формирование векторов корней неизвестных  по выходным таблицам корней решателя (например, вектора х1 по таблице rs3all(x1,y1) из рис. 6.7, Г).

– Щелчком ЛКМ на номере столбца таблицы rs3all(x1,y1)0 выделите все записи в ней, и командой Copy Selection контекстного меню этой таблицы запишите их в буфер.

– Введите присваивание для х1, и в правое место ввода вставьте содержимое буфера.

– Повторите первые два пункта для остальных неизвестных системы.

·  Ввод выражений компонент невязки системы: в правые части присваиваний компонент копируются функции уравнений системы из действующего решателя, а затем – векторизуются аналогично этапу 3.