This function inverse a Vandermonde Matrix B.
Matrix B is a n-by-n matrix, its (i,j) entry is i^(j-1),
where i,j = 1,2,...,n
for example, n = 4
1 1 1 1
1 2 4 8
1 3 9 27
1 4 16 64
This routine uses a Stirling polynomial(the first kind) coefficients
For fast operation, a C Stirling coefficient function has
been posted with name: mStirling.c. The C-version of this
inverse function is also available upon request.
Seems to work well for low dimension case. But gets worse quickly as n grows.
You can use the code to invert an Vandermonde matrix with n=20, and multiply it back and compare the result with the identity matrix. The error is huge.
First sentence should read: "This function inverts..." instead of "This function inverse..."