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

l = length(sp);
s = 0;

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

res = l^(-2)*s;




end
