m=m+1; it=it+1;

xout=[xout;x1'];

dxout=[dxout;dx'];

mout=[mout;m];

end;

disp('Количество итераций равно'); disp(it);

2.4 Используемые технические средства

ЭВМ, совместимые с IBM PC. Процессор не ниже Pentium1. ОП не меньше 32 Мб, мышь, стандартная клавиатура, видеокарта.

2.5 Вызов и загрузка

Для запуска программы из среды MatLab необходимо вызвать головную программу, набрав в командной строке ее имя main. При этом необходимо перейти в тот каталог, в котором находится данный пакет программ. Затем производятся все предусмотренные программой операции. После вывода на экран каждого графика или сообщения система ожидает нажатия любой клавиши.

Программа со всеми сопутствующими файлами, занимает 4,3 Кб.

2.6 Входные данные

Правые части системы в файле funf.m, матрица Якоби в файле dif.m, система дифференциальных уравнений, составленная по исходной системе нелинейных САУ (в файле dif.m), начальное приближение, значение шага, интервал времени для интегрирования, допустимая ошибка.

2.7 Выходные данные

На выходе программы на экран выводятся приближенные решения, полученные по методам Рунге-Кутта, значения шага и времени на каждом шаге интегрирования, количество шагов интегрирования, графики приближенных решений по методам Рунге-Кутта, количество итераций, необходимых для нахождения уточненного решения дискретным методом Ньютона, уточненное решение при каждой итерации, графики уточненных значений по методу Ньютона и ошибки для каждой составляющей решения.

3. ОПИСАНИЕ ТЕСТОВЫХ ЗАДАЧ

Решение системы нелинейных САУ.

Для интегрирования возьмем систему:

x² +y² -4=0

xy – 1=0

Тогда при запуске программы на экране появляются следующие сообщения:

Метод Рунге - Кутта 1го порядка

t =

0

h =

0.1000

y =

2 0

…

t =