Вычислить действительный корень многочлена на заданном интервале скомбинировав методы хорд и касательных

from math import log10 
 
def y(x): 
  return x**4 + 3*x**3 - x**2 - 3*x + 1
    
def dy(x):
  return 4*x**3 + 9*x**2 - 2*x - 3 
  
def ddy(x): 
  return 2*(6*x**2 + 9*x-1)
 
def combi(a, b, eps=.001): 
  cnt=0 
  while abs (a - b) > eps: 
      if y(a)*ddy(a) < 0: 
          a -= y(a)*(a-b)/(y(a)-y(b)) 
      else:
          a -= y(a)/dy(a)
        
      if y(b)*ddy(b) < 0:
          b -= y(b)*(b-a)/(y(b)-y(a))
      else:
          b -= y(b)/dy(b)
      cnt -= 1
  x = (a+b)/2
  # Проверка
  k = -int(log10(eps)) 
  print(f'y({round(x, k)}) ={round(y(x), k)}') 
    
  return f'x = {round(x, k)}\nЧисло итераций: {cnt}' 
    
print(combi(-2,0))

from math import log10 
 
def y(x): 
    return x**4 + 3*x**3 - x**2 - 3*x + 1
 
def dy(x):
    return 4*x**3 + 9*x**2 - 2*x - 3 
 
def ddy(x): 
    return 2*(6*x**2 + 9*x-1)
 
def combi(a, b, eps=.001): 
    cnt=0 
    while abs (a - b) > eps: 
        if y(a)*ddy(a) < 0: 
            a -= y(a)*(a-b)/(y(a)-y(b)) 
        else:
            a -= y(a)/dy(a)
 
        if y(b)*ddy(b) < 0:
            b -= y(b)*(b-a)/(y(b)-y(a))
        else:
            b -= y(b)/dy(b)
        cnt += 1
        x = (a+b)/2
    # Проверка
    k = -int(log10(eps)) 
    print(f'y({round(x, k)}) ={round(y(x), k)}') 
    
    return f'x = {round(x, k)}\nЧисло итераций: {cnt}\n'
 
 
 
print(combi(-4,-2))
print(combi(-2,-1))
print(combi(-1,2))

print(combi(-4,-2))
print(combi(-2,-1))
print(combi(-1,4))

print(combi(-4,-2.5))
print(combi(-2.5,-1))
print(combi(-0.5, 0.5))
print(combi(0.5,2))

print(combi(-4,-2.5))
print(combi(-2.5,-1))
print(combi(-0.5, 0.5))
print(combi(0.5,2))
print('\n=====================\n')
print(combi(-2,1))
print(combi(-2,0))