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

Функция:  SF67R (модуль SF67R_p)

Назначение

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

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

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

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

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

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

procedure SF67R(XOCH :Real; ROCH :Real; var T :Real);


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

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

Unit TSF67R_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF67R_p;

function TSF67R: String;

implementation

function TSF67R: String;
var
T :Real;
const
ХОСН :Real = 4.0;
RОСН :Real = 3.0;
begin
Result := '';  { результат функции }
SF67R(XOCH,ROCH,T);
Result := Result + Format('     %20.16f ',[T]) + #$0D#$0A;
UtRes('TSF67R',Result);  { вывод результатов в файл TSF67R.res }
exit;
end;

end.

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