clc
clear
b=0.1;
sum=0;
N=input('请输入N的值:');
for i=1:N
    x=2.0;
    y=0.0;
    r=sqrt(x*x+y*y);
    while(r>1&&r<3)
        deg = rand(1);
        if(deg<0.25)
            x=x+b;
        elseif(deg<0.5)
            x=x-b;
        elseif(deg<0.75)
            y=y+b;
        else
            y=y-b;
        end
        r=sqrt(x*x+y*y);
        if(r<=1)
            sum=sum+4;break;
        end
        if(r>=3)
            sum=sum+6;break;
        end
    end
end

我对程序基本不太懂，论文里有个小程序，前半部分能看懂，后半部分不太明白。

 返回小木虫查看更多