Текст подпрограммы и версий
sf21r_p.zip 		Тексты тестовых примеров
tsf21r_p.zip

Функция:  SF21R (модуль SF21R_p)

Назначение

Вычисление значений гамма-функции вещественного аргумента.

Математическое описание

Функция SF21R вычисляет значения гамма-функции 

                   ∞
     Г (x)  =  ∫ t x - 1 e - t dt
                   0

для вещественных значений аргумента  x.

J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.

Использование 

function SF21R(X :Real; var IERR :Integer): Real;


   Параметры
X - заданное значение аргумента x (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда значение аргумента pавно нулю или целому отрицательному числу; значение функции полагается равным 3.4E38;
IERR=66 - если значение аргумента больше чем 21.43; значение функции полагается равным 3.4E38. 
   Версии:  нет 
   Вызываемые подпрограммы
UTSF10 - подпрограмма выдачи диагностических сообщений при работе функции SF21R. 
   Замечания по использованию:  нет

Пример использования 

Unit tsf21r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF21R_p;

function tsf21r: String;

implementation

function tsf21r: String;
var
IERR :Integer;
X,Y :Real;
begin
Result := '';  { результат функции }
X := 1.25;
Y := SF21R(X,IERR);
Result := Result + Format(' %20.16f ',[Y]) + #$0D#$0A;
Result := Result + Format(' %5d ',[IERR]) + #$0D#$0A;
UtRes('tsf21r',Result);  { вывод результатов в файл tsf21r.res }
exit;
end;

end.


Результаты:

       Y  =  0.906402477054,    IERR  =  0