function res = momexp( t,j,pp )


l = length(pp);
som = 0;

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

res = l^(-1)*som;



end


