Текст подпрограммы и версий ( Фортран )
sf62r.zip  		Тексты тестовых примеров ( Фортран )
tsf62r.zip 
Текст подпрограммы и версий ( Си )
sf62r_c.zip  		Тексты тестовых примеров ( Си )
tsf62r_c.zip 
Текст подпрограммы и версий ( Паскаль )
sf62r_p.zip 		Тексты тестовых примеров ( Паскаль )
tsf62r_p.zip

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

Назначение

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

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

SF62R производит вычисление экспоненциальных интегралов 

                    ∞
     En(x)  =  ∫ t -n e -xt dt    ,
                  1

для  n = 0, 1, 2, ... и x > 0, а также выражения вида e x En(x).

I.A.Stegun and R.Zucker, Automatic Computing Methods for Special Functions. The Exponential Integral En (x), Journal of Research of the National Bureau of Standards - B. Mathematical Sciences, Vol. 78 B, No. 4, October - December 1974.

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

    SUBROUTINE  SF62R (N, X, ENX, EXPENX) 
   Параметры
N - заданное значение параметра  n (тип: целый);
X - заданное значение аргумента  x (тип: вещественный);
ENX - вещественная переменная, значение которой в pезультате работы подпрограммы полагается pавной En (x);
EXPENX - вещественная переменная, значение которой в pезультате работы подпрограммы полагается pавной e x En (x). 
   Версии:  нет 
   Вызываемые подпрограммы:  нет 
   Замечания по использованию
 

Подпрограмма проводит вычисления для абсолютных значений N и X.

Подпрограмма выделяет несколько специальных случаев:
EN(0) = 9*1018,   если N ≤ 1;
EN(0) = 1 / (N - 1),   если N > 1;
E0(x) = e - x / x,   если x > 9*10 - 18;
E0(x) = 9*10 18,   если x ≤ 9*10 - 18 .

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

         N = 1
         X = 0.1
         CALL  SF62R (N, X, ENX, EXPENX)

Результаты:

       ENX         =  1.822923958416
       EXPENX  =  2.014642544706