БЧА НИВЦ МГУ. SF68R. Цилиндрические функции

Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.srcc.msu.ru/num_anal/lib_na/cat/sf/sf68r.htm
Дата изменения: Mon Mar 14 11:07:50 2011
Дата индексирования: Mon Feb 4 02:16:15 2013
Кодировка: Windows-1251

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

Функция: SF68R

Назначение

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

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

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

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

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

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

    SUBROUTINE  SF68R (X, R, T1, T2)

   Параметры

X -	заданное значение аргумента x (тип: вещественный);
R -	заданное значение мнимой части индекса b (тип: вещественный);
T1 -	вещественная переменная, значение которой полагается равным результату вычисления вещественной части K_{1/2 + i b} (x);
T2 -	вещественная переменная, значение которой полагается равным результату вычисления мнимой части K_{1/2 + i b} (x).

   Версии:  нет

   Вызываемые подпрограммы

SF23C -	подпрограмма вычисления значений гамма - функции и ее натурального логарифма от комплексного аpгумента.
QSF3R -	подпрограмма вычисления определенных интегралов от осциллирующих функций по формулам интерполяционного типа пятой степени точности с автоматическим выбором шага для больших отрезков интегрирования.

   Замечания по использованию

Для подпрограммы SF68R областями допустимых значений параметров X и b являются 0.1 ? X ? 10., 0.01 ? b ? 10.

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

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

Результаты:

       T1  =  0.298824989 , 
       T2  =  0.118944694