|
> 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)
Sorry, war noch ein kleiner Bug drin:
factor[x_, {i_}] := (x - #)/(X[[i]] - #) &
p[x_] := Plus@@ MapIndexed[ #1 Times@@ factor[x, #2]/@ Drop[X, #2] &, Y ]
TH
|