function res = momfibothtronc( t,k,phit,psi,l,eps )



som = 0;

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

res = l^(-1)*som;



end


