uinform

%產生一組[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為要產生的隨機數個數



riddling

%產生一個N(0,1)正態分布隨機數
%采用篩選法，精度較高
%其他方法可參閱《現代應用數學手冊--概率統計與隨機過程卷》清華大學出版社馬振華主編
%function y=riddling()



normal

%產生n個N(a,b)正態分布隨機數
%其中a為均值，b為方差
%function x=normal(a,b,n)
function x=normal(a,b,n)