R3=200;
F=zeros(2*R3+1,2*R3+1);
sigma=R3/3;
r=0;
for x=1:2*R3+1
    for y=1:2*R3+1
        r=(x-R3)*(x-R3)+(y-R3)*(y-R3);
        F(x,y)=exp(-r/(2*sigma*sigma));
        K=K+F(x,y);
    end
end
posted on 2015-10-20 18:00  gyearth  阅读(255)  评论(0)    收藏  举报