Integration by substitution

Integration by substitution
Topics in Calculus
Fundamental theorem
Limits of functions
Mean value theorem

In calculus, integration by substitution is a method for finding antiderivatives and integrals. Using the fundamental theorem of calculus often requires finding an antiderivative. For this and other reasons, integration by substitution is an important tool for mathematicians. It is the counterpart to the chain rule of differentiation.

Let   I\subseteq {\mathbb{R}}   be an interval and   g : [a,b] \to I   be a continuously differentiable function. Suppose that f : I \to \mathbb{R} is a continuous function. Then

\int_a^b f(g(t))g'(t)\, dt = \int_{g(a)}^{g(b)} f(x)\,dx.

Using Leibniz notation: the substitution x = g(t) yields dx / dt = g'(t) and thus, formally, dx = g'(t)\,dt, which is the required substitution for dx. (One could view the method of integration by substitution as a major justification of Leibniz's notation for integrals and derivatives.)

The formula is used to transform one integral into another integral that is easier to compute. Thus, the formula can be used from left to right or from right to left in order to simplify a given integral. When used in the former manner, it is sometimes known as u-substitution.


Relation to the fundamental theorem of calculus

Integration by substitution can be derived from the fundamental theorem of calculus as follows. Let ƒ and g be two functions satisfying the above hypothesis that ƒ is continuous on I and g'\, is continuous on the closed interval [a,b]. Then the function f(g(t))g'(t) is also continuous on [a,b]. Hence the integrals

\int_{g(a)}^{g(b)} f(x)\,dx


\int_a^b f(g(t))g'(t)\,dt

in fact exist, and it remains to show that they are equal.

Since ƒ is continuous, it possesses an antiderivative F. The composite function F\circ g is then defined. Since F and g are differentiable, the chain rule gives

(F \circ g)'(t) = F'(g(t))g'(t) = f(g(t))g'(t).

Applying the fundamental theorem of calculus twice gives

\int_a^b f(g(t))g'(t)\,dt & {} = (F \circ g)(b) - (F \circ g)(a) \\
& {} = F(g(b)) - F(g(a)) \\
& {} = \int_{g(a)}^{g(b)} f(x)\,dx,

which is the substitution rule.


Consider the integral

\int_{0}^2 x \cos(x^2+1) \,dx

If we make the substitution u = x2 + 1, we obtain du = 2x dx and

\int_{x=0}^{x=2} x \cos(x^2+1) \,dx & {} = \frac{1}{2} \int_{u=1}^{u=5}\cos(u)\,du \\
& {} = \frac{1}{2}(\sin(5)-\sin(1)).

Here we substituted from right to left. It is important to note that since the lower limit x = 0 was replaced with u = 02 + 1 = 1, and the upper limit x = 2 replaced with u = 22 + 1 = 5, a transformation back into terms of x was unnecessary.

For the integral

\int_0^1 \sqrt{1-x^2}\; dx

the formula needs to be used from left to right: the substitution x = sin(u), dx = cos(udu is useful, because \sqrt{(1-\sin^2(u))} = \cos(u):

\int_0^1 \sqrt{1-x^2}\; dx = \int_0^\frac{\pi}{2} \sqrt{1-\sin^2(u)} \cos(u)\;du = \int_0^\frac{\pi}{2} \cos^2(u)\;du=\frac{\pi}{4}

The resulting integral can be computed using integration by parts or a double angle formula followed by one more substitution. One can also note that the function being integrated is the upper right quarter of a circle with a radius of one, and hence integrating the upper right quarter from zero to one is the geometric equivalent to the area of one quarter of the unit circle, or pi over 4.


Substitution can be used to determine antiderivatives. One chooses a relation between x and u, determines the corresponding relation between dx and du by differentiating, and performs the substitutions. An antiderivative for the substituted function can hopefully be determined; the original substitution between u and x is then undone.

Similar to our first example above, we can determine the following antiderivative with this method:

& {} \quad \int x \cos(x^2+1) \,dx = \frac{1}{2} \int 2x \cos(x^2+1) \,dx \\
& {} = \frac{1}{2} \int\cos u\,du = \frac{1}{2}\sin u + C = \frac{1}{2}\sin(x^2+1) + C

where C is an arbitrary constant of integration.

Note that there were no integral boundaries to transform, but in the last step we had to revert the original substitution u = x2 + 1.

Substitution for multiple variables

One may also use substitution when integrating functions of several variables. Here the substitution function (v1,...,vn) = φ(u1, ..., un ) needs to be one-to-one and continuously differentiable, and the differentials transform as

dv_1\cdots dv_n = |\det(\operatorname{D}\varphi)(u_1, \ldots, u_n)| \, du_1\cdots du_n

where det(Dφ)(u1, ..., un ) denotes the determinant of the Jacobian matrix containing the partial derivatives of φ . This formula expresses the fact that the absolute value of the determinant of given vectors equals the volume of the spanned parallelotope.

More precisely, the change of variables formula is stated in the next theorem:

Theorem. Let U be an open set in Rn and φ : URn an injective differentiable function with continuous partial derivatives, the Jacobian of which is nonzero for every x in U. Then for any real-valued, compactly supported, continuous function f, with support contained in φ(U),

 \int_{\varphi(U)} f(\mathbf{v})\, d \mathbf{v} = \int_U f(\varphi(\mathbf{u})) \left|\det(\operatorname{D}\varphi)(\mathbf{u})\right| \,d \mathbf{u}.

The conditions on the theorem can be weakened in various ways. First, the requirement that φ be continuously differentiable can be replaced by the weaker assumption that φ be merely differentiable and have a continuous inverse (Rudin 1987, Theorem 7.26). This is guaranteed to hold if φ is continuously differentiable by the inverse function theorem. Alternatively, the requirement that Det(Dφ)≠0 can be eliminated by applying Sard's theorem (Spivak 1965).

For Lebesgue measurable functions, the theorem can be stated in the following form (Fremlin 2010, Theorem 263D):

Theorem. Let U be a measurable subset of Rn and φ : URn an injective function, and suppose for every x in U there exists φ'(x) in Rn,n such that φ(y) = φ(x) + φ'(x) (yx) + o(||yx||) as yx. Then φ(U) is measurable, and for any real-valued function f defined on φ(U),

 \int_{\varphi(U)} f(v)\, dv \;=\; \int_U f(\varphi(u)) \; \left|\det \varphi'(u)\right| \,du

in the sense that if either integral exists (or is properly infinite), then so does the other one, and they have the same value.

Another very general version in measure theory is the following (Hewitt & Stromberg 1965, Theorem 20.3):

Theorem. Let X be a locally compact Hausdorff space equipped with a finite Radon measure μ, and let Y be a σ-compact Hausdorff space with a σ-finite Radon measure ρ. Let φ : X → Y be a continuous and absolutely continuous function (where the latter means that ρ(φ(E)) = 0 whenever μ(E) = 0). Then there exists a real-valued Borel measurable function w on X such that for every Lebesgue integrable function f : Y → R, the function (f °φ)w is Lebesgue integrable on X, and

\int_Y f(y)\,d\rho(y) = \int_X f\circ\varphi(x)w(x)\,d\mu(x).

Furthermore, it is possible to write

w(x) = g\circ\varphi(x)

for some Borel measurable function g on Y.

In geometric measure theory, integration by substitution is used with Lipschitz functions. A bi-Lipschitz function is a Lipschitz function T : URn which is one-to-one, and such that its inverse function T-1 T(U) → U is also Lipschitz. By Rademacher's theorem a bi-Lipschitz mapping is differentiable almost everywhere. In particular, the Jacobian determinant of a bi-Lipschitz mapping det DT is well-defined almost everywhere. The following result then holds:

Theorem. Let U be an open subset of Rn and T : URn be a bi-Lipschitz mapping. Let f : T(U) → R be measurable. Then

\int_U (f\circ T)|\det DT| = \int_{T(U)}f

in the sense that if either integral exists (or is properly infinite), then so does the other one, and they have the same value.

The above theorem was first proposed by Euler when he developed the notion of double integrals in 1769. Although generalized to triple integrals by Lagrange in 1773, and used by Legendre, Laplace, Gauss, and first generalized to n variables by Mikhail Ostrogradski in 1836, it resisted a fully rigorous formal proof for a surprisingly long time, and was first satisfactorily resolved 125 years later, by Élie Cartan in a series of papers beginning in the mid-1890s (Katz 1982; Ferzola 1994).

Application in probability

Substitution can be used to answer the following important question in probability: given a random variable X with probability density px and another random variable Y related to X by the equation y = Φ(x), what is the probability density for Y?

It is easiest to answer this question by first answering a slightly different question: what is the probability that Y takes a value in some particular subset S? Denote this probability P(Y \in S). Of course, if Y has probability density py then the answer is

P(Y \in S) = \int_S p_y(y)\,dy,

but this isn't really useful because we don't know py; it's what we're trying to find in the first place. We can make progress by considering the problem in the variable X. Y takes a value in S whenever X takes a value in Φ − 1(S), so

 P(Y \in S) = \int_{\Phi^{-1}(S)} p_x(x)\,dx.

Changing from variable x to y gives

P(Y \in S) = \int_{\Phi^{-1}(S)} p_x(x)~dx = \int_S p_x(\Phi^{-1}(y)) ~ \left|\frac{d\Phi^{-1}}{dy}\right|~dy.

Combining this with our first equation gives

\int_S p_y(y)~dy = \int_S p_x(\Phi^{-1}(y)) ~ \left|\frac{d\Phi^{-1}}{dy}\right|~dy


p_y(y) = p_x(\Phi^{-1}(y)) ~ \left|\frac{d\Phi^{-1}}{dy}\right|.

In the case where X and Y depend on several uncorrelated variables, i.e. p_x=p_x(x_1\ldots x_n), and y = Φ(x), py can be found by substitution in several variables discussed above. The result is

p_y(y) = p_x(\Phi^{-1}(y)) ~ \left|\det \left[ D\Phi ^{-1}(y) \right] \right|.

See also


Wikimedia Foundation. 2010.

Look at other dictionaries:

  • Integration durch Substitution — Die Integration durch Substitution oder Substitutionsregel ist eine wichtige Methode in der Integralrechnung, um Stammfunktionen und Integrale zu berechnen. Durch Einführung einer neuen Integrationsvariablen wird ein Teil des Integranden ersetzt …   Deutsch Wikipedia

  • Integration by reduction formulae — can be used when we want to integrate a function raised to the power n. If we have such an integral we can establish a reduction formula which can be used to calculate the integral for any value of n. How to find the reduction formula The… …   Wikipedia

  • Substitution — (von lat. substituere ‚ersetzen‘) steht für: Mathematik/Informatik Substitution (Mathematik), das Ersetzen eines Ausdrucks Integration durch Substitution, eine Methode zur Lösung von Integralen Substitution (Logik), in der Logik die Ersetzung… …   Deutsch Wikipedia

  • Integration by parts — Topics in Calculus Fundamental theorem Limits of functions Continuity Mean value theorem Differential calculus  Derivative Change of variables Implicit differentiation Taylor s theorem Related rates …   Wikipedia

  • Substitution of variables — In mathematics, substitution of variables (also called variable substitution or coordinate transformation) refers to the substitution of certain variables with other variables. Though the study of how variable substitutions affect a certain… …   Wikipedia

  • Substitution (Mathematik) — Unter Substitution versteht man allgemein das Ersetzen eines Terms durch einen anderen. Dies wird üblicherweise angewandt um den Ausdruck, der den Term enthält, zu vereinfachen beziehungsweise in eine Standardform zu überführen, wie zum Beispiel… …   Deutsch Wikipedia

  • Integration horizontale — Intégration horizontale En économie, l intégration horizontale (ou concentration horizontale) consiste pour une entreprise à étendre son réseau, en acquérant ou développant des activités économiques au même niveau de la chaîne de valeur que ses… …   Wikipédia en Français

  • Intégration Horizontale — En économie, l intégration horizontale (ou concentration horizontale) consiste pour une entreprise à étendre son réseau, en acquérant ou développant des activités économiques au même niveau de la chaîne de valeur que ses produits. Les… …   Wikipédia en Français

  • Intégration horizontale — Pour les articles homonymes, voir Intégration. En management, l intégration horizontale (ou concentration horizontale) consiste pour une entreprise à étendre son réseau, en acquérant ou développant des activités économiques au même niveau de la… …   Wikipédia en Français

  • Integration of the normal density function — Main article: Normal distribution The probability density function for the normal distribution is given by:f(x;mu,sigma)=frac{1}{sigmasqrt{2pi , exp left( frac{(x mu)^2}{2sigma^2} ight),where mu is the mean and sigma the standard deviation.By the …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”

We are using cookies for the best presentation of our site. Continuing to use this site, you agree with this.