Ich habe das zu Hause gemacht. Und zwar mit einem Computerprogramm namens wxMaxima:
pts:[[1, 3200],[2, 760],[3, 1056],[4, 784],[5, 296],[6, 248],[7, 192],[8, 88],[9, 64],[10, 48],[11, 24],[12, 16],[13, 8],[14, 8],[15, 8],[16, 0],[17, 0]];
tpl(x) := sum(a[i]*x^i, i, 0, length(pts)-1);
eqs: [];
for i:1 step 1 thru length(pts) do eqs: append(eqs, [tpl(pts[i][1]) = pts[i][2]]);
s: solve(eqs, listofvars(eqs));
f(x) := ''(at(tpl(x), s[1]));