Some other values might be checked. E.g., since Gamma(z) Gamma(1-z) == pi/sin(pi z), setting z = 1/2 + i gives |Gamma(1/2 + i)|^2 = 2pi /(e^pi + e^-pi). --Dan David Makin wrote: << Thanks Warren, will take a look to verify I got the UF version at least roughly correct - of course I had no simple way of checking it beyond the obvious positive integer values - which I seem to remember it got correct to a minimum of around 10 or 11 sf, not that great but good enough to play with.

...

________________________________________________________________________________________ It goes without saying that .