# Integration by substitution

Last updated

In calculus, integration by substitution, also known as u-substitution or change of variables, [1] is a method for evaluating integrals and antiderivatives. It is the counterpart to the chain rule for differentiation, and can loosely be thought of as using the chain rule "backwards".

## Substitution for a single variable

### Introduction

Before stating the result rigorously, consider a simple case using indefinite integrals.

Compute ${\displaystyle \textstyle \int (2x^{3}+1)^{7}(x^{2})\,dx}$. [2]

Set ${\displaystyle u=2x^{3}+1}$. This means ${\displaystyle \textstyle {\frac {du}{dx}}=6x^{2}}$, or in differential form, ${\displaystyle du=6x^{2}\,dx}$. Now

${\displaystyle \int (2x^{3}+1)^{7}(x^{2})\,dx={\frac {1}{6}}\int \underbrace {(2x^{3}+1)^{7}} _{u^{7}}\underbrace {(6x^{2})\,dx} _{du}={\frac {1}{6}}\int u^{7}\,du={\frac {1}{6}}\left({\frac {1}{8}}u^{8}\right)+C={\frac {1}{48}}(2x^{3}+1)^{8}+C}$,

where ${\displaystyle C}$ is an arbitrary constant of integration.

This procedure is frequently used, but not all integrals are of a form that permits its use. In any event, the result should be verified by differentiating and comparing to the original integrand.

${\displaystyle {\frac {d}{dx}}\left[{\frac {1}{48}}(2x^{3}+1)^{8}+C\right]={\frac {1}{6}}(2x^{3}+1)^{7}(6x^{2})=(2x^{3}+1)^{7}(x^{2}).}$

For definite integrals, the limits of integration must also be adjusted, but the procedure is mostly the same.

### Definite integrals

Let φ : [a, b] → I be a differentiable function with a continuous derivative, where IR is an interval. Suppose that f : IR is a continuous function. Then [3]

