(12C) ~Γ(x+1)

[Gerson W. Barbosa]

(12-11-2017 01:58 PM)Dieter Wrote:  

(12-11-2017 02:19 AM)Gerson W. Barbosa Wrote:  This is based on formula 27 here, slightly modified to improve accuracy.

Maybe you should change the final steps where x^x is divided by e^x. Since both terms are calculated separately especially the first one may cause an overflow much earlier than neccessary. What about (x/e)^x? Or maybe you have an even better idea?!

Dieter

Currently run out of good ideas here. Whatever you may suggest to improve accuracy or extended the range is much appreciated. Rewriting of the 12C program for speed is also welcome.

Thank you for your interest!

Gerson.