Stirling numbers of the first kind

This c-mex function obtains the Stirling numbers of the first kind.

Description

The Stirling numbers of the first kind are defined as the coefficients of powers of x in the polynomials:

Q(x)=(x-1)(x-2)...(x-n). For example,

Q0(x)=1;
Q1(x)=x-1; %
Q2(x)=(x-1)(x-2)=x^2-3x+2;
Q3(x)=(x-1)(x-2)(x-3)=x^3-6x^2+11x-6;
...

This function calculate n>=2 case(n=0 and 1 are trivial case).

To use:
a = mStirling(4)
returns
a = 1 -10 35 -50 24

the coefficients are listed in ascending order of x.

• MATLAB release: MATLAB 6.5 (R13)
• Other requirements: VC6++

[Björn Andersson]

The numbers calculated doesn't seem to be correct, at least not for n=100.

[Andrew]