function res = cumul4(k,t,sp,v)

l = length(sp);
s = 0;

for m = 1:l
    s = s +abs(v(1,m))^2*sp(m)^k*exp(-i*t*sp(m));
end

res = l^(-1)*s;




end
