def bisec(fun, a, b, eps=0.0001):
    while
    c = (a+b)/2

    if np.sign(fun(a)) != np.sign(fun(c)):
        b = c
    else:
        a = c