VF_polyfitEven, VF_polyfitEvenwW:
Das Polynom besteht ausschließlich aus geraden Termen,
P_i = a₀ + a₂X_i² ... a_2nX_i²ⁿ
A enthät deg/2+1 Elemente. Die Koeffizienten der ungeraden Terme (alle = 0) werden nicht gespeichert.

VF_polyfitOdd, VF_polyfitOddwW:
Das Polynom besteht ausschließlich aus ungeraden Termen,
P_i = a₁X_i + a₃X_i³ ... a_2n+1X_i²ⁿ⁺¹
A enthät (deg+1)/2 Elemente. Die Koeffizienten der geraden Terme (alle = 0) werden nicht gespeichert.

Die Koeffizienten werden so berechnet, dass die Abweichung zwischen den "theoretischen" (durch Auswertung des Polynoms berechneten) P_i-Werten einerseits und den tatsächlichen Y_i-Werten minimiert wird. Genauer gesagt, wird die Anpassungstestfunktion chi-Quadrat,
c² = Summe( 1/s_i² * (P_i - Y_i)² );
minimiert. In den Einzelpunkt-gewichteten Varianten VF_polyfitwW etc. müssen die Werte 1/s_i² als Vektor InvVar bekannt sein. In den nicht-gewichteten Varianten werden alle s_i willkürlich als 1.0 angenommen.

Argumente:

A	Vektor der Länge deg+1; gibt die berechneten Koeffizienten zurück
Covar	Matrix der Dimensionen [deg+1, deg+1]; gibt die Kovarianzen der Koeffizenten zurück Falls die Kovarianzen nicht benötigt werden, rufe man die Funktion mit Covar=NULL / nil.
deg	Grad des Fit-Polynoms
X, Y, InvVar	Vektoren der Länge sizex mit den Eingabe-Daten

Es empfiehlt sich, die X-Achse so zu wählen, dass sie den Nullpunkt einschließt. Weit entfernt vom Nullpunkt wird die Anpassung nämlich immer schlechter. Dies liegt an der Form der für die Anpassung verwendeten Polynome, die nur maximal deg-1 Maxima bzw. Minima haben können. Falls Ihre X-Achse beispielsweise von 9 bis 10 reicht, sollten Sie sie lieber um den Betrag -9,5 verschieben, so dass x' = 0 in der Mitte des Bereiches liegt. Eine Verschiebung um -9,0 mit x' =0 am Beginn des Bereiches wäre ebenfalls gut.

In dem sehr seltenen Fall, dass diese Funktion keine Lösung finden kann, gibt sie 1 (TRUE) zurück und setzt alle A[i] = 0.