Текст подпрограммы и версий ( Фортран )
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 

Подпрограмма:  SFA2R

Назначение

Вычисление функций  Риккати - Бесселя первого рода  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