function res = cumulmoy(k,j,t,sp)

l = length(sp);
s = 0;

for m = 1:l
    for mm = 1:l
        s = s +sp(m)^k*sp(mm)^j*sinc(t*(sp(m)-sp(mm)));
    end
end

res = l^(-2)*s;




end

