function res = momthouless( t,j,N,pp,th )


l = length(pp);
som = 0;

for k = 1:l
    som = som + pp(k)^(j)*exp(i*t*pp(k))*exp(-N*th(k));
end

res = l^(-1)*som;



end

