liti9.m

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

mm=10000
Randnum1=unifrnd(-1,1,1,2*mm);
Randnum2=unifrnd(-1,1,1,2*mm);
xRandnum=zeros(1,mm);
yRandnum=zeros(1,mm);
s=Randnum1.^2+Randnum2.^2
ii=1;jj=1;
while ii<mm
    if  s(1,jj)<=1;
       xRandnum(1,ii)=Randnum1(1,jj);
       yRandnum(1,ii)=Randnum2(1,jj);
       ii=ii+1;
    end
    jj=jj+1;
end
[Y,JJ] =sort(xRandnum);
plot(xRandnum(JJ),yRandnum(JJ),'.')