function res = momthprod3( N,k,t,p,th )

s = 0;

for j = 1:(k)
    for jj = 1:(k)
        s = s +momthouless3(t,j,N,p,th)*conj(momthouless3(t,jj,N,p,th));
    end;
end;

res = s;
end


