function res = momfibo( t,k,psi,eps,k0 )


l = length(eps);
som = 0;

for j = 1:l
    som = som + psi(k0,j)*psi(k,j)*exp(i*t*eps(j));
end

res = l^(-1)*som;



end
