Unleashing Mathematical Mysteries

Mathematics
Author

Abhirup Moitra

Published

April 15, 2023

“Nothing takes place in the world whose meaning is not that of some maximum or minimum.”― Leonhard Euler

A Brief Introduction

Esteemed folks of mathematical sciences are quite familiar with the name Leonhard Euler, was one of the greatest mathematicians in history: not only did he produce outstanding mathematics, but he also produced it at an outrageous rate, publishing more than any other single mathematician before or after him. If the quality of output multiplied by the quantity of output is equal to greatness, then Euler is the greatest mathematician ever. His discoveries greatly influenced the fields of mathematics and physics. Perhaps the best-known of Euler’s findings is the Euler identity, which shows the relationship between fundamental mathematical constants and is often called the most beautiful equation in mathematics. He also introduced a notation for writing mathematical functions that are widely used today.

Story Behind Famous Eulerian Integral

Imagine it is 1743 and you are a mathematician at the peak of your powers. You are confronted by the following integrals which have never before been encountered by mathematicians:

\[ \text{I}_1 = \int_0^{\infty} \text{sin}(s^2)ds \ \; \text{and} \;\;\ \text{I}_2= \int_0^{\infty} \text{cos}(s^2)ds. \]

It is not at all obvious just what to do. While primarily a mathematician, Euler was often stimulated in his mathematics by questions raised during the study of physical problems. In this case, the above two integrals occurred in an analysis of the physics of a coiled spring. Much later, in 1815, the French mathematicians Cauchy and Poisson developed equivalent integrals in certain hydro-dynamical applications. 3 years later, in 1818, the French scientist Augustin Jean Fresnel (1788-1827) also ran up against \(\text{I}_1\) and \(\text{I}_2\) in his study of the diffraction of light, and today they are usually called the Fresnel Integrals. One does still see them also called the Euler integrals, however, and it was Euler who first evaluated them. But it wasn’t easy to do these evaluations, even for a man with Euler’s powerful intellect. As he wrote of them in 1743,

“We must admit that analysis will make no small gain should anyone find a method whereby, approximately at least, the value of [these integrals] would be determined…. This problem does not seem to be unworthy of the best strength of geometers.”

The best that even Euler could do in 1743, however, was to find convergent infinite series for \(\text{I}_1\), and \(\text{I}_2\) that allowed him to make numerical calculations; nearly 40 years later he solved the problem exactly using complex quantities. What I will show you next is essentially what Euler reported to the St. Petersburg Academy on April 30, 1781, but which was published only years after his death. The key idea to evaluating \(\text{I}_1\), and \(\text{I}_1\) is to start with yet another of Euler’s creations, with what today is called the gamma function. It is defined as,

\[ \Gamma(n)=\int_0^{\infty} e^{-x}x^{n-1}dx \;\;\text{; where}\;\ n\gt 0 \]

Dealing Gamma Integral As Improper Integral

While you learn integration you’ll work with two main types of integrals, definite integrals, and indefinite integrals. The difference between the two is those definite integrals have limits of integration, and indefinite integrals do not. When you first start working with definite integrals you’ll use simple real numbers for your limits of integration, but you’ll soon see that more is possible. It’s very common to see upper limits of integration set to infinity, and lower limits set to negative infinity.

Indefinite Integral Proper Definite Integral Improper Definite Integral
\(\int e^{-x}dx\) \(\int_0^1 e^{-x}dx\) \(\int_0^{\infty} e^{-x}dx\)

Any definite integral that has one or more infinite limits of integration, or an integrand that approaches infinity within its limits of integration is known as an improper integral. Improper integrals pop up all the time in calculus and other higher-level math courses. Here we’ll look at the properties of one famous function defined by an improper integral known as the gamma function.

The Necessity of Improper Integral

The purpose of using improper integrals is that one is often able to compute values for improper integrals, even when the function is not integrable in the conventional sense (as a Riemann integral, for instance) because of a singularity in the function as an integrand or because one of the bounds of integration is infinite.

Convergence of Gamma Integral

As it’s already claimed that gamma integral is an improper integral, so before going to compute \(\int_0^{\infty} e^{-x}x^{n-1}dx\) we need to check whether this improper integral is convergent or not.

\(\text{Here}\ \; \text{f}(x)=x^{n-1}e^{-x}\). We write the given integral as,

\[ \int_0^{\infty} x^{n-1}e^{-x}dx=\int_0^1 x^{n-1}e^{-x}dx\; + \int^{\infty}_1 x^{n-1}e^{-x}dx \hspace{0.75 cm}...(*) \]

We consider the two integrals separately.

First, consider \(\int_0^{\infty} x^{n-1}e^{-x}dx\). When \(n\ge 1\), this is a proper integral, hence convergent.

When \(m \lt 1\), the integrand is unbounded at \(x=0\). Then we have,

\[ \lim_{x \;\to \;0+0}\;(x-0)^{\mu}f(x)= \lim_{x \;\to\; 0+0} x^{1-n}.x^{n-1}e^{-x}=1 \]

which is a finite and non-zero limit. So, the integral will be convergent when \(0\lt 1-n\lt 1\), i.e.,if ,\(0\lt n\lt1\) and divergent if \(n \le0\).

