Improper integral

The original definition of the Riemann integral does not apply to a function such as ${\displaystyle 1/{x^{2}}}$ on the interval [1, ∞), because in this case the domain of integration is unbounded. However, the Riemann integral can often be extended by continuity, by defining the improper integral instead as a limit

${\displaystyle \int _{1}^{\infty }{\frac {dx}{x^{2}}}=\lim _{b\to \infty }\int _{1}^{b}{\frac {dx}{x^{2}}}=\lim _{b\to \infty }\left(-{\frac {1}{b}}+{\frac {1}{1}}\right)=1.}$

The narrow definition of the Riemann integral also does not cover the function ${\textstyle 1/{\sqrt {x}}}$ on the interval [0, 1]. The problem here is that the integrand is unbounded in the domain of integration (the definition requires that both the domain of integration and the integrand be bounded). However, the improper integral does exist if understood as the limit

${\displaystyle \int _{0}^{1}{\frac {dx}{\sqrt {x}}}=\lim _{a\to 0^{+}}\int _{a}^{1}{\frac {dx}{\sqrt {x}}}=\lim _{a\to 0^{+}}\left(2-2{\sqrt {a}}\right)=2.}$

The improper integral
${\displaystyle \int _{0}^{\infty }{\frac {dx}{(x+1){\sqrt {x}}}}=\pi }$
has unbounded intervals for both domain and range.

Sometimes integrals may have two singularities where they are improper. Consider, for example, the function 1/((x + 1)√x) integrated from 0 to ∞ (shown right). At the lower bound, as x goes to 0 the function goes to ∞, and the upper bound is itself ∞, though the function goes to 0. Thus this is a doubly improper integral. Integrated, say, from 1 to 3, an ordinary Riemann sum suffices to produce a result of π/6. To integrate from 1 to ∞, a Riemann sum is not possible. However, any finite upper bound, say t (with t > 1), gives a well-defined result, 2 arctan(√t) − π/2. This has a finite limit as t goes to infinity, namely π/2. Similarly, the integral from 1/3 to 1 allows a Riemann sum as well, coincidentally again producing π/6. Replacing 1/3 by an arbitrary positive value s (with s < 1) is equally safe, giving π/2 − 2 arctan(√s). This, too, has a finite limit as s goes to zero, namely π/2. Combining the limits of the two fragments, the result of this improper integral is

${\displaystyle {\begin{aligned}\int _{0}^{\infty }{\frac {dx}{(1+x){\sqrt {x}}}}&{}=\lim _{s\to 0^{+}}\int _{s}^{1}{\frac {dx}{(1+x){\sqrt {x}}}}+\lim _{t\to \infty }\int _{1}^{t}{\frac {dx}{(1+x){\sqrt {x}}}}\\&{}=\lim _{s\to 0^{+}}\left({\frac {\pi }{2}}-2\arctan {\sqrt {s}}\right)+\lim _{t\to \infty }\left(2\arctan {\sqrt {t}}-{\frac {\pi }{2}}\right)\\&{}={\frac {\pi }{2}}+\left(\pi -{\frac {\pi }{2}}\right)\\&{}=\pi .\end{aligned}}}$

This process does not guarantee success; a limit might fail to exist, or might be infinite. For example, over the bounded interval from 0 to 1 the integral of 1/x does not converge; and over the unbounded interval from 1 to ∞ the integral of 1/√x does not converge.

The improper integral
${\displaystyle \int _{-1}^{1}{\frac {dx}{\sqrt[{3}]{x^{2}}}}=6}$
converges, since both left and right limits exist, though the integrand is unbounded near an interior point.

It might also happen that an integrand is unbounded near an interior point, in which case the integral must be split at that point. For the integral as a whole to converge, the limit integrals on both sides must exist and must be bounded. For example:

