Abstract
We give an apparently new proof of Stirling’s original asymptotic formula for the behavior of for large z. Stirling’s original formula is not the formula widely known as “Stirling’s formula”, which was actually due to De Moivre. We also show by experiment that this old formula is quite effective for numerical evaluation of over when coupled with the sequence acceleration method known as Levin’s u-transform. As an homage to Stirling, who apparently used inverse symbolic computation to identify the constant term in his formula, we do the same in our proof.
Notes
1 Of course, there is no hope of changing the popular meaning of the name “Stirling’s formula”.
2 https://isc.carma.newcastle.edu.au. Remark: The ISC is currently down because a security flaw was found. Discussion is under way as to how or if this can be resolved.
3 Correctly, in the sense of Euler summation, taking even if by redefining what the infinite sum actually means: see e.g. [CitationHardy 00], for more classical work on making sense of divergent series.
4 Except of course for rounding error. We do not attempt a numerical analysis here, which appears involved. The main difficulty is predicting the number of arithmetic operations.