(*Liniowy generator kongruencyjny liczb pseudolosowych o rozkładzie jednostajnym (0,1)*)

a=7^5
m=2^31-1
c=0;
seed=1;
L=10000;

(*Sprawdź poniższe wartości z najbardziej niesławnego generatora w historii (RANDU IBM)*)
(*Na wykresie korelacji pojawią się wyraźne płaszczyzny*)

(*a=2^16+3*)
(*m=2^31*)

Ran[x_]:=Mod[a*x+c,m];

X=RecurrenceTable[{z[i]==Ran[z[i-1]],z[0]==seed},z,{i,1,L}];

Y=Table[{X[[i]],X[[i+1]],X[[i+2]]},{i,1,L-2}];

X[[10000]]

ListPlot[X]

(*Wykres korelacji między punktami*)
(*Nie powinny występować żadne widoczne struktury, chyba, że dla RANDU*) 
ListPointPlot3D[Y]