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

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

Назначение

Вычисление логарифмической производной ψ (1 + ix) гамма - функции.

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

sf65r_c вычисляет логарифмическую производную гамма - функции (пси - функцию) ψ (1 + ix).

Подпрограмма  sf65r_c представлена в Библиотеку Вычислительным центром AH CCCP (группа Л.Н.Кармазиной) в виде процедуры на АЛГОЛе - 60. Перевод этой процедуры на стандарт языка ФОРТРАН (а затем на Си) осуществлен в НИВЦ МГУ. Алгоритм взят из кандидатской диссертации:

Г.Неметх, "Разложение обобщенных гипергеометрических функций по полиномам Чебышева", ОИЯИ, Дубна, 1971.

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

    void sf65r_c (complex * z, real *x)

   Параметры
x - заданное значение аргумента  x (тип: вещественный);
z - результат вычисления (тип: комплексный). 
   Версии:  нет 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

int main(void)
{
    /* System generated locals */
    complex q__1;

    /* Local variables */
    extern /* Complex */ void sf65r_c(complex *, float *);
    static complex z;
    static float x;

    x = 2.f;
    sf65r_c(&q__1, &x);
    z.r = q__1.r, z.i = q__1.i;
/* l1: */
    printf("\n %16.7e \n",x);
    printf("\n z = %16.7e %16.7e \n",z.r, z.i);
    return 0;
} /* main */


Результат:    z = (0.7145915153687, 1.320807282642)