Переделать задачу без numpy

import numpy as np
import matplotlib.pyplot as plt

import numpy as np
import matplotlib.pyplot as plt
def func():
    z = 0
    coefficients = []
    try:
        n = int(input("Введите степень полинома:\n"))+1
    except:
        print("Степень полинома должна быть ЦЕЛЫМ числом!")
        func()
    print('Последовательно вводите коэффициенты полинома через Enter.')
    while len(coefficients) != n:
        try:
            i = float(input())
            coefficients.append(i)
        except:
            print("Степень полинома должна быть числом.")
    return n,coefficients
def f(x):
    s = 0
    for i in range(n):
        s +=coefficients[i]*x**(n-i)
    return s
def dichotomy (a,b, eps): # отрезок от a до b делим на n частей, погрешность eps
 
    assert a!=0,  'a равно 0'
    assert b!=0, 'b равно 0'
    # сначала отделим корни
    n = 10001
    setka=np.linspace(a, b, n)
    # далее уточним корни
    for x,y in zip(setka, setka[1:]):
        if f(x) * f(y) > 0: # если на отрезке нет корня, смотрим следующий
            continue
        root = None
        while ( abs(f(y)-f(x)) )>eps:     # пока отрезок больше заданной погрешности, выполняем нижестоящие операции:
            mid = (y+x)/2                   # получаем середину отрезка
            if f(mid) == 0 or f(mid)<eps:    # если функция в середине отрезка равну нулю или меньше погрешности:
                root = mid                  # корень равень серединному значению
                break
            elif (f(mid) * f(x)) < 0:       # иначе если произведение функции в середине отрезка на функцию в т. а <0
                y = mid                     # серединой становится точка b
            else:
                x = mid                     #в другом случае - точка а
        if root:
            return root
if __name__ == '__main__':
    n,coefficients = func()
    a=float(input("Введите границу интервала a:\n"))
    b=float(input("Введите границу интервала b:\n"))
    eps=float(input("Введите границу интервала eps:\n"))
    res = (dichotomy(a,b,eps))
    print('Корень по методу дихотомии находится в точке:')
    print(res)