@lockie Как-то плохо у меня получаетмя с maxima.
(%i32) expand(0+0.0*(1*(x-2.0)*(x-4.0)*(x-8.0)*(x-16.0)*(x-32.0))/(1*(1.0-2.0)*(1.0-4.0)*(1.0-8.0)*(1.0-16.0)*(1.0-32.0))+1.0*(1*(x-1.0)*(x-4.0)*(x-8.0)*(x-16.0)*(x-32.0))/(1*(2.0-1.0)*(2.0-4.0)*(2.0-8.0)*(2.0-16.0)*(2.0-32.0))+4.0*(1*(x-1.0)*(x-2.0)*(x-8.0)*(x-16.0)*(x-32.0))/(1*(4.0-1.0)*(4.0-2.0)*(4.0-8.0)*(4.0-16.0)*(4.0-32.0))+-79.0*(1*(x-1.0)*(x-2.0)*(x-4.0)*(x-16.0)*(x-32.0))/(1*(8.0-1.0)*(8.0-2.0)*(8.0-4.0)*(8.0-16.0)*(8.0-32.0))+1.0*(1*(x-1.0)*(x-2.0)*(x-4.0)*(x-8.0)*(x-32.0))/(1*(16.0-1.0)*(16.0-2.0)*(16.0-4.0)*(16.0-8.0)*(16.0-32.0))+31.0*(1*(x-1.0)*(x-2.0)*(x-4.0)*(x-8.0)*(x-16.0))/(1*(32.0-1.0)*(32.0-2.0)*(32.0-4.0)*(32.0-8.0)*(32.0-16.0)));
5 4 3
(%o32) - .0027467757936507934 x + .1519159226190476 x - 2.40656001984127 x
2
+ 12.121279761904763 x - 20.71150793650794 x + 10.847619047619048