Вторая пачка часть 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
Индекс
Элементарные функции
Линейные уравнения
Нелинейные уравнения
Случайные числа
