Текст подпрограммы и версий ( Фортран )
gsi1r.zip
Тексты тестовых примеров ( Фортран )
tgsi1r.zip
Текст подпрограммы и версий ( Си )
gsi1r_c.zip
Тексты тестовых примеров ( Си )
tgsi1r_c.zip
Текст подпрограммы и версий ( Паскаль )
gsi1r_p.zip
Тексты тестовых примеров ( Паскаль )
tgsi1r_p.zip

Подпрограмма:  GSI1R

Назначение

Генерация одного псевдослучайного числа, имеющего биномиальное распределение.

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

Генерируется псевдослучайное число R такое, что

     P{R=K}  =  CNK  pK (1 - p) (N-K) ,     K = 0, 1, ..., N ,    0 ≤ p ≤ 1 , 

где N и p - параметры биномиального распределения.

Если N < 15, то число R получается как сумма N независимых псевдослучайных чисел, принимающих значения 1 и 0 с вероятностями  p и (1 - p) соответственно.

Для N ≥ 15 используется нормальная аппроксимация биномиального распределения.

A.D.Relles, A simple algorithm for generating binomial random variables when N is large, Journal of the American Statistical Association, 339(67), 1972, 612-613.

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

    REAL FUNCTION  GSI1R (ISEED, N, P) 

Параметры

ISEED - целая переменная, значение которой перед обращением к подпрограмме может быть любым целым числом в пределах [1,2147483646]; по окончании работы ей присваивается новое значение, котоpое может быть использовано при последующем вхождении в подпрограмму;
N - параметр биномиального распределения (N ≥ 0) (тип: целый);
P - параметр биномиального распределения (0 ≤ P ≤ 1) (тип: вещественный).

Версии : нет

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

GSU2R - генерация одного псевдослучайного числа, pавномеpно распределенного в интервале (0, 1) ;
GSN2R - генерация одного псевдослучайного числа, ноpмально распределенного с нулевым средним значением и единичной дисперсией;
SF35R - вычисление значений дополнительной функции ошибок (дополнительного интеграла вероятностей).

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

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

       N = 2
       P = 0.01
       ISEED = 123456
       R = GSI1R (ISEED, N, P)

Результаты:    R = 0.0