function res = momfibotronc( t,k,psi,l,eps )


som = 0;

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

res = l^(-1)*som;



end


