A384406 Consecutive internal states of the IMSL pseudo-random number generator RNUN when started with ISEED=1 and RNOPT=3.

1, 397204094, 2083249653, 858616159, 557054349, 1979126465, 2081507258, 1166038895, 1141799280, 106931857, 142950581, 1759473232, 1125003378, 1832650327, 144277780, 2055193084, 638219178, 585429359, 1481612600, 2097586569, 486421192, 1477976737, 886403653

Offset: 1

Author

Sean A. Irvine, May 27 2025

Keywords

• nonn
• easy

Comments

Periodic with period 2^31-2.

Also used by SAS and GPSS/PC.

Links

• Alois P. Heinz, Table of n, a(n) for n = 1..10000
• Visual Numerics, IMSL Fortran Numerical Library: User's Guide, Stat Library, Version 7.0 (see p. 1430).
• Index entries for sequences related to pseudo-random numbers.

Crossrefs

Cf. A096550.

Programs

Maple

a:= proc(n) option remember; `if`(n<2, n,
      irem(397204094*a(n-1), 2^31-1))
    end:
seq(a(n), n=1..23);  # Alois P. Heinz, May 29 2025

Mathematica

NestList[Mod[397204094*#, 2^31 - 1] &, 1, 50] (* Paolo Xausa, May 30 2025 *)

Formula

a(n) = 397204094 * a(n-1) mod (2^31-1).