, L: K' j# B4 wfor i=1:M . a$ @2 F) x) \# l% r0 b. t2 Y* e( \* E2 _/ y1 {
for j=1:N, k6 M/ B8 B' \! J! P) x
7 j" c5 B' O4 R; F% u d=sqrt((i-n1)^2+(j-n2)^2);- E6 @9 w9 l% ] @ [+ N R
5 @; a! M3 U& u+ W5 C+ p( }
h=1/(1+(d/D0)^(2*n));$ L; S u, U/ h2 M3 @& q4 Y. v
3 p& A; F) F0 b% K9 K. B7 N result(i,j)=h*g(i,j);0 [3 F+ F% N3 [! g* U" Y2 f
4 w# ^8 O" X) f0 o+ C" ?- n
end4 X/ q9 y5 W' P5 T