${\displaystyle {\begin{aligned}\int _{-1}^{1}{\frac {dx}{\sqrt[{3}]{x^{2}}}}&{}=\lim _{s\to 0^{-}}\int _{-1}^{s}{\frac {dx}{\sqrt[{3}]{x^{2}}}}+\lim _{t\to 0^{+}}\int _{t}^{1}{\frac {dx}{\sqrt[{3}]{x^{2}}}}\\&{}=\lim _{s\to 0^{-}}3\left(1-{\sqrt[{3}]{s}}\right)+\lim _{t\to 0^{+}}3\left(1-{\sqrt[{3}]{t}}\right)\\&{}=3+3\\&{}=6.\end{aligned}}}$

But the similar integral

${\displaystyle \int _{-1}^{1}{\frac {dx}{x}}}$

cannot be assigned a value in this way, as the integrals above and below zero do not independently converge. (However, see Cauchy principal value.)

An improper integral converges if the limit defining it exists. Thus for example one says that the improper integral

${\displaystyle \lim _{t\to \infty }\int _{a}^{t}f(x)\ dx}$

exists and is equal to L if the integrals under the limit exist for all sufficiently large t, and the value of the limit is equal to L.

It is also possible for an improper integral to diverge to infinity. In that case, one may assign the value of ∞ (or −∞) to the integral. For instance

${\displaystyle \lim _{b\to \infty }\int _{1}^{b}{\frac {dx}{x}}=\infty .}$

However, other improper integrals may simply diverge in no particular direction, such as

${\displaystyle \lim _{b\to \infty }\int _{1}^{b}x\sin(x)\,dx,}$

which does not exist, even as an extended real number. This is called divergence by oscillation.

A limitation of the technique of improper integration is that the limit must be taken with respect to one endpoint at a time. Thus, for instance, an improper integral of the form

${\displaystyle \int _{-\infty }^{\infty }f(x)\,dx}$

can be defined by taking two separate limits; to wit

${\displaystyle \int _{-\infty }^{\infty }f(x)\,dx=\lim _{a\to -\infty }\lim _{b\to \infty }\int _{a}^{b}f(x)\,dx}$

provided the double limit is finite. It can also be defined as a pair of distinct improper integrals of the first kind:

${\displaystyle \lim _{a\to -\infty }\int _{a}^{c}f(x)\,dx+\lim _{b\to \infty }\int _{c}^{b}f(x)\,dx}$

where c is any convenient point at which to start the integration. This definition also applies when one of these integrals is infinite, or both if they have the same sign.

An example of an improper integral where both endpoints are infinite is the Gaussian integral ${\textstyle \int _{-\infty }^{\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}}$. An example which evaluates to infinity is ${\textstyle \int _{-\infty }^{\infty }e^{x}\,dx}$. But one cannot even define other integrals of this kind unambiguously, such as ${\textstyle \int _{-\infty }^{\infty }x\,dx}$, since the double limit is infinite and the two-integral method

${\displaystyle \lim _{a\to -\infty }\int _{a}^{c}x\,dx+\lim _{b\to \infty }\int _{c}^{b}x\,dx}$

yields ${\displaystyle \infty -\infty }$. In this case, one can however define an improper integral in the sense of Cauchy principal value:

${\displaystyle \operatorname {p.v.} \int _{-\infty }^{\infty }x\,dx=\lim _{b\to \infty }\int _{-b}^{b}x\,dx=0.}$

The questions one must address in determining an improper integral are:

• Does the limit exist?
• Can the limit be computed?

The first question is an issue of mathematical analysis. The second one can be addressed by calculus techniques, but also in some cases by contour integration, Fourier transforms and other more advanced methods.

There is more than one theory of integration. From the point of view of calculus, the Riemann integral theory is usually assumed as the default theory. In using improper integrals, it can matter which integration theory is in play.

