Inszwischen habe ich eine weitere Lösung (Werner Sans, Praktische Numerik mit Mathematica, Teubner Verlag) gefunden:
L[i_,x_]:=Product[IF[i==j,1,(x-xw[i])]/IF[i==j,1,(xw[i]-xw[j])],{j,0,n}]
pn[x_]:=Sum[yw[i]*L[i,x],{i,0,n}]
-> Wobei die Funktionen xw[i oder j],yw[i] noch definiert werden müssen...
Mit freundlichen Grüßen aus Düsseldorf
Winfried Bilgic
______________________________________________________________________________
Viren? Wir wissen nicht was Ihr Arzt empfiehlt. Wir empfehlen den
Virencheck für Ihre E-Mail-Anhänge! http://freemail.web.de/features/?mc=021159
DMUG-Archiv, http://www.mathematica.ch/archiv.html