MODULE ObxRandom;
(**
   project   = "BlackBox"
   organization   = "www.oberon.ch"
   contributors   = "Oberon microsystems"
   version   = "System/Rsrc/About"
   references   = "Martin Reiser, Niklaus Wirth, Programming In Oberon, ISBN 0201565439"
   changes   = ""
   issues   = ""

**)

   VAR z: INTEGER;   (* global variable *)

   PROCEDURE Uniform* (): REAL;
      CONST a = 16807; m = 2147483647; q = m DIV a; r = m MOD a;
      VAR gamma: INTEGER;
   BEGIN
      gamma := a * (z MOD q) - r * (z DIV q);
      IF gamma > 0 THEN
         z := gamma
      ELSE
         z := gamma + m
      END;
      RETURN z * (1.0 / m)   (* value of the function *)
   END Uniform;

   PROCEDURE InitSeed* (seed: INTEGER);
   BEGIN
      z := seed
   END InitSeed;

BEGIN
   z := 314159   (* initial value of seed *)
END ObxRandom.

Arr[i]:= SHORT(ENTIER(ObxRandom.Uniform()));

MODULE RandomNumbers;
(** Детали алгоритма см. в Park S.K. and Miller K.W. (1988)
  "Random number generators, good ones are hard to find".
  Communications of the ACM, 31, 1192-1201. *)

CONST modulo* = 2147483647;  (* =2^31-1 *)
VAR z* : LONGINT;

(** Инициализация генератора целым числом: *)
   
PROCEDURE PutSeed*(seed: LONGINT);
BEGIN
    seed := seed MOD modulo;
    IF (seed = 0) THEN z := 1 ELSE z := seed END
END PutSeed;

PROCEDURE NextRND;
CONST
    a = 16807;
    q = 127773;      (* m div a *)
    r = 2836;        (* m mod a *)
VAR lo, hi, test: LONGINT;
BEGIN
    hi := z DIV q;
    lo := z MOD q;
    test := a * lo - r * hi;
    IF (test > 0) THEN z := test ELSE z := test + modulo END
END NextRND;

(** Целые случайные числа в диапазоне 0 <= x < range.
   Range должен быть в интервале 1 .. 2^31-2 *)
PROCEDURE RND*(range : LONGINT): LONGINT;
BEGIN
    NextRND;
    RETURN z MOD range
END RND;

(** Вещественные случайные числа в диапазоне 0.0 <= x < 1.0 *)
PROCEDURE Random*(): REAL;
BEGIN
    NextRND;
    RETURN (z-1)*(1 / (modulo-1))
END Random;

BEGIN
  z := 314053
END RandomNumbers.

ENTIER(x)   real type   LONGINT   largest integer not greater
         than x

Arr[i]:= SHORT(ENTIER(ObxRandom.Uniform()));