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