Название: Вычислительная математика (В.Т. Кононов,Г.П. Чикильдин)

Жанр: Педагогика

Просмотров: 1574

2.4. описание подпрограмм

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

N1YDUA  (N, NK, A, B, YO, DT, KON, X, Y)

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

Входные параметры подпрограммы:

N – порядок дифференциального уравнения;

NK – количество коэффициентов в левой части дифференциального уравнения, NK=N+1;

A(I) – массив коэффициентов   в левой части дифференциального уравнения, I = 1, NK;

B – коэффициент  в правой части дифференциального уравнения;

YO(I) – массив начальных условий   

I = 1,NK;

DT – шаг дискретизации  решения дифференциального уравнения;

KON – количество дискретных значений решения на интервале ;

X(K) – массив реализации функции  в правой части дифференциального уравнения, K=1,KON.

Выходные параметры подпрограммы:

Y(I,K) – массив решения  и всех его производных на интервале решения, I=1,NK, K=1,KON.

В подпрограмме N1YDUA осуществляется обращение к подпрограммам N1YRKC и N1YADC, реализующим алгоритмы Рунге-Кутта четвертого порядка точности и Адамса четвертого порядка точности с прогнозом и коррекцией решения, соответственно.

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

N1YRKC  (NM, NK, DT, A, B, X1, X2, Y)

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

Входные параметры подпрограммы:

NM – порядок  дифференциального уравнения;

NK, DT, A(I), B – см. описание подпрограммы N1YDUA;

X1 – значение функции  в -й момент времени;

X2 – значение функции  в -й момент времени;

Y(J) – массив решения  и всех его производных в -й момент времени J=1,NK.

Выходные параметры подпрограммы:

Y(J) – массив решения и всех его производных в -й момент времени, J=1,NK.

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

N1YADC  (NM, NK, DT, A, B, X, Y, YY, YYY)

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

Входные параметры подпрограммы:

NM, NK, DT, A, B – см. описание подпрограмм N1YDUA, N1YRKC;

X – значение функции  в -й момент времени;

Y(J) – массив решения  и всех его производных в -й момент времени J=1,NK;

YY(J) – массив решения  и всех его производных в -й момент времени J=1,NK;

YYY(J) – массив решения  и всех его производных в -й момент времени J=1,NK.

Выходные параметры подпрограммы:

Y(J) – массив решения  и всех его производных в -й момент времени, J=1,NK.