2 lines
132 B
Text
2 lines
132 B
Text
lsf(X,Y,n)=my(M=matrix(#X,n+1,i,j,X[i]^(j-1))); Polrev(matsolve(M~*M,M~*Y~))
|
|
lsf([0..10], [1,6,17,34,57,86,121,162,209,262,321], 2)
|