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