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

s = 0;

for j = 1:(k)
    
        s = s +abs(momthouless3(t,j,N,p,th))^2;
end;
res = s;
end
