A Vandermonde matrix inversion

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
B =
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.