|
Текст подпрограммы и версий sf79r_p.zip |
Тексты тестовых примеров tsf79r_p.zip |
Вычисление обратной функции распределения вероятностей Коши.
SF79R вычисляет значение x такое, что F (X, λ ,A) = P, где F (Х, λ, А) - функция распределения Коши:
x
F(X, λ, A) = 1/πλ ∫ dz / (λ2 + (A - Z)2) , λ > 0
-∞
Справочник по специальным функциям. Под ред. M.Абрамовица, И.Стеган. M.: "Hаука", 1979.
procedure SF79R(P :Real; RL :Real; A :Real; var X :Real;
var IERR :Integer);
Параметры
| P - | заданное значение аргумента (тип: вещественный); |
| RL, A - | параметры λ, A распределения Коши (тип: вещественный); |
| X - | вещественная переменная, содержащая вычисленное значение функции; |
| IERR - | целая переменная, служащая для сообщения об ошибках, при этом: |
| IERR=65 - | когда λ ≤ 0; |
| IERR=66 - | когда P < 0 либо P > 1; значение x в обоих случаях полагается равным нулю. |
Версии: нет
Вызываемые подпрограммы
| UTSF16 - | подпрограмма выдачи диагностических сообщений при работе подпрограммы SF79R. |
Замечания по использованию: нет
Unit tsf79r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF79R_p;
function tsf79r: String;
implementation
function tsf79r: String;
var
IER :Integer;
FX :Real;
begin
Result := ''; { результат функции }
SF79R(0.5,1.0,1.0,FX,IER);
Result := Result + Format(' %16.7f %10d ',[FX,IER]) + #$0D#$0A;
UtRes('tsf79r',Result); { вывод результатов в файл tsf79r.res }
exit;
end;
end.
Результаты: FX = 1.000 000 000 005 , IER = 0