Lagrange Interpolation in Python

Ich möchte ein Polynom mit der Lagrange-Methode interpolieren, aber dieser Code funktioniert nicht:

def interpolate(x_values, y_values):
    def _basis(j):
        p = [(x - x_values[m])/(x_values[j] - x_values[m]) for m in xrange(k + 1) if m != j]
        return reduce(operator.mul, p)

    assert len(x_values) != 0 and (len(x_values) == len(y_values)), 'x and y cannot be empty and must have the same length'

    k = len(x_values)
    return sum(_basis(j) for j in xrange(k))

Ich folgte Wikipedia, aber wenn ich es ausführe, erhalte ich einen IndexError in Zeile 3!

Vielen Dan

Antworten auf die Frage(4)

Ihre Antwort auf die Frage