def find_bounds(f, y):
    """
        given a function f,
        return lv & hv such that lv <= t(y) <= hv
    """
    x = 1
    while f(x) < y:
        x = x * 2
    lv = 0 if x == 1 else x / 2
    return lv, x