Программа выдает псевдослучайные числа, равномерно распределенные на интервале (0,1), не включая конечные точки.
Предлагается 2 различных генератора последовательности случайных чисел.
Выбор генератора задается при обращении к функции.

Структура:

Обращение:

X=D_Random(N), где:

RandomInitiate(i,k) - инициализация генератора целыми числами i,k.

N=1

Этот генератор случайных чисел впервые опубликован в статье:

Он имеет период 2^144. Полностью переносим (дает идентичную последовательность результатов на всех машинах с мантиссой не менее 48 бит в представлении с плавающей точкой).

Для получения одного случайного числа выполняется 5 операций сложения с плавающей запятой и 3 операции сложения с целыми числами.

Подпрограмма RandomInitiate(i,k) инициализирует генератор двумя целыми числами i и k
( 0 <= i <= 31328, 0 <= k <= 30081).
Это позволяет получить около 900000000 различных серий случайных чисел.
Инициализация по умолчанию: i=1802, k=9373.
RandomInitiate определяет только серию, числа серии начинаются с начала.
Полное состояние этого генератора содержится в COMMON-блоке:
COMMON /RanGen_Parameters/ istat(4),rstat(100)
real*8 rstat

D_Random(-999) выводит istat(1) и istat(2), соответствующие значениям i и k текущей серии.

N = 2

Переносимый генератор случайных чисел, описанный в статье:

Подпрограмма RandomInitiate(i,k) инициализирует генератор
двумя целыми числами i > 0 и k > 0.
Инициализация по умолчанию: i=1802, k=9373.
Эти два числа полностью определяют состояние генератора.
D_Random(-999) выдает значения i и k.

      . . .
      implicit real*8 (a-h,o-z)
      . . .
      do i=1,5
        write(*,'(2f18.14)') D_Random(1),D_Random(2)
      enddo
      . . .

     .47444060942140   .85597053361232
     .80073601515072   .38790535493994
     .51635387886227   .14669151631475
     .30258604127470   .19100434993394
     .01901665275539   .71819570289829