function res = psik(k,t,sp,ep0,vep)

l = length(sp);
s = 0;

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

res = s;
end

