Текст подпрограммы и версий ( Фортран )
gsb1r.zip
Тексты тестовых примеров ( Фортран )
tgsb1r.zip
Текст подпрограммы и версий ( Си )
gsb1r_c.zip
Тексты тестовых примеров ( Си )
tgsb1r_c.zip
Текст подпрограммы и версий ( Паскаль )
gsb1r_p.zip
Тексты тестовых примеров ( Паскаль )
tgsb1r_p.zip

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

Назначение

Генерация массива псевдослучайных чисел, имеющих бета-распределение.

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

Генерируются псевдослучайные числа, с плотностью распределения

     f(x) = x(P-1)*(1-x)(Q-1)*(Г(P+Q)/Г(P)*(Q)) ,  0≤x≤1 ,   P>0,  Q>0.

Здесь:
P и Q - параметры бета-распределения,
Г - гамма функция.

1.

Johnk,M.D., "Erzung von betaverteiler und gammaverteiler zufallzahlen", Metrica, 8(2), 1964.

2.

Phillips, Don T., and Beightler, Charles S., "Procedure for generating gamma variates with non-integer parameter sets", Journal of Statistical Computation and Simulation, 1, 1972, 203.

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

    SUBROUTINE  GSB1R (ISEED, P, Q, N, R) 

Параметры

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

Версии:  нет

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

GSU1R - генерация массива псевдослучайных чисел, pавномеpно распределенных в интервале (0,1);

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

  Время работы подпрограммы pастет с pостом P и Q.

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

       DIMENSION R(3)
       ISEED = 1432518450
       P = 0.5
       Q = 0.8
       N = 3
       CALL  GSB1R (ISEED, P, Q, N, R)

Результаты:
       R(1) = 0.364229459422
       R(2) = 0.959182891531
       R(3) = 0.132878061951