uniform.m

基于MATLAB

function U=uniform(s,k)
X(1,1)=s;
for i=2:k
    X(i,1)=mod(16807*X(i-1,1),(2^31-1));
end
U=X/(2^31-1);
hist(U);
title('uniform variable U by congruent');