Series Solutions & Special Functions

the Gamma function

/ GAM-uh /

The factorial n! is defined only for whole numbers — you can multiply 5 times 4 times 3 times 2 times 1, but what could 'half factorial' possibly mean? The Gamma function is the answer: a smooth, continuous function that agrees with the factorial at the integers but is defined for almost all numbers in between, filling in the gaps with a single elegant curve. It is the quiet workhorse behind nearly every special-function formula, because series coefficients constantly call for factorials of non-integer order.

It is defined by an integral, Gamma(x) = integral from 0 to infinity of t^(x-1) e^(-t) dt, which converges for x > 0 (and extends to most other values). Its single most useful property is the recursion Gamma(x + 1) = x Gamma(x) — the continuous echo of the rule n! = n times (n-1)!. Combined with Gamma(1) = 1, this recursion forces Gamma(n + 1) = n! exactly at the integers, so Gamma really is the factorial in disguise, shifted by one: factorial of n equals Gamma of (n+1). The celebrated special value Gamma(1/2) = sqrt(pi) is what makes half-integer factorials possible, and it is exactly the constant that lets the Bessel function series make sense for fractional orders.

Why it belongs to this field: the moment you allow Frobenius exponents to be non-integers, you need factorials of non-integers to write the coefficients, and Gamma supplies them. It appears in the Bessel series, in normalisation constants for orthogonal polynomials, in probability distributions, and across analysis. One honest caveat: Gamma is not defined at zero or the negative integers — it shoots off to infinity there — so it generalises the factorial everywhere except where the factorial itself would already be undefined for negative whole numbers.

Use the recursion to find Gamma(7/2): Gamma(7/2) = (5/2) Gamma(5/2) = (5/2)(3/2) Gamma(3/2) = (5/2)(3/2)(1/2) Gamma(1/2) = (15/8) sqrt(pi). This is the value of '(5/2) factorial', and such half-integer factorials are exactly what Bessel and other series coefficients demand.

The recursion Gamma(x+1) = x Gamma(x) plus Gamma(1/2) = sqrt(pi) lets you compute any half-integer factorial by hand.

Watch the shift by one: Gamma(n) = (n-1)!, not n!. The factorial of n is Gamma(n+1) — a constant source of off-by-one errors in series coefficients.

Also called
Gamma(x)generalized factorial伽瑪函數廣義階乘