
n=500;
m=(2^13)-1;
a=17;
x=100;
for i=1:n,
    x=mod(a*x,m);
    u(i)=x/m;
    i=i+1;
end
hist(u)

plot(u(1:n-1),u(2:n),'*b')