import math

#x = [0, 2, 3, 4, 6]
x = [-4, -2, 0, 2, 4]
#y = [1, 3, 2, 5, 7]
y = [-734, -66, 2, -2, -318]

x_zadany = 1

Buffor = 0

il_roz_matrix = [[0] * (len(x)-1) for i in range(len(x)-1)]

def roz_il():
    for xi in range(len(il_roz_matrix[0])):
        il_roz_matrix[xi][0] = (y[xi+1] - y[xi])/(x[xi+1]-x[xi])
        print(il_roz_matrix[xi][0])
    for xz in range(1, len(il_roz_matrix[0])):
        for xr in range(len(il_roz_matrix[0])-xz):
            il_roz_matrix[xr][xz] = (il_roz_matrix[xr+1][xz-1] - il_roz_matrix[xr][xz-1])/((x[xr+xz+1]-x[xr]))


def el_mult(st):
    buff = 1
    buff *= il_roz_matrix[0][st]
    for xi in range(1, st+1):
        buff *= (x_zadany - x[xi])
    return buff



def w_roz_il():

    Buffor = x_zadany
    for xi in range(len(x)-1):
        Buffor += el_mult(xi)
    return Buffor



roz_il()

for i in range(len(il_roz_matrix)):
    print(il_roz_matrix[i])


print("Wynik: ", w_roz_il())