Returns the Log Gamma function of «X». Without numeric overflow, this function is exactly equivalent to
Because the gamma function grows so rapidly, it is often much more convenient to use LGamma() to avoid numeric overflow.
LGamma(a) + LGamma(b) - LGamma(a+b)
which is less susceptible to numeric underflow for really large values of