|
Текст подпрограммы и версий ( Фортран ) gsi1r.zip |
Тексты тестовых примеров ( Фортран ) tgsi1r.zip |
|
Текст подпрограммы и версий ( Си ) gsi1r_c.zip |
Тексты тестовых примеров ( Си ) tgsi1r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) gsi1r_p.zip |
Тексты тестовых примеров ( Паскаль ) tgsi1r_p.zip |
Генерация одного псевдослучайного числа, имеющего биномиальное распределение.
Генерируется псевдослучайное число 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