${\displaystyle \int _{a}^{b}f(\varphi (x))\varphi '(x)\,dx=\int _{\varphi (a)}^{\varphi (b)}f(u)\,du.}$

In Leibniz notation, the substitution u = φ(x) yields

${\displaystyle {\frac {du}{dx}}=\varphi '(x).}$

Working heuristically with infinitesimals yields the equation

${\displaystyle du=\varphi '(x)\,dx,}$

which suggests the substitution formula above. (This equation may be put on a rigorous foundation by interpreting it as a statement about differential forms.) One may view the method of integration by substitution as a partial 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 read 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 or w-substitution in which a new variable is defined to be a function of the original variable found inside the composite function multiplied by the derivative of the inner function. The latter manner is commonly used in trigonometric substitution, replacing the original variable with a trigonometric function of a new variable and the original differential with the differential of the trigonometric function.

### Proof

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

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

and

${\displaystyle \int _{\varphi (a)}^{\varphi (b)}f(u)\,du}$

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

Since f is continuous, it has an antiderivative F. The composite function Fφ is then defined. Since φ is differentiable, combining the chain rule and the definition of an antiderivative gives

${\displaystyle (F\circ \varphi )'(x)=F'(\varphi (x))\varphi '(x)=f(\varphi (x))\varphi '(x).}$

Applying the fundamental theorem of calculus twice gives

{\displaystyle {\begin{aligned}\int _{a}^{b}f(\varphi (x))\varphi '(x)\,dx&=\int _{a}^{b}(F\circ \varphi )'(x)\,dx\\&=(F\circ \varphi )(b)-(F\circ \varphi )(a)\\&=F(\varphi (b))-F(\varphi (a))\\&=\int _{\varphi (a)}^{\varphi (b)}f(u)\,du,\end{aligned}}}

which is the substitution rule.

### Examples

#### Example 1

Consider the integral

${\displaystyle \int _{0}^{2}x\cos(x^{2}+1)dx.}$

Make the substitution ${\displaystyle u=x^{2}+1}$ to obtain ${\displaystyle du=2xdx}$, meaning ${\displaystyle \textstyle xdx={\frac {1}{2}}du}$. Therefore,

{\displaystyle {\begin{aligned}\int _{x=0}^{x=2}x\cos(x^{2}+1)dx&={\frac {1}{2}}\int _{u=1}^{u=5}\cos(u)\,du\\[6pt]&={\frac {1}{2}}(\sin(5)-\sin(1)).\end{aligned}}}

Since the lower limit ${\displaystyle x=0}$ was replaced with ${\displaystyle u=1}$, and the upper limit ${\displaystyle x=2}$ with ${\displaystyle 2^{2}+1=5}$, a transformation back into terms of ${\displaystyle x}$ was unnecessary.

Alternatively, one may fully evaluate the indefinite integral (see below) first then apply the boundary conditions. This becomes especially handy when multiple substitutions are used.

#### Example 2

For the integral

${\displaystyle \int _{0}^{1}{\sqrt {1-x^{2}}}\,dx,}$

a variation of the above procedure is needed. The substitution ${\displaystyle x=\sin u}$ implying ${\displaystyle dx=\cos u\,du}$ is useful because ${\displaystyle {\sqrt {1-\sin ^{2}u}}=\cos(u)}$. We thus have

{\displaystyle {\begin{aligned}\int _{0}^{1}{\sqrt {1-x^{2}}}\,dx&=\int _{0}^{\pi /2}{\sqrt {1-\sin ^{2}u}}\cos(u)\,du\\[6pt]&=\int _{0}^{\pi /2}\cos ^{2}u\,du\\[6pt]&=\left[{\frac {u}{2}}+{\frac {\sin(2u)}{4}}\right]{\Biggl |}_{0}^{\pi /2}\\[6pt]&={\frac {\pi }{4}}+0={\frac {\pi }{4}}.\end{aligned}}}

The resulting integral can be computed using integration by parts or a double angle formula, ${\displaystyle 2\cos ^{2}u=1+\cos(2u)}$, 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 ${\displaystyle \pi /4}$.

### Antiderivatives

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

Similar to example 1 above, the following antiderivative can be obtained with this method:

{\displaystyle {\begin{aligned}\int x\cos(x^{2}+1)\,dx&={\frac {1}{2}}\int 2x\cos(x^{2}+1)\,dx\\[6pt]&={\frac {1}{2}}\int \cos u\,du\\[6pt]&={\frac {1}{2}}\sin u+C={\frac {1}{2}}\sin(x^{2}+1)+C,\end{aligned}}}

where ${\displaystyle C}$ is an arbitrary constant of integration.

There were no integral boundaries to transform, but in the last step reverting the original substitution ${\displaystyle u=x^{2}+1}$ was necessary. When evaluating definite integrals by substitution, one may calculate the antiderivative fully first, then apply the boundary conditions. In that case, there is no need to transform the boundary terms.

The tangent function can be integrated using substitution by expressing it in terms of the sine and cosine:

${\displaystyle \int \tan x\,dx=\int {\frac {\sin x}{\cos x}}\,dx}$

Using the substitution ${\displaystyle u=\cos x}$ gives ${\displaystyle du=-\sin x\,dx}$ and

{\displaystyle {\begin{aligned}\int \tan x\,dx&=\int {\frac {\sin x}{\cos x}}\,dx\\&=\int -{\frac {du}{u}}\\&=-\ln |u|+C\\&=-\ln |\cos x|+C=\ln |\sec x|+C.\end{aligned}}}

## 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 injective and continuously differentiable, and the differentials transform as

${\displaystyle dv_{1}\cdots dv_{n}=\left|\det(D\varphi )(u_{1},\ldots ,u_{n})\right|\,du_{1}\cdots du_{n},}$

where det()(u1, ..., un) denotes the determinant of the Jacobian matrix of partial derivatives of φ at the point (u1, ..., un). This formula expresses the fact that the absolute value of the determinant of a matrix equals the volume of the parallelotope spanned by its columns or rows.

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),

${\displaystyle \int _{\varphi (U)}f(\mathbf {v} )\,d\mathbf {v} =\int _{U}f(\varphi (\mathbf {u} ))\left|\det(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. [4] This is guaranteed to hold if φ is continuously differentiable by the inverse function theorem. Alternatively, the requirement that det() ≠ 0 can be eliminated by applying Sard's theorem. [5]

For Lebesgue measurable functions, the theorem can be stated in the following form: [6]

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 (here o is little-o notation). Then φ(U) is measurable, and for any real-valued function f defined on φ(U),

${\displaystyle \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 (including the possibility of being properly infinite), then so does the other one, and they have the same value.

Another very general version in measure theory is the following: [7]

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 φ : XY 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 : YR, the function (fφ) ⋅ w is Lebesgue integrable on X, and

${\displaystyle \int _{Y}f(y)\,d\rho (y)=\int _{X}(f\circ \varphi )(x)\,w(x)\,d\mu (x).}$

Furthermore, it is possible to write

${\displaystyle 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 φ : URn which is injective and whose inverse function φ1 : φ(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 is well-defined almost everywhere. The following result then holds:

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

${\displaystyle \int _{U}(f\circ \varphi )(x)|\det D\varphi (x)|\,dx=\int _{\varphi (U)}f(x)\,dx}$

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. [8] [9]

## Application in probability

Substitution can be used to answer the following important question in probability: given a random variable ${\displaystyle X}$ with probability density ${\displaystyle p_{X}}$ and another random variable ${\displaystyle Y}$ such that ${\displaystyle Y=\phi (X)}$, what is the probability density for ${\displaystyle Y}$?

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

${\displaystyle P(Y\in S)=\int _{S}p_{Y}(y)\,dy,}$

but this isn't really useful because we don't know ${\displaystyle p_{Y}}$; it's what we're trying to find. We can make progress by considering the problem in the variable ${\displaystyle X}$. ${\displaystyle Y}$ takes a value in ${\displaystyle S}$ whenever ${\displaystyle X}$ takes a value in ${\displaystyle \phi ^{-1}(S)}$, so

${\displaystyle P(Y\in S)=\int _{\phi ^{-1}(S)}p_{X}(x)\,dx.}$

Changing from variable ${\displaystyle x}$ to ${\displaystyle y}$ gives

${\displaystyle 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

${\displaystyle \int _{S}p_{Y}(y)\,dy=\int _{S}p_{X}(\phi ^{-1}(y))\left|{\frac {d\phi ^{-1}}{dy}}\right|\,dy,}$

so

${\displaystyle p_{Y}(y)=p_{X}(\phi ^{-1}(y))\left|{\frac {d\phi ^{-1}}{dy}}\right|.}$

In the case where ${\displaystyle X}$ and ${\displaystyle Y}$ depend on several uncorrelated variables, i.e. ${\displaystyle p_{X}=p_{X}(x_{1},\ldots ,x_{n})}$ and ${\displaystyle y=\phi (x)}$, ${\displaystyle p_{Y}}$ can be found by substitution in several variables discussed above. The result is

${\displaystyle p_{Y}(y)=p_{X}(\phi ^{-1}(y))\left|\det D\phi ^{-1}(y)\right|.}$

## Notes

1. Swokowski 1983 , p. 257
2. Swokowsi 1983 , p. 258
3. Briggs & Cochran 2011 , pg.361
4. Rudin 1987 , Theorem 7.26
5. Spivak 1965 , p. 72
6. Fremlin 2010 , Theorem 263D
7. Hewitt & Stromberg 1965 , Theorem 20.3

## Related Research Articles

In mathematics, the Laplace transform, named after its inventor Pierre-Simon Laplace, is an integral transform that converts a function of a real variable to a function of a complex variable . The transform has many applications in science and engineering because it is a tool for solving differential equations. In particular, it transforms linear differential equations into algebraic equations and convolution into multiplication.

In mathematics, the polar coordinate system is a two-dimensional coordinate system in which each point on a plane is determined by a distance from a reference point and an angle from a reference direction. The reference point is called the pole, and the ray from the pole in the reference direction is the polar axis. The distance from the pole is called the radial coordinate, radial distance or simply radius, and the angle is called the angular coordinate, polar angle, or azimuth. Angles in polar notation are generally expressed in either degrees or radians.

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace who first studied its properties. This is often written as

In mathematics, the Dirac delta function, also known as the unit impulse symbol, is a generalized function or distribution over the real numbers, whose value is zero everywhere except at zero, and whose integral over the entire real line is equal to one. It can also be interpreted as a linear functional that maps every function to its value at zero, or as the weak limit of a sequence of bump functions, which are zero over most of the real line, with a tall spike at the origin. Bump functions are thus sometimes called "approximate" or "nascent" delta functions.

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.

In vector calculus, the Jacobian matrix of a vector-valued function of several variables is the matrix of all its first-order partial derivatives. When this matrix is square, that is, when the function takes the same number of variables as input as the number of vector components of its output, its determinant is referred to as the Jacobian determinant. Both the matrix and the determinant are often referred to simply as the Jacobian in literature.

In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.

In mathematics, the Mellin transform is an integral transform that may be regarded as the multiplicative version of the two-sided Laplace transform. This integral transform is closely connected to the theory of Dirichlet series, and is often used in number theory, mathematical statistics, and the theory of asymptotic expansions; it is closely related to the Laplace transform and the Fourier transform, and the theory of the gamma function and allied special functions.

The path integral formulation is a description in quantum mechanics that generalizes the action principle of classical mechanics. It replaces the classical notion of a single, unique classical trajectory for a system with a sum, or functional integral, over an infinity of quantum-mechanically possible trajectories to compute a quantum amplitude.

In probability theory, the Borel–Kolmogorov paradox is a paradox relating to conditional probability with respect to an event of probability zero. It is named after Émile Borel and Andrey Kolmogorov.

Arc length is the distance between two points along a section of a curve.

In mathematics, a volume element provides a means for integrating a function with respect to volume in various coordinate systems such as spherical coordinates and cylindrical coordinates. Thus a volume element is an expression of the form

In mathematics, a change of variables is a basic technique used to simplify problems in which the original variables are replaced with functions of other variables. The intent is that when expressed in new variables, the problem may become simpler, or equivalent to a better understood problem.

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

In mathematics, a multiple integral is a definite integral of a function of several real variables, for instance, f(x, y) or f(x, y, z). Integrals of a function of two variables over a region in are called double integrals, and integrals of a function of three variables over a region in are called triple integrals. For multiple integrals of a single-variable function, see the Cauchy formula for repeated integration.

The gradient theorem, also known as the fundamental theorem of calculus for line integrals, says that a line integral through a gradient field can be evaluated by evaluating the original scalar field at the endpoints of the curve. The theorem is a generalization of the fundamental theorem of calculus to any curve in a plane or space rather than just the real line.

In physics, the Green's function for Laplace's equation in three variables is used to describe the response of a particular type of physical system to a point source. In particular, this Green's function arises in systems that can be described by Poisson's equation, a partial differential equation (PDE) of the form

In differential calculus, there is no single uniform notation for differentiation. Instead, several different notations for the derivative of a function or variable have been proposed by different mathematicians. The usefulness of each notation varies with the context, and it is sometimes advantageous to use more than one notation in a given context. The most common notations for differentiation are listed below.

Common integrals in quantum field theory are all variations and generalizations of Gaussian integrals to the complex plane and to multiple dimensions. Other integrals can be approximated by versions of the Gaussian integral. Fourier integrals are also considered.

In calculus, the integral of the secant function can be evaluated using a variety of methods and there are multiple ways of expressing the antiderivative, all of which can be shown to be equivalent via trigonometric identities,

## References

• Briggs, William; Cochran, Lyle (2011), Calculus /Early Transcendentals (Single Variable ed.), Addison-Wesley, ISBN   978-0-321-66414-3
• Ferzola, Anthony P. (1994), "Euler and differentials", The College Mathematics Journal , 25 (2): 102–111, doi:10.2307/2687130
• Fremlin, D.H. (2010), Measure Theory, Volume 2, Torres Fremlin, ISBN   978-0-9538129-7-4 .
• Hewitt, Edwin; Stromberg, Karl (1965), Real and Abstract Analysis, Springer-Verlag, ISBN   978-0-387-04559-7 .
• Katz, V. (1982), "Change of variables in multiple integrals: Euler to Cartan", Mathematics Magazine , 55 (1): 3–11, doi:10.2307/2689856
• Rudin, Walter (1987), Real and Complex Analysis, McGraw-Hill, ISBN   978-0-07-054234-1 .
• Swokowski, Earl W. (1983), Calculus with analytic geometry (alternate ed.), Prindle, Weber & Schmidt, ISBN   0-87150-341-7
• Spivak, Michael (1965), Calculus on Manifolds, Westview Press, ISBN   978-0-8053-9021-6 .