def ln(x):
    assert x>0
    n = 0
    while not(0.999<x<1.001):
        x **= 0.5
        n += 1
    return (x-1)*2**n