liti1.m

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

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