function res = momthprod( k,t,p,th )

s = 0;

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

res = s;
end

