liti31.m

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

function proguji=liti31(10,mm)
frq=0;
randnum=unidrnd(10,mm,3);proguji=0;
for i=1:mm
    a=(randnum(i,1)+1)*(randnum(i,2)+1)*(randnum(i,3)+1)
    if  mod(a,2)==1
        frq=frq+1
    end
end; 
proguji=frq/mm