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

3), реализован в процедуре-функции GAMF1, зна-чения формальных параметров в которой те же, что и в процедуре-функции GAMF. FUNCTION GAMF1(X: REAL;K: INTEGER) : REAL; BEGIN GAMF1:=1; IF (K<> -1) AND INT(X)=-ABS(X) THEN BEGIN WRITE (“ОШИБКА В ИСХОДНЫХ ДАННЫХ ‘); EXIT; END; WHILE X<15 DO BEGIN GAMF1:=GAMF1/X; X:=X+1; END; GAMF1:=GAMF1*SQRT(2*3.141592655/X) *EXP(-X+X*LN(X)); X:=1/X; GAMF1:=GAMF1*(1+X*(1/12+X*(1/288- X*(1339/51840+X*571/2488320)))); IF K <>-1 THEN EXIT; GAMF1:=1/GAMF1 END {****GAMF1****}. Процедура-функция GAMF1 была получена путем пе¬ре¬работки и перевода Бейсик-программы вы¬чис¬ле¬ния гам¬ма-функции, приведенной в работе Гринчишина и др. (1988), на язык FORTRAN [Белашов, 1997], а затем на язык PAS¬CAL и протестирована на IBM PC/AT-286. По¬лу¬чен¬ные при этом результаты совпадают с точ¬нос¬тью 10-8 с ре¬зуль¬татами вычислений процедуры-функ¬ции GAMF, при¬ве-денными в табл. 5.2. Замечание. Процедуру-функцию GAMF1 ре-ко¬мен¬ду¬ет¬¬ся использовать для значений аргумента х ? 15, од¬на¬ко ее можно с успехом применять и для малых x, так как в этом случае предусматривается приведение зна¬че¬ния ар¬¬-гумента к x = 15 в соответствие с со¬от¬но¬ше¬ни¬ем (5.1), на что, правда, расходуется дополнительное про¬¬цес¬сор¬¬ное время. Вычисление логарифмической производной гамма-функ¬ции ? (х) с использованием формул (5.5) и асим¬п¬то¬¬ти¬чес¬ко¬го разложения (5.6) реализовано в процедуре-фун¬к¬ции GLOGD. Исключение полюсов функции ? (х) при х = 0, -1, -2, ... производится так же, как в про¬це¬ду¬рах-функциях GAMF и GAMF1. Формальные параметры процедуры. Входной: x (тип do¬¬uble) - значение аргумента x. Выходной: glogd (тип do¬u¬b¬le) - значение функции ? (x). Как и в процедурах-функ¬ци¬ях GAMF и GAMF1, при x = 0, -1, -2, ... , когда ? (x) имеет по¬люсы, осуществляется выход к внешней процедуре об¬ра¬ботки ошибки. FUNCTION GLOGD(X: DOUBLE) : DOUBLE; VAR PI : REAL; BEGIN IF INT(X)=-ABS(X) THEN BEGIN WRITE (“ОШИБКА В ИСХОДНЫХ ДАННЫХ ‘); EXIT; END; IF X=1 THEN BEGIN GLOGD:=-0

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