%產生一組[0，1]分布的隨機數，采用余數法
%從文獻中得知下列參數組較為適用：
%y的初值為1，w=7,模M=10^10(偽隨機數周期為5*10^7)
%y的初值為1，w=5^13,模M=2^36(偽隨機數周期為2^34,約2*10^10)
%y的初值為1，w=5^17,模M=2^42(偽隨機數周期為2^40,約10^12)
%--------------------------------------------------------
%function x=uniform(y,w,M,n)
%n為要產生的隨機數個數
function x=uniform(y,w,M,n)
for i=1:n
    x(i)=y./M;
    y=rem(w*y,M);%rem()為求余函數
end