[Maxima] Implementation of the Logarithm of the Gamma function
drdieterkaiser at web.de
Wed Oct 8 14:51:59 CDT 2008
> Note that using logs reduces the accuracy of the result, especially for
> large arguments that otherwise would have fit in a float.
> Don't know if this really matters or not.
Perhaps using the Logarithm of the Gamma function is only interesting for double
float values. Here we get an overflow for the Gamma function for values greater
than about 170.0. With the help of gamma_log we no longer have the problem of an
I have done some tests to compare exp(gamma_log(z)) with gamma(z). We get an
increasing absolute error for the difference of the two result with increasing
z, but the relative error is almost constant and is within the expected error of
double float calculation.
More information about the Maxima