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

Функция:  sf67r_c

Назначение

Вычисление модифицированной функции Бесселя второго рода чисто мнимого индекса и вещественного аргумента   Ki τ (x).

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

sf67r_c вычисляет модифицированную функцию Бесселя мнимого индекса  Кi τ (x).

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

1.  Журина М.И., Кармазина Л.Н. Таблицы модифицированных функций Бесселя мнимого индекса. M., ВЦ AH CCCP, 1967.
2.  Кармазина Л.Н., Раппопорт Ю.М. Вычисление модифицированной функции Бесселя  Ki β (x), "Алгоритмы и программы", информ. бюлл. ВНТИЦентра, 1979., N 3, П003594.

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

    int sf67r_c (real *xoch, real *roch, real *t)

   Параметры
xoch - заданное значение аргумента  x (тип: вещественный);
roch - заданное значение индекса  τ (тип: вещественный);
t - вещественная переменная, значение которой полагается равным результату вычисления   Ki τ (x). 
   Версии:  нет 
   Вызываемые подпрограммы
sf23c_c - подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
qsf3r_c - подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования. 
   Замечания по использованию
  Для подпрограммы sf67r_c областями допустимых значений параметров  x и  τ являются 0.1 ≤ x ≤ 10,   0.01 ≤ τ ≤ 10.

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

int main(void)
{
    /* Initialized data */
    static float xoch = 4.f;
    static float roch = 3.f;

    /* Local variables */
    extern int sf67r_c(float *, float *, float *);
    static float t;

    sf67r_c(&xoch, &roch, &t);

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


Результат:    t  =  0.0039263859478