Вторая пачка часть 221
4 и n = 2, x = 1.4, eps = 1e-8. Полученные при этом результаты
0.1775967712, 0.5419477138,
0.1775967712, -0.5419477138,
0.2073558994, 0.2073558994
совпадают с табличными [Справочник..., 1979] вплоть до девятого десятичного знака.
В случаях когда v = 0, 1, для вычисления функции Бес¬селя I рода вместо формулы (5.23) может быть ис¬поль¬зовано ее разложение в ряд [Справочник ..., 1979]:
a) при |x| ? 3
; (5.27)
б) при x > 3
(5.28)
где коэффициенты ai, bi, ci при v = 0 и v = 1 вы-би¬ра¬ют сответственно из табл. 5.9a и 5.9б
для v = 0
|?(x)| < 5?10-8, |?1(x)| < 7?10-8, |?2(x)| < 5?10-8;
для v = 1
|?(x)| < 1.3?10-8, |?1(x)| < 9?10-8, |?2(x)| < 4?10-8.
Таблица 5.9a
i ai bi ci
0 -0.78539816 0.79788456
1 -2.2499997 -0.04166397 -7.7e-7
2 1.2656208 -0.00003954 -0.00552740
3 -0.3163866 0.00262573 -0.00009512
4 0.0444479 -0.00054125 0.00137237
5 -0.0039444 -0.00029333 -0.00072805
6 0.0002100 0.00013558 0.00014476
Таблица 5.9б
i ai bi ci
0 -2.35619449 0.79788456
1 -0.56249985 0.12499612 0.00000156
2 0.21093573 0.00005650 0.01659667
3 -0.03954289 -0.00637879 0.00017105
4 0.00443319 0.00074348 -0.00249511
5 -0.00031761 0.00079824 0.00113653
6 0.00001109 -0.00029166 -0.00020033
В приведенной далее процедуре-функции JNX01 для вы¬¬числения функции Бесселя при v = 0; v = 1 ис¬поль¬¬¬¬зуются формулы (5.27), (5.28). При этом значения ко¬эф¬¬ф謬циентов ai, bi, ci передаются в процедуру в виде мас¬си¬вов из вызывающей программы.
Формальные параметры процедуры. Входные: n (тип in¬¬¬¬teger) - значение порядка v функции (0 или 1); x (тип re¬¬al) - значение аргумента; a[1:6], b[0:6], c[0:6] (тип do¬u¬ble) - зна-чения коэффициентов разложений (5.27), (5.28). Вы¬¬ход¬ной: jnx01 (идентификатор процедуры-функ¬ции, тип do¬ub¬le) - вычисленное значение функции Бес¬селя или .
FUNCTION JNX01(N : INTEGER; X : DOUBLE;
A,B,C : ARRAY OF DOUBLE) : DOUBLE;
VAR Y,S1,S2 : DOUBLE;
BEGIN
X1:=DOUBLE(X); S1:=0
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
