|
Текст подпрограммы и версий (Паскаль) de44e_p.zip |
Тексты тестовых примеров (Паскаль) tde44e1_p.zip , tde44e2_p.zip , tde44e3_p.zip , tde44e4_p.zip |
Выполнение одного шага приближенного интегрирования системы обыкновенных дифференциальных уравнений первого порядка методом рядов Чебышёва.
Решается задача Коши для системы M обыкновенных дифференциальных уравнений первого порядка
(1) Y ' = F (x, Y) ,
Y = ( y1, ... , yM ) , F = ( f1 (x, y1, ... , yM), ... , fM (x, y1, ... , yM) ) , XN ≤ x ≤ XK
с начальными условиями, заданными в точке XN :
(2) Y(XN) = YN , YN = ( y10, ... , yM0 ) ,
при условии, что правая часть системы (1) имеет непрерывные ограниченные частные производные по переменным x, Y. Предполагается, что на отрезке [XN, XK] задача (1), (2) имеет единственное решение. Тогда решение задачи Коши и его производная Y'(XN + αΔ) = F(XN + αΔ, Y(XN + αΔ)) = Φ(α) разлагаются на промежутке интегрирования [XN, XK] в равномерно сходящиеся ряды по смещенным многочленам Чебышёва первого рода
∞ 1
(3) Y(XN + αΔ) = ∑ ' ai*[Y] Ti* (α), 0 ≤ α ≤ 1, ai*[Y] = 2/π ∫ Y(XN + αΔ) Ti* (α) / √ (α (1 - α) ) dα ,
i=0 0
∞ 1
(4) Φ(α) = ∑ ' ai*[Φ] Ti* (α), 0 ≤ α ≤ 1, ai*[Φ] = 2/π ∫ Φ(α) Ti* (α) / √ (α (1 - α) ) dα .
i=0 0
Здесь: штрих у знака суммы означает, что слагаемое с индексом 0 берется с дополнительным множителем 1/2; Ti* (α) - смещенный многочлен Чебышёва первого рода на [0, 1]: Ti* (α) = Ti (2α-1) ; Ti (t) - многочлен Чебышёва первого рода i-го порядка на [-1, 1]; Δ = XK - XN. Если ряд Чебышёва (3) (и ряд (4)) является быстросходящимся, то его сумма на [XN, XK] (и сумма ряда (4)) хорошо приближается частичной суммой некоторого порядка. Эта частичная сумма принимается в качестве приближенного аналитического решения задачи (1), (2) на промежутке [XN, XK]. В противном случае, т.е. при медленной сходимости ряда (3) на интервале [XN, XK], получение аналитического решения в виде одной частичной суммы на всем отрезке интегрирования [XN, XK] может быть затруднительным. Поэтому целесообразно использовать разбиение промежутка интегрирования [XN, XK] на такие элементарные сегменты некоторой длины H: [xs, xs+H], x0 = XN, s = 0, 1, ... , на каждом из которых ряды Чебышёва для решения Y(x) и его производной Y'(x) будут сходиться значительно быстрее. На каждом подобном сегменте решение исходной задачи Коши приближенно представляется в виде (K + 1) - й частичной суммы смещенного ряда Чебышёва
K+1 1
(5) Y (xs + αH) ≈ ∑ ' ai*[Y] Ti*(α) , 0 ≤ α ≤ 1 , ai*[Y] = 2/π ∫ Y(xs + αH) Ti* (α) / √ (α (1 - α) ) dα ,
i=0 0
а его производная - в виде частичной суммы K-го порядка
K 1
(6) Y' (xs + αH) ≈ ∑ ' ai*[Y ' ] Ti*(α) , 0 ≤ α ≤ 1 , ai*[Y ' ] = 2/π ∫ Y ' (xs + αH) Ti* (α) / √ (α (1 - α) ) dα .
i=0 0
В этом случае аналитическое решение задачи (1), (2) будет состоять из совокупности частичных сумм рядов Чебышёва, построенных на этих элементарных сегментах. Порядок частичных сумм и элементарные сегменты задаются пользователем при обращении к подпрограмме.
При обращении к подпрограмме DE44E задаются начало элементарного сегмента и его длина. По заданному значению решения в начале X = xs элементарного сегмента [xs, xs + H] подпрограмма DE44E вычисляет значение решения в конце элементарного сегмента, т.е. в узле X + H = xs + H. Одновременно вычисляются коэффициенты Чебышёва ai*[Y] (i = 0, 1, ... , K + 1) на элементарном сегменте [X, X + H] = [xs, xs + H] для решения задачи Коши Y(X + αH), 0 ≤ α ≤ 1, и коэффициенты Чебышёва ai*[Y'] (i = 0, 1, ... , K) его производной Y'(X + αH). Значение H может изменяться от сегмента к сегменту, в общем случае H = Hs.
При разбиении промежутка интегрирования на элементарные сегменты решение задачи на [XN, XK] сводится к определению нескольких наборов коэффициентов ai*[Y], i = 0, 1, ... , K + 1. Коэффициенты ai*[Y] ряда Чебышёва для решения на сегменте [X, X + H] = [xs, xs + H] выражаются через коэффициенты ai*[Φ] ряда Чебышёва его производной Φ(α) = F(xs + αH, Y(xs + αH)), 0 ≤ α ≤ 1, на [X, X + H] = [xs, xs + H], которые, в свою очередь, вычисляются приближенно итерационным способом, исходя из некоторого начального приближения. Вычисления выполняются с помощью квадратурной формулы Маркова на [X, X + H] с K + 1 узлом. При этом один из узлов квадратурной формулы совпадает с X, а остальные K узлов лежат внутри интервала (X, X + H). Количество итераций, которое предписывается выполнить в этом итерационном процессе, задается при обращении к подпрограмме и может меняться от сегмента к сегменту. Если при выбранном H ряды Чебышёва для решения Y(X + αH), 0 ≤ α ≤ 1 , и его производной Y'(X + αH) на элементарном сегменте [X, X + H] = [xs, xs + H] быстро сходятся, то для того, чтобы приближенное решение в конце такого сегмента имело максимальный порядок точности относительно H, необходимо выполнить не менее K итераций; при этом погрешность приближенного решения в конце элементарного сегмента является величиной порядка O(HK + 2) при H --> 0. Если H подобрано достаточно малым, то хорошая точность приближенного решения может быть получена и при меньшем числе итераций. Вообще, число итераций зависит от K и H. С увеличением H число итераций может также возрастать.
Начальное приближение коэффициентов ai*[Φ] ряда Чебышёва для производной на сегменте [X, X + H] = [xs, xs + H] может вычисляться двумя способами. В первом способе начальное приближение определяется только с использованием значения решения в точке X = xs. При этом погрешность начального приближения для всех коэффициентов a0*[Φ], a1*[Φ], ..., aK*[Φ] является величиной O(H2) при H --> 0. Во втором способе начальное приближение определяется через коэффициенты ряда Чебышёва производной решения Φ(α) на предыдущем элементарном сегменте [X - H', X] = [xs - 1, xs], где H' - длина предыдущего элементарного сегмента. В этом случае погрешности начального приближения для коэффициентов a0*[Φ], a1*[Φ], ... , aK*[Φ] имеют, соответственно, порядки O(H), O(H2), ..., O(HK + 1). Заметим, что длина H текущего элементарного сегмента [X, X + H] может быть больше или меньше длины H' предыдущего элементарного сегмента [X - H', X] или равна ей. Второй способ определения начального приближения в некоторых случаях может привести к более быстрой сходимости итерационного процесса и, тем самым, к меньшему числу выполняемых итераций. Второй способ может быть применен только начиная со второго элементарного сегмента [x0 + H0, x1 + H], H0 -длина начального (первого) элементарного сегмента, x1 = x0 + H0. На начальном элементарном сегменте [x0, x0 + H0] всегда применяется исключительно первый способ. Способ выбора начального приближения задается пользователем при обращении к подпрограмме и может меняться от сегмента к сегменту..
В дальнейшем при описании параметров подпрограммы коэффициенты ряда Чебышёва будем называть коэффициентами Чебышёва.
О.Б.Арушанян, С.Ф.Залеткин. Приближенное решение задачи Коши для обыкновенных дифференциальных уравнений методом рядов Чебышёва. Вычислительные методы и программирование: Новые вычислительные технологии (Электронный научный журнал) (17), 121 - 131, 2016.
О.Б.Арушанян, С.Ф.Залеткин. Использование рядов Чебышёва для приближенного аналитического решения обыкновенных дифференциальных уравнений. Вестник Московского университета. Серия 1: Математика. Механика. 5, 52 - 56, 2016.
О.Б.Арушанян, Н.И.Волченскова, С.Ф.Залеткин. Вычисление коэффициентов разложения решения задачи Коши в ряд по многочленам Чебышёва. Вестник Московского университета. Серия 1: Математика. Механика. 5, 24 - 30, 2012.
procedure DE44E(F :Proc_F70E; var M :Integer; var K :Integer;
INIAPR :Integer; IMAX :Integer; var JSTART :Integer;
var YX :Array of Extended; var X :Extended;
var XJ0 :Array of Extended; var AU :Array of Extended;
var AJK :Array of Extended; var H :Extended;
var HD4I :Array of Extended; var RABC :Array of Extended;
var ALP :Array of Extended; var ALPN :Array of Extended;
var CMAR :Array of Extended; var P :Array of Extended;
var S :Array of Extended; var XP :Extended;
var YP :Array of Extended; var AJKP :Array of Extended;
var U :Array of Extended; var ZFI :Array of Extended);
Параметры
| F - | имя подпрограммы вычисления значений правой части дифференциального уравнения. Первый оператор подпрограммы должен соответствовать процедурному типу: |
Procedure (X :Extended; var Y :Array of Extended; var Z :Array of Extended; M :Integer);
Здесь: X, Y - значения независимой и зависимой переменных, соответственно. Вычисленное значение правой части должно быть помещено в Z. Для скалярного случая при M = 1 параметры Y и Z представляют собой массивы длины 1. В случае системы уравнений, т.е. когда M ≠ 1 , параметры Y и Z представляют массивы длины M (тип параметров X, Y и Z: Extended); |
| M - | количество уравнений в системе (тип: целый); |
| K - | порядок частичной суммы смещенного ряда Чебышёва, с помощью которой аппроксимируется производная решения задачи Коши на элементарном сегменте [X, X + H] разбиения интервала интегрирования; при этом само решение задачи Коши приближается на элементарном сегменте [X, X + H] частичной суммой (K + 1) - го порядка; K ≥ 2. Если обращение к подпрограмме DE44E осуществляется со значением параметра JSTART = 1 (см. ниже), то значение параметра K должно совпадать с его значением при предыдущем обращении к подпрограмме. Если же при очередном обращении к подпрограмме DE44E значение параметра K необходимо изменить, то это обращение должно выполняться при нулевом значении параметра JSTART(JSTART = 0) (см. "Математическое описание", "Замечания по использованию" и "Примеры"; тип: целый); |
| INIAZR - |
целый указатель способа выбора начального
приближения коэффициентов Чебышёва для производной
решения на элементарном сегменте [X, X + H]: |
| INIAZR=1 - |
для первого способа, когда начальное
приближение определяется только с
использованием значения решения в
начале X элементарного сегмента; |
| INIAZR=2 - | для второго способа, когда начальное приближение коэффициентов Чебышёва на текущем элементарном сегменте [X, X + H] (начиная со второго) определяется через коэффициенты Чебышёва, вычисленные на предыдущем элементарном сегменте, т.е. путем экстраполяции коэффициентов с предыдущего сегмента на следующий (см. "Математическое описание"). |
| Значение этого параметра может меняться от сегмента к сегменту. | |
| IMAX - | целая переменная, задающая количество итераций, которое предполагается выполнить в итерационном процессе вычисления коэффициентов Чебышёва для производной решения задачи Коши на элементарном сегменте [X, X + H], исходя из некоторого начального приближения, способ определения которого задается параметром INIAZR; IMAX ≥ 1. Для получения максимального порядка точности приближенного решения необходимо выполнить не менее K итераций. Значение IMAX может изменяться от сегмента к сегменту (см. "Математическое описание", "Замечания по использованию" и "Примеры"); |
| JSTART - |
целый указатель режима использования подпрограммы, имеющий следующие значения: |
| 0 - |
первое обращение к подпрограмме должно быть исполнено с нулевым значением JSTART;
выполнить первый (начальный) шаг интегрирования для значений независимой и зависимой переменных и
шага интегрирования (длины элементарного сегмента), заданных параметрами X, YX, H,
соответственно. При данном значении параметра JSTART будет применен исключительно первый способ определения
начального приближения для коэффициентов Чебышёва производной решения независимо от значения параметра INIAZR.
Нулевое значение параметра JSTART также может означать, что необходимо выполнить очередной шаг
интегрирования с измененным значением параметра K; |
| 1 - | выполнить следующий (очередной) шаг интегрирования системы дифференциальных уравнений для значений независимой и зависимой переменных и шага интегрирования (длины элементарного сегмента), заданных параметрами X, YX, H, соответственно. При данном значении параметра JSTART способ определения начального приближения для коэффициентов Чебышёва производной решения определяется параметром INIAZR. При обращении к подпрограмме DE44E со значением JSTART = 1 значение параметра K (см. выше) должно совпадать с его значением при предыдущем обращении к подпрограмме. |
| На выходе из подпрограммы параметр JSTART всегда принимает значение, равное 1 (см. "Примеры"); | |
| X, YX - | начальное значение аргумента и решения (начало элементарного сегмента X и решение в нем YX); в результате работы подпрограммы в X получается новое значение аргумента, равное X + H (конец элементарного сегмента), а в YX - соответствующее значение решения. В случае системы уравнений, т.е. когда M ≠ 1, YX задается одномерным массивом длины M (тип параметров X, YX: с расширенной (Extended) точностью); |
| XJ0 - | одномерный рабочий массив длины K (тип: с расширенной (Extended) точностью); |
| AU - | двумерный массив размера M * (K + 2). На выходе из подпрограммы содержит коэффициенты Чебышёва ai*[Y] для решения Y(X + αH), 0 ≤ α ≤ 1, на элементарном сегменте [X, X + H]. При этом переменная с индексом AU(N, I + 1) представляет I-й коэффициент Чебышёва N-й компоненты решения yN(x) (I = 0, 1, ... , K + 1) (тип: с расширенной (Extended) точностью); |
| AJK - | двумерный массив размера M * (K + 1). На выходе из подпрограммы содержит коэффициенты Чебышёва ai*[Y'] производной решения Y'(X + αH), 0 ≤ α ≤ 1, на элементарном сегменте [X, X + H]. При этом переменная с индексом AJK(N, I + 1) представляет I-й коэффициент Чебышёва N-й компоненты производной решения y'N(x) (I = 0, 1, ... , K). Если при JSTART = 1 (т.е. при повторных обращениях к подпрограмме) параметр INIAZR = 2, то содержащиеся в массиве AJK при входе в подпрограмму коэффициенты Чебышёва производной решения, вычисленные на предыдущем элементарном сегменте [X - H', X] во время предыдущего обращения к подпрограмме DE44E, используются в подпрограмме при вычислении коэффициентов Чебышёва производной решения на текущем элементарном сегменте [X, X + H]. Заметим, что длина H текущего элементарного сегмента может быть больше или меньше длины H' предыдущего сегмента [X - H', X] или равна ей. Независимо от значения параметра INIAZR при повторных обращениях к подпрограмме значения массива AJK при входе в подпрограмму обязательно запоминаются в массиве AJKP и могут быть снова доступны на выходе из подпрограммы (тип: с расширенной (Extended) точностью); |
| H - | длина текущего элементарного сегмента [X, X + H]. Значение H может быть переменной величиной и изменяться от сегмента к сегменту (тип: с расширенной (Extended) точностью); |
| HD4I - | одномерный рабочий массив длины K + 1 (тип: с расширенной (Extended) точностью); |
|
RABC, - ALP, ALPN | одномерные рабочие массивы длины K (тип: с расширенной (Extended) точностью); |
| CMAR - | двумерный рабочий массив размера K * K (тип: с расширенной (Extended) точностью); |
| P - | одномерный рабочий массив размера (K + 2) * (K - 1) / 2 (тип: с расширенной (Extended) точностью); |
| S - | одномерный рабочий массив размера (K + 3) * K / 2 (тип: с расширенной (Extended) точностью); |
| XP, YP - | переменная и одномерный массив длины M, соответственно. На выходе из подпрограммы содержат начальные значения аргумента и решения, т.е. те значения, которые имели параметры X, YX на входе в подпрограмму. Таким образом, на выходе из подпрограммы DE44E границы элементарного сегмента, к которому относятся вычисленные подпрограммой коэффициенты Чебышёва решения и коэффициенты Чебышёва его производной и содержащиеся в массивах AU и AJK соответственно, показываются параметрами XP и X, а именно: параметр XP содержит начало сегмента, а параметр X содержит конец этого сегмента (тип: с расширенной (Extended) точностью); |
| AJKP - | двумерный массив размера M * (K + 1). Если обращение к подпрограмме было выполнено со значением JSTART = 1 (т.е. при повторных обращениях к подпрограмме), то на выходе из подпрограммы содержит коэффициенты Чебышёва производной решения, относящиеся к предыдущему элементарному сегменту [X - H', X]; тем самым значение этого параметра определено только при повторных обращениях к подпрограмме (тип: с расширенной (Extended) точностью); |
| U - | двумерный рабочий массив размера M * K (тип: с расширенной (Extended) точностью); |
| ZFI - | двумерный рабочий массив размера M * (K + 3) (тип: с расширенной (Extended) точностью). |
Версии: нет
Вызываемые подпрограммы: нет
|
DE44E использует рабочие подпрограммы DE70EK, DE70EH, DE70E0, DE70EI, DE70EF, DE70EQ, DE71EE, DE70EP, DE71ET, DE71EP, DE71EI, DE71EF, DE71ES, DE70EA, DE70EC. |
Замечания по использованию
|
Разбиение промежутка интегрирования на элементарные сегменты (шаги) выполняется для того, чтобы на каждом таком сегменте ряды Чебышёва для решения и его производной были быстросходящимися рядами. Другими словами, длина элементарных сегментов, задаваемая параметром H, подбирается таким образом, чтобы убывание коэффициентов этих рядов Чебышёва на элементарном сегменте [X, X + H] происходило достаточно быстро, вследствие чего можно было бы считать частичные суммы этих рядов близкими к многочленам наилучшего равномерного приближения на элементарном сегменте [X, X + H] для решения и его производной. Порядок этих частичных сумм задается параметром K. Если начальное приближение для коэффициентов Чебышёва функции Φ(α) определяется первым способом (т.е. при INIAZR = 1), то для получения максимального порядка точности приближенного решения в конце элементарного сегмента [X, X + H] необходимо выполнить в итерационном процессе не менее K итераций; тогда IMAX≥K. Если начальное приближение коэффициентов Чебышёва функции Φ(α) определяется вторым способом (т.е. при INIAZR = 2), то для получения максимального порядка точности приближенного решения необходимо выполнить в итерационном процессе не менее K + 1 итераций; в этом случае IMAX ≥ K + 1. Однако в некоторых случаях при втором способе определения начального приближения итерационный процесс может сойтись за значительно меньшее число итераций. Если длина H элементарного сегмента [X, X + H] подобрана достаточно малой, то хорошая точность приближенного решения может быть получена и с существенно меньшим числом итераций при любом способе выбора начального приближения. Вообще, число итераций зависит от K и H. С увеличением H число итераций может также возрастать. Если правая часть дифференциального уравнения не зависит от переменной Y, т.е. дифференциальное уравнение имеет вид Y'= F(X), то число итераций можно положить равным 1 при любых H и K, удовлетворяющих описанным выше условиям. В этом случае параметр IMAX = 1. Как следует из вышеописанного, управлять точностью приближенного решения задачи Коши можно с помощью четырех параметров H, K, IMAX, INIAZR, подбирая для любой конкретной задачи и для каждого элементарного сегмента наиболее подходящий набор их значений. При работе подпрограммы значения параметров M, K, INIAZR, IMAX, H сохраняются. При многократном использовании подпрограммы DE44E для вычисления коэффициентов Чебышёва решения задачи Коши (1), (2) и его производной на последовательности элементарных сегментов, образующей промежуток интегрирования [XN, XK] системы дифференциальных уравнений, значения параметров M, K, YX, X, AU, AJK, HD4I, RABC, ALP, ALPN, CMAR, P, S не должны изменяться в вызывающей программе между последовательными обращениями к подпрограмме. При интегрировании системы уравнений с помощью подпрограммы DE44E используется глобальная запись (структура данных) с именем _COM70D из модуля Lstruct для хранения промежуточных значений типа: type COM70D = record elm1: Extended; // WC1 elm2: Extended; // WC2 elm3: Extended; // WC3 elm4: Integer; // LASN elm5: Extended; // HD2 elm6: Extended; // HD4 elm7: Extended; // HOLD end; var _COM70D : COM70D; Эта структура данных используется и пользователь не должен портить её элементы, если вызывающая (главная) программа содержит хотя бы одно обращение к подпрограмме DE44E со значением параметра JSTART = 1 (и, следовательно, число обращений к DE44E более одного). Если вызывающая программа содержит одно обращение к подпрограмме DE44E со значением параметра JSTART = 0 и при этом других обращений к DE44E нет (как в первом примере), или если вызывающая программа содержит несколько обращений к DE44E и при этом все обращения осуществляются с нулевым значением параметра JSTART (как в четвертом примере), то указанная структура данных не задействована (см. "Примеры использования"). |
Использование подпрограммы DE44E иллюстрируется на четырех примерах. Во всех примерах решаются задачи Коши для одной и той же системы обыкновенных дифференциальных уравнений.
(7) y'1 = -2q (y2-1) + (1 -e1-y1+cos q(2x-1) ) / (x+1),
y'2 = 2q (y1-1) + ( 1 -e1-y2+sin q(2x-1) ) / (x+1), q = 0,5.
Частное решение системы (7) имеет вид
y1(x) = 1 + cos q(2x - 1), y2(x) = 1 + sin q(2x - 1).
Функции y1(x), y2(x) разлагаются на отрезке [0, 1] в смещенные ряды Чебышёва, коэффициенты которых выражаются через цилиндрические функции:
∞
(8) y1(x) = 1 + J0 (q) + 2 ∑ (-1)i J2i (q) T2i* (x),
i=1
∞
(9) y2(x) = 1 + 2 ∑ (-1)i J2i+1 (q) T2i+1* (x).
i=0
Задачи Коши отличаются начальными условиями. Вычисления во всех четырех примерах на Паскале проводились с 19-20 значащими цифрами.
1) Начальное условие задачи Коши для системы (7)
XN = 0, y1(0) = 1 + cos q, y2(0) = 1 - sin q.
Рассматривается один элементарный сегмент [0, 1]. Выполняется одно обращение к подпрограмме DE44E из начальной точки X = 0 с шагом H = 1. Вычисляются решение YX в конце данного сегмента, т.е. при x = 1, коэффициенты Чебышёва AU решения на данном сегменте и коэффициенты Чебышёва AJK его производной. Приводятся подпрограмма вычисления правой части системы (7), фрагмент вызывающей программы, результаты счета, включая точное значение решения YT в точке x = 1, абсолютную погрешность DELY приближенного решения YX, вычисленного в точке x = 1. Кроме вышеперечисленного приводятся: значения параметров XP, X на выходе из подпрограммы, которые представляют элементарный сегмент [0, 1]; значения параметров H, K, INIAPR, IMAX, при которых были вычислены приближенное решение YX и коэффициенты Чебышёва. Даются также точные значения коэффициентов Чебышёва на данном сегменте для компонент решения y1(x), y2(x) в (8), (9) и абсолютные погрешности их приближенных значений, вычисленных подпрограммой DE44E.
Unit Fsc_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc;
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
implementation
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
var
R,Q :Extended;
begin
Q := 0.5e0;
R := 2.e0*X-1.e0;
Z[0] := -2.e0*Q*(Y[1]-1.e0)+(1.e0-Exp(1.e0-Y[0]+Cos(Q*R)))/(1.e0+X);
Z[1] := 2.e0*Q*(Y[0]-1.e0)+(1.e0-Exp(1.e0-Y[1]+Sin(Q*R)))/(1.e0+X);
end;
end.
function tde44e1: String;
var
M,IMAX,JSTART,INIAPR,_i,NX,KP2,KP1,L,J :Integer;
H,XK,Q,X,ХР :Extended;
DELY :Array [0..1] of Extended;
YT :Array [0..1] of Extended;
YX :Array [0..1] of Extended;
XJ0 :Array [0..10] of Extended;
AU :Array [0..25] of Extended;
AJK :Array [0..23] of Extended;
HD4I :Array [0..11] of Extended;
RАВС :Array [0..10] of Extended;
ALP :Array [0..10] of Extended;
ALPN :Array [0..10] of Extended;
CMAR :Array [0..120] of Extended;
P :Array [0..64] of Extended;
S :Array [0..76] of Extended;
YP :Array [0..1] of Extended;
AJКР :Array [0..23] of Extended;
U :Array [0..21] of Extended;
ZFI :Array [0..27] of Extended;
const
K :Integer = 11;
label
_20;
begin
M := 2;
Q := 0.5e0;
X := 0.0;
YX[0] := Cos(Q)+1.e0;
YX[1] := -Sin(Q)+1.e0;
H := 1.0e0;
ХК := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IМАХ := 16;
JSTART := 0;
INIAPR := 1;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
end;
Результаты:
----------------------------------------------------------------
XP= 0,00000000000000000E+000 X= 1,00000000000000000E+000
H= 1,0000000000000000 K=11 INIAPR=1 IMAX=16 JSTART=1
YX= 1,87758256189037256E+000 1,47942553860420291E+000
YT= 1,87758256189037272E+000 1,47942553860420300E+000
DELTA Y = 1,54932490448178584E-016 9,13982431405280238E-017
COEFFICIENTS AU AND AJK ON 1 SEGMENT
------------------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,87693961448162577E+000 -7,16583097113560052E-016
1 -4,05644078440084432E-017 -4,84536915349748445E-001
2 -6,12080469173653167E-002 -5,54325465737526280E-016
3 -2,53084974310112398E-017 5,12745998917408848E-003
4 3,21472952728558635E-004 -2,50623496565391402E-016
5 -9,47372470186044757E-018 -1,61072544828496729E-005
6 -6,72136925728463897E-007 -6,11490025281824501E-017
7 -1,96923059765949767E-018 2,40317346334606051E-008
8 7,51644630286200614E-010 -6,01054579371651521E-018
9 -1,62818555416659954E-020 -2,08935356978145765E-011
10 -5,22635220826872873E-013 -5,42439899421653937E-018
11 -1,23281795323103168E-019 1,18731352603384382E-014
12 2,47356984590384130E-016
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 1,99999999999999998E+000 1,87693961448162531E+000
1 4,84536915349747752E-001 -4,75788570868107552E-016
2 -1,86398070372781333E-017 -6,12080469173656969E-002
3 -5,12745998917450329E-003 -3,26670114569882486E-016
4 -1,11774467719677473E-017 3,21472952728342590E-004
5 1,61072544827075208E-005 -1,47830966218398530E-016
6 -4,44861703897958538E-018 -6,72136925807826225E-007
7 -2,40317346578953275E-008 -4,10641572828884804E-017
8 -1,04523865691180662E-018 7,51644613242944803E-010
9 2,08935345562023931E-011 -7,61652026171066865E-018
10 5,81860812786869075E-018 -5,22630780341350187E-013
11 -1,18779722804852315E-014 -2,40360845376458299E-016
12 -5,00751761200954789E-018
**************************************************************************************
Number of Chebyshev coefficients for 1 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 3,87693961448162577E+000 3,87693961448162581E+000 4,27175655959288747E-017
1 -4,05644078440084432E-017 0,00000000000000000E+000 4,05644078440084432E-017
2 -6,12080469173653167E-002 -6,12080469173652826E-002 3,40913820610910800E-017
3 -2,53084974310112398E-017 0,00000000000000000E+000 2,53084974310112398E-017
4 3,21472952728558635E-004 3,21472952728575194E-004 1,65588588441111309E-017
5 -9,47372470186044757E-018 0,00000000000000000E+000 9,47372470186044757E-018
6 -6,72136925728463897E-007 -6,72136925723769759E-007 4,69413768647027798E-018
7 -1,96923059765949767E-018 0,00000000000000000E+000 1,96923059765949767E-018
8 7,51644630286200614E-010 7,51644630959521991E-010 6,73321377438174419E-019
9 -1,62818555416659954E-020 0,00000000000000000E+000 1,62818555416659954E-020
10 -5,22635220826872873E-013 -5,22635472164560617E-013 -2,51337687743974838E-019
11 -1,23281795323103168E-019 0,00000000000000000E+000 1,23281795323103168E-019
12 2,47356984590384130E-016 2,47676511895986538E-016 3,19527305602407805E-019
-----------------------------------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 1,99999999999999998E+000 2,00000000000000000E+000 2,24429849704499418E-017
1 4,84536915349747752E-001 4,84536915349747773E-001 2,11148373091551989E-017
2 -1,86398070372781333E-017 0,00000000000000000E+000 1,86398070372781333E-017
3 -5,12745998917450329E-003 -5,12745998917448815E-003 1,51432550310123815E-017
4 -1,11774467719677473E-017 0,00000000000000000E+000 1,11774467719677473E-017
5 1,61072544827075208E-005 1,61072544827149482E-005 7,42737217976061802E-018
6 -4,44861703897958538E-018 0,00000000000000000E+000 4,44861703897958538E-018
7 -2,40317346578953275E-008 -2,40317346555260458E-008 2,36928174945624626E-018
8 -1,04523865691180662E-018 0,00000000000000000E+000 1,04523865691180662E-018
9 2,08935345562023931E-011 2,08935351786579598E-011 6,22455566630061421E-019
10 5,81860812786869075E-018 0,00000000000000000E+000 -5,81860812786869075E-018
11 -1,18779722804852315E-014 -1,18837079244649228E-014 -5,73564397969128273E-018
12 -5,00751761200954789E-018 0,00000000000000000E+000 5,00751761200954789E-018
------------------------------------------------------------------------------------------------------------
2) Начальное условие задачи Коши для системы (7)
XN = - 0,5, y1(XN) = 1 + cos q (2 * XN - 1), y2(XN) = 1 + sin q (2 * XN - 1).
Рассматриваются два элементарных сегмента: [- 0,5, 0] и [0, 1]. Выполняются два обращения к подпрограмме DE44E. Первое обращение осуществляется из начальной точки X = - 0,5 с шагом H = 0,5. Второе обращение осуществляется из точки X = 0 с шагом H = 1. Приводятся фрагмент вызывающей программы и результаты счета. В качестве выходных данных, полученных после каждого обращения к подпрограмме DE44E, приводятся: значения параметров XP, X (которые последовательно представляют границы каждого из двух элементарных сегментов); значения параметров H, K, INIAPR, IMAX, при которых было вычислено приближенное значение решения YX в конце каждого элементарного сегмента; приближенное значение решения YX в конце каждого элементарного сегмента; точное значение решения YT в конце каждого элементарного сегмента и абсолютная погрешность DELY приближенного значения YX. После этого для каждого элементарного сегмента приводятся приближенные значения коэффициентов Чебышёва AU для компонент решения y1(x), y2(x) и приближенные значения коэффициентов Чебышёва AJK их производных y'1(x), y'2(x). Даются также точные значения коэффициентов Чебышёва решения в (8), (9) и абсолютные погрешности приближенных значений AU, вычисленных за два обращения к подпрограмме DE44E.
Unit Fsc_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc;
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
implementation
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
var
R,Q :Extended;
begin
Q := 0.5e0;
R := 2.e0*X-1.e0;
Z[0] := -2.e0*Q*(Y[1]-1.e0)+(1.e0-Exp(1.e0-Y[0]+Cos(Q*R)))/(1.e0+X);
Z[1] := 2.e0*Q*(Y[0]-1.e0)+(1.e0-Exp(1.e0-Y[1]+Sin(Q*R)))/(1.e0+X);
end;
end.
function tde44e2: String;
var
M,IMAX,JSTART,INIAPR,_i,NX,KP2,KP1,L,J,IERR :Integer;
H,XK,Q,V0,V1,W1,W2,X,XP :Extended;
DELY :Array [0..1] of Extended;
JK :Array [0..31] of Extended;
EXACT :Array [0..31] of Extended;
ERROR :Array [0..41] of Extended;
RAB :Array [0..12] of Extended;
YT :Array [0..1] of Extended;
YX :Array [0..1] of Extended;
XJ0 :Array [0..10] of Extended;
AU :Array [0..25] of Extended;
AJK :Array [0..23] of Extended;
HD4I :Array [0..11] of Extended;
RABC :Array [0..10] of Extended;
ALP :Array [0..10] of Extended;
ALPN :Array [0..10] of Extended;
CMAR :Array [0..120] of Extended;
P :Array [0..64] of Extended;
S :Array [0..76] of Extended;
YP :Array [0..1] of Extended;
AJKP :Array [0..23] of Extended;
U :Array [0..21] of Extended;
ZFI :Array [0..27] of Extended;
const
K :Integer = 11;
label
_1,_15,_20,_2,_8,_3,_4,_5,_6,_800,_300,_400,_500,_600;
begin
M := 2;
Q := 0.5e0;
X := -0.5e0;
YX[0] := Cos(Q*(2.e0*X-1.e0))+1.e0;
YX[1] := Sin(Q*(2.e0*X-1.e0))+1.e0;
H := 0.5e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 14;
JSTART := 0;
INIAPR := 1;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY;
. . . . . . . . . . . . . . . . . . . . . . .
H := 1.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 16;
INIAPR := 2;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
end;
Результаты:
----------------------------------------------------------------
RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
XP= -5,00000000000000000E-001 X= 0,00000000000000000E+000
H= 0,5000000000000000 K=11 INIAPR=1 IMAX=14 JSTART=1
YX= 1,87758256189037268E+000 5,20574461395797029E-001
YT= 1,87758256189037272E+000 5,20574461395797000E-001
DELTA Y = 3,17671236538252799E-017 -2,94360889829814454E-017
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 1 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,44060162317046235E+000 1,34205937233528711E+000
1 1,69081826785891476E-001 -1,81496854116472721E-001
2 -1,13732097613172041E-002 -1,05952419518446992E-002
3 -4,42044443620042861E-004 4,74502064602543738E-004
4 1,48397891478254549E-005 1,38246950363294697E-005
5 3,45797545843563867E-007 -3,71188127870819074E-007
6 -7,73596476359903190E-009 -7,20679741308494302E-009
7 -1,28728736309408262E-010 1,38180781934457188E-010
8 2,15954346014731081E-012 2,01182024191963807E-012
9 2,79467987621283430E-014 -2,99995149707042719E-014
10 -3,75047554319364596E-016 -3,49268953602627219E-016
11 -3,96896538184803658E-018 4,28937484489577692E-018
12 4,46809879676643429E-020
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 6,57940627664712603E-001 1,44060162317046259E+000
1 1,81496854116472453E-001 1,69081826785891700E-001
2 1,05952419518444771E-002 -1,13732097613170320E-002
3 -4,74502064602704752E-004 -4,42044443619932880E-004
4 -1,38246950364312117E-005 1,48397891478820657E-005
5 3,71188127816013586E-007 3,45797545865894704E-007
6 7,20679738754471149E-009 -7,73596475847772953E-009
7 -1,38180791086587680E-010 -1,28728736251447007E-010
8 -2,01182316307157535E-012 2,15954237118057800E-012
9 2,99988530615133786E-014 2,79461851338154432E-014
10 3,49376780896812691E-016 -3,75049248385259104E-016
11 -4,26192327710521709E-018 -3,95733792957209118E-018
12 -4,12222700997092832E-020
----------------------------------------------------------------
RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
XP= 0,00000000000000000E+000 X= 1,00000000000000000E+000
H= 1,0000000000000000 K=11 INIAPR=1 IMAX=16 JSTART=1
YX= 1,87758256189037248E+000 1,47942553860420289E+000
YT= 1,87758256189037272E+000 1,47942553860420300E+000
DELTA Y = 2,38741318381308076E-016 1,13082286590238112E-016
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 2 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,87693961448162564E+000 -8,15099805142813216E-016
1 -6,73399663396641313E-017 -4,84536915349748425E-001
2 -6,12080469173653141E-002 -5,45739939784156691E-016
3 -2,45862606713190731E-017 5,12745998917408839E-003
4 3,21472952728558631E-004 -2,50704811728327814E-016
5 -9,47812927318616993E-018 -1,61072544828497068E-005
6 -6,72136925728466156E-007 -6,11422262646044157E-017
7 -1,96777854117849030E-018 2,40317346334809338E-008
8 7,51644630286306493E-010 -6,04442711160668722E-018
9 -1,62818555416659954E-020 -2,08935356808739175E-011
10 -5,22635220742169579E-013 -5,45828031210671138E-018
11 -1,24051825275152531E-019 1,18731488128655943E-014
12 2,47357266934699882E-016
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 2,00000000000000000E+000 1,87693961448162521E+000
1 4,84536915349747726E-001 -5,23080114379209649E-016
2 -2,46689875558342431E-017 -6,12080469173656950E-002
3 -5,12745998917450314E-003 -3,25728213932535704E-016
4 -1,11198485315544548E-017 3,21472952728342583E-004
5 1,61072544827075198E-005 -1,47810637427664426E-016
6 -4,44579359582207105E-018 -6,72136925807812673E-007
7 -2,40317346578958115E-008 -4,11115911279347213E-017
8 -1,04883854693763740E-018 7,51644613270049857E-010
9 2,08935345578964590E-011 -7,54875762593032462E-018
10 5,82131863329990451E-018 -5,22630814222668078E-013
11 -1,18779730505151836E-014 -2,40401502957926505E-016
12 -5,00836464495680219E-018
**************************************************************************************
Number of Chebyshev coefficients for 1 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 3,87693961448162564E+000 3,87693961448162581E+000 1,63497687610814069E-016
1 -6,73399663396641313E-017 0,00000000000000000E+000 6,73399663396641313E-017
2 -6,12080469173653141E-002 -6,12080469173652826E-002 3,15231781650160414E-017
3 -2,45862606713190731E-017 0,00000000000000000E+000 2,45862606713190731E-017
4 3,21472952728558631E-004 3,21472952728575194E-004 1,65626704923737753E-017
5 -9,47812927318616993E-018 0,00000000000000000E+000 9,47812927318616993E-018
6 -6,72136925728466156E-007 -6,72136925723769759E-007 4,69639645822921888E-018
7 -1,96777854117849030E-018 0,00000000000000000E+000 1,96777854117849030E-018
8 7,51644630286306493E-010 7,51644630959521991E-010 6,73215498319767631E-019
9 -1,62818555416659954E-020 0,00000000000000000E+000 1,62818555416659954E-020
10 -5,22635220742169579E-013 -5,22635472164560617E-013 -2,51422391038729850E-019
11 -1,24051825275152531E-019 0,00000000000000000E+000 1,24051825275152531E-019
12 2,47357266934699882E-016 2,47676511895986538E-016 3,19244961286656364E-019
-----------------------------------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 2,00000000000000000E+000 2,00000000000000000E+000 2,71050543121376108E-018
1 4,84536915349747726E-001 4,84536915349747773E-001 4,68104287970616539E-017
2 -2,46689875558342431E-017 0,00000000000000000E+000 2,46689875558342431E-017
3 -5,12745998917450314E-003 -5,12745998917448815E-003 1,49840128369285730E-017
4 -1,11198485315544548E-017 0,00000000000000000E+000 1,11198485315544548E-017
5 1,61072544827075198E-005 1,61072544827149482E-005 7,42838795755283314E-018
6 -4,44579359582207105E-018 0,00000000000000000E+000 4,44579359582207105E-018
7 -2,40317346578958115E-008 -2,40317346555260458E-008 2,36976576805245055E-018
8 -1,04883854693763740E-018 0,00000000000000000E+000 1,04883854693763740E-018
9 2,08935345578964590E-011 2,08935351786579598E-011 6,20761500735552820E-019
10 5,82131863329990451E-018 0,00000000000000000E+000 -5,82131863329990451E-018
11 -1,18779730505151836E-014 -1,18837079244649228E-014 -5,73487394973923330E-018
12 -5,00836464495680219E-018 0,00000000000000000E+000 5,00836464495680219E-018
------------------------------------------------------------------------------------------------------------
3) Начальное условие задачи Коши для системы (7)
XN = - 4/9, y1(XN) = 1 + cos q (2 * XN - 1), y2(XN) = 1 + sin q (2 * XN - 1).
Рассматриваются три элементарных сегмента: [- 4/9, - 1/3], [- 1/3, 0] и [0, 1]. Выполняются три обращения к подпрограмме DE44E. Первое обращение осуществляется из начальной точки X = - 4/9 с шагом H = 1/9. Второе обращение осуществляется из точки X = - 1/3 с шагом H = 1/3. Третье обращение осуществляется из точки X = 0 с шагом H = 1. Приводятся фрагмент вызывающей программы и результаты счета, аналогичные результатам из примера 2. А именно, после каждого обращения к подпрограмме DE44E приводятся: значения параметров XP, X (которые последовательно представляют границы каждого из трех элементарных сегментов); значения параметров H, K, INIAPR, IMAX, при которых было вычислено приближенное значение решения YX в конце каждого элементарного сегмента; приближенное значение решения YX в конце каждого элементарного сегмента; точное значение решения YT в конце каждого элементарного сегмента и абсолютная погрешность DELY приближенного значения YX. После этого для каждого элементарного сегмента приводятся приближенные значения коэффициентов Чебышёва AU для компонент решения y1(x) и y2(x) и приближенные значения коэффициентов Чебышёва AJK их производных y'1(x), y'2(x). Даются также, как и в примерах 1 и 2, точные значения коэффициентов Чебышёва для компонент решения y1(x), y2(x) в (8), (9) и абсолютные погрешности приближенных значений AU, вычисленных за три обращения к подпрограмме DE44E.
Unit Fsc_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc;
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
implementation
procedure Fsc(X :Extended; var Y :Array of Extended; var Z :Array of Extended;
M :Integer);
var
R,Q :Extended;
begin
Q := 0.5e0;
R := 2.e0*X-1.e0;
Z[0] := -2.e0*Q*(Y[1]-1.e0)+(1.e0-Exp(1.e0-Y[0]+Cos(Q*R)))/(1.e0+X);
Z[1] := 2.e0*Q*(Y[0]-1.e0)+(1.e0-Exp(1.e0-Y[1]+Sin(Q*R)))/(1.e0+X);
end;
end.
function tde44e3: String;
var
M,IMAX,JSTART,INIAPR,_i,NX,KP2,KP1,L,J,IERR :Integer;
H,XK,Q,V0,V1,W1,W2,X,XP :Extended;
DELY :Array [0..1] of Extended;
JK :Array [0..31] of Extended;
EXACT :Array [0..31] of Extended;
ERROR :Array [0..41] of Extended;
RAB :Array [0..12] of Extended;
YT :Array [0..1] of Extended;
YX :Array [0..1] of Extended;
XJ0 :Array [0..10] of Extended;
AU :Array [0..25] of Extended;
AJK :Array [0..23] of Extended;
HD4I :Array [0..11] of Extended;
RABC :Array [0..10] of Extended;
ALP :Array [0..10] of Extended;
ALPN :Array [0..10] of Extended;
CMAR :Array [0..120] of Extended;
P :Array [0..64] of Extended;
S :Array [0..76] of Extended;
YP :Array [0..1] of Extended;
AJKP :Array [0..23] of Extended;
U :Array [0..21] of Extended;
ZFI :Array [0..27] of Extended;
const
K :Integer = 11;
label
_1,_2,_15,_20,_7,_8,_3,_4,_5,_6,_800,_300,_400,_500,_600;
begin
Result := ''; { результат функции }
M := 2;
Q := 0.5e0;
X := -4.e0/9.e0;
YX[0] := Cos(Q*(2.e0*X-1.e0))+1.e0;
YX[1] := Sin(Q*(2.e0*X-1.e0))+1.e0;
H := 1.e0/9.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 10;
JSTART := 0;
INIAPR := 1;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
H := 1.e0/3.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
INIAPR := 2;
IMAX := 14;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
H := 1.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
INIAPR := 2;
IMAX := 17;
DE44E(Fsc,M,K,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE THIRD CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
end;
Результаты:
----------------------------------------------------------------
RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
XP= -4,44444444444444444E-001 X= -3,33333333333333333E-001
H= 0,1111111111111111 K=11 INIAPR=1 IMAX=10 JSTART=1
YX= 1,67241224408305669E+000 2,59823146803962944E-001
YT= 1,67241224408305669E+000 2,59823146803962944E-001
DELTA Y = -1,08420217248550443E-019 5,42101086242752217E-020
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 1 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,25957764277102818E+000 1,55154596888058375E+000
1 4,31151352007385917E-002 -3,50017733622691490E-002
2 -4,86198270661226205E-004 -5,98898346005547975E-004
3 -5,54571166090718438E-006 4,50212533913777602E-006
4 3,12659655792634606E-008 3,85133724279375042E-008
5 2,13968683567506418E-010 -1,73704276162304615E-010
6 -8,04201239859050574E-013 -9,90614213651052483E-013
7 -3,93106208018041328E-015 3,19164725030851582E-015
8 1,10827908401120792E-017 1,34305544116641862E-017
9 4,08039822245219744E-020 -1,96511643762997679E-019
10 -6,77626357803440271E-022 2,10064170919066484E-019
11 5,30465078078450717E-022 4,74338450462408190E-020
12 1,09800567236668562E-022
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 4,48454031119416248E-001 1,25957764277102818E+000
1 3,50017733622691502E-002 4,31151352007385895E-002
2 5,98898346005548580E-004 -4,86198270661228160E-004
3 -4,50212533913709962E-006 -5,54571166090821801E-006
4 -3,85133724277203108E-008 3,12659655785981632E-008
5 1,73704276556476103E-010 2,13968683506745921E-010
6 9,90613955384432547E-013 -8,04201567535314799E-013
7 -3,19131942166799357E-015 -3,93085629150913274E-015
8 -1,36484654016918446E-017 1,09267250195804744E-017
9 3,32538860773910503E-020 -9,82558218814988393E-020
10 -5,08219768352580203E-022 1,52465930505774061E-019
11 3,85014976024681972E-022 8,47032947254300339E-020
12 1,96072441494051004E-022
----------------------------------------------------------------
RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
XP= -3,33333333333333333E-001 X= 2,71050543121376108E-020
H= 0,3333333333333333 K=11 INIAPR=1 IMAX=14 JSTART=1
YX= 1,87758256189037272E+000 5,20574461395797000E-001
YT= 1,87758256189037272E+000 5,20574461395797000E-001
DELTA Y = -2,16840434497100887E-019 -5,42101086242752217E-020
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 2 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,56087835592218985E+000 1,22816603568964620E+000
1 1,02704194891647743E-001 -1,30526940340578536E-001
2 -5,44492816475257603E-003 -4,28430301012671280E-003
3 -1,19077351393356962E-004 1,51335613483289228E-004
4 3,15392064583584882E-006 2,48164003413782946E-006
5 4,13702452651637928E-008 -5,25775168315156219E-008
6 -7,30364065563845305E-010 -5,74681771998107865E-010
7 -6,84229824554752899E-012 8,69588908124002344E-012
8 9,05909168390864790E-014 7,12806278845706798E-014
9 6,60055318116053231E-016 -8,38935312278549228E-016
10 -6,99180522867904701E-018 -5,34647196306914374E-018
11 -4,05035754777965435E-020 8,13151629364128326E-020
12 5,64688631502866893E-022
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 7,71833964310353801E-001 1,56087835592218985E+000
1 1,30526940340578536E-001 1,02704194891647743E-001
2 4,28430301012671249E-003 -5,44492816475257619E-003
3 -1,51335613483289221E-004 -1,19077351393356926E-004
4 -2,48164003413796464E-006 3,15392064583576976E-006
5 5,25775168316943459E-008 4,13702452653767934E-008
6 5,74681771712563647E-010 -7,30364065890988252E-010
7 -8,69588876949156395E-012 -6,84229792778916906E-012
8 -7,12808110907385630E-014 9,05907463031197652E-014
9 8,38866451637096517E-016 6,59936921732981463E-016
10 5,49882495584861723E-018 -6,83047368665867793E-018
11 -5,17460127777172571E-020 7,79270311473956312E-020
12 5,41159938523580772E-022
----------------------------------------------------------------
RESULTS AFTER THE THIRD CALL SUBROUTINE DE44E
XP= 2,71050543121376108E-020 X= 1,00000000000000000E+000
H= 1,0000000000000000 K=11 INIAPR=2 IMAX=17 JSTART=1
YX= 1,87758256189037272E+000 1,47942553860420300E+000
YT= 1,87758256189037272E+000 1,47942553860420300E+000
DELTA Y = -4,33680868994201774E-019 -2,16840434497100887E-018
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 3 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,87693961448162581E+000 2,30392961653169692E-019
1 8,13151629364128326E-020 -4,84536915349747773E-001
2 -6,12080469173652827E-002 -9,48676900924816380E-020
3 -1,69406589450860068E-021 5,12745998917448825E-003
4 3,21472952728575194E-004 -7,45388993583784298E-020
5 2,37169225231204095E-021 -1,61072544827148524E-005
6 -6,72136925723773170E-007 -1,21972744404619249E-019
7 -2,05708001476044368E-021 2,40317346557036903E-008
8 7,51644630957897741E-010 -6,43745039913268258E-020
9 1,58206931592719874E-019 -2,08935349490374511E-011
10 -5,22635199989862371E-013 -5,75982404132924231E-018
11 -1,30905091848391871E-019 1,18730505570437128E-014
12 2,47355219938410684E-016
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 2,00000000000000000E+000 1,87693961448162581E+000
1 4,84536915349747773E-001 3,72694496791892149E-019
2 4,06575814682064163E-020 -6,12080469173652828E-002
3 -5,12745998917448816E-003 4,74338450462408190E-020
4 -1,69406589450860068E-020 3,21472952728575097E-004
5 1,61072544827149649E-005 3,18484388167616928E-019
6 -5,64688631502866893E-022 -6,72136925724200356E-007
7 -2,40317346555415439E-008 3,32036915323685733E-019
8 4,48927462044779180E-020 7,51644630962874059E-010
9 2,08935349957183780E-011 -1,10453096321960764E-018
10 5,97090465178501395E-018 -5,22628882987548338E-013
11 -1,18779291588079168E-014 -2,39940717034620166E-016
12 -4,99876493822125345E-018
**************************************************************************************
Number of Chebyshev coefficients for 1 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 3,87693961448162581E+000 3,87693961448162581E+000 -8,67361737988403547E-019
1 8,13151629364128326E-020 0,00000000000000000E+000 -8,13151629364128326E-020
2 -6,12080469173652827E-002 -6,12080469173652826E-002 4,06575814682064163E-020
3 -1,69406589450860068E-021 0,00000000000000000E+000 1,69406589450860068E-021
4 3,21472952728575194E-004 3,21472952728575194E-004 -2,11758236813575085E-022
5 2,37169225231204095E-021 0,00000000000000000E+000 -2,37169225231204095E-021
6 -6,72136925723773170E-007 -6,72136925723769759E-007 3,41041396676784854E-021
7 -2,05708001476044368E-021 0,00000000000000000E+000 2,05708001476044368E-021
8 7,51644630957897741E-010 7,51644630959521991E-010 1,62425026551425001E-021
9 1,58206931592719874E-019 0,00000000000000000E+000 -1,58206931592719874E-019
10 -5,22635199989862371E-013 -5,22635472164560617E-013 -2,72174698246460209E-019
11 -1,30905091848391871E-019 0,00000000000000000E+000 1,30905091848391871E-019
12 2,47355219938410684E-016 2,47676511895986538E-016 3,21291957575854272E-019
-----------------------------------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 2,00000000000000000E+000 2,00000000000000000E+000 -2,16840434497100887E-019
1 4,84536915349747773E-001 4,84536915349747773E-001 -2,43945488809238498E-019
2 4,06575814682064163E-020 0,00000000000000000E+000 -4,06575814682064163E-020
3 -5,12745998917448816E-003 -5,12745998917448815E-003 1,01643953670516041E-020
4 -1,69406589450860068E-020 0,00000000000000000E+000 1,69406589450860068E-020
5 1,61072544827149649E-005 1,61072544827149482E-005 -1,66925047613200985E-020
6 -5,64688631502866893E-022 0,00000000000000000E+000 5,64688631502866893E-022
7 -2,40317346555415439E-008 -2,40317346555260458E-008 1,54981932816999074E-020
8 4,48927462044779180E-020 0,00000000000000000E+000 -4,48927462044779180E-020
9 2,08935349957183780E-011 2,08935351786579598E-011 1,82939581777697899E-019
10 5,97090465178501395E-018 0,00000000000000000E+000 -5,97090465178501395E-018
11 -1,18779291588079168E-014 -1,18837079244649228E-014 -5,77876565700604718E-018
12 -4,99876493822125345E-018 0,00000000000000000E+000 4,99876493822125345E-018
-------------------------------------------------------------------------------------------------------------
Заметим, что во втором и третьем примерах последовательные обращения к подпрограмме DE44E выполнялись при одном и том же значении K (K = 11). При этом изменялись только параметры H, IMAX, INIAPR.
4) Четвертый пример служит иллюстрацией того, как последовательные обращения к подпрограмме DE44E можно выполнять не только с разными значениями параметров H, INIAPR, IMAX, но и с разными значениями параметра K.
Рассматриваются начальные условия и элементарные сегменты из примера 3. Выполняются также, как и в примере 3, три обращения к подпрограмме. Первое обращение осуществляется из начальной точки X = - 4/9 с шагом H = 1/9 со значением параметра K = 6. Второе обращение осуществляется из точки X = - 1/3 с шагом H = 1/3 и со значением параметра K = 8. Третье обращение осуществляется из точки X = 0 с шагом H = 1 и со значением K = 11. Дополнительным и существенным отличием является только то, что каждое последующее обращение к подпрограмме DE44E с новым значением параметра K (отличным от значения K при предыдущем обращении) выполняется с нулевым значением параметра JSTART. Приводится фрагмент вызывающей программы и результаты счета, полностью аналогичные результатам из примера 3. Обратим внимание на то, что в примере 4 дается значение параметра JSTART на входе в подпрограмму DE44E при каждом обращении к ней, в то время как в первых трех примерах приведено значение параметра JSTART на выходе из подпрограммы DE44E, которое всегда полагается равным 1.
function tde44e4: String;
var
M,IMAX,JSTARN,JSTART,INIAPR,K0,_i,NX,KP2,KP1,L,J :Integer;
H,XK,Q,X,XP :Extended;
DELY :Array [0..1] of Extended;
YT :Array [0..1] of Extended;
YX :Array [0..1] of Extended;
XJ0 :Array [0..13] of Extended;
AU :Array [0..31] of Extended;
AJK :Array [0..29] of Extended;
HD4I :Array [0..14] of Extended;
RABC :Array [0..13] of Extended;
ALP :Array [0..13] of Extended;
ALPN :Array [0..13] of Extended;
CMAR :Array [0..195] of Extended;
P :Array [0..103] of Extended;
S :Array [0..118] of Extended;
YP :Array [0..1] of Extended;
AJKP :Array [0..29] of Extended;
U :Array [0..27] of Extended;
ZFI :Array [0..33] of Extended;
const
K :Integer = 14;
label
_1,_2,_15,_20,_7;
begin
Result := ''; { результат функции }
{ }
{ take into consideraitoin: all arraies dependent on K }
{ must be described for maximal value K }
{ }
M := 2;
Q := 0.5e0;
X := -4.e0/9.e0;
YX[0] := Cos(Q*(2.e0*X-1.e0))+1.e0;
YX[1] := Sin(Q*(2.e0*X-1.e0))+1.e0;
H := 1.e0/9.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 10;
JSTARN := 0;
JSTART := JSTARN;
INIAPR := 1;
K0 := 6;
DE44E(Fsc,M,K0,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
H := 1.e0/3.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 14;
{ }
{ take into consideraitoin: parameter JSTART must be equal to zero }
{ }
JSTART := JSTARN;
K0 := 8;
DE44E(Fsc,M,K0,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
H := 1.e0;
XK := X+H;
YT[0] := Cos(Q*(2.e0*XK-1.e0))+1.e0;
YT[1] := Sin(Q*(2.e0*XK-1.e0))+1.e0;
IMAX := 17;
{ }
{ take into consideraitoin: parameter JSTART must be equal to zero }
{ }
JSTART := JSTARN;
K0 := 11;
DE44E(Fsc,M,K0,INIAPR,IMAX,JSTART,YX,X,XJ0,AU,AJK,H,
HD4I,RABC,ALP,ALPN,CMAR,P,S,XP,YP,AJKP,U,ZFI);
DELY[0] := YT[0]-YX[0];
DELY[1] := YT[1]-YX[1];
// RESULTS AFTER THE THIRD CALL SUBROUTINE DE44E
// Операторы вывода на печать: H, K, INIAPR, IMAX,JSTART,YX,YT, DELY,AU,AJK
. . . . . . . . . . . . . . . . . . . . . . .
end;
Результаты:
----------------------------------------------------------------
RESULTS AFTER THE FIRST CALL SUBROUTINE DE44E
XP= -4,44444444444444444E-001 X= -3,33333333333333333E-001
H= 0,1111111111111111 K= 6 INIAPR=1 IMAX=10 JSTART(INITIAL)=0
YX= 1,67241224408305670E+000 2,59823146803962931E-001
YT= 1,67241224408305669E+000 2,59823146803962944E-001
DELTA Y = -1,06251812903579434E-017 1,26309553094561267E-017
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 1 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,25957764277102819E+000 1,55154596888058376E+000
1 4,31151352007385919E-002 -3,50017733622691488E-002
2 -4,86198270661226204E-004 -5,98898346005547704E-004
3 -5,54571166090715382E-006 4,50212533913792510E-006
4 3,12659655788080392E-008 3,85133724249085143E-008
5 2,13968700840541091E-010 -1,73704210432547908E-010
6 -8,04186159409944021E-013 -9,93726388882117812E-013
7 -3,94335868604015005E-015
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 4,48454031119416244E-001 1,25957764277102818E+000
1 3,50017733622691501E-002 4,31151352007385898E-002
2 5,98898346005548586E-004 -4,86198270661227916E-004
3 -4,50212533913712104E-006 -5,54571166090838064E-006
4 -3,85133724274543189E-008 3,12659655811562028E-008
5 1,73704254886418109E-010 2,13968645041285720E-010
6 9,90595578894841298E-013 -8,00298399056785514E-013
7 -3,17578729784438696E-015
----------------------------------------------------------------
RESULTS AFTER THE SECOND CALL SUBROUTINE DE44E
XP= -3,33333333333333333E-001 X= 2,71050543121376108E-020
H= 0,3333333333333333 K= 8 INIAPR=1 IMAX=14 JSTART(INITIAL)=0
YX= 1,87758256189037273E+000 5,20574461395796990E-001
YT= 1,87758256189037272E+000 5,20574461395797000E-001
DELTA Y = -1,63714528045311170E-017 9,26992857475106291E-018
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 2 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,56087835592218988E+000 1,22816603568964626E+000
1 1,02704194891647748E-001 -1,30526940340578533E-001
2 -5,44492816475257593E-003 -4,28430301012671326E-003
3 -1,19077351393356975E-004 1,51335613483289241E-004
4 3,15392064583585983E-006 2,48164003413782269E-006
5 4,13702452651266363E-008 -5,25775168320306179E-008
6 -7,30364065207526778E-010 -5,74681769775493411E-010
7 -6,84230781548710668E-012 8,69586291131008507E-012
8 9,05819053261467195E-014 7,20867254235500743E-014
9 6,67469679847685873E-016
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 7,71833964310353776E-001 1,56087835592218985E+000
1 1,30526940340578535E-001 1,02704194891647751E-001
2 4,28430301012671283E-003 -5,44492816475257603E-003
3 -1,51335613483289215E-004 -1,19077351393356895E-004
4 -2,48164003413796259E-006 3,15392064583569184E-006
5 5,25775168317145052E-008 4,13702452653090308E-008
6 5,74681771645271585E-010 -7,30364067178478332E-010
7 -8,69588092326403392E-012 -6,84229315052334655E-012
8 -7,12738869846181932E-014 8,99303757005177895E-014
9 8,32688663893683236E-016
----------------------------------------------------------------
RESULTS AFTER THE THIRD CALL SUBROUTINE DE44E
XP= 2,71050543121376108E-020 X= 1,00000000000000000E+000
H= 1,0000000000000000 K=11 INIAPR=1 IMAX=17 JSTART(INITIAL)=0
YX= 1,87758256189037275E+000 1,47942553860420301E+000
YT= 1,87758256189037272E+000 1,47942553860420300E+000
DELTA Y = -3,27429056090622339E-017 -8,34835672813838414E-018
----------------------------------------------------------------
COEFFICIENTS AU AND AJK ON 3 SEGMENT
----------------------------------------------------------------
Number of Chebyshev coefficients for 1 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 3,87693961448162586E+000 3,18653794757067788E-017
1 8,48642309854083510E-018 -4,84536915349747786E-001
2 -6,12080469173652846E-002 -2,08031291845656163E-018
3 -4,22669440679895869E-019 5,12745998917449080E-003
4 3,21472952728575173E-004 2,99172036970218880E-018
5 4,45539330255761978E-020 -1,61072544827119691E-005
6 -6,72136925723714724E-007 2,10064170919066484E-018
7 5,07009721285074060E-020 2,40317346571843038E-008
8 7,51644630993579004E-010 6,81014489592457473E-019
9 1,74394672362468725E-019 -2,08935346102242722E-011
10 -5,22635190587796656E-013 -5,59719371545641664E-018
11 -1,27208948078554924E-019 1,18730132875940336E-014
12 2,47354443491542367E-016
----------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
---------------------------------------------------
coefficient for Y for Y'
----------------------------------------------------------------
0 2,00000000000000000E+000 1,87693961448162581E+000
1 4,84536915349747784E-001 -3,31494814237442981E-017
2 9,31736241979730373E-021 -6,12080469173653280E-002
3 -5,12745998917449008E-003 -3,32240203231026765E-017
4 -1,27309051972321341E-018 3,21472952728553013E-004
5 1,61072544827142049E-005 -1,28545720075312620E-017
6 -4,20975374785387269E-019 -6,72136925731085040E-007
7 -2,40317346557457999E-008 -2,75116301268196750E-018
8 -3,81164826264435153E-020 7,51644629797356724E-010
9 2,08935349659781101E-011 -1,53143556863577501E-018
10 5,96209550913356923E-018 -5,22628977855238430E-013
11 -1,18779313148917825E-014 -2,40015255933978544E-016
12 -5,00031783195788634E-018
**************************************************************************************
Number of Chebyshev coefficients for 1 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 3,87693961448162586E+000 3,87693961448162581E+000 -5,29090660172926164E-017
1 8,48642309854083510E-018 0,00000000000000000E+000 -8,48642309854083510E-018
2 -6,12080469173652846E-002 -6,12080469173652826E-002 1,94817577868489078E-018
3 -4,22669440679895869E-019 0,00000000000000000E+000 4,22669440679895869E-019
4 3,21472952728575173E-004 3,21472952728575194E-004 2,03287907341032081E-020
5 4,45539330255761978E-020 0,00000000000000000E+000 -4,45539330255761978E-020
6 -6,72136925723714724E-007 -6,72136925723769759E-007 -5,50348593937788749E-020
7 5,07009721285074060E-020 0,00000000000000000E+000 -5,07009721285074060E-020
8 7,51644630993579004E-010 7,51644630959521991E-010 -3,40570126375731518E-020
9 1,74394672362468725E-019 0,00000000000000000E+000 -1,74394672362468725E-019
10 -5,22635190587796656E-013 -5,22635472164560617E-013 -2,81576763960963221E-019
11 -1,27208948078554924E-019 0,00000000000000000E+000 1,27208948078554924E-019
12 2,47354443491542367E-016 2,47676511895986538E-016 3,22068404444170706E-019
-----------------------------------------------------------------------------------------
Number of Chebyshev coefficients for 2 component
----------------------------------------------------------------------------------------------
coefficient approximate exact absolute error
------------------------------------------------------------------------------------------------------------------
0 2,00000000000000000E+000 2,00000000000000000E+000 -6,50521303491302660E-019
1 4,84536915349747784E-001 4,84536915349747773E-001 -1,09233368877914572E-017
2 9,31736241979730373E-021 0,00000000000000000E+000 -9,31736241979730373E-021
3 -5,12745998917449008E-003 -5,12745998917448815E-003 1,93165863621343192E-018
4 -1,27309051972321341E-018 0,00000000000000000E+000 1,27309051972321341E-018
5 1,61072544827142049E-005 1,61072544827149482E-005 7,43266448131973229E-019
6 -4,20975374785387269E-019 0,00000000000000000E+000 4,20975374785387269E-019
7 -2,40317346557457999E-008 -2,40317346555260458E-008 2,19754137815140579E-019
8 -3,81164826264435153E-020 0,00000000000000000E+000 3,81164826264435153E-020
9 2,08935349659781101E-011 2,08935351786579598E-011 2,12679849702639044E-019
10 5,96209550913356923E-018 0,00000000000000000E+000 -5,96209550913356923E-018
11 -1,18779313148917825E-014 -1,18837079244649228E-014 -5,77660957314030861E-018
12 -5,00031783195788634E-018 0,00000000000000000E+000 5,00031783195788634E-018
-----------------------------------------------------------------------------------------------------------
В процессе интегрирования дифференциальных уравнений (1), (2) с помощью подпрограммы DE44E нужно следить за тем, чтобы параметры H, K, INIAPR, IMAX удовлетворяли приведенным выше условиям, сформулированным в "Математическом описании" и в "Замечаниях по использованию". Особенно это относится к длине элементарного сегмента H и параметру K, определяющему порядок используемых частичных сумм рядов Чебышёва для решения и его производной. Необходимо заботиться о том, чтобы указанные ряды Чебышёва на каждом сегменте [X, X + H] длиной H быстро сходились, при этом сумма ряда Чебышёва для решения должна хорошо приближаться частичной суммой порядка K + 1, а сумма ряда Чебышёва для производной должна хорошо приближаться частичной суммой K-го порядка.