Текст подпрограммы и версий
sf41r_p.zip , sf41e_p.zip
Тексты тестовых примеров
tsf41r_p.zip , tsf41e_p.zip

Подпрограмма:  SF41R (модуль SF41R_p)

Назначение

Вычисление биномиальной функции распределения вероятностей BI (n, p).

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

Пусть мы имеем  n  испытаний Бернулли с вероятностью  p  успеха каждого такого испытания. Тогда SF41R вычисляет вероятность  РК  того, что будет получено pовно  K  успешных испытаний, и вероятность  PS  того, что число успешных испытаний будет меньше или pавно  K.

W.H.Beyer, Handbook of Tables for Probability and Statistics, The Chemical Rubber Co., Cleveland, Ohio, 1966.

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

procedure SF41R(N :Integer; K :Integer; P :Real; var PK :Real;
                var PS :Real; var IERR :Integer);

Параметры

N - заданное число испытаний Бернулли  n (тип: целый);
K - заданное число успешных испытаний (тип: целый);
P - заданная вероятность успеха каждого испытания Бернулли (тип: вещественный);
PK - вещественная переменная, значение которой в pезультате работы функции полагается равным веpоятности того, что будет получено pовно  K успешных испытаний;
PS - вещественная переменная, значение которой в pезультате работы функции полагается равным веpоятности того, что число успешных испытаний будет меньше или pавно  K;
IERR - целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом:
IERR=65 - когда число успешных испытаний Бернулли задано меньшим 0 или большим общего числа испытаний; значения выходных вероятностей положены равными 3.4E38;
IERR=66 - когда заданная вероятность успешного испытания Бернулли не принадлежит отрезку [0, 1]; значения выходных вероятностей положены равными 3.4E38.

Версии

SF41E - вычисление биномиальной функции распределения вероятностей Bi (n, p) с расширенной (Extended) точностью.

Вызываемые подпрограммы

UTSF12 - подпрограмма выдачи диагностических сообщений при работе подпрограммы SF41R.
UTSF13 - подпрограмма выдачи диагностических сообщений при работе подпрограммы SF41E.

Замечания по использованию

  Для функции SF41E параметры P, PK и PS должны иметь тип Extended и при  IERR ≠ 0 значения выходных вероятностей полагаются равными 1.7E308.

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

Unit tsf41r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF41R_p;

function tsf41r: String;

implementation

function tsf41r: String;
var
N,K,IERR :Integer;
P,PK,PS :Real;
begin
Result := '';  { результат функции }
N := 5;
K := 0;
P := 0.05;
SF41R(N,K,P,PK,PS,IERR);
Result := Result + Format(' %20.16f  %20.16f  %6d  %6d  %20.16f %6d ',
 [PK,PS,N,K,P,IERR]) + #$0D#$0A;
UtRes('tsf41r',Result);  { вывод результатов в файл tsf41r.res }
exit;
end;

end.


Результаты:

       PK  =  0.773780937
       PS  =  0.773780937
       IERR  =  0