Текст подпрограммы и версий ( Фортран )
sf68r.zip  		Тексты тестовых примеров ( Фортран )
tsf68r.zip 
Текст подпрограммы и версий ( Си )
sf68r_c.zip  		Тексты тестовых примеров ( Си )
tsf68r_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf68r_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf68r_p.zip

Функция:  SF68R

Назначение

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

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

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

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

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

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

    SUBROUTINE  SF68R (X, R, T1, T2) 
   Параметры
X - заданное значение аргумента  x (тип: вещественный);
R - заданное значение мнимой части индекса  β (тип: вещественный);
T1 - вещественная переменная, значение которой полагается равным результату вычисления вещественной части  K1/2  + i β (x);
T2 - вещественная переменная, значение которой полагается равным результату вычисления мнимой части   K1/2  + i β (x). 
   Версии:  нет 
   Вызываемые подпрограммы
SF23C - подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
QSF3R - подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования. 
   Замечания по использованию
  Для подпрограммы SF68R областями допустимых значений параметров  X и  β являются  0.1 ≤ X ≤ 10.,   0.01 ≤ β ≤ 10.

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

         X = 1.
         R = 1.
         CALL  SF68R (X, R, T1, T2)

Результаты:

       T1  =  0.298824989 , 
       T2  =  0.118944694