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

                                     5
Задание 2.
Решение задачи Коши с апостериорной оценкой глобальной погрешности,
максимальное допустимое значение которой задано.
Назначение.
Интегрирование обыкновенного дифференциального уравнения

                        y′ = f(x,y),     x∈[А,В]
с начальным условием
                                y(c)=yc,
где точка c совпадает либо с началом, либо с концом отрезка интегрирования.
Описание параметров.
 data – имя файла исходных данных;
    f – имя процедуры – функции с двумя параметрами, которая должна
         быть описана в программе (функция f – вычисляет значение правой
         части уравнения (1));
  rez – имя файла выходных данных;
Icod – код завершения.
Замечание о структуре файла исходных данных.
      1. Первая строка – значения A, B, C, yc;
      2. Вторая строка – начальное значение H шага интегрирования;
         максимально допустимое значение ε абсолютной погрешности в
         конечной точке интегрирования.
Замечание о структуре выходного файла.
      1. Первая строка – вычисленное значение εr абсолютной погрешности в
         конечной точке интегрирования; шаг интегрирования, с которым
         получена погрешность εr; значение Icod – индикатор ошибки,
         принимающий следующие значения:
           • Icod = 0 – завершение в соответствии с назначением (εr ≤ε);
           • Icod = 1 – процесс решения прекращен, т.к. с уменьшением шага
              погрешность не уменьшается;
           • Icod = 2 – процесс решения прекращен, т.к. значение шага стало
              недопустимо малым;
           • Icod = 3 – процесс решения прекращен, т.к. дальнейшее
              применение метода невозможно (в случае, когда реализуется
              расчетная схема 2 или 3);
           • Icod = 4 – решение не получено, двухсторонний метод Рунге –
              Кутта с данным начальным шагом неприменим (в случае, когда
              реализуются расчетные схемы 2 или 3).
      2. Вторая строка – X–координата конца отрезка интегрирования;
         полученное в конце отрезка интегрирования значение решения.
          Если Icod = 4, то выходной файл во второй строке содержит только
          X – координату конца отрезка интегрирования.