|
Текст подпрограммы и версий sf34r_p.zip , sf34e_p.zip |
Тексты тестовых примеров tsf34r_p.zip , tsf34e_p.zip |
Вычисление значений функции ошибок (интеграла вероятностей).
Функция SF34R вычисляет значения функции ошибок (интеграла вероятностей)
x
erf(x) = 2/√π ∫ e - t*t dt
0
для вещественных значений аргумента x.
J.F.Hart, E.W.Cheney, C.L.Lawson, Computer Approximations, Wiley, New York, 1968.
function SF34R(X :Real): Real;
Параметр
| X - | заданное значение аргумента x (тип: вещественный). |
Версии
| SF34D - | вычисление функции ошибок c расширенной (Extended) точностью. |
Вызываемые подпрограммы: нет
Замечания по использованию
| Для X ≥ 6.5 функция SF34R полагается равной 1, а для X ≥ 8 SF34D = 1. Для SF34D параметр X должен иметь тип Extended . |
Unit tsf34r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF34R_p;
function tsf34r: String;
implementation
function tsf34r: String;
var
X,Y :Real;
begin
Result := ''; { результат функции }
X := 1.2;
Y := SF34R(X);
Result := Result + Format(' %20.16f %20.16f ',[X,Y]) + #$0D#$0A;
UtRes('tsf34r',Result); { вывод результатов в файл tsf34r.res }
exit;
end;
end.
Результат: Y = 0.910313978229