[ENH] added stirling numbers ( + tests)

I added functions for calculating the first and second order Stirling numbers, respectively.