Алгоритм бинарного поиска для решения кубического уравнения с одним корнем

a, b, c, d = map(int, input().split())
l = -1000000
r = 10001000
for i in range(100):
    m = (r + l)/2
    if a * m**3 + b*m**2 + c*m + d > 0:
        r = m
    else:
        l = m
print(f"{m:.{4}f}")

>>> '{:.4f}'.format(m)
'0.1235'

#
# Уточнение корня уравнения f(x)=0 на отрезке [a,b] с точностью eps
#
def root(f,a,b,eps=1.0e-8):
    fa=f(a)
    fb=f(b)
    if fa*fb>0:
        return None
    while(True):
        c=(a+b)*0.5
        if abs(a-b)<eps:
            return c
        fc=f(c)    
        if abs(fc)<eps:
            return c
        if fc*fa<0:
            b=c
            fb=fc
        else:
            a=c
            fa=fc
 
#
# Поиск всех корней уравнения f(x)=0 на большом отрезке [a,b] с шагом h и точностью eps
#
def all_roots(f,a,b,h,eps=1.0e-8):
    res=[]
    x1=a
    f1=f(x1)
    while(x1<b):
        x2=x1+h
        f2=f(x2)
        if f1*f2<0:
            r=root(f,x1,x2,eps)
            res.append(r)
        x1=x2
        f1=f2
    return res
 
r=all_roots(lambda x: (x-5)*(x+7)*(x-12)*(x+5.6789)*(x-34)*(x+99),-100,100,0.1,1.0e-15)
print(r)