## The Gamma Function, Beta Function, and Duplication Formula

The title might as well continue — *because I constantly forget them and hope that writing about them will make me remember.* At least afterwards I’ll have a centralized repository for my preferred proofs, regardless.

In this note, we will play with the Gamma and Beta functions and eventually get to Legendre’s Duplication formula for the Gamma function. This is part reference, so I first will write the results themselves.

**1. Results **

We define the Gamma function for $latex {s > 0}$ by $$ \Gamma(s) := \int_0^\infty t^s e^{t} \frac{dt}{t}. \tag{1}$$

Similarly, we define the Beta function by $$ B(a,b) := \int_0^1 t^{a-1}(1-t)^{b-1}dt \tag{2}$$

for $latex {a, b > 0}$.

From these defininitions, it is not so obvious that these two functions are intimately related – but they are! In fact, $$ \frac{\Gamma(x)\Gamma(y)}{\Gamma(x + y)} = B(x,y) \tag{3}$$

Evaluating the Gamma function at integers is easy. We can use the relation with the Beta function to evalate it at half-integers too. $$ \Gamma(\frac{1}{2}) = \sqrt \pi \tag{4}$$

Finally, we can relate the values at half-integers and integers in an intimate way. $$ \Gamma(z)\Gamma(z + 1/2)=2^{1-2z}\sqrt{\pi}\Gamma(2z) \tag{5}$$

for $latex {\text{Re}(z) > 0}$.

** 1.1. Proof of 3 **

We begin by writing down a different representation of the Beta function.

$$ B(a,b) = \int_0^\infty \frac{u^a}{(1+u)^{a+b}}\frac{du}{u}, $$

which is in terms of the Haar measure and is generally more agreeable. *Proof:* Consider the (un-inspired) substitution $latex {u = \frac{t}{1-t}}$, or equivalently $latex {t = \frac{u}{1+u}}$. Then the bounds $latex {0 \mapsto 0}$ and $latex {1 \mapsto \infty}$, and the integrand transforms exactly into the form in the proposition. $latex \Box$

We will also want a different representation of the Gamma function. $$ \int_0^\infty e^{-pt} t^z \frac{dt}{t} = \frac{\Gamma(z)}{p^z}. $$

*Proof:* This comes quite quickly. Performing the change of variables $latex {s = pt}$ in the integral definition of the Gamma function pops out the extra $latex {p^z}$ factor and gives this form of the integral. $latex \Box$

We can now put these together. Rearranging the lemma above gives $$ \frac{1}{p^z} = \frac{1}{\Gamma(z)}\int_0^\infty e^{-pt}t^{z-1}dt. $$

Thinking of $latex {p = 1+u}$ and $latex {z = a+b}$, we can substitute this expression inside the lemma-given integral expression for the Beta function.

$$\begin{align} B(a,b) &= \frac{1}{\Gamma(a+b)}\int_0^\infty e^{-t}t^{a+b-1}dt\int_0^\infty e^{-ut}u^{a-1}du \\

&= \frac{\Gamma(a)}{\Gamma(a+b)}\int_0^\infty e^{-t}t^{b-1}dt \\

&= \frac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)}, \end{align}$$

where the first Gamma factor pulled out $latex {a}$ factors of $latex {t}$ from the first integral. This completes the proof of Prop 3.

** 1.2. Proof of Prop 4 **

We begin with another integral representation of the Beta function. $$ B(a,b) = 2\int_0^{\pi/2}(\cos u)^{2a-1}(\sin u)^{2b – 1}du $$

for $latex {a,b}$ with positive real part. *Proof:* This comes immediately from the change of variables $latex {t = \cos^2 u}$ in the integral definition of the Beta function. It’s necessary to flip the bounds of integration to cancel the negative sign from the sign of the change of variables. $latex \Box$

In this form, it is particularly easy to see that $latex {B(\frac{1}{2}, \frac{1}{2}) = \pi}$, since we integrate the constant function $latex {1}$ from $latex {0}$ to $latex {\pi/2}$ and multiply the result by $latex {2}$. And from Prop 1, we know that $latex {B(\frac{1}{2}, \frac{1}{2}) = (\Gamma(\frac{1}{2})^2}$ (as $latex {\Gamma(1) = 1}$).

Thus $latex {\Gamma(\frac{1}{2}) = \sqrt \pi}$, and we know it’s the positive square root because $latex {\Gamma(\frac{1}{2})}$ is clearly positive. This completes the proof.

This is my favorite proof, as it uses neither complex analysis nor multivariable integration – both of which are dear to my heart, but separate from the pleasant theory of the Gamma function.

** 1.3. Proof of Theorem 5**

Start from $$ \frac{\Gamma(z)\Gamma(z)}{\Gamma(2z)} = B(z,z) = \int_0^1 u^{z-1}(1-u)^{z-1}du. $$

Perform the substitution $latex {u = \frac{1+x}{2}}$, so that $latex {du = dx/2}$. This transforms the above integral into $$ 2^{1-2z}\cdot 2\int_0^1 (1-x^2)^{z-1}dx. $$

$$ B(m,n) = 2\int_0^1 x^{2m – 1}(1-x^2)^{n-1}dx $$

*Proof:* This is immediate upon the change of variables $latex {t = x^2}$ in the defining integral for the Beta function. $latex \Box$

This allows us to recognize the integral above $$ 2^{1-2z}\cdot 2\int_0^1 (1-x^2)^{z-1}dx = 2^{1-2z}B(\frac{1}{2}, z). $$

Rewriting in terms of Gammas, $$ B(\frac{1}{2}, z) = \frac{\Gamma(\frac{1}{2})\Gamma(z)}{\Gamma(z + \frac{1}{2})}. $$

In total, we have that $$ \frac{\Gamma(z)\Gamma(z)}{\Gamma(2z)} = 2^{1-2z}\frac{\Gamma(\frac{1}{2})\Gamma(z)}{\Gamma(z + \frac{1}{2})}. $$

Rearranging, and using that $latex {\Gamma(1/2)= \sqrt \pi}$ as above, we see that $$ \Gamma(2z) = \frac{2^{2z-1}}{\sqrt \pi} \Gamma(z)\Gamma(z + \frac{1}{2}), $$

which is what we wanted to show.