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

35619449019D0); F1 := DSIN(N1*2.35619449019D0); END; N1 := 1; FOR I := 1 TO NA DO N1 := N1*I; X1 := X1*X1; F1 := F1/N1; F2 := F2/N1/(NA+1.); N1 := 1; I := 2; BENX := F1+K*F2; REPEAT F1 := -X1*F1/(I*(I-1)*(NA+I-1)*(NA+I)); F2 := -X1*F2/(I*(I+1)*(NA+I+1)*(NA+I)); F := BENX+F1+K*F2; IF (F<>BENX)THEN BEGIN BENX := F; INC (I, 2); END; UNTIL F = BENX; BENX := BENX*EXP (LN(X0/2.)*NA); IF(N < 0) THEN IF NA MOD 2 = 1 THEN BENX := (-1)*BENX; END. Процедура-функция BENX была получена путем не¬ко¬то¬¬рой модификации и перевода вначале на язык FORTRAN [Белашов, 1997], а затем на PASCAL Бей¬сик-программ вычисления bern(x) и bein(x), опуб¬ли¬ко¬ван¬ных в работе Гринчишина и др. (1988), и про¬тес¬ти¬ро¬ва¬на на IBM PC/AT-286 для n = 0, 2; x = 3. Полученные ре¬зультаты ber0 (3) = -0.221380250, ber2 (3) = 0.808368465, bei0 (3) = 1.937586785, bei2 (3) = -0.891022363 совпадают с табличными значениями [Справочник ..., 1979] и результатами контрольных примеров, при¬ве¬ден¬ны¬ми в работе Гринчишина и др. (1988). В процедуре-функции BE0X в зависимости от спо¬со¬ба об¬¬ращения реализовано вычисление функций ber(x) или bei(x) на основании их разложений в ряды (5.45) с точ-нос¬тью ?, задаваемой при вызове про¬це¬ду¬ры. Вы¬чис¬ления пре¬кращаются при выполнении не-ра¬вен¬ства , где m = 2k для функции ber(x) и m = 2k + 1 - для функ¬ции bei(x). Формальные параметры процедуры. Входные: x (тип double) - значение аргумента; eps (тип double) - за¬да¬ва¬е¬мая точность вычислений; n (тип integer) - параметр, по зна-чению которого выбирается вычисляемая функция: при n = 0 вычисляется функция ber(x), при n = 1 - bei(x). Выходной: be0x (идентификатор процедуры-функции, тип double) - вычисленное значение функции. FUNCTION BE0X (X1,EPS : DOUBLE;N : INTEGER) : DOUBLE; VAR F,F1,X : DOUBLE; BEGIN X := DOUBLE(X1*X1/4.); BE0X := 1.; I := 1+N; K := 1; F := 1.; REPEAT F := F/(I*I*(I+1)*(I+1)); F1 := EXP ( LN(X)*(2*K))*F; IF K MOD 2 = 1 THEN F1 := - F1; BE0X := BE0X+F1; I := I+2; K := K+1; UNTIL (ABS(F1) <=EPS); IF(N

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