|
Текст подпрограммы и версий ( Фортран ) gse1r.zip |
Тексты тестовых примеров ( Фортран ) tgse1r.zip |
|
Текст подпрограммы и версий ( Си ) gse1r_c.zip |
Тексты тестовых примеров ( Си ) tgse1r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) gse1r_p.zip |
Тексты тестовых примеров ( Паскаль ) tgse1r_p.zip |
Генерация массива псевдослучайных чисел, имеющих экспоненциальное распределение со средним значением 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