How to Integrate Using the Gamma Function

Last Updated: March 22, 2023

The Gamma function is a special function that extends the factorial function into the real and complex plane. It is widely encountered in physics and engineering, partially because of its use in integration. In this article, we show how to use the Gamma function to aid in doing integrals that cannot be done using the techniques of elementary calculus.

Preliminaries

The Gamma function is defined by the integral below for $\operatorname {Re} (z)>0.$ $\operatorname {Re}(z)>0.$ The Greek letter $\Gamma$ $\Gamma$ is used to denote this function.
- $\Gamma (z)=\int _{0}^{\infty }x^{z-1}e^{-x}\mathrm {d} x$
For positive integers $z=n,$ $z=n,$ the Gamma function is equal to the factorial function with its argument shifted by 1.
- $\Gamma (n)=(n-1)!$
Because the Gamma function extends the factorial function, it satisfies a recursion relation. This recursion relation is important because an answer that is written in terms of the Gamma function should have its argument between 0 and 1.
- $z\Gamma (z)=\Gamma (z+1)$
The Gamma function also satisfies Euler's reflection formula. It is from here that we can continue the function into the entire complex plane, minus the poles at the negative real numbers. Using the reflection formula, we also obtain the famous $\Gamma (1/2)={\sqrt {\pi }}.$ $\Gamma (1/2)={\sqrt {\pi }}.$ Alternatively, we can use the u-sub $u={\sqrt {x}}$ $u={\sqrt {x}}$ into the definition of the Gamma function, resulting in a Gaussian function.
- $\Gamma (z)\Gamma (1-z)={\frac {\pi }{\sin(\pi z)}}$
Below is a plot of the Gamma function along the real axis, showing the locations of the poles. This function grows faster than any exponential function.

Part 1

Part 1 of 4:

Example 1

1
Evaluate the integral below. The most important thing to check before doing the integral is to check that the integral actually converges. This integral certainly converges because the exponential decay term dominates for large $x.$ $x.$ This integral is one example of a more general integral which always converges, which we will evaluate next.
- $\int _{0}^{\infty }x^{1/3}e^{-3x}\mathrm {d} x$
- Notice that no amount of integration by parts will solve this integral.
2
Make the u-sub $u=3x$ . This allows the integral to be written with a $e^{-u}$ $e^{{-u}}$ term, which is what the Gamma function demands. It doesn't matter what the exponent on the power term is. Each time we u-sub, we also have to back-sub in order to rewrite the power term in terms of $u.$ $u.$
- $\int _{0}^{\infty }x^{1/3}e^{-3x}\mathrm {d} x={\frac {1}{3^{4/3}}}\int _{0}^{\infty }u^{1/3}e^{-u}\mathrm {d} u$
Advertisement
3
Evaluate the integral. Instead of directly evaluating, we use the Gamma function to write our answer in terms of that function. Since the argument is shifted by 1, the integral will equal $\Gamma (4/3).$ $\Gamma (4/3).$
- ${\frac {1}{3^{4/3}}}\int _{0}^{\infty }u^{1/3}e^{-u}\mathrm {d} u={\frac {\Gamma (4/3)}{3^{4/3}}}$
4
Use the recursion relation to rewrite the answer in terms of an argument between 0 and 1. It may seem pointless to write our answer in terms of this function, when we do not have a way of determining the actual value. However, there are methods to do so via other definitions. It is for this reason that we simplify our answer this way, so that we can allow computers to determine these specific values to extreme accuracy. The specific value $\Gamma (1/3)$ $\Gamma (1/3)$ has been proven to be transcendental, so there is no way to write this number algebraically.
- $\Gamma \left({\frac {4}{3}}\right)={\frac {1}{3}}\Gamma \left({\frac {1}{3}}\right)$
- ${\frac {\Gamma (4/3)}{3^{4/3}}}={\frac {\Gamma (1/3)}{3^{7/3}}}$
5
Consider the generalized integral. We assume that $a,$ $a,$ $m,$ $m,$ and $n$ $n$ are real numbers. Because this is a generalization, we have to be careful for which values that the integral fails to converge.
- $\int _{0}^{\infty }x^{m}e^{-ax^{n}}\mathrm {d} x$
6
Make the u-sub $ax^{n}$ . We can use the same technique used to evaluate the previous integral.
- $\int _{0}^{\infty }x^{m}e^{-ax^{n}}\mathrm {d} x={\frac {1}{an}}\int _{0}^{\infty }\left({\frac {u^{1/n}}{a^{1/n}}}\right)^{m-n+1}e^{-u}\mathrm {d} u$
7
Evaluate the integral in terms of the Gamma function. Of course, we pull out constants. For our answer to be consistent with where the Gamma function converges, we must put the qualifier that ${\frac {m+1}{n}}>0.$ ${\frac {m+1}{n}}>0.$
- $\int _{0}^{\infty }x^{m}e^{-ax^{n}}\mathrm {d} x={\frac {1}{na^{1+{\frac {m}{n}}-1+{\frac {1}{n}}}}}\Gamma \left({\frac {m}{n}}+{\frac {1}{n}}\right)={\frac {\Gamma \left({\frac {m+1}{n}}\right)}{na^{\frac {m+1}{n}}}}$

