# Вторая пачка часть 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

Индекс
Элементарные функции Линейные уравнения Нелинейные уравнения Случайные числа