Равномерная и чебышевская сетка для полинома Ньютона

from numpy import linspace
import math
def F(x):                       #Исходная функция
    return abs(3*x - 5)
 
def NN(X):
    Y = []
    for i in range(0, len(X)):
        Y.append(F(X[i]))
    return Y
 
def n(j, x, X):                # (x - x1)(x - x2)*....*(x-xn)
    n = 1
    for i in range(j):
        n *= (x - X[i])
    return n
 
def a(j, l, X, Y):             # Раздельные разности
    if j == 0:
        return Y[0]
    elif j - l == 1:
        return (Y[j] - Y[l]) / (X[j] - X[l])
    else:
        return (a(j, l + 1, X, Y) - a(j - 1, l, X, Y)) / (X[j] - X[l])
 
def N(xc, X, Y):
    N = 0
    for j in range(len(X)):
        N += a(j, 0, X, Y) * n(j, xc, X)
    return N
 
def w(x, X):
    w = 1
    for i in range(len(X)):
        w *= (x - X[i])
    return w
 
def M(X, Y):
    B = []
    j = 1
    for j in range(len(X)):
        B.append(abs(a(j, 0, X, Y)))
    return max(B)
 
def X_Chebushev(a, b, m):
    X1 = []
    for i in range(0, m):
        x = (b+a)/2 + ((b-a)/2)*((math.cos((2*i + 1)*math.pi/(2*(m+1)))))
        X1.append(x)
    return X1
 
def M_Chebushev(X, Y):
    C = []
    j = 1
    for j in range(len(X)):
        C.append(abs(a(j, 0, X, Y)))
    return max(C)
 
def M_Re(X):
    B = []
    for i in range(len(X) - 1):
        B.append((X[i] + X[i + 1])/2)
    return B
 
a1 = -2                                 #Cut boundaries
b1 = 5
 
print("Заданная функция: F(x) = |3*x - 5|")
print("Входное число: ", end = '')
x = float(input())
 
m = [10,15,25]                              # размерность сетки
for i in range(len(m)):
    m1 = m[i]                               #Chebyshev grid
    x2 = X_Chebushev(a1, b1, m1)            #Nodal values X
    y2 = NN(x2)                             #Nodal values Y
    i1 = M_Chebushev(x2,y2)                 #M(n+1)
    g = (i1 * (b1 - a1)**(m1 + 1)) / (math.factorial(m1 + 1) * 2**(2 * m1 + 1))
    print('Погрешность аппроксимации функции при использовании чебышёвской сетки: |{}|'.format(m1), g)
print('\t')
for i in range(len(m)):
    m2 = m[i]
    X = linspace(a1, b1, m2)
    x3 = M_Re(X)
    y3 = NN(x3)
    p3 = (w(x, x3))
    i3 = (M(x3,y3))
    l = abs((p3 * i3)/math.factorial(m2 + 1))
    print('Погрешность аппроксимации функции при использовании равномерной сетки: |{}|'.format(m2), l)