N=5;  
for m=1:5
    for n=1:5 
     r0=max(m,n);
     for r=r0:N
    i(m,n)=((-1)^(N-r)*(2*N-1)/(4*(N+r-1)))*nchoosek(N+r-1,2*r-1)*nchoosek(2*r-1,r-m)*nchoosek(2*r-1,r-n)*(10558/2)^(2*r-1)
    end
    end
end
式中有三个参量,反正这个很长的式子,想实现i(m,n)=sum(i(m,n)),但一直没有成功,求达人解之