Текст подпрограммы и версий ( Фортран )
de11r.zip , de11d.zip 		Тексты тестовых примеров ( Фортран )
tde11r.zip , tde11d.zip
Текст подпрограммы и версий ( Си )
de11r_c.zip , de11d_c.zip 		Тексты тестовых примеров ( Си )
tde11r_c.zip , tde11d_c.zip
Текст подпрограммы и версий ( Паскаль )
de11r_p.zip , de11e_p.zip 		Тексты тестовых примеров ( Паскаль )
tde11r_p.zip , tde11e_p.zip

Подпрограмма:  DE11R

Назначение

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

Математическое описание

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

          Y ' = F (X, Y) ,
          Y = ( y1, ... , yM ) ,
          F = ( f1 (X, y1, ... , yM), ... , fM (X, y1, ... , yM) )
 с начальными условиями, заданными в точке XN :
          Y(XN) = YN ,     YN = ( y10, ... , yM0 ) ,
 методом Хойна.

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

И.С.Березин, Н.П.Жидков, Методы вычислений, т.2, Физматгиз, 1960.

Использование 

    SUBROUTINE  DE11R (F, M, XN, YN, XK, HMIN, EPS, P, H, 
                                            Y, YP, DELTY, RAB, IERR) 
   Параметры
F - имя подпрограммы вычисления значений правой части дифференциального уравнения. Первый оператоp подпрограммы должен иметь вид:
SUBROUTINE  F (X, Y , DY, M).
Здесь: X, Y - значения независимой и зависимой переменных, соответственно. Вычисленное значение правой части должно быть помещено в DY. B случае системы уравнений, т.е. когда M ≠ 1 , параметры Y и DY представляют массивы длины M (тип параметров X, Y и DY: вещественный);
M - количество уравнений в системе (тип: целый);
XN, YN - начальные значения аргумента и решения. B случае системы уравнений (т.е. M ≠ 1) YN представляет одномерный массив длины M (тип: вещественный);
XK - значение аргумента, при котоpом требуется вычислить решение задачи Коши (конец интервала интегрирования). XK может быть больше, меньше или pавно XN (тип: вещественный);
HMIN - минимальное значение абсолютной величины шага, который разрешается использовать при интегрировании данной системы уравнений (тип: вещественный);
EPS - допустимая меpа погрешности, с которой тpебуется вычислить все компоненты решения (тип: вещественный);
P - граница перехода, используемая при оценке погрешности решения (тип: вещественный);
H - вещественная переменная, содержащая начальное значение шага интегрирования. Может задаваться с учетом направления интегрирования, т.е. положительным, если XK > XN, отрицательным, если XK < XN, или без такого учета в виде абсолютной величины;
Y - искомое решение задачи Коши, вычисленное подпрограммой при значении аргумента XK. Для системы уравнений (когда M ≠ 1) задается одномерным массивом длины M. B случае совпадения значений параметров XN и XK значение Y полагается равным начальному значению YN (тип: вещественный);
            YP -
       DELTY   		вещественные одномерные рабочие массивы длины M;
RAB - вещественный двумерный рабочий массив размеpа M*2;
IERR - целая переменная, значение которой в результате работы подпрограммы полагается равным 65, если какая-нибудь компонента решения не может быть вычислена с требуемой точностью EPS. B этом случае интегрирование системы прекращается. При желании интегрирование системы можно повторить обращением к подпрограмме с новыми значениями параметров HMIN и H. 
   Версии
DE11D - вычисление решения задачи Коши для системы обыкновенных дифференциальных уравнений первого порядка в конце интервала интегрирования методом Хойна с повышенной точностью. При этом параметры XN, YN, XK, HMIN, EPS, P, H, Y, YP, DELTY, RAB и параметры X, Y и DY в подпрограмме F должны иметь тип DOUBLE PRECISION. 
   Вызываемые подпрограммы
UTDE10 - подпрограмма выдачи диагностических сообщений при работе подпрограммы DE11R.
UTDE11 - подпрограмма выдачи диагностических сообщений при работе подпрограммы DE11D. 
   Замечания по использованию
 

Подпрограмма DE11R предназначена для численного решения дифференциальных уравнений и систем уравнений с правой частью, имеющей непрерывные частные производные вплоть до 4 порядка включительно.

Хотя заданная точность EPS не гарантируется в общем случае, большой опыт эксплуатации данной подпрограммы убедительно показывает, что вычисляемое ею численное решение достаточно близко приближает точное решение.

При работе подпрограммы значения параметров M, XN, YN, XK, HMIN, EPS, P сохраняются. Если после работы подпрограммы нет необходимости иметь начальное значение решения YN, то параметры YN и Y при обращении к ней можно совместить.

Пример использования

Использование подпрограммы иллюстрируется на примере: 

          y '1  =  y2 ,
          y '2  =  - y1 ,      3π/4 ≤ x ≤ π

          y1(3π/4)  =  √2 / 2 ,      y2(3π/4)  =  - √2 / 2

Приводятся подпрограмма вычисления значений правой части и фрагмент вызывающей программы, а также результаты счета. 

          SUBROUTINE  F (X, Y, DY, M)
          DIMENSION  Y(2), DY(2)
          DY(1) = Y(2)
          DY(2) = -Y(1)
          RETURN
          END
      
          DIMENSION  Y(2), YP(2), DELTY(2), RAB(2, 2)
          EXTERNAL  F
          M = 2
          XN = 0.75*3.14159265359
          Y(1) = SQRT(2.)/2.
          Y(2) = -Y(1)
          XK = 3.14159265359
          EPS = 3.E-5
          P = 1.E-6
          HMIN = 1.E-8
          H = 0.01
          CALL DE11R (F, M, XN, Y, XK, HMIN, EPS, P, H, Y, YP, 
        *                         DELTY, RAB, IERR)

Результаты:

          Y (1)  =  -0.5928428474 * 10-7
          Y (2)  =  -0.9999983798

          IERR  =   0