# Вторая пачка часть 138

; IF (X =0) AND (N=0) THEN EXIT; NA:=ABS(N); XA:=ABS(X); IF (XA > 3) THEN JANX:=JNX01(NA,XA,A,B,C) ELSE IF (X =0) THEN BEGIN JANX:=0.; EXIT; END ELSE JANX:=JNX(NA,XA,EPS) END; IF(N=0) THEN EXIT; IF(N < 0) THEN IF (NA MOD 2 = 1) THEN JANX := -JANX; IF(X < 0.) THEN IF (NA MOD 2 = 1) THEN JANX := -JANX; END. Тестирование процедуры-функции JANX на IBM PC/AT-386 при тех же значениях n и x, что и для про¬це¬дур JNX и JNX01, дало аналогичные значения ; ре¬зуль¬та¬ты при отрицательных значениях порядка и ар¬гу¬мента пред¬ставлены в табл. 5.11 (ср. с табл. 5.10). Таблица 5.11 x n = -1 n = 1 2.9 -0.3754275162-3.44e-8 0.3754275162+3.44e-8 -2.9 0.3754275162+3.44e-8 -0.3754275162-3.44e-8 4 0.0660433224+0.56e-8 -0.0660433224-0.56e-8 -4 -0.0660433224-0.56e-8 0.0660433224+0.56e-8 15 -0.2051040490+1.04e-8 0.2051040490-1.04e-8 -15 0.2051040490-1.04e-8 -0.2051040490+1.04e-8 Вычисление функции Бесселя II рода для любого це¬ло¬го v= = n и вещественного аргумента x может быть вы¬пол¬нено на основании ее разложения [Справочник ..., 1979] (5.29) где C = 0.57721566490153 - постоянная Эйлера. Такой алгоритм реализован в приведенной далее процедуре-функции ynx, использующей константы 2C= =1.1544... и 1/? = 0.318... . При этом, если n < 0, вначале вы-числяется функция , а затем используется соотношение . Формальные параметры процедуры. Входные: n (тип in¬¬teger) - порядок n функции; x (тип real) - значение ар¬гу¬¬¬мен¬та. Выходной: ynx (идентификатор процедуры-функ¬¬ции, тип double) - значение . При x = 0, ког¬да функ-ция Неймана (Вебера) не определена, осу¬щес¬тв-ля¬¬¬ет¬ся вы¬ход к программе обработки соответствующей ошиб¬ки. FUNCTION YNX (N : INTEGER; X : DOUBLE) : DOUBLE; VAR X1,X2,A,R,B,S,D,P,T : DOUBLE; BEGIN IF(X = 0.) THEN EXIT; YNX := 0.; NA := IABS(N); X1 := DOUBLE(X/2.); X2 := X1*X1; A := 1.; R := 1.; B := 0.; S := 0.; IF(N <>0)THEN FOR I : := 1 TO NA DO BEGIN R := R*I; S := S+1

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