Article 3J7CK Approximating gamma ratios

Approximating gamma ratios

by
John
from John D. Cook on (#3J7CK)

Ratios of gamma functions come up often in applications. If the two gamma function arguments differ by an integer, then it's easy to calculate their ratio exactly by using (repeatedly if necessary) the fact at I(x + 1) = x I(x).

If the arguments differ by 1/2, there is no closed formula, but the there are useful approximations. I've needed something like this a few times lately.

The simplest approximation is

gamma_ratio_approx1.svg

You could motivate or interpret this as saying I(x + 1/2) is approximately the geometric mean between I(x + 1) and I(x). As we'll see in the plot below, this approximation is good to a couple significant figures for moderate values of x.

There is another approximation that is a little more complicated but much more accurate.

gamma_ratio_approx2.svg

The following plot shows the relative error in both approximations.

gamma_ratio_approx_plot.svg

By the way, the first approximation above is a special case of the more general approximation

gamma_ratio_approx3.svg

Source: J. S. Frame. An Approximation to the Quotient of Gamma Function. The American Mathematical Monthly, Vol. 56, No. 8 (Oct., 1949), pp. 529-535

hwEuFyeED8o
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments