Численное решение задачи Коши для обыкновенных дифференциальных уравнений методами типа Рунге-Кутта. Часть 2. Корзунина В.В - 13 стр.

                                                13
        коэффициентов,     которые      после интегрирования всех r
        задач Коши дадут коэффициенты интерполяционных полиномов
        откорректированного решения.
Варианты Задания 6
                      Вариант                        Порядок точности r
                                  1                           4
                                  2                           5
                                  3                           6
                                  4                           7

Задание 7
Решение задачи Коши для системы обыкновенных дифференциальных
уравнений первого порядка с апостериорной оценкой погрешности в
контрольных точках.
Назначение.
Интегрирование системы обыкновенных дифференциальных уравнений
y′ =f (x, y),x∈[A, B] с начальным условием y(c) =yc (точка                  C совпадает либо с

началом, либо с концом                отрезка интегрирования), где          y =(y1, y2,...,yM)T ,
yc =(yc1, yc2,...,ycM )T , f =( f 1(x, y1,...,yM ),...,f M (x, y1,...,yM ))T , M   - размерность
системы уравнений.
Описание параметров.
 data - имя файла исходных данных;
    f -    имя            подпрограммы,      вычисляющей      значения        функций
         f i (x, y1,...,yM ), i =1, M . Список её параметров: (M, x,Y1, R1) , где М -
        описан выше;
   х - входная переменная х, определяющая значение аргумента х;
  Y1 - входной массив размера M , содержащий значения аргументов
           y1 , y 2 ,..., y M ;
 R1 -                                          M,
          выходной массив размера                      содержащий значения             функций
        f i ( x, y1 ,..., y M ) , i =1, M
 rez - имя файла выходных данных;
Icod - выходная переменная – код завершения подпрограммы,
       принимающей следующие значения:
                                   Icod=0 – нет ошибки, решение получено;
                                   Icod=4 – в четвертой контрольной точке решение не
                                   получено из-за невозможности дальнейшего
                                   применения двухстороннего метода Рунге-Кутта;
                                   Icod=1 – ошибка входных данных.