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

5772156649; EXIT; END; GLOGD:=0; PI:=3.1415926535; IF X<-2 THEN BEGIN GLOGD:=-PI*COS(PI*X)/SIN(PI*X); X:=1-X END; WHILE X<15 DO BEGIN GLOGD:=GLOGD-1/X; X:=X+1; END; X:=1/(X-1); GLOGD:=GLOGD-LN(X)+X/2; X:=X*2; GLOGD:=GLOGD-X*((X/21-0.1)*X+1)/12 END {****GLOGD****}. Процедура-функция GLOGD была получена путем пе¬ре¬¬ра¬¬ботки и перевода на язык PASCAL FOR¬TRAN-про¬¬грам¬¬мы вычисления ?(x) [Белашов, 1997] и про¬тес¬ти¬рована на IBM PC/AT-286. При этом ре¬зуль¬та¬ты, по¬¬лученные для x= = 0.5, 1, ?(0.5) = -1.96351003; ?(1) = -0.577215665 совпадают с точностью 10-8 с результатами, при-ве¬ден¬ны¬¬ми в работе Гринчишина и др. (1988). Замечание. Для процедуры GLOGD справедливо то же замечание, что и для процедуры GAMF1. Вычисление неполной гамма-функции Гх(а) с ис¬поль¬зованием ее разложения в степенной ряд (5.8) реализовано в процедуре-функции UGAMF. При этом вначале исключаются точки a = 0, -1, -2, ... , в которых функция не определена. Точность вычисления Гх(а) за¬да¬ется при вызове функции. Формальные параметры процедуры. Входные: x, a (тип real) - амплитуда и аргумент функции Гх(а); eps (тип real) - точность вычисления. Выходной: ugamf (тип double) - значение неполной гамма-функции с точ¬нос¬тью eps. FINCTION UGAMF (X,A,EPS : REAL) ; REAL; VAR S, S1: REAL I : INTEGER; LABEL A1; BEGIN IF INT(A)=-ABS(A) THEN BEGIN WRITE (“ОШИБКА В ИСХОДНЫХ ДАННЫХ ‘); EXIT; END; UGAMF:=1/A; I:=1; S:=1; REPEAT S:=-S*X/I; S1:=S/(A+I); I:=I+1; UGAMF:=UGAMF+S1; UNTIL ABS(S1) <= EPS*ABS(UGAMF); UGAMF:=UGAMF*X A END { *** UGAMF *** }. Процедура-функция UGAMF получена путем пе¬¬¬ре¬ра¬¬¬ботки и перевода на язык PASCAL про¬грам¬мы вы¬чис¬ления неполной гамма-функции Гх(а) [Бе¬ла¬шов, 1997] и про¬тестирована на IBM PC/AT-286 для тех же зна¬че¬ний x и a. Результат вычисления Г2(3) = 0.646647167 при за¬да¬вав¬шей¬ся точности ? = 1e-8 совпадает с конт¬роль¬ным с точ¬нос¬тью до вось¬¬ми десятичных знаков

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