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 Γ(*x* + 1) = *x* Γ(*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

You could motivate or interpret this as saying Γ(*x* + 1/2) is approximately the geometric mean between Γ(*x* + 1) and Γ(*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.

The following plot shows the relative error in both approximations.

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

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

The next one in that sequence is a bit more accurate still for large x, but gives a poorer approximation near x=0

(x^4 + x^3 + (1/2)*x^2 + (1/8)*x)^(1/8).

The function

(2*x+1)/(4*x+3)^(1/2)

is an overestimate with a relative error at most 0.023 for all x >= 0. The relative error decreases for large x and for x=100 it is about 10^(-6).

The function

((128*x^4+384*x^3+416*x^2+188*x+31)/(128*x^3+352*x^2+324*x+97))^(1/2)

has a relative error at most 0.002 for all x>=0. The relative error decreases for large x and for x=100 is about 10^(-11).