• For the Riemann integral (or the Darboux integral, which is equivalent to it), improper integration is necessary both for unbounded intervals (since one cannot divide the interval into finitely many subintervals of finite length) and for unbounded functions with finite integral (since, supposing it is unbounded above, then the upper integral will be infinite, but the lower integral will be finite).
• The Lebesgue integral deals differently with unbounded domains and unbounded functions, so that often an integral which only exists as an improper Riemann integral will exist as a (proper) Lebesgue integral, such as ${\textstyle \int _{1}^{\infty }{\frac {dx}{x^{2}}}}$. On the other hand, there are also integrals that have an improper Riemann integral but do not have a (proper) Lebesgue integral, such as ${\textstyle \int _{0}^{\infty }{\frac {\sin x}{x}}\,dx}$. The Lebesgue theory does not see this as a deficiency: from the point of view of measure theory, ${\textstyle \int _{0}^{\infty }{\frac {\sin x}{x}}\,dx=\infty -\infty }$ and cannot be defined satisfactorily. In some situations, however, it may be convenient to employ improper Lebesgue integrals as is the case, for instance, when defining the Cauchy principal value. The Lebesgue integral is more or less essential in the theoretical treatment of the Fourier transform, with pervasive use of integrals over the whole real line.
• For the Henstock–Kurzweil integral, improper integration is not necessary, and this is seen as a strength of the theory: it encompasses all Lebesgue integrable and improper Riemann integrable functions.

Figure 1

Figure 2

In some cases, the integral

${\displaystyle \int _{a}^{c}f(x)\ dx}$

can be defined as an integral (a Lebesgue integral, for instance) without reference to the limit

${\displaystyle \lim _{b\to c^{-}}\int _{a}^{b}f(x)\,dx}$

but cannot otherwise be conveniently computed. This often happens when the function f being integrated from a to c has a vertical asymptote at c, or if c = ∞ (see Figures 1 and 2). In such cases, the improper Riemann integral allows one to calculate the Lebesgue integral of the function. Specifically, the following theorem holds (Apostol 1974, Theorem 10.33):

• If a function f is Riemann integrable on [a,b] for every b ≥ a, and the partial integrals

${\displaystyle \int _{a}^{b}|f(x)|\,dx}$

are bounded as b → ∞, then the improper Riemann integrals

${\displaystyle \int _{a}^{\infty }f(x)\,dx,\quad {\mbox{and }}\int _{a}^{\infty }|f(x)|\,dx}$

both exist. Furthermore, f is Lebesgue integrable on [a, ∞), and its Lebesgue integral is equal to its improper Riemann integral.

For example, the integral

${\displaystyle \int _{0}^{\infty }{\frac {dx}{1+x^{2}}}}$

can be interpreted alternatively as the improper integral

${\displaystyle \lim _{b\to \infty }\int _{0}^{b}{\frac {dx}{1+x^{2}}}=\lim _{b\to \infty }\arctan {b}={\frac {\pi }{2}},}$

or it may be interpreted instead as a Lebesgue integral over the set (0, ∞). Since both of these kinds of integral agree, one is free to choose the first method to calculate the value of the integral, even if one ultimately wishes to regard it as a Lebesgue integral. Thus improper integrals are clearly useful tools for obtaining the actual values of integrals.

In other cases, however, a Lebesgue integral between finite endpoints may not even be defined, because the integrals of the positive and negative parts of f are both infinite, but the improper Riemann integral may still exist. Such cases are "properly improper" integrals, i.e. their values cannot be defined except as such limits. For example,

${\displaystyle \int _{0}^{\infty }{\frac {\sin(x)}{x}}\,dx}$

cannot be interpreted as a Lebesgue integral, since

${\displaystyle \int _{0}^{\infty }\left|{\frac {\sin(x)}{x}}\right|\,dx=\infty .}$

But ${\displaystyle f(x)={\frac {\sin(x)}{x}}}$ is nevertheless integrable between any two finite endpoints, and its integral between 0 and ∞ is usually understood as the limit of the integral:

${\displaystyle \int _{0}^{\infty }{\frac {\sin(x)}{x}}\,dx=\lim _{b\to \infty }\int _{0}^{b}{\frac {\sin(x)}{x}}\,dx={\frac {\pi }{2}}.}$

