def Ln(x,eps):
         assert x>0
         n = 0
         while abs(x-1)>eps:
             x **= 0.5
             n += 1
         return (x-1)*2**n