def find_root_bounds(x, power):
    """x a float, power a positive int
       return low, high such that low**power <=x and high**power >= x"""
    low = min(-1, x)
    high = max(1, x)
    return low, high