Spouge's approximation

In mathematics, Spouge's approximation is a formula for the gamma function due to John L. Spouge. The formula is a modification of Stirling's approximation, and has the form

$\Gamma (z+1)=(z+a)^{z+1/2}e^{-(z+a)}{\sqrt {2\pi }}\left[c_{0}+\sum _{k=1}^{a-1}{\frac {c_{k}}{z+k}}+\epsilon _{a}(z)\right]$

where a is an arbitrary positive integer and the coefficients are given by

$c_{0}=1\,$

$c_{k}={\frac {1}{\sqrt {2\pi }}}{\frac {(-1)^{k-1}}{(k-1)!}}(-k+a)^{k-1/2}e^{-k+a}\quad k=1,2,\dots ,a-1.$

Spouge has proved that, if Re(z) > 0 and a > 2, the relative error is bounded by

$\epsilon _{a}(z)\leq a^{-1/2}(2\pi )^{-(a+1/2)}.$

The formula is similar to the Lanczos approximation, but has some distinct features. Whereas the Lanczos formula exhibits faster convergence, Spouge's coefficients are much easier to calculate and the error can be set arbitrarily low. The formula is therefore feasible for arbitrary-precision evaluation of the gamma function.

References

• Spouge, John L. "Computation of the gamma, digamma, and trigamma functions", SIAM Journal on Numerical Analysis 31 (1994), no. 3, 931-944.

External links

• Gamma function with Spouge's formula - Mathematica implementation at LiteratePrograms