function res = pk(k,t,sp,vep)

l = length(sp);
s = 0;


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




res = s;
end
