function res = cumulmoyap(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*heav(abs(sp(m)-sp(mm)),1/t);
    end
end

res = l^(-2)*s;




end
