Текст подпрограммы и версий ( Фортран )
sf67r.zip  		Тексты тестовых примеров ( Фортран )
tsf67r.zip 
Текст подпрограммы и версий ( Си )
sf67r_c.zip  		Тексты тестовых примеров ( Си )
tsf67r_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf67r_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf67r_p.zip

Функция:  SF67R

Назначение

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

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

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

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

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

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

    SUBROUTINE  SF67R (XOCH, ROCH, T) 
   Параметры
XOCH - заданное значение аргумента  x (тип: вещественный);
ROCH - заданное значение индекса  τ (тип: вещественный);
T - вещественная переменная, значение которой полагается равным результату вычисления   Ki τ (x). 
   Версии:  нет 
   Вызываемые подпрограммы
SF23C - подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
QSF3R - подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования. 
   Замечания по использованию
  Для подпрограммы SF67R областями допустимых значений параметров  x и  τ являются 0.1 ≤ x ≤ 10,   0.01 ≤ τ ≤ 10.

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

         XOCH = 4.0
         ROCH = 3.0
         CALL  SF67R (XOCH, ROCH, T)

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