Part 2

Part 2 of 4:

Example 2

1
Evaluate the integral below. The integral is a product of three functions which also converges because the exponential decay term still dominates. The way we integrate this is to use Euler's formula and then to take the real part of our result.
- $\int _{0}^{\infty }x^{2}{\sqrt {x}}e^{-x}\cos x\mathrm {d} x$
2
Use Euler's formula and make a u-sub. Our u-sub will be $u=(1-i)x$ $u=(1-i)x$ from the way we have set up our integral. Every complex number should be rewritten in polar form in order to simplify the algebra.
- $\int _{0}^{\infty }x^{2}{\sqrt {x}}e^{-(1-i)x}\mathrm {d} x={\frac {1}{(1-i)^{7/2}}}\int _{0}^{\infty }u^{5/2}e^{-u}\mathrm {d} u$
3
Evaluate the integral in terms of the Gamma function. We then use the recursion relation to get the argument between 0 and 1. After simplifying further, we multiply by $(-1)e^{i\pi },$ $(-1)e^{{i\pi }},$ or 1, in order to get the angle in the exponent to something more manageable.
- ${\frac {1}{(1-i)^{7/2}}}\int _{0}^{\infty }u^{5/2}e^{-u}\mathrm {d} u={\frac {\Gamma (7/2)}{({\sqrt {2}}e^{-i\pi /4})^{7/2}}}=-{\frac {15{\sqrt {\pi }}}{8\cdot 2^{7/4}}}e^{-i\pi /8}$
4
Take the real part of the result. We can evaluate $\cos {\frac {\pi }{8}}$ $\cos {\frac {\pi }{8}}$ using the half-angle identity.
- $\int _{0}^{\infty }x^{2}{\sqrt {x}}e^{-x}\cos x\mathrm {d} x=-{\frac {15{\sqrt {\pi }}}{8\cdot 2^{7/4}}}\cos {\frac {\pi }{8}}=-{\frac {15{\sqrt {2\pi }}}{64}}{\sqrt {{\sqrt {2}}+1}}$
- We can take the imaginary part as well to obtain the sine integral for free. This is the benefit to working with trigonometric functions.
- $\int _{0}^{\infty }x^{2}{\sqrt {x}}e^{-x}\sin x\mathrm {d} x={\frac {15{\sqrt {2\pi }}}{64}}{\sqrt {{\sqrt {2}}-1}}$

Part 3

Part 3 of 4:

Example 3

