R3=;
F=zeros(*R3+,*R3+);
sigma=R3/;
r=;
for x=:*R3+
for y=:*R3+
r=(x-R3)*(x-R3)+(y-R3)*(y-R3);
F(x,y)=exp(-r/(*sigma*sigma));
K=K+F(x,y);
end
end
R3=;
F=zeros(*R3+,*R3+);
sigma=R3/;
r=;
for x=:*R3+
for y=:*R3+
r=(x-R3)*(x-R3)+(y-R3)*(y-R3);
F(x,y)=exp(-r/(*sigma*sigma));
K=K+F(x,y);
end
end