liti5.m

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

function proguji=liti5(mm)
%mm 是随机实验次数
frq=0;
randnum1=unifrnd(0,60,mm,1);
randnum2=unifrnd(0,60,mm,1);
randnum=randnum1-randnum2;
proguji=0;
for ii=1:mm
    if  abs(randnum(ii,1))<=20
        frq=frq+1;
    end
end
proguji=frq/mm