Ich würde dir gerne noch einen zweiten Ansatz vorschlagen (v.a. da ich finde dass man bei Tschaka ebenfalls ein Gleichungssystem lösen muss [auch wenn es sehr einfach ist, da die Koeffizienentenmatrix Dreiecksgestalt hat])
Der Trick ist spezielle Polynome einzuführen. Wir betrachten die Stützstellen
$$ x_1 = -3, x_2 = -2, x_3 = -1, x_4 = 0, x_5 = 1, x_6 = 2 $$
und definieren die sogenannte Lagrange-Polynome:
$$ l_i := \prod_{j=1\\j\neq i}^6 \frac{x - x_j}{x_i - x_j} $$
für \( i = 2 \) also z.B.
$$ l_2 = \frac{x - x_1}{x_2 - x_1}\cdot\frac{x - x_3}{x_2 - x_3}\cdot\frac{x - x_4}{x_2 - x_4}\cdot\frac{x - x_5}{x_2 - x_5}\cdot\frac{x - x_6}{x_2 - x_6} $$
Diese Polynome haben die nette Eigenschaft, dass
$$ l_i(x_j) = \delta_{ij} = \begin{cases}1 & i = j \\ 0 & i \neq j\end{cases} $$
Überlege dir das vielleicht mal oben an \( l_2 \). Der "Witz" ist nun, dass man mit diesen Polynomen das Interpolationspolynom auf einfachstem Weg angeben kann:
$$ f = \sum_{i=1}^6 f_i l_i $$
Wobei mit \( f_i \) der Funktionswert an der Stelle \( x_i \) gemeint ist, d.h.
$$ f_1= -196, f_2 = -22, f_3 = 4, f_4 = 8, f_5 = 20, f_6 = 94 $$
Das gesuchte Polynom ist damit
$$ f = -196 l_1 -22 l_2 + 4 l_3 + 8l_4 + 20 l_5 + 94 l_6 $$