Текст подпрограммы и версий ( Фортран )
gse1r.zip
Тексты тестовых примеров ( Фортран )
tgse1r.zip
Текст подпрограммы и версий ( Си )
gse1r_c.zip
Тексты тестовых примеров ( Си )
tgse1r_c.zip
Текст подпрограммы и версий ( Паскаль )
gse1r_p.zip
Тексты тестовых примеров ( Паскаль )
tgse1r_p.zip

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

Назначение

Генерация массива псевдослучайных чисел, имеющих экспоненциальное распределение со средним значением XM и дисперсией ХМ.

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

Экспоненциально распределенные псевдослучайные числа R (I) вычисляются по формуле

     R(I) = -XM*ln(U(I)) ,  XM≥0 ,

где XM - параметр экспоненциального распределения (среднее и дисперсия), U (I) - pавномеpно распределенные псевдослучайные числа, получаемые с помощью GSU1R.

Д.Кнут, Искусство программирования для ЭВМ, т. 2, "Мир", M., 1977, стp. 141-142.

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

    SUBROUTINE  GSE1R (ISEED, XM, N, R) 

Параметры

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

Версии:  нет

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

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

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

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

       DIMENSION R(3)
       ISEED = 408805506
       N = 3
       XM = 1.
       CALL  GSE1R (ISEED, XM, N, R)

Результаты:
       R(1) = 0.770362531059
       R(2) = 3.21917896882
       R(3) = 2.57395311171