liti7.m

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

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