Wallis' integrals

Last updated

In mathematics, and more precisely in analysis, the Wallis integrals constitute a family of integrals introduced by John Wallis.

Contents

Definition, basic properties

The Wallis integrals are the terms of the sequence defined by

or equivalently,

The first few terms of this sequence are:

...
...

The sequence is decreasing and has positive terms. In fact, for all

Since the sequence is decreasing and bounded below by 0, it converges to a non-negative limit. Indeed, the limit is zero (see below).

Recurrence relation

By means of integration by parts, a reduction formula can be obtained. Using the identity , we have for all ,

Integrating the second integral by parts, with:

  • , whose anti-derivative is
  • , whose derivative is

we have:

Substituting this result into equation (1) gives

and thus

for all

This is a recurrence relation giving in terms of . This, together with the values of and give us two sets of formulae for the terms in the sequence , depending on whether is odd or even:

Another relation to evaluate the Wallis' integrals

Wallis's integrals can be evaluated by using Euler integrals:

  1. Euler integral of the first kind: the Beta function:
    for Re(x), Re(y) > 0
  2. Euler integral of the second kind: the Gamma function:
    for Re(z) > 0.

If we make the following substitution inside the Beta function:
we obtain:

so this gives us the following relation to evaluate the Wallis integrals:

So, for odd , writing , we have:

whereas for even , writing and knowing that , we get :

Equivalence

(equivalence of two sequences).
Indeed, for all  :
(since the sequence is decreasing)
(since )
(by equation ).
By the sandwich theorem, we conclude that , and hence .
(and consequently ).
Proof

For all , let .

It turns out that, because of equation . In other words is a constant.

It follows that for all , .

Now, since and , we have, by the product rules of equivalents, .

Thus, , from which the desired result follows (noting that ).


Deducing Stirling's formula

