liti2.m

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

function liti2(n,p,mm)
pro=zeros(1,mm);
randnum = binornd(n,p,2,mm);a=0;
 for i=1:mm
 a=a+randnum(1,i)*randnum(2,i);   pro(i)=a/i;
end    
pro=pro,num=1:mm;plot(num,pro)