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

Функция:  SF68R (модуль SF68R_p)

Назначение

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

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

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

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

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

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

procedure SF68R(X :Real; R :Real; var T1 :Real; var T2 :Real);


   Параметры
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.

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

Unit tsf68r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF68R_p;

function tsf68r: String;

implementation

function tsf68r: String;
var
T1,T2 :Real;
const
X :Real = 1.0;
R :Real = 1.0;
begin
Result := '';  { результат функции }
SF68R(X,R,T1,T2);
Result := Result + Format(' %20.16f %20.16f ',[T1,T2]) + #$0D#$0A;
UtRes('tsf68r',Result);  { вывод результатов в файл tsf68r.res }
exit;
end;

end.


Результаты:

       T1  =  0.298824989 , 
       T2  =  0.118944694