Вторая пачка часть 177
/X;
G := SQRT(1.+(JINX1*X)* (JINX1*X))/X;
END;
IF(N <>0)THEN
FOR I := 1 TO N DO
BEGIN
C := JINX1;
JINX1 := K*(G-JINX1*(2*I-1)/X);
G := C;
END;
END.
Процедура-функция JINX1 тестировалась на IBM PC/AT-286 для n = 2, 3; x = 1, 3. Полученные ре-зуль¬та¬ты, при¬веденные в табл. 5.18, совпадают с табличными зна¬че¬ни¬ями функций [Справочник ..., 1979] и ре¬зуль¬та¬та¬ми конт¬рольных примеров, приведенными в работах [Бе¬ла¬шов, 1997; Гринчишин и др., 1988].
Таблица 5.18
n jn (x) in (x)
x = 1 x = 3 x = 1 x = 3
2 6.2035052e-22 0.29863749 7.1562870e-22 1.09650152
3 9.0065812e-32 0.15205166 1.0065091e-22 0.41528758
В процедуре-функцие YINX1 в зависимости от спо¬со¬ба ее вызова выполняется вычисление сфе¬ри¬чес¬кой или модифицированной сферической функ¬ций Бес¬селя II рода в соответствии с со¬от¬но¬ше¬ни¬ем (5.50) для произвольного целого n ? 0. Значения фор¬маль¬¬ных параметров и структура процедуры те же, что и для JINX1.
FUNCTION YINX1(N,K : INTEGER; X1 : DOUBLE) : DOUBLE;
VAR X,G,C : DOUBLE;
BEGIN IF(X1 = 0.) THEN EXIT;
X := DOUBLE(X1);
IF(K = -1) THEN
BEGIN
YINX1 := -COS(X)/X; G := SIN(X)/X;
END
ELSE
BEGIN
YINX1 := (EXP(X)+EXP(-X))/2./X ;
G := SQRT((YINX1*X)* (YINX1*X)-1)/X
END;
IF(N <> 0)THEN
FOR I := 1 TO N DO
BEGIN C := YINX1;
YINX1 := K*(G-YINX1*(2*I-1)/X);
G := C ;
END ;
END.
Процедура-функция YINX1 тестировалась на IBM PC/AT-286 для n = 2, 3; x = 2, 3. Полученные ре-зуль¬та¬ты, приведенные в табл. 5.19, совпадают с табличными зна¬чениями функций [Справочник ..., 1979].
Таблица 5.19
x yn (x) i-n (x)
n=2 n=3 n=2 n=3
2 -0.73399142 -1.48436655 0.57177592 -0.55655852
3 -0.26703834 -0.50802305 1.13522480 0.32862120
С помощью процедуры-функции JINX2 вычисляется в за¬висимости от способа ее вызова значение функции или в соответствии с их разложениями (5.51) для .
Формальные параметры про¬це¬ду¬ры. Входные: n (тип in¬¬teger) - порядок; x (тип real) - аргумент функции; k (тип in¬teger) - параметр, значение которого определяет вы¬¬бор вы¬числяемой функции: при k = -1 вычисляется , при k = 1 -
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
