Текст подпрограммы и версий 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