This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: sf_gamma questions
- From: Brian Gough <bjg at network-theory dot co dot uk>
- To: jonathan at leto dot net
- Cc: gsl-discuss at sources dot redhat dot com
- Date: Mon, 3 Dec 2001 22:59:58 +0000 (GMT)
- Subject: Re: sf_gamma questions
- References: <20011201214553.A31518@leto.net>
Jonathan Leto writes:
> Why is the error so large for gamma(x>18) ?
It's computed with an approximation so the relative error is the
appropriate measure. Note that 1/Gamma(19) is the first value which is
less than DBL_EPSILON.
> Would it be possible to just return factorial(x-1) if gamma is
> given an integer argument?
Sounds reasonable to me.