One can speak of the singularities of an improper integral, meaning those points of the extended real number line at which limits are used.

Consider the difference in values of two limits:

${\displaystyle \lim _{a\to 0^{+}}\left(\int _{-1}^{-a}{\frac {dx}{x}}+\int _{a}^{1}{\frac {dx}{x}}\right)=0,}$

${\displaystyle \lim _{a\to 0^{+}}\left(\int _{-1}^{-a}{\frac {dx}{x}}+\int _{2a}^{1}{\frac {dx}{x}}\right)=-\ln 2.}$

The former is the Cauchy principal value of the otherwise ill-defined expression

${\displaystyle \int _{-1}^{1}{\frac {dx}{x}}{\ }\left({\mbox{which}}\ {\mbox{gives}}\ -\infty +\infty \right).}$

Similarly, we have

${\displaystyle \lim _{a\to \infty }\int _{-a}^{a}{\frac {2x\,dx}{x^{2}+1}}=0,}$

but

${\displaystyle \lim _{a\to \infty }\int _{-2a}^{a}{\frac {2x\,dx}{x^{2}+1}}=-\ln 4.}$

The former is the principal value of the otherwise ill-defined expression

${\displaystyle \int _{-\infty }^{\infty }{\frac {2x\,dx}{x^{2}+1}}{\ }\left({\mbox{which}}\ {\mbox{gives}}\ -\infty +\infty \right).}$

All of the above limits are cases of the indeterminate form ∞ − ∞.

These pathologies do not affect "Lebesgue-integrable" functions, that is, functions the integrals of whose absolute values are finite.

An improper integral may diverge in the sense that the limit defining it may not exist. In this case, there are more sophisticated definitions of the limit which can produce a convergent value for the improper integral. These are called summability methods.

One summability method, popular in Fourier analysis, is that of Cesàro summation. The integral

${\displaystyle \int _{0}^{\infty }f(x)\,dx}$

is Cesàro summable (C, α) if

${\displaystyle \lim _{\lambda \to \infty }\int _{0}^{\lambda }\left(1-{\frac {x}{\lambda }}\right)^{\alpha }f(x)\ dx}$

exists and is finite (Titchmarsh 1948, §1.15). The value of this limit, should it exist, is the (C, α) sum of the integral.

An integral is (C, 0) summable precisely when it exists as an improper integral. However, there are integrals which are (C, α) summable for α > 0 which fail to converge as improper integrals (in the sense of Riemann or Lebesgue). One example is the integral

${\displaystyle \int _{0}^{\infty }\sin x\,dx}$

which fails to exist as an improper integral, but is (C,α) summable for every α > 0. This is an integral version of Grandi's series.

The improper integral can also be defined for functions of several variables. The definition is slightly different, depending on whether one requires integrating over an unbounded domain, such as ${\displaystyle \mathbb {R} ^{2}}$, or is integrating a function with singularities, like ${\displaystyle f(x,y)=\log \left(x^{2}+y^{2}\right)}$.

• Apostol, T (1974), Mathematical analysis, Addison-Wesley, ISBN 978-0-201-00288-1.
• Apostol, T (1967), Calculus, Vol. 1 (2nd ed.), Jon Wiley & Sons.
• Autar Kaw, Egwu Kalu (2008), Numerical Methods with Applications (1st ed.), autarkaw.com
• Titchmarsh, E (1948), Introduction to the theory of Fourier integrals (2nd ed.), New York, N.Y.: Chelsea Pub. Co. (published 1986), ISBN 978-0-8284-0324-5.
• Cooper, Jeffery (2005), Working analysis, Gulf Professional
• Ghorpade, Sudhir; Limaye, Balmohan (2010), A course in multivariable calculus and analysis, Springer

• Numerical Methods to Solve Improper Integrals at Holistic Numerical Methods Institute