from pylab import plot,show N=100 a=1664525 c=101390422 m=4294967296 x=1 results=[] for i in range(N): x=(a*x+c)%m results.append(x) plot(results,"o") show()