Next, we consider \(\int_1^{\infty} x^{n-1}e^{-x}dx\). We see that,

\(e^x \gt \frac{x^k}{k!}\gt \frac{x^{n+1}}{k!},\text{where}\; x\gt 1\), \(k\) is a positive integer such that \(k\gt n+1\).

\[ \therefore\; e^{-x} \lt \frac{k!}{x^{n+1}},\text{or} \ \; e^{-x} x^{n-1} \lt \frac{k!}{x^{n+1}}.x^{n+1}, \]

\[ e^{-x}x^{n-1} \lt \frac{k!}{x^2} \]

\[ \int_1^{\infty} x^{n-1}e^{-x}dx \lt \int_1^{\infty} \frac{k!}{x^2} < k! \; \;\text{ a finite number} \]

\(\text{So},\; \int_1^{\infty} x^{n-1}e^{-x}dx\) is convergent for every value of \(\sqrt{n}\). Hence, both the integrals in \((*)\) are convergent if \(n \gt 0\).

Hence, the given integral is convergent when \(n \gt 0\).

Now, we will evaluate the \(\int_0^{\infty} x^{n-1}e^{-x}\) and then we will be unraveling the mystery behind this math miracle.

Evaluation of Gamma Integral

We have \(\Gamma(n) = \int_0^{\infty} x^{n-1}e^{-x}\). We will solve it by applying integration by parts, keeping \(x^{n-1}\) as first function,

\[ \Gamma(n) =\;\big[ -e^{-x}x^{n-1}\big]_0^{\infty}+\;(n-1) \int_0^{\infty}x^{n-2}e^{-x}dx \]

\[ = (n-1) \int_0^{\infty} x^{n-2}e^{-x}dx \]

\[ \text{since},\lim_{x \;\to \;\infty} \dfrac{x^{n-1}}{e^x}= \lim_{x \;\to \;\infty} \dfrac{x^{n-1}}{1+\frac{x}{1!}+\frac{x^2}{2!}+\ldots+\frac{x^n}{n!}+\ldots}=0 \]

\[ \text{i.e.,} \Gamma(n)=(n-1)\Gamma(n-1). \]

\(\text{Applying this result again}, \Gamma(n)=(n-1)(n-2)\Gamma(n-2).\)

\(\text{Hence, if}\;n\; \text{is positive integer then repeating the process above,}\)

\(\Gamma(n)=(n-1)(n-2)\ldots3.2.1\;\Gamma(1)\). \(\text{But}\; \int_0^{\infty}x^0e^{-x}=\ \big[-e^x\big]_0^{\infty}=1.\)

\[ \text{So}, \boxed{\Gamma(n)=(n-1)!} \]

Remarks :

  • \(\Gamma(n)=(n-1)\Gamma(n-1)\; \forall n \in \mathbb{Z^{+}}\)

  • \(\Gamma(n)=(n-1)!\; \forall n\in\mathbb{Z^+}\)

  • \(\Gamma(1)=1,\Gamma(0)={\infty},\Gamma(-n)={-\infty}\; \text{where}\ n\ \text{is positive integer}\)

Corollary

We have,

\[ I_n=\int_0^{\infty} e^{-x}x^ndx = \big[-x^ne^{-x}\big]_0^{\infty}\ +\ n\int_0^{\infty} e^{-x}x^{n-1} \]

\[ \lim_{x \;\to \;\infty} x^ne^{-x}= \lim_{x \;\to \;\infty}\dfrac{x^n}{e^x} =\lim_{x \;\to \;\infty}\dfrac{nx^{n-1}}{e^x} \]

\[ = \lim_{x \;\to \;\infty} \bigg\{\dfrac{n(n-1)(n-2)\ldots 1}{e^x}x^0 \bigg\} \]

\[ = \dfrac{n(n-1)(n-2)\ldots 1}{\infty}=0 \]

\(\therefore \ I_n = n\ I_{n-1}= n(n-1) I_{n-2} =\ldots \ldots =n(n-1)(n-2)\ldots 1\ . I_0\)

\[ I_0=\int_0^{\infty} e^{-x}x^0dx=\big[-e^{-x}\big]_0^\infty=1 \]

\[ I_n=n(n-1)(n-2)\ldots 1=n! \]

\[ \boxed{\Gamma(n)=\int_0^\infty e^{-x}x^ndx= n!} \]

Euler found a continuous function of the real variable \(x\) which reduces to \(x!\) when \(x\) is a positive integer. Gamma function can be defined in a number of ways but it is convenient to define it is as a definite integral. Mathematicians and enthusiasts love this function because it generalizes the factorial. It is often used in probability and statistics, as it shows up in the normalizing constants of important probability distributions such as the Chi-square and the Gamma. The motivation behind the gamma function was to find a smooth curve that connects the points \((x, y)\) given by \(y = (x − 1)!\) at the positive integer values for \(x\).

This article is a short introduction about Eulerian Integral. Moreover, we will definitely explore lots of fascinating mathematical insights and problems which is highly motivate us for core understanding of improper integrals, probability distributions and their applications, and lots more insightful applications of mathematics.

References