cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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

This page as a plain text file.
%I A384406 #14 May 30 2025 11:26:54
%S A384406 1,397204094,2083249653,858616159,557054349,1979126465,2081507258,
%T A384406 1166038895,1141799280,106931857,142950581,1759473232,1125003378,
%U A384406 1832650327,144277780,2055193084,638219178,585429359,1481612600,2097586569,486421192,1477976737,886403653
%N A384406 Consecutive internal states of the IMSL pseudo-random number generator RNUN when started with ISEED=1 and RNOPT=3.
%C A384406 Periodic with period 2^31-2.
%C A384406 Also used by SAS and GPSS/PC.
%H A384406 Alois P. Heinz, <a href="/A384406/b384406.txt">Table of n, a(n) for n = 1..10000</a>
%H A384406 Visual Numerics, <a href="https://www.afs.enea.it/software/imsl/help/stat/stat.pdf">IMSL Fortran Numerical Library: User's Guide, Stat Library</a>, Version 7.0 (see p. 1430).
%H A384406 <a href="/index/Ps#PRN">Index entries for sequences related to pseudo-random numbers</a>.
%F A384406 a(n) = 397204094 * a(n-1) mod (2^31-1).
%p A384406 a:= proc(n) option remember; `if`(n<2, n,
%p A384406       irem(397204094*a(n-1), 2^31-1))
%p A384406     end:
%p A384406 seq(a(n), n=1..23);  # _Alois P. Heinz_, May 29 2025
%t A384406 NestList[Mod[397204094*#, 2^31 - 1] &, 1, 50] (* _Paolo Xausa_, May 30 2025 *)
%Y A384406 Cf. A096550.
%K A384406 nonn,easy
%O A384406 1,2
%A A384406 _Sean A. Irvine_, May 27 2025