Gamma Function by Stieltjes Continued Fraction

[Dieter]

(09-04-2015 06:32 PM)lcwright1964 Wrote:  Dieter, I like your rewrite so much, using registers and some other tweaks, that I have adopted it in toto for all of the units I am trying to program.

Ah, that's nice.

(09-04-2015 06:32 PM)lcwright1964 Wrote:  Here is the HP41 version, where, again inspired by JMP, I use the synthetic registers M, N, and O in place of your 0, 1, and 2.

The HP41 version can be done with a single register, similar to your original version. But if you want to use M, N and O be sure to add a CLA at the end to clean up the Alpha register so that no weird characters are left there. ;-)

(09-04-2015 06:32 PM)lcwright1964 Wrote:  It seems to run just as fast and works great:

It should deliver about nine digits, and I did not observe any cases where the error exceeded 6 units in the last place. But maybe someone else will. ;-)

Dieter