>>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]);
スタックオーバーフローで計算できんかったぜ