Interpolationspolynom

Wertetabellen können über das Clipboard von einem normalen Textverarbeitungssystem oder Editor nach folgender Syntax eingegeben werden:

Die Zahlen müssen reihenweise in der Datei stehen und durch Leerzeichen getrennt sein. Andere Trennzeichen sind nicht erlaubt. In Fließkommazahlen wird der Dezimalpunkt verwendet. Exponenten müssen in der Form ''123e456`` eingegeben werden. Achtung! Nach der letzten Zahl darf kein Leerzeichen oder Leerzeile folgen.

Das Einlesen von Datendateien wird mittels readdata(`dateiname`, zahlenformat, spaltenzahl) erreicht. Als Zahlenformat sind integer und float möglich.

readlib(readdata):
s:=readdata(`dateiname`,zahlenformat,spaltenzahl);
p1:=plot(s,style=point):
s1:=convert(s,array); #`(* Umwandlung der Datenliste in eine Matrix *)`:
with(linalg):
xa:=col(s1,1);# `(* Umwandlung der ersten Spalte in den x-Vektor *)`:
ya:=col(s1,2);# `(* Umwandlung der zweiten Spalte in den y-Vektor *)`:
f:=interp(xa,ya,x);# `(* Interpolationsfunktion *)`:
p2:=plot(f,x=0..11):
with(plots):
display([p1,p2]);
#`(* Graphische Darstellung des Polynoms und der Datenpunkte *)`: