liti3.asv

来自「此文件包括蒙特卡罗方法的基本思想及应用」· ASV 代码 · 共 21 行

function proguji=liti3(nn,num,mm)
%nn 是每盒中的火柴数
%num 是剩余的火柴数
%mm 是随机实验次数
frq=0;
randnum=binornd(1,0.5,mm,2*nn);proguji=0;
for i=1:mm
    a1=0;a2=0;j=1;
    while (a1<20)&(a2<20) 
        if  randnum(i,j)==1
            a1=a1+1
        else a2=a2+1
        end    
       j=j+1; 
    end    
    if abs(a1-a2)>=5  
        frq=frq+1
    end
pause
end    
proguji=frq/mm