> >>671の式を出したくて、Maximaに放り込んでみたが
> solve([f(1)=1,f(2)=3,f(3)=10,f(5)=23,f(6)=33,f(7)=44,f(8)=100],[a,b,c,d,e,f,g]);
> スタックオーバーフローで計算できんかったぜ

Maxima は使ったことがないけど、多項式「f(x)」の係数に同じ文字「f」を使ってもOKなの?