liti6.m

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

function liti6(p,mm)
efreq=zeros(1,mm);
randnum1 = binornd(1,p,1,mm);
randnum2 = unidrnd(6,1,mm);k1=0;k2=0;k3=0;
for i=1:mm
    if randnum1(i)==0         k1=k1+1;
    else
        if randnum2(i)<=3            k1=k1+1;
        elseif  randnum2(i)==6              k3=k3+1;
        else                k2=k2+1;
        end
    end
   efreq(i)=(k2+k3)/i;    
end    
num=1:mm;plot(num,efreq)