Текст подпрограммы и версий ( Фортран )
sfa5r.zip  sfa5d.zip  		Тексты тестовых примеров ( Фортран )
tsfa5r.zip  tsfa5d.zip 
Текст подпрограммы и версий ( Си )
sfa5r_c.zip  sfa5d_c.zip  		Тексты тестовых примеров ( Си )
tsfa5r_c.zip  tsfa5d_c.zip 
Текст подпрограммы и версий ( Паскаль )
sfa5r_p.zip  sfa5e_p.zip  		Тексты тестовых примеров ( Паскаль )
tsfa5r_p.zip  tsfa5e_p.zip 

Функция:  SFA5R

Назначение

Вычисление экспоненциального интеграла Bn (a).

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

Подпрограмма - функция SFA5R вычисляет значение интеграла 

                                   1
          Bn(a)     =        ∫    t n   e-at   dt ,
                              -1

для  n = 0, 1, 2,...  и вещественного параметра  a . Если  | a | < 1, то используется разложение  e- at  в ряд Тейлора. В противном случае используется рекуррентная формула 

           Bn(a)     =     ( (-1)n  ea   -   e-a   +   n Bn-1(a) ) / a     ,

где      B0(a)     =     ( ea   -   e-a) / a

Абрамовиц М., Стиган И. Справочник по специальным функциям. М.: Наука, 1979.

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

    FUNCTION  SFA5R (N, A) 
   Параметры
N - заданное значение параметра  n  (тип: целый);
A - заданное значение параметра  a  (тип: вещественный). 
   Версии
SFA5D - вычисление экспоненциального интеграла Bn (a) в режиме удвоенной точности; при этом параметр A должен иметь тип DOUBLE PRECISION. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

       N = 3 
       A = 0.25 
       RES = SFA5R (N, A) 

Результат:        RES   =   -0.1007459