Текст подпрограммы и версий
sf66r_p.zip , sf66e_p.zip 		Тексты тестовых примеров
tsf66r_p.zip , tsf66e_p.zip

Подпрограмма:  SF66R (модуль SF66R_p)

Назначение

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

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

SF66R вычисляет значение  ln Г(x), где 

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

Справочник по специальным функциям. Под ред. M.Абрамовица, И.Стеган. M.: "Hаука", 1979.

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

procedure SF66R(X :Real; var ALNGX :Real; var IERR :Integer);


   Параметры
X - заданное значение аргумента  x (тип: вещественный);
ALNGX - вещественная переменная, содержащая вычисленное значение натурального логарифма  Г(x);
IERR - целая переменная, служащая для сообщения об ошибках; при этом
IERR=67 - когда  x ≤ 0;  значение  ln Г(x) полагается равным 1. 
   Версии
SF66E - вычисление натурального логарифма гамма - функции вещественного аргумента с расширенной (Extended) точностью. 
   Вызываемые подпрограммы
UTSF16 - подпрограмма выдачи диагностических сообщений при работе подпрограммы SF66R. 
   Замечания по использованию:  нет

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

Unit tsf66r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF66R_p;

function tsf66r: String;

implementation

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

end.


Результаты:

       Y  =  0.572 364 942 9228
        I  =  0