function [ ans ] = guo( D )
%GUO Summary of this function goes here
%   Detailed explanation goes here
N=100;
M=100;
m=1.47;
ans=0;
for a=1:N
    for b=1:M
        ans=(((N+1-2*a)^2+(M+1-2*b)^2)*((D^2)/4+1)^(-(m+6)/2)*(1-((m+3)*(N+1-2*a)^2-(M+1-2*b)^2)*((D^2)/4)))+ans;
    end
end


end

