def lagrange(XX, YY): n = len(XX) P = 0 for i in range(0,n ): P = P + YY[i] * coefficient(XX, i) return expand(P) print(lagrange([1, -3, 8], [6, -8, -4]))