Suppose that we have the following equivalence (known as Stirling's formula):

for some constant that we wish to determine. From above, we have

(equation (3))

Expanding and using the formula above for the factorials, we get

From (3) and (4), we obtain by transitivity:

Solving for gives In other words,

Deducing the Double Factorial Ratio

Similarly, from above, we have:

Expanding and using the formula above for double factorials, we get:

Simplifying, we obtain:

or

Evaluating the Gaussian Integral

The Gaussian integral can be evaluated through the use of Wallis' integrals.

We first prove the following inequalities:

In fact, letting , the first inequality (in which ) is equivalent to ; whereas the second inequality reduces to , which becomes . These 2 latter inequalities follow from the convexity of the exponential function (or from an analysis of the function ).

Letting and making use of the basic properties of improper integrals (the convergence of the integrals is obvious), we obtain the inequalities:

for use with the sandwich theorem (as ).

The first and last integrals can be evaluated easily using Wallis' integrals. For the first one, let (t varying from 0 to ). Then, the integral becomes . For the last integral, let (t varying from to ). Then, it becomes .

As we have shown before, . So, it follows that .

Remark: There are other methods of evaluating the Gaussian integral. Some of them are more direct.

Note

The same properties lead to Wallis product, which expresses (see ) in the form of an infinite product.

Related Research Articles

<span class="mw-page-title-main">Bessel function</span> Families of solutions to related differential equations

Bessel functions, first defined by the mathematician Daniel Bernoulli and then generalized by Friedrich Bessel, are canonical solutions y(x) of Bessel's differential equation

<span class="mw-page-title-main">Gamma function</span> Extension of the factorial function

In mathematics, the gamma function is one commonly used extension of the factorial function to complex numbers. The gamma function is defined for all complex numbers except the non-positive integers. For every positive integer n,

In calculus, and more generally in mathematical analysis, integration by parts or partial integration is a process that finds the integral of a product of functions in terms of the integral of the product of their derivative and antiderivative. It is frequently used to transform the antiderivative of a product of functions into an antiderivative for which a solution can be more easily found. The rule can be thought of as an integral version of the product rule of differentiation.

Integration is the basic operation in integral calculus. While differentiation has straightforward rules by which the derivative of a complicated function can be found by differentiating its simpler component functions, integration does not, so tables of known integrals are often useful. This page lists some of the most common antiderivatives.

<span class="mw-page-title-main">Trigonometric integral</span> Special function defined by an integral

In mathematics, trigonometric integrals are a family of integrals involving trigonometric functions.

<span class="mw-page-title-main">Beta function</span> Mathematical function

In mathematics, the beta function, also called the Euler integral of the first kind, is a special function that is closely related to the gamma function and to binomial coefficients. It is defined by the integral

<span class="mw-page-title-main">Fresnel integral</span> Special function defined by an integral

The Fresnel integralsS(x) and C(x) are two transcendental functions named after Augustin-Jean Fresnel that are used in optics and are closely related to the error function (erf). They arise in the description of near-field Fresnel diffraction phenomena and are defined through the following integral representations:

<span class="mw-page-title-main">Trigonometric substitution</span> Technique of integral evaluation

In mathematics, trigonometric substitution is the replacement of trigonometric functions for other expressions. In calculus, trigonometric substitution is a technique for evaluating integrals. Moreover, one may use the trigonometric identities to simplify certain integrals containing radical expressions. Like other methods of integration by substitution, when evaluating a definite integral, it may be simpler to completely deduce the antiderivative before applying the boundaries of integration.

<span class="mw-page-title-main">Airy function</span> Special function in the physical sciences

In the physical sciences, the Airy function (or Airy function of the first kind) Ai(x) is a special function named after the British astronomer George Biddell Airy (1801–1892). The function Ai(x) and the related function Bi(x), are linearly independent solutions to the differential equation

<span class="mw-page-title-main">Digamma function</span> Mathematical function

In mathematics, the digamma function is defined as the logarithmic derivative of the gamma function:

<span class="mw-page-title-main">Gaussian integral</span> Integral of the Gaussian function, equal to sqrt(π)

The Gaussian integral, also known as the Euler–Poisson integral, is the integral of the Gaussian function over the entire real line. Named after the German mathematician Carl Friedrich Gauss, the integral is

<span class="mw-page-title-main">Sinc function</span> Special mathematical function defined as sin(x)/x

In mathematics, physics and engineering, the sinc function, denoted by sinc(x), has two forms, normalized and unnormalized.

In the mathematical field of complex analysis, contour integration is a method of evaluating certain integrals along paths in the complex plane.

<span class="mw-page-title-main">Stieltjes constants</span>

In mathematics, the Stieltjes constants are the numbers that occur in the Laurent series expansion of the Riemann zeta function:

In calculus, the Leibniz integral rule for differentiation under the integral sign states that for an integral of the form

<span class="mw-page-title-main">Carl Johan Malmsten</span>

Carl Johan Malmsten was a Swedish mathematician and politician. He is notable for early research into the theory of functions of a complex variable, for the evaluation of several important logarithmic integrals and series, for his studies in the theory of Zeta-function related series and integrals, as well as for helping Mittag-Leffler start the journal Acta Mathematica. Malmsten became Docent in 1840, and then, Professor of mathematics at the Uppsala University in 1842. He was elected a member of the Royal Swedish Academy of Sciences in 1844. He was also a minister without portfolio in 1859–1866 and Governor of Skaraborg County in 1866–1879.

Volume of an <i>n</i>-ball Size of a mathematical ball

In geometry, a ball is a region in a space comprising all points within a fixed distance, called the radius, from a given point; that is, it is the region enclosed by a sphere or hypersphere. An n-ball is a ball in an n-dimensional Euclidean space. The volume of a n-ball is the Lebesgue measure of this ball, which generalizes to any dimension the usual volume of a ball in 3-dimensional space. The volume of a n-ball of radius R is where is the volume of the unit n-ball, the n-ball of radius 1.

In mathematics, a Borwein integral is an integral whose unusual properties were first presented by mathematicians David Borwein and Jonathan Borwein in 2001. Borwein integrals involve products of , where the sinc function is given by for not equal to 0, and .