1
Evaluate the integral below. We cannot directly use the Gamma function because our bounds are from 0 to 1 and there exists a logarithm inside a square root.
- $\int _{0}^{1}x^{3}{\sqrt {x\ln {\frac {1}{x}}}}\mathrm {d} x$
2
Use the u-sub $u=\ln {\frac {1}{x}}$ . This has the effect of changing the bounds, which are then negated because of the differential $\mathrm {d} u=-{\frac {1}{x}}\mathrm {d} x.$ ${\mathrm {d}}u=-{\frac {1}{x}}{\mathrm {d}}x.$ It works out nicely that the back-sub puts the exponential function into the integrand, allowing the Gamma function to do its work.
- $\int _{0}^{1}x^{3}{\sqrt {x\ln {\frac {1}{x}}}}\mathrm {d} x=\int _{0}^{\infty }u^{1/2}e^{-9u/2}\mathrm {d} u$
3
Evaluate the integral in terms of the Gamma function. Another u-sub should be used. The value $\Gamma \left({\frac {3}{2}}\right)={\frac {\sqrt {\pi }}{2}}$ $\Gamma \left({\frac {3}{2}}\right)={\frac {{\sqrt {\pi }}}{2}}$ occurs often enough that you may as well memorize it. Otherwise, going back to the recursion relation is a good way to check your work. As standard, if you can write the value in terms of constants, do so. Otherwise, just leave it in terms of the Gamma function.
- $\int _{0}^{\infty }u^{1/2}e^{-9u/2}\mathrm {d} u=\left({\frac {2}{9}}\right)^{3/2}\Gamma \left({\frac {3}{2}}\right)={\frac {\sqrt {2\pi }}{27}}$

Part 4

Part 4 of 4:

Example 4: Regularization

1
Evaluate the integral below. The integral below is divergent. You can verify this by using the u-sub $u={\sqrt {x}}.$ $u={\sqrt {x}}.$ However, there is a method by which we can assign a value to this integral in a way that makes sense. This is called regularization. The standard method is by introducing a term $e^{-\epsilon f(x)},$ $e^{{-\epsilon f(x)}},$ where $f(x)$ $f(x)$ is a positive function on the interval $[0,\infty ).$ $[0,\infty ).$
- $\int _{0}^{\infty }\cos {\sqrt {x}}\mathrm {d} x$
2
Multiply the integrand by $e^{-\epsilon {\sqrt {x}}}$ . The integral changes to taking the limit as $\epsilon \to 0^{+}.$ $\epsilon \to 0^{{+}}.$ Because this is an exponential term, it does not matter what function we choose in the exponent, as long as it is a positive function. We simply choose ${\sqrt {x}}$ ${\sqrt {x}}$ for convenience.
- $\int _{0}^{\infty }\cos {\sqrt {x}}\mathrm {d} x=\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }e^{-\epsilon {\sqrt {x}}}\cos {\sqrt {x}}\mathrm {d} x$
3
U-sub $u={\sqrt {x}}$ and rewrite the integral in terms of the complex exponential. This allows us to rewrite the integral in terms of the Gamma function.
- $\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }e^{-\epsilon {\sqrt {x}}}\cos {\sqrt {x}}\mathrm {d} x=2\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }ue^{-\epsilon u}\cos u\mathrm {d} u$
- $\int _{0}^{\infty }ue^{-(\epsilon -i)u}\mathrm {d} u$
4
Evaluate the integral in terms of the Gamma function. Remember to set $\epsilon =0$ $\epsilon =0$ at the earliest convenient time.
- $2\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }ue^{-(\epsilon -i)u}\mathrm {d} u=2\lim _{\epsilon \to 0^{+}}{\frac {1}{(\epsilon -i)^{2}}}=-2$
- Finally, we take the real part of our answer. The handling of these integrals must be done very carefully because of the divergence.
  - $\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }e^{-\epsilon {\sqrt {x}}}\cos {\sqrt {x}}\mathrm {d} x=-2$
- We can also figure out the corresponding sine integral simply by taking the imaginary part of our result.
  - $\lim _{\epsilon \to 0^{+}}\int _{0}^{\infty }e^{-\epsilon {\sqrt {x}}}\sin {\sqrt {x}}\mathrm {d} x=0$

Community Q&A

Question

Was this function used in general relativity?

Community Answer

Perhaps it may be seen when computing certain integrals as a result of solving the geodesic equation, or a similar calculation. You may instead be looking for the Lorentz factor, which is seen in special relativity and denoted by a lowercase gamma. The Gamma function is more commonly seen in quantum mechanics and quantum field theory.