Skip to main content
2 of 3
added 1962 characters in body

Does anyone know of any rapidly converging series for the reciprocal gamma function?

I have seen rapidly converging infinite series for Pi and the such but none for either the gamma or reciprocal gamma function. Rapidly converging algorithms are also welcome as are recursive solutions. The information being requested is required for use as a comparison.

So far some time has been spent with Borwein's gamma function (can't get the link to post, help would be appreciated!) with one specific example being, $\Gamma(1/2)= \frac{1}{\sqrt[4]{2}}\frac{AG[1]}{\sqrt{\frac{1}{\sqrt{2}}\Sigma[1]}}$

By taking advantage of AGM and making use of summations that take the difference of bn^2 and an^2 it would seem this function would be rapidly convergent however it is only good for a countable number of values of gamma without the introduction of some rather profound mathematical gymnastics.

Overall the scope of his algorithm may be limited but for the values it can calculate the decimal approximation for accuracy appears to be astounding and grows radically with each step although calculations still need to be run in order to verify this claim and to see to what degree and to how many terms this holds true.

The algorithm for this purpose can be simplified by removing the constants leaving us with,

$\frac{1}{\sum_{n=0}^{\infty} (an)^2-(bn)^2}-\left({\sum_{n=0}^{\infty} (an)^2-(bn)^2}\right)^\frac{1}{2}$

(My apologies for the badly written Latex, I will work on improving it)

From here the next step is to build a function that shows the rate of convergency for each additional calculated step for the approximation of Borwein's algorithm and that should wrap this portion up (not asking for help, just outlining steps). One method would be to use MatLab but pen and paper are more enjoyable so will likely prevail.

Are there any other links to fast converging series for either the gamma or reciprocal gamma function preferably with full functionality for the variable? Recursive identities are welcome as are exceptions for extremely rapidly converging algorithms such as Borwein's.

Thank you all for the helpful comments and suggestions. If you see any errors please point them out immediately!