Algorithm 291: logarithm of gamma function | Synapse