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

                                     27
с начальным условием
                                 y(c)=yc,
где точка c совпадает либо с началом, либо с концом отрезка интегрирования.
Описание параметров.
 data – имя файла исходных данных;
    f – имя процедуры – функции с двумя параметрами, которая должна
         быть описана в программе (функция f – вычисляет значение правой
         части уравнения (1));
  rez - имя файла выходных данных;
Icod - выходная переменная – код завершения.
Замечание о структуре файла исходных данных.
      Первая строка – значения A, B, C, yc.
      Вторая строка – начальное значение Н шага интегрирования;
      m – число верных знаков решения; hmin – наименьший допустимый шаг
      интегрирования.
Замечание о структуре выходного файла.
      Первая и последующие строки - х–координата точки интегрирования;
      полученное приближенное решение в этой точке, число верных знаков
      приближенного решения.
      Последняя строка – значение Icod – индикатор ошибки, принимающий
      следующие значения:
           Icod=0 –завершение в соответствии с назначением
                    (приближенное решение с заданным числом
                     верных знаков получено);
           Icod=L – в L точках требуемая точность не достигается;
           Icod=1 – ошибка входных данных;
      значение реального наименьшего шага интегрирования;
      значение наибольшего шага интегрирования.
Метод. (См. раздел 4.4. Вложенные методы локальной погрешности;
        Раздел 5.1. Метод удвоения и деления шага пополам; раздел 4.5.
        Мера погрешности приближенного решения).
   1. Для получения приближенных значений решения и оценки их
      локальных погрешностей используются вложенные методы, конкретный
      вид которых зависит от номера Вашего варианта.
   2. Для достижения заданной точности решения (для обеспечения m
      верных знаков в приближенных значениях решения) шаг в каждой точке
      интегрирования выбирается методом удвоения и деления шага пополам.
      При этом метод удвоения и деления шага пополам должен быть
      реализован с учетом Замечания 1 и Замечания 5 из раздела 5.1. Кроме
      того, необходимо выполнить следующее требование к алгоритму выбора
      шага. Если текущее значение шага, большее hmin., не обеспечивает
      требуемую точность, а при делении текущего шага пополам получается