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

Функция:  SF56R (модуль SF56R_p)

Назначение

Вычисление обратной функции t - распределения вероятностей Стьюдента.

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

По заданному значению вероятности  P SF56R вычисляет такое значение  X, что абсолютное значение случайной величины, подчиненной закону t - распределения вероятностей Стьюдента со степенью свободы DF (DF ≥ 1), больше  X.

G.W.Hill, Algorithm 396, Student's t - quantilies, CACM, 13 (10), 1970.

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

function SF56R(P :Real; DF :Real; var IERR :Integer): Real;


   Параметры
P - заданное значение вероятности (тип: вещественный);
DF - заданное значение степени свободы t - распределения вероятностей Стьюдента (тип: вещественный);
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда заданная степень свободы t - распределения Стьюдента меньше 1; значение функции полагается равным 3.4E38;
IERR=66 - когда заданное значение вероятности не принадлежит интервалу (0, 1); значение функции полагается равным 3.4E38. 
   Версии:  нет 
   Вызываемые подпрограммы
SF51R - вычисление обратной функции нормального pаспределения вероятностей.
UTSF14 - подпрограмма выдачи диагностических сообщений в ходе работы функции SF56R. 
   Замечания по использованию:  нет

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

Unit tsf56r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF56R_p;

function tsf56r: String;

implementation

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

end.


Результаты:    X  =  0.325000000004 ,   IERR  =  0