|
Текст подпрограммы и версий sf48r_p.zip |
Тексты тестовых примеров tsf48r_p.zip |
Вычисление обратной функции χ2 - распределения вероятностей.
SF48R вычисляет значение X такое, что
χ2 (X; AN) = P ,
где P - заданное значение вероятности (0 < P < 1) ,
AN - заданное значение степени свободы χ2 - распределения
вероятностей (0.5 < AN < 200000).
M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, Dover Publications, INC., 1965.
function SF48R(P :Real; AN :Real; var IERR :Integer): Real;
Параметры
| P - | заданное значение вероятности (тип: вещественный); |
| AN - | заданное значение степени свободы (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда значение вероятности P не принадлежит интервалу (0, 1); значение функции полагается равным 3.4E38; |
| IERR=66 - | когда заданное значение степени свободы меньше 0.5; значение функции полагается равным 3.4E38. |
Версии: нет
Вызываемые подпрограммы
| SF47R - | вычисление функции χ2 - распределения вероятностей с вещественными степенями свободы. |
| SF51R - | вычисление обратной функции нормального pаспределения вероятностей. |
| UTSF12 - | подпрограмма выдачи диагностических сообщений в ходе работы функции SF48R. |
Замечания по использованию: нет
Unit tsf48r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF48R_p;
function tsf48r: String;
implementation
function tsf48r: String;
var
IERR :Integer;
AN,P,Y :Real;
begin
Result := ''; { результат функции }
AN := 2.0;
P := 0.99;
Y := SF48R(P,AN,IERR);
Result := Result + Format(' %20.16f %20.16f %20.16f %5d ',
[P,AN,Y,IERR]) + #$0D#$0A;
UtRes('tsf48r',Result); { вывод результатов в файл tsf48r.res }
exit;
end;
end.
Результаты:
Y = 9.210340350925
IERR = 0