Текст подпрограммы и версий ( Фортран )
sf75r.zip  		Тексты тестовых примеров ( Фортран )
tsf75r.zip 
Текст подпрограммы и версий ( Си )
sf75r_c.zip  		Тексты тестовых примеров ( Си )
tsf75r_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf75r_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf75r_p.zip

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

Назначение

Вычисление модифицированных функций Бесселя второго рода   Кn (x) для последовательности целых индексов.

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

SF75R вычисляет значения модифицированных функций Бесселя второго рода  Kp(x) для действительного аргумента  x и последовательности целых индексов  p от 0 до  n. Расчеты осуществляются по рекуррентной формуле: 

     Kp+1(x) = (2p/x)*Kp(x) + Kp-1(x) .

Ошибка вычисления  Kp (x) не pастет в pекуppентном процессе, примененном для возрастающих значений  p. Поэтому сначала вычисляются с заданной точностью значения  K0 (x) и  K1 (x) с помощью разложений по полиномам Чебышева, а затем последовательно применяется рекуррентная формула для  p = 2, 3, ..., n.

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

    SUBROUTINE  SF75R (X, N1, S, T) 
   Параметры
X - заданное значение аргумента  x (тип: вещественный);
N1 - максимальный порядок последовательности функций   Kp (x), увеличенный на единицу (тип: целый);
S - логическая переменная, задающая режим работы подпрограммы; при этом:
если S=.TRUE. , то вычисляются функции  Kp (x),
если S=.FALSE. , то вычисляются функции  ex Kp (x).
T - вещественный вектоp длины N1 вычисленных значений функции Бесселя порядков от  0 до  n;    Kj (x) = T (j + 1),  j = 0, ..., N.  (N = N1 - 1). 
   Версии:  нет 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

      DIMENSION  T(21)
      LOGICAL S
      X = 5.
      N1 = 21
      S = .TRUE.
      CALL  SF75R (X, N1, S, T, )

Результаты:

       T(6)    =  3.2706273712E-02
       T(11)  =  9.7585628291E+00
       T(16)  =  3.0169766300E+04
       T(21)  =  4.8270005205E+08