liti7.asv

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

mm=100
Randnum=unifrnd(0,4,1,100);
for ii=1:10000
    if Randnum(1,ii)<=1
       xRandnum(1,ii)=0;
       yRandnum=Randnum(1,ii);
   else if  Randnum(1,ii)<=2
       xRandnum(1,ii)=Randnum(1,ii)-1;
       yRandnum=1;
   else if  Randnum(1,ii)<=3
       xRandnum(1,ii)=1;
       yRandnum=Randnum(1,ii)-2;
   else     
       xRandnum(1,ii)=1-(Randnum(1,ii)-3);
       yRandnum=0;
   end    
end       
[Y,II] =sort(xRandnum)
plot(xRandnum(II),yRandnum(II),'.')