Текст подпрограммы и версий ( Си )
sf66r_c.zip  sf66d_c.zip  		Тексты тестовых примеров ( Си )
tsf66r_c.zip  tsf66d_c.zip 

Подпрограмма:  sf66r_c

Назначение

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

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

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

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

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

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

    int sf66r_c (real *x, real *alngx, integer *ierr)

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

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

int main(void)
{
    /* Local variables */
    extern int sf66r_c(float *, float *, int *);
    static int ierr;
    static float y;
    static float c_b1 = .5f;

    sf66r_c(&c_b1, &y, &ierr);

    printf("\n %16.7e %5i \n",y,ierr);
    return 0;
} /* main */


Результаты:

       y  =  0.572 364 942 9228
       ierr  =  0