|
Текст подпрограммы и версий ( Фортран ) sf32r.zip sf32d.zip |
Тексты тестовых примеров ( Фортран ) tsf32r.zip tsf32d.zip |
|
Текст подпрограммы и версий ( Си ) sf32r_c.zip sf32d_c.zip |
Тексты тестовых примеров ( Си ) tsf32r_c.zip tsf32d_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) sf32r_p.zip , sf32e_p.zip |
Тексты тестовых примеров ( Паскаль ) tsf32r_p.zip , tsf32e_p.zip |
Вычисление сферической функции Бесселя первого рода целых порядков jn (x).
SF32R вычисляет сферическую функцию Бесселя первого рода целых порядков для вещественных положительных значений аргумента
jn (x) = (π / 2x) 1/2 Jn + 1/2 (x) ,
где Jn + 1/2 (x) - функция Бесселя первого рода и 0 ≤ n ≤ 29.
M.A.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, Nat. Bur. Stand. Washington D.C., 1964.
SUBROUTINE SF32R (X, N, Y, IERR)
Параметры
| X - | заданное значение аргумента x (тип: вещественный); |
| N - | заданное значение порядка n (тип: целый); |
| Y - | вещественный вектоp длины N + 1, в который в pезультате работы подпрограммы помещаются вычисленные значения j0 (X), ..., jN (X); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы подпрограммы; при этом: |
| IERR=65 - | когда значение аргумента отрицательно; |
| IERR=66 - | когда заданный порядок сферической функции Бесселя не принадлежит отрезку (0, 29). |
Версии
| SF32D - | вычисление сферической функции Бесселя первого рода целых порядков jn (x) с повышенной точностью. |
Вызываемые подпрограммы
| UTSF10 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы SF32R. |
| UTSF11 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы SF32D. |
Замечания по использованию
| Для подпрограммы SF32D параметры X и Y должны иметь тип DOUBLE PRECISION. |
DIMENSION Y(3)
X = 0.1
N = 2
CALL SF32R (X, N, Y, IERR)
Результаты:
Y(1) = 0.998334166468,
Y(2) = 0.033300011902,
Y(3) = 0.000666190608,
IERR = 0