|
Текст подпрограммы и версий sfa4r_p.zip sfa4e_p.zip |
Тексты тестовых примеров tsfa4r_p.zip tsfa4e_p.zip |
Вычисление экспоненциального интеграла An(b) .
Подпрограмма - функция SFA4R вычисляет значение интеграла
∞
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.
function SFA4R(N :Integer; B :Real): Real;
Параметры
| N - | заданное значение параметра n (тип: целый); |
| B - | заданное значение параметра b > 0 (тип: вещественный). |
Версии
| SFA4E - | вычисление экспоненциального интеграла An (b) в режиме расширенной (Extended) точности; при этом параметр B должен иметь тип Extended. |
Вызываемые подпрограммы: нет
Замечания по использованию: нет
Unit tsfa4r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SFA4R_p;
function tsfa4r: String;
implementation
function tsfa4r: String;
var
N :Integer;
B,RES :Real;
begin
Result := ''; { результат функции }
N := 3;
B := 0.25;
RES := SFA4R(N,B);
Result := Result + Format('%15.8f ',[RES]) + #$0D#$0A;
UtRes('tsfa4r',Result); { вывод результатов в файл tsfa4r.res }
exit;
end;
end.
Результат: RES = 0.1535795*104