Post Reply 
Gamma Function by Stieltjes Continued Fraction
09-05-2015, 05:49 AM
Post: #7
RE: Gamma Function by Stieltjes Continued Fraction
(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. ;-)

Find all posts by this user
Quote this message in a reply
Post Reply 

Messages In This Thread
RE: Gamma Function by Stieltjes Continued Fraction - Dieter - 09-05-2015 05:49 AM

User(s) browsing this thread: 1 Guest(s)