Вторая пачка часть 139
/I;
END;
END;
D := 0.;
IF(N <>0) THEN
D := R/NA
R := 1./R;
P := EXP (NA * LN(X2));
T := LN (X2)+1.1544313298031D0;
I := 0;
WHILE (I > NA ) AND (B = YNX) DO
BEGIN
B := YNX;
YNX := YNX+A*R*(T-S);
IF (I < NA) THEN
YNX := YNX-A*D/P;
I1 := I+1;
A := A*X2/I1;
R := -R/(I1+NA);
S := S+1./I1+1./(I1+NA);
IF (I1 < NA) THEN
D := D/(NA-I1);
I := I1;
END;
P := EXP ( NA * LN(X1));
YNX := 0.318309886184D0*YNX*P;
IF(N < 0) THEN
YNX := (-1)**NA*YNX;
END.
Процедура-функция YNX тестировалась на IBM PC/AT-386 при следующих значениях входных па¬ра¬мет¬¬¬ров: n = 1, x = ?8 и n = ?7, x = 4. Полученные при этом результаты
N1 (8) = -0.1580554, N1 (-8) = 0.1580554,
N7 (4) = -3.706224, N-7 (4) = 3.706224
совпадают с табличными [Справочник ..., 1979] с точ¬нос¬тью до шести - восьми десятичных знаков.
Вычисление функции Вебера Yv(x) в част-ных слу¬ча¬ях v=n= 0, 1 выполняется в процедуре-функ¬ции YNX01 на основании следующих полиномиальных ап¬прок¬си¬ма¬ций:
a) при 0 ? x ? 3
, (5.30)
где и y определяются выражениями (5.27), зна¬че¬ния ai представлены в табл. 5.12; |? (x) | < 1.4?10-8 - при n=0 и |? (x) | << 11?10-7 при n = 1;
б) при x > 3
, (5.31)
где ? и S определяются так же, как и в вы¬ра¬же-нии (5.28).
Таблица 5.12
i ai
n = 0 n = 1
0 0.36746691 -0.6366198
1 0.60559366 0.2212091
2 -0.74350384 2.1682709
3 0.25300117 -1.3164827
4 -0.04261214 0.3123951
5 0.00427916 -0.0400976
6 -0.00024846 0.0027873
Формальные параметры процедуры. Входные: n (тип in¬teger) - значение порядка n функции (0 или 1); x1 (тип re¬al) - значение аргумента; a[0:6], b[0:6], c[0:6] (тип do¬u¬b¬le) - значения коэффициентов разложений (5.28), (5.30). Вы¬ходной: ynx01 (идентификатор процедуры-функции, тип do¬uble) - вычисленное значение функции В嬬бера или . В процедуре YNX01 используется внеш¬няя про-цедура-функция JNX.
FUNCTION YNX01(N : INTEGER; X1 : DOUBLE;
A,B,C : ARRAY OF DOUBLE) : DOUBLE;
VAR X,PI,Y,S1,S2 : DOUBLE;
BEGIN
{ *** ИСПОЛЬЗУЕТСЯ ВНЕШНЯЯ ФУНКЦИЯ JNX ***}
PI := 3
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа