Полиномиальное уравнение

import numpy as np
coeff = [1. ,-2. , -12., -6. , 4.]
res = np.roots(coeff)
print(res)

def fun(x):
    return 2 - x*(4 + x*(4 - x))
 
def dfun(x):
    return 3*x*x - 8*x - 4
 
def y_fun(x):
    return x**4 - 2*x**3 - 12*x*x - 6*x + 4
 
 
x0 = -4
xk = 10
h = 1
xo = [x0]
while x0 < xk:
    if dfun(x0)*dfun(xo[-1]) < 0:
        xo.append(x0)
    x0 += h 
xo.append(xk)
n = len(xo)
x = []
for i in range(1, n):
    c = 1
    L = xo[i-1]
    R = xo[i]
    eps = 1e-7
    while R - L > eps:
        M = (R + L) / 2
        if fun(M)*fun(L) > 0:
            L = M
        else:
            R = M
        c += 1
    x.append(M)
    
    print(f'x = {M:.3f}  y = {fun(M):.4f}  итераций = {c}')

def delit(a) :
    res = set()
    i = 1
    while i * i < a + 1 :
        if a % i == 0 :
            res.update({i, a//i, -i, - a//i})
        i += 1
    return res 
 
def fun(x):
    return 2 - x*(4 + x*(4 - x))
 
def dfun(x):
    return 3*x*x - 8*x - 4
 
def y_fun(x):
    return x**4 - 2*x**3 - 12*x*x - 6*x + 4
 
x = []
dd = delit(4)
for d in dd:
    if y_fun(d) == 0:
        x.append([d, 0])
x0 = -4
xk = 10
h = 1
xo = [x0]
while x0 < xk:
    if dfun(x0)*dfun(xo[-1]) < 0:
        xo.append(x0)
    x0 += h 
xo.append(xk)
n = len(xo)
for i in range(1, n):
    c = 1
    L = xo[i-1]
    R = xo[i]
    eps = 1e-7
    while R - L > eps:
        M = (R + L) / 2
        if fun(M)*fun(L) > 0:
            L = M
        else:
            R = M
        c += 1
    x.append([M, c])
for i in sorted(x):
    print(f'x = {i[0]:.3f}  y = {y_fun(i[0]):.4f}  итераций = {i[1]}')