|
> Lagrange'sche Interpolationsformel:
> p(x)=Sum (von i=0 über n) fi * Product(von k ungleich i UND k=0 über n)
> (x-xk) / (xi-xk)
Wenn X das Array der Stützstellen und Y das dazugehörige Array
der Ordinatenwerte darstellt:
factor[x_, {i_}] := (x - #)/(X[[i]] - #)
p[x_] := MapIndexed[ #1 Times@@ factor[x, #2]/@ Drop[X, #2] &, Y ]
Gruß,
Thomas Hahn
|