function res = phith(polys, psi,k,ep)

s = 0;

if k ==1
    s =  psi(1,ep);
else
    s = polys(k-1,ep)*psi(1,ep);
end



res = s;

end

