Текст подпрограммы и версий
sfa4r_c.zip  sfa4d_c.zip  		Тексты тестовых примеров
tsfa4r_c.zip  tsfa4d_c.zip 

Функция:  sfa4r_c

Назначение

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

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

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

                         ∞
        An(b)  =   ∫   t n  e-bt  dt
                       1

для  n = 0, 1, 2, ... и положительного вещественного параметра  b по рекуррентной формуле 

        An(b)  =  A0(b)  +  (n / b)  An-1(b)  ,
где  A0 (b)  =  (e - b) / b

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

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

    real sfa4r_c (integer *n, real *b)

   Параметры
n - заданное значение параметра  n  (тип: целый);
b - заданное значение параметра  b > 0 (тип: вещественный). 
   Версии
sfa4d_c - вычисление экспоненциального интеграла An (b) в режиме удвоенной точности; при этом параметр b должен иметь тип double. 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию:  нет

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

int main(void)
{
    /* Local variables */
    extern float sfa4r_c(int *, float *);
    static float b;
    static int n;
    static float res;

    n = 3;
    b = .25f;
    res = (float)sfa4r_c(&n, &b);

    printf("\n %16.7e \n", res);
    return 0;
} /* main */
                                      

Результат:   res  =  0.1535795*104