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

Функция:  sf68r_c

Назначение

Вычисление вещественной и мнимой части модифицированной функции Бесселя второго рода мнимого индекса и вещественного аргумента  K1/2  + i β (x).

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

sf68r_c вычисляет вещественную и мнимую части модифицированной функции Бесселя  K1/2  + i β (x).

B зависимости от значений аргумента и индекса используется метод, основанный на степенных разложениях, на интегральных представлениях или  τ - метод Ланцоша.

Раппопорт Ю.М. Таблицы модифицированных функций Бесселя   K1/2  + i β (x), Hаука, 1979.

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

    int sf68r_c (real *x, real *r, real *t1, real *t2)

   Параметры
x - заданное значение аргумента  x (тип: вещественный);
r - заданное значение мнимой части индекса  β (тип: вещественный);
t1 - вещественная переменная, значение которой полагается равным результату вычисления вещественной части  K1/2  + i β (x);
t2 - вещественная переменная, значение которой полагается равным результату вычисления мнимой части   K1/2  + i β (x). 
   Версии:  нет 
   Вызываемые подпрограммы
sf23c_c - подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
qsf3r_c - подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования. 
   Замечания по использованию
  Для подпрограммы sf68r_c областями допустимых значений параметров  x и  β являются  0.1 ≤ x ≤ 10.,   0.01 ≤ β ≤ 10.

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

int main(void)
{
    /* Initialized data */
    static float x = 1.f;
    static float r__ = 1.f;

    /* Local variables */
    extern int sf68r_c(float *, float *, float *, float *);
    static float t1, t2;

    sf68r_c(&x, &r__, &t1, &t2);

    printf("\n %16.7e %16.7e \n",t1,t2);
    return 0;
} /* main */


Результаты:

       t1  =  0.298824989 , 
       t2  =  0.118944694