Текст подпрограммы и версий
sf47r_p.zip 		Тексты тестовых примеров
tsf47r_p.zip

Функция:  SF47R (модуль SF47R_p)

Назначение

Вычисление функции χ2 - распределения вероятностей с вещественными степенями свободы.

Математическое описание

SF47R вычисляет вероятность того, что случайная величина, подчиненная χ2 - распределению с вещественным параметром  AN  (AN ≥ 0.5), меньше или pавна  x (x ≥ 0).

M.Abramowitz, I.A.Stegun, Handbook of Mathematical Functions, New York, Dover Publications, INC., 1965.

Использование 

function SF47R(X :Real; AN :Real; var IERR :Integer): Real;


   Параметры
X - заданное значение аргумента  x (тип: вещественный);
AN - заданное значение степени свободы χ2 - распределения (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда заданная степень свободы меньше 0.5; значение функции полагается равным 3.4E38;
IERR=66 - когда верхний предел интегрирования меньше 0; значение функции полагается равным 3.4E38. 
   Версии:  нет 
   Вызываемые подпрограммы
SF23C - вычисление гамма - функции  Г (x) комплексного аpгумента.
SF49R - вычисление функции нормального распределения вероятностей.
UTSF12 - подпрограмма выдачи диагностических сообщений в ходе работы функции SF47R. 
   Замечания по использованию:  нет

Пример использования 

Unit tsf47r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF47R_p;

function tsf47r: String;

implementation

function tsf47r: String;
var
IERR :Integer;
X,AN,P :Real;
begin
Result := '';  { результат функции }
X := 9.2205;
AN := 2.0;
P := SF47R(X,AN,IERR);
Result := Result + Format('%s',['  X=']);
Result := Result + Format('%20.16f ',[X]);
Result := Result + Format('%s',[' AN=']);
Result := Result + Format('%20.16f ',[AN]);
Result := Result + Format('%s',['  IERR=']);
Result := Result + Format('%3d ',[IERR]);
Result := Result + Format('%s',[' P=']);
Result := Result + Format('%20.16f ',[P]) + #$0D#$0A;
UtRes('tsf47r',Result);  { вывод результатов в файл tsf47r.res }
exit;
end;

end.


Результаты:

       P  =  0.9900506693175
       IERR  =  0