m = 2^31-1;
a = 48271; 
N = 10000;
zs=zeros(N,1);
zs(1)=26778143;
for i=2:N
    zs(i) = mod( a*zs(i-1), m ) ;
end
xs=zs/m;
hist(xs,20)