liti9.asv

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

mm=10
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,ii);
       yRandnum(1,ii)=Randnum2(1,ii);
       ii=ii+1;
    end
        jj=jj+1;
end
         
xRandnum
yRandnum
 [Y,JJ] =sort(xRandnum)
plot(xRandnum(JJ),yRandnum(JJ),'.')