%產(chǎn)生一個(gè)N(0,1)正態(tài)分布隨機(jī)數(shù)
%采用篩選法，精度較高
%其他方法可參閱《現(xiàn)代應(yīng)用數(shù)學(xué)手冊(cè)--概率統(tǒng)計(jì)與隨機(jī)過程卷》清華大學(xué)出版社馬振華主編
%function y=riddling()
function y=riddling()
sign=0;
while 1
    x=rand(1,2);
    v1=2*x(1)-1;v2=2*x(2)-1;
    s=v1^2+v2^2;
    if s<=1
        a=sqrt(-2*log(s)/s);
        y=v1*a;
        sign=1;
    end
    if sign==1
        break;
    end
end