Текст подпрограммы и версий ( Фортран ) sfa2r.zip sfa2d.zip |
Тексты тестовых примеров ( Фортран ) tsfa2r.zip tsfa2d.zip |
Текст подпрограммы и версий ( Си ) sfa2r_c.zip sfa2d_c.zip |
Тексты тестовых примеров ( Си ) tsfa2r_c.zip tsfa2d_c.zip |
Текст подпрограммы и версий ( Паскаль ) sfa2r_p.zip sfa2e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsfa2r_p.zip tsfa2e_p.zip |
Вычисление функций Риккати - Бесселя первого рода Sn (x) и второго рода Cn (x) от вещественных аргументов.
Подпрограмма SFA2R вычисляет для вещественных аргументов x функции Риккати - Бесселя первого рода
Sn(x) = (πx/2)1/2 Jn+1/2(x) = x jn(x) , n = 0, 1, 2, ..., N и второго рода Cn(x) = (πx/2)1/2 Yn+1/2(x) = x yn(x) , n = 0, 1, 2, ..., N , где Jn + 1/2 и Yn + 1/2 - функции Бесселя первого и второго рода, jn и yn - сферические функции Бесселя первого и второго рода
Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979
SUBROUTINE SFA2R (X, N, S, C)
Параметры
X - | заданное значение аргумента (тип: вещественный); |
N - | заданное значение максимального порядка, для которого необходимо вычислить функции Sn (x) и Cn (x), n = 0, 1, ..., N (тип: целый); |
S, C - | вещественные одномерные массивы длины N + 1, в которых размещаются вычисленные значения S0 (x), ..., SN (x) и C0 (x), ..., CN (x) соответственно. |
Версии
SFA2D - | вычисление функций Риккати - Бесселя первого рода Sn (x) и второго рода Cn (x) от вещественных аргументов в режиме удвоенной точности; при этом параметры X, S и C должны иметь тип DOUBLE PRECISION. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
DIMENSION S(3), C(3) X = 3 N = 2 CALL SFA2R (X, N, S, C) Результаты: S = 0.141120, 1.03703, 0.895913 C = -0.989992, -0.188878, 0.801115