FUNCTION pella_tau, U

RETURN, (U*ALOG(U)/(U-1.)) - 1.
END