Текст подпрограммы и версий ( Фортран ) sfa6r.zip sfa6d.zip |
Тексты тестовых примеров ( Фортран ) tsfa6r.zip tsfa6d.zip |
Текст подпрограммы и версий ( Си ) sfa6r_c.zip sfa6d_c.zip |
Тексты тестовых примеров ( Си ) tsfa6r_c.zip tsfa6d_c.zip |
Текст подпрограммы и версий ( Паскаль ) sfa6r_p.zip sfa6e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsfa6r_p.zip tsfa6e_p.zip |
Вычисление функции Бесселя третьего рода (функции Ханкеля) Hn (x) целых порядков от вещественных аргументов.
Подпрограмма SFA6R вычисляет комплексно - значную функцию Бесселя третьего рода (функцию Ханкеля)
Hn(x) = Jn(x) ± i Nn(x) = H1 ± i H2
для n ≥ 0 и вещественных значений x . Здесь Jn (x) и Nn (x) - функции Бесселя первого и второго родов соответственно. Для вычисления Hn (x) используется разложение этой функции в ряд.
Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979.
SUBROUTINE SFA6R (N, X, H1, H2)
Параметры
N - | заданный порядок функции Hn (x), N ≥ 0 (тип: целый); |
X - | заданное значение аргумента x (тип: вещественный); |
H1, H2 - | вещественные переменные, значения которых на выходе полагаются равными вещественной и мнимой части функции Hn (x) соответственно. |
Версии
SFA6D - | вычисление функции Бесселя третьего рода (функции Ханкеля) Hn (x) целых порядков от вещественных аргументов в режиме удвоенной точности; при этом параметры X, H1 и H2 должны иметь тип DOUBLE PRECISION. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
N = 3 X = 4.0 CALL SFA6R (N, X, H1, H2) Результаты: H1 = 0.4301715, H2 = -0.1820221