Текст подпрограммы и версий ( Фортран )
sf28r.zip  sf28d.zip  		Тексты тестовых примеров ( Фортран )
tsf28r.zip  tsf28d.zip 
Текст подпрограммы и версий ( Си )
sf28r_c.zip  sf28d_c.zip  		Тексты тестовых примеров ( Си )
tsf28r_c.zip  tsf28d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf28r_p.zip , sf28e_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf28r_p.zip , tsf28e_p.zip

Функция:  SF28R

Назначение

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

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

SF28R вычисляет модифицированную функцию Бесселя первого рода нулевого порядка  I0 (x) для вещественных  x.

J.F.Hart, E.W.Cheney, C.L.Lawson et al, Computer Approximations, Wiley, New York, 1968.

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

    REAL  FUNCTION  SF28R (X, IERR) 
   Параметры
X - заданное значение аргумента  x (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=1 - когда значение аргумента меньше 0; функция вычисляется для абсолютного значения аргумента;
IERR=65 - когда значение аргумента больше 67.35; значение функции полагается равным 3.4E38. 
   Версии
SF28D - вычисление модифицированной функции Бесселя первого рода нулевого порядка  I0 (x) с повышенной точностью. 
   Вызываемые подпрограммы
UTSF10 - подпрограмма выдачи диагностических сообщений при работе функции SF28R.
UTSF11 - подпрограмма выдачи диагностических сообщений при работе функции SF28D. 
   Замечания по использованию
 

Для функции SF28D параметр X имеет тип DOUBLE PRECISION и его областью допустимых значений является полуинтервал  0 ≤ X < 2070. При IERR = 65 значение функции SF28D полагается равным 1.7D308.

Если значение X меньше нуля, то вычисляется  I0 ( |X| ).

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

         X = 1.
         Y = SF28R (X, IERR)

Результаты:

       Y  =  1.26606587768
       IERR  =  0