Reference sheet

Integrals reference sheet

The main integrals, methods, theorems and conditions of the integrals course on one page, unit by unit, with links to the chapter pages that explain each topic in full.

This sheet follows the seven units of the course in order. Each unit starts with the big idea in a few sentences, then gives the antiderivatives, methods and theorems you will look up most often, each with the conditions it needs. Most topics end with the integral of the matching chapter worked as an example, and the links under each topic lead to the chapter pages with the full explanation.

Angles are in radians throughout. Definite integrals show their limits, as in ∫01x dx\int_0^1 x\,dx, and an antiderivative carries + C+\,C for an arbitrary constant on each interval where it holds. Improper integrals are defined as limits of ordinary integrals. A function is continuously differentiable on an interval when its derivative exists and is continuous there, with one-sided derivatives at the endpoints. Vectors are bold, as in F\mathbf{F} and r\mathbf{r}. A formula holds wherever both sides are defined unless a condition says otherwise.

Accumulation and foundations

A definite integral adds up rate times change in input across an interval. For a positive integrand it is the area under the graph. A limit of Riemann sums defines it, the fundamental theorem of calculus evaluates it with an antiderivative, and linearity splits a sum into separate signed parts.

The definite integral

In ∫abf(x) dx\int_a^b f(x)\,dx the function ff is the integrand, the numbers aa and bb are the limits of integration and xx is the variable of integration. The integral is the signed area between the graph and the xx-axis, where area above the axis counts as positive and area below as negative. Its units are the units of f(x)f(x) times the units of xx, so a rate in liters per second integrated over seconds gives liters.

For ff integrable on an interval that contains every limit in the row,

RuleFormula
Constant integrand∫abc dx=c (b−a)\int_a^b c\,dx = c\,(b - a)
Equal limits∫aaf(x) dx=0\int_a^a f(x)\,dx = 0
Reversed limits∫baf(x) dx=−∫abf(x) dx\int_b^a f(x)\,dx = -\int_a^b f(x)\,dx
Adjacent intervals∫abf(x) dx+∫bcf(x) dx=∫acf(x) dx\int_a^b f(x)\,dx + \int_b^c f(x)\,dx = \int_a^c f(x)\,dx
Positivity∫abf(x) dx≥0\int_a^b f(x)\,dx \ge 0 when f≥0f \ge 0 on [a,b][a, b] and a<ba < b

A flow of 22 liters per second for 33 seconds accumulates

∫032 dx=2 (3−0)=6\int_0^3 2\,dx = 2\,(3 - 0) = 6

liters, the area of a rectangle of height 22 and width 33.

Lessons What an integral measures

Riemann sums

A regular partition splits [a,b][a, b] into nn subintervals of equal width Δx=b−an\Delta x = \frac{b - a}{n}, with points xk=a+k Δxx_k = a + k\,\Delta x. The mesh width of a partition is the width of its widest subinterval, here Δx\Delta x. The right-endpoint and left-endpoint sums are

Rn=∑k=1nf(xk) Δx,Ln=∑k=0n−1f(xk) Δx.\begin{aligned} R_n &= \sum_{k=1}^{n} f(x_k)\,\Delta x, \\[0.4em] L_n &= \sum_{k=0}^{n-1} f(x_k)\,\Delta x. \end{aligned}

A function ff is integrable on [a,b][a, b] when its Riemann sums, with any sample point in each subinterval, approach one limit as the mesh width tends to 00. That limit is ∫abf(x) dx\int_a^b f(x)\,dx. Continuous functions are integrable, and so are bounded functions with finitely many discontinuities. For an increasing ff the integral lies between LnL_n and RnR_n, and the gap is Rn−Ln=(f(b)−f(a)) ΔxR_n - L_n = \big(f(b) - f(a)\big)\,\Delta x.

SumClosed form
1+2+⋯+n1 + 2 + \cdots + nn(n+1)2\dfrac{n(n + 1)}{2}
12+22+⋯+n21^2 + 2^2 + \cdots + n^2n(n+1)(2n+1)6\dfrac{n(n + 1)(2n + 1)}{6}
13+23+⋯+n31^3 + 2^3 + \cdots + n^3(n(n+1)2)2\left(\dfrac{n(n + 1)}{2}\right)^2

For f(x)=xf(x) = x on [0,1][0, 1] the points are xk=knx_k = \frac{k}{n}, and

Rn=1n2∑k=1nk=12+12n,\begin{aligned} R_n &= \frac{1}{n^2}\sum_{k=1}^{n} k \\[0.4em] &= \frac{1}{2} + \frac{1}{2n}, \end{aligned}

so R4=0.625R_4 = 0.625, R10=0.55R_{10} = 0.55 and R100=0.505R_{100} = 0.505. The left sums are Ln=12−12nL_n = \frac{1}{2} - \frac{1}{2n}, both tend to the same limit, and

∫01x dx=lim⁡n→∞Rn=12.\int_0^1 x\,dx = \lim_{n\to\infty} R_n = \frac{1}{2}.

Lessons The Riemann-sum limit

The fundamental theorem of calculus

An antiderivative of ff on an interval is a function FF with F′=fF' = f there. Any two antiderivatives on the same interval differ by a constant, and the indefinite integral ∫f(x) dx=F(x)+C\int f(x)\,dx = F(x) + C names the whole family. Let ff be continuous on [a,b][a, b].

Part 1. The accumulation function A(x)=∫axf(t) dtA(x) = \int_a^x f(t)\,dt is continuous on [a,b][a, b] and differentiable on (a,b)(a, b), with A′(x)=f(x)A'(x) = f(x).

Part 2. If FF is an antiderivative of ff on [a,b][a, b], then

∫abf(x) dx=[F(x)]ab=F(b)−F(a).\begin{aligned} \int_a^b f(x)\,dx &= \big[F(x)\big]_a^b \\[0.4em] &= F(b) - F(a). \end{aligned}

The constant of integration cancels in the difference. The reverse power rule reads the power rule backward,

∫xr dx=xr+1r+1+C,\int x^r\,dx = \frac{x^{r+1}}{r + 1} + C,

for r≠−1r \ne -1 on each interval where xrx^r is defined. For example F(x)=x33F(x) = \frac{x^3}{3} has F′(x)=x2F'(x) = x^2, so

∫02x2 dx=233−0=83.\int_0^2 x^2\,dx = \frac{2^3}{3} - 0 = \frac{8}{3}.

The theorem needs ff continuous on all of [a,b][a, b]. See the improper integrals topic for what goes wrong at a point where ff blows up.

Lessons The fundamental theorem

Standard antiderivatives

Each entry holds on every interval inside the stated set, and adding + C+\,C gives all antiderivatives there. Differentiating an answer checks it.

IntegrandAntiderivativeHolds for
xrx^rxr+1r+1\dfrac{x^{r+1}}{r + 1}r≠−1r \ne -1, where xrx^r is defined
1x\dfrac{1}{x}ln⁡∣x∣\ln\lvert x \rvertx>0x > 0 or x<0x < 0
ekxe^{kx}ekxk\dfrac{e^{kx}}{k}k≠0k \ne 0
axa^xaxln⁡a\dfrac{a^x}{\ln a}a>0a > 0, a≠1a \ne 1
sin⁡x\sin x−cos⁡x-\cos xall xx
cos⁡x\cos xsin⁡x\sin xall xx
sec⁡2x\sec^2 xtan⁡x\tan xcos⁡x≠0\cos x \ne 0
sec⁡xtan⁡x\sec x \tan xsec⁡x\sec xcos⁡x≠0\cos x \ne 0
tan⁡x\tan xln⁡∣sec⁡x∣\ln\lvert \sec x \rvertcos⁡x≠0\cos x \ne 0
sec⁡x\sec xln⁡∣sec⁡x+tan⁡x∣\ln\lvert \sec x + \tan x \rvertcos⁡x≠0\cos x \ne 0
ln⁡x\ln xxln⁡x−xx\ln x - xx>0x > 0
1a2+x2\dfrac{1}{a^2 + x^2}1aarctan⁡xa\dfrac{1}{a}\arctan\dfrac{x}{a}a>0a > 0
1a2−x2\dfrac{1}{\sqrt{a^2 - x^2}}arcsin⁡xa\arcsin\dfrac{x}{a}a>0a > 0, ∣x∣<a\lvert x \rvert < a
1x2+a2\dfrac{1}{\sqrt{x^2 + a^2}}ln⁡(x+x2+a2)\ln\big(x + \sqrt{x^2 + a^2}\big)a>0a > 0

Lessons The fundamental theorem, Integration by parts, Trigonometric substitution, Completing the square

Linearity and average value

For ff and gg integrable on [a,b][a, b] and constants α\alpha and β\beta,

∫ab(αf+βg) dx=α∫abf dx+β∫abg dx.\begin{aligned} &\int_a^b \big(\alpha f + \beta g\big)\,dx \\[0.4em] &\quad = \alpha\int_a^b f\,dx \\[0.4em] &\qquad + \beta\int_a^b g\,dx. \end{aligned}

The sum rule is the case α=β=1\alpha = \beta = 1 and the constant multiple rule the case β=0\beta = 0. Each term keeps its own sign, so a term can make a negative signed contribution even when the whole integrand is positive. The integral of a product is not the product of the integrals. For a<ba < b,

QuantityFormulaConditions
Average valuefˉ=1b−a∫abf(x) dx\bar{f} = \dfrac{1}{b - a}\displaystyle\int_a^b f(x)\,dxff integrable on [a,b][a, b]
Mean value theoremf(c)=fˉf(c) = \bar{f} for some cc in (a,b)(a, b)ff continuous on [a,b][a, b]
Boundsm(b−a)≤∫abf(x) dx≤M(b−a)m(b - a) \le \displaystyle\int_a^b f(x)\,dx \le M(b - a)m≤f≤Mm \le f \le M on [a,b][a, b]

For example the three terms of 3x2−4x+23x^2 - 4x + 2 contribute 88, −8-8 and 44 on [0,2][0, 2], so

∫02(3x2−4x+2) dx=4,\int_0^2 \big(3x^2 - 4x + 2\big)\,dx = 4,

although the integrand 3(x−23)2+233\left(x - \frac{2}{3}\right)^2 + \frac{2}{3} is positive everywhere. One antiderivative F(x)=x3−2x2+2xF(x) = x^3 - 2x^2 + 2x confirms F(2)−F(0)=4F(2) - F(0) = 4. The average value is 22, which the integrand takes at c=43c = \frac{4}{3}.

Lessons Linearity and signed contributions

Recognizing symbolic methods

Each method rewrites an integrand until it matches an entry of the table of standard antiderivatives. The work is mostly recognition, choosing the method from the shape of the integrand, and every answer can be checked by differentiating it.

Substitution

Substitution rule. If gg is continuously differentiable on [a,b][a, b] and ff is continuous on the set of values g([a,b])g([a, b]), then

∫abf(g(x)) g′(x) dx=∫g(a)g(b)f(u) du.\begin{aligned} &\int_a^b f\big(g(x)\big)\,g'(x)\,dx \\[0.4em] &\quad = \int_{g(a)}^{g(b)} f(u)\,du. \end{aligned}

The indefinite form is ∫f(g(x)) g′(x) dx=F(g(x))+C\int f\big(g(x)\big)\,g'(x)\,dx = F\big(g(x)\big) + C for an antiderivative FF of ff, the chain rule in reverse. Look for an inner function u=g(x)u = g(x) whose derivative also appears as a factor, write the differential du=g′(x) dxdu = g'(x)\,dx, and either change the bounds to g(a)g(a) and g(b)g(b) or substitute back before using the original bounds. When leftover factors of xx must be written in terms of uu, the substitution must be one-to-one on the interval. For example u=x2u = x^2 is one-to-one for x≥0x \ge 0, where x=ux = \sqrt{u}.

With u=x2u = x^2 and du=2x dxdu = 2x\,dx, the bounds 00 and π/2\sqrt{\pi/2} become 00 and π2\frac{\pi}{2}, and

∫0π/22xcos⁡(x2) dx=∫0π/2cos⁡u du=1.\begin{aligned} &\int_0^{\sqrt{\pi/2}} 2x\cos(x^2)\,dx \\[0.4em] &\quad = \int_0^{\pi/2} \cos u\,du = 1. \end{aligned}

The check is ddxsin⁡(x2)=2xcos⁡(x2)\frac{d}{dx}\sin(x^2) = 2x\cos(x^2).

Lessons Substitution

Integration by parts

Integration by parts. If uu and vv are continuously differentiable on [a,b][a, b], then

∫abu v′ dx=[u v]ab−∫abu′ v dx.\begin{aligned} \int_a^b u\,v'\,dx &= \big[u\,v\big]_a^b \\[0.4em] &\quad - \int_a^b u'\,v\,dx. \end{aligned}

The indefinite form is ∫u dv=uv−∫v du\int u\,dv = uv - \int v\,du. It is the product rule (uv)′=u′v+uv′(uv)' = u'v + uv' integrated and rearranged, and the boundary term [uv]ab[uv]_a^b must be kept. Choose uu as the factor that gets simpler when differentiated, such as a polynomial or ln⁡x\ln x, and dvdv as a factor you can integrate, such as ex dxe^x\,dx or sin⁡x dx\sin x\,dx.

IntegralChoiceResult
∫xex dx\int x e^x\,dxu=xu = x, dv=ex dxdv = e^x\,dx(x−1) ex+C(x - 1)\,e^x + C
∫ln⁡x dx\int \ln x\,dxu=ln⁡xu = \ln x, dv=dxdv = dxxln⁡x−x+Cx\ln x - x + C for x>0x > 0
∫0πxsin⁡x dx\int_0^\pi x\sin x\,dxu=xu = x, dv=sin⁡x dxdv = \sin x\,dxπ\pi
∫exsin⁡x dx\int e^x\sin x\,dxparts twice, then solve for the integral12ex(sin⁡x−cos⁡x)+C\frac{1}{2}e^x(\sin x - \cos x) + C

With u=xu = x and v=exv = e^x,

∫01xex dx=[xex]01−∫01ex dx=e−(e−1)=1.\begin{aligned} &\int_0^1 x e^x\,dx \\[0.4em] &\quad = \big[x e^x\big]_0^1 - \int_0^1 e^x\,dx \\[0.4em] &\quad = e - (e - 1) = 1. \end{aligned}

The same number is a complementary area. The rectangle [0,1]×[0,e][0, 1] \times [0, e] has area ee, the region under y=exy = e^x has area e−1e - 1, and the region between the curve and the yy-axis has area ∫1eln⁡y dy=1\int_1^e \ln y\,dy = 1. In general, for ff continuously differentiable and strictly increasing on [a,b][a, b],

∫abf(x) dx+∫f(a)f(b)f−1(y) dy=b f(b)−a f(a).\begin{aligned} &\int_a^b f(x)\,dx \\[0.4em] &\quad + \int_{f(a)}^{f(b)} f^{-1}(y)\,dy \\[0.4em] &\quad = b\,f(b) - a\,f(a). \end{aligned}

Lessons Integration by parts

Powers of sine and cosine

The cosine of a doubled angle, cos⁡2x=1−2sin⁡2x=2cos⁡2x−1\cos 2x = 1 - 2\sin^2 x = 2\cos^2 x - 1, gives the half-angle or power-reducing identities

sin⁡2x=1−cos⁡2x2,cos⁡2x=1+cos⁡2x2,\begin{aligned} \sin^2 x &= \frac{1 - \cos 2x}{2}, \\[0.6em] \cos^2 x &= \frac{1 + \cos 2x}{2}, \end{aligned}

and the antiderivatives that follow from them,

∫sin⁡2x dx=x2−sin⁡2x4+C,∫cos⁡2x dx=x2+sin⁡2x4+C.\begin{aligned} &\int \sin^2 x\,dx \\[0.4em] &\quad = \frac{x}{2} - \frac{\sin 2x}{4} + C, \\[0.8em] &\int \cos^2 x\,dx \\[0.4em] &\quad = \frac{x}{2} + \frac{\sin 2x}{4} + C. \end{aligned}

Here sin⁡2x\sin^2 x means (sin⁡x)2(\sin x)^2, not sin⁡(x2)\sin(x^2). Over a full period the oscillation cancels. For a nonzero integer kk,

∫02πcos⁡kx dx=0,∫02πsin⁡kx dx=0,∫02πcos⁡2kx dx=π,∫02πsin⁡2kx dx=π.\begin{aligned} \int_0^{2\pi} \cos kx\,dx &= 0, \\[0.4em] \int_0^{2\pi} \sin kx\,dx &= 0, \\[0.4em] \int_0^{2\pi} \cos^2 kx\,dx &= \pi, \\[0.4em] \int_0^{2\pi} \sin^2 kx\,dx &= \pi. \end{aligned}

For ∫sin⁡mxcos⁡nx dx\int \sin^m x \cos^n x\,dx with whole numbers mm and nn,

CaseStepThen
nn oddkeep one cos⁡x\cos x, rewrite the rest with cos⁡2x=1−sin⁡2x\cos^2 x = 1 - \sin^2 xu=sin⁡xu = \sin x
mm oddkeep one sin⁡x\sin x, rewrite the rest with sin⁡2x=1−cos⁡2x\sin^2 x = 1 - \cos^2 xu=cos⁡xu = \cos x
mm and nn evenuse the half-angle identities, repeatedly if neededlower powers of cos⁡2x\cos 2x

For example ∫sin⁡3x dx=∫(1−cos⁡2x)sin⁡x dx=−cos⁡x+13cos⁡3x+C\int \sin^3 x\,dx = \int (1 - \cos^2 x)\sin x\,dx = -\cos x + \frac{1}{3}\cos^3 x + C. As xx runs from 00 to π\pi the doubled angle 2x2x runs through one full period of cosine, so

∫0πsin⁡2x dx=π2−12∫0πcos⁡2x dx=π2.\begin{aligned} &\int_0^\pi \sin^2 x\,dx \\[0.4em] &\quad = \frac{\pi}{2} - \frac{1}{2}\int_0^\pi \cos 2x\,dx \\[0.4em] &\quad = \frac{\pi}{2}. \end{aligned}

Lessons The half-angle identity

Trigonometric substitution

For a>0a > 0, a square root of a2−x2a^2 - x^2, a2+x2a^2 + x^2 or x2−a2x^2 - a^2 becomes a single trigonometric function after the substitution below. The angle range makes the substitution one-to-one and makes the trigonometric function in the last column nonnegative, so the square root needs no absolute value.

RootSubstitutionAngle rangeIdentityRoot becomes
a2−x2\sqrt{a^2 - x^2}x=asin⁡θx = a\sin\theta−π2≤θ≤π2-\frac{\pi}{2} \le \theta \le \frac{\pi}{2}1−sin⁡2θ=cos⁡2θ1 - \sin^2\theta = \cos^2\thetaacos⁡θa\cos\theta
a2+x2\sqrt{a^2 + x^2}x=atan⁡θx = a\tan\theta−π2<θ<π2-\frac{\pi}{2} < \theta < \frac{\pi}{2}1+tan⁡2θ=sec⁡2θ1 + \tan^2\theta = \sec^2\thetaasec⁡θa\sec\theta
x2−a2\sqrt{x^2 - a^2}x=asec⁡θx = a\sec\theta0≤θ<π20 \le \theta < \frac{\pi}{2} for x≥ax \ge a, π≤θ<3π2\pi \le \theta < \frac{3\pi}{2} for x≤−ax \le -asec⁡2θ−1=tan⁡2θ\sec^2\theta - 1 = \tan^2\thetaatan⁡θa\tan\theta

The differentials are dx=acos⁡θ dθdx = a\cos\theta\,d\theta, asec⁡2θ dθa\sec^2\theta\,d\theta and asec⁡θtan⁡θ dθa\sec\theta\tan\theta\,d\theta. Change the bounds with the angle, or return to xx with a right triangle, as in θ=arcsin⁡xa\theta = \arcsin\frac{x}{a} and cos⁡θ=a2−x2a\cos\theta = \frac{\sqrt{a^2 - x^2}}{a}. The same substitution gives the quarter circle ∫0aa2−x2 dx=πa24\int_0^a \sqrt{a^2 - x^2}\,dx = \frac{\pi a^2}{4}.

With x=2sin⁡θx = 2\sin\theta and 0≤θ≤π60 \le \theta \le \frac{\pi}{6}, which covers 0≤x≤10 \le x \le 1, the cosine is positive, so 4−x2=2cos⁡θ\sqrt{4 - x^2} = 2\cos\theta and dx=2cos⁡θ dθdx = 2\cos\theta\,d\theta. Then

∫01dx4−x2=∫0π/6dθ=π6,\begin{aligned} \int_0^1 \frac{dx}{\sqrt{4 - x^2}} &= \int_0^{\pi/6} d\theta \\[0.4em] &= \frac{\pi}{6}, \end{aligned}

which agrees with arcsin⁡12=π6\arcsin\frac{1}{2} = \frac{\pi}{6}.

Lessons Trigonometric substitution

Partial fractions

A rational function PQ\frac{P}{Q}, a ratio of polynomials, is proper when the degree of PP is less than the degree of QQ. If it is not proper, divide first to write it as a polynomial plus a proper fraction. Factor QQ over the real numbers into linear factors and irreducible quadratics, those with a negative discriminant. A proper fraction is then a sum of terms, with one group for each distinct factor of QQ, where mm is the power to which that factor appears.

Factor of QQTerms
x−rx - rAx−r\dfrac{A}{x - r}
(x−r)m(x - r)^mA1x−r+A2(x−r)2+⋯+Am(x−r)m\dfrac{A_1}{x - r} + \dfrac{A_2}{(x - r)^2} + \cdots + \dfrac{A_m}{(x - r)^m}
x2+px+qx^2 + px + q with p2<4qp^2 < 4qBx+Dx2+px+q\dfrac{Bx + D}{x^2 + px + q}
(x2+px+q)m(x^2 + px + q)^m with p2<4qp^2 < 4qB1x+D1x2+px+q+⋯+Bmx+Dm(x2+px+q)m\dfrac{B_1x + D_1}{x^2 + px + q} + \cdots + \dfrac{B_mx + D_m}{(x^2 + px + q)^m}

Find the constants by clearing denominators and comparing coefficients, or by putting the roots of QQ into the cleared identity. That identity is between polynomials, so it holds at every xx, even where the original fraction is undefined. The pieces integrate to logarithms, powers and, after completing the square, arctangents,

∫dxx−r=ln⁡∣x−r∣+C,∫dx(x−r)k=(x−r)1−k1−k+C,\begin{aligned} &\int \frac{dx}{x - r} = \ln\lvert x - r \rvert + C, \\[0.8em] &\int \frac{dx}{(x - r)^k} \\[0.4em] &\quad = \frac{(x - r)^{1-k}}{1 - k} + C, \end{aligned}

the second for k≥2k \ge 2.

FractionDecomposition
x+5(x−1)2(x+2)\dfrac{x + 5}{(x - 1)^2(x + 2)}2(x−1)2−13(x−1)+13(x+2)\dfrac{2}{(x - 1)^2} - \dfrac{1}{3(x - 1)} + \dfrac{1}{3(x + 2)}
2x2+x+1(x−1)(x2+1)\dfrac{2x^2 + x + 1}{(x - 1)(x^2 + 1)}2x−1+1x2+1\dfrac{2}{x - 1} + \dfrac{1}{x^2 + 1}

The poles of 1x2−1\frac{1}{x^2 - 1} are x=±1x = \pm 1, and neither lies in [2,3][2, 3], where both x−1x - 1 and x+1x + 1 are positive. From 1x2−1=12(x−1)−12(x+1)\frac{1}{x^2 - 1} = \frac{1}{2(x - 1)} - \frac{1}{2(x + 1)},

∫23dxx2−1=12[ln⁡x−1x+1]23=12ln⁡32≈0.2027.\begin{aligned} &\int_2^3 \frac{dx}{x^2 - 1} \\[0.4em] &\quad = \frac{1}{2}\left[\ln\frac{x - 1}{x + 1}\right]_2^3 \\[0.4em] &\quad = \frac{1}{2}\ln\frac{3}{2} \approx 0.2027. \end{aligned}

With a pole inside the interval the integral is improper. The integral of 1x2−1\frac{1}{x^2 - 1} over [0,2][0, 2] diverges, because near x=1x = 1 the integrand behaves like 12(x−1)\frac{1}{2(x - 1)}.

Lessons Partial fractions

Completing the square

For a quadratic x2+bx+cx^2 + bx + c,

x2+bx+c=(x+b2)2+c−b24.\begin{aligned} &x^2 + bx + c \\[0.4em] &\quad = \left(x + \frac{b}{2}\right)^2 + c - \frac{b^2}{4}. \end{aligned}

When the discriminant b2−4cb^2 - 4c is negative, the quadratic is irreducible over the real numbers and c−b24=k2c - \frac{b^2}{4} = k^2 with k>0k > 0. With h=−b2h = -\frac{b}{2}, the shift and scale substitution u=x−hku = \frac{x - h}{k} leads to ∫du1+u2\int \frac{du}{1 + u^2}, and

∫dx(x−h)2+k2=1karctan⁡x−hk+C,∫dxk2−(x−h)2=arcsin⁡x−hk+C,\begin{aligned} &\int \frac{dx}{(x - h)^2 + k^2} \\[0.4em] &\quad = \frac{1}{k}\arctan\frac{x - h}{k} + C, \\[0.8em] &\int \frac{dx}{\sqrt{k^2 - (x - h)^2}} \\[0.4em] &\quad = \arcsin\frac{x - h}{k} + C, \end{aligned}

the second for ∣x−h∣<k\lvert x - h \rvert < k. The inverse tangent takes values in (−π2,π2)\left(-\frac{\pi}{2}, \frac{\pi}{2}\right). A linear numerator is split into a multiple of the derivative 2x+b2x + b plus a constant, and ∫2x+bx2+bx+c dx=ln⁡(x2+bx+c)+C\int \frac{2x + b}{x^2 + bx + c}\,dx = \ln(x^2 + bx + c) + C when the quadratic is positive.

The quadratic x2+4x+8=(x+2)2+4x^2 + 4x + 8 = (x + 2)^2 + 4 has discriminant 16−32=−1616 - 32 = -16. With u=x+22u = \frac{x + 2}{2} and dx=2 dudx = 2\,du the bounds become 11 and 22, and

∫02dxx2+4x+8=12∫12du1+u2=12(arctan⁡2−π4)≈0.160875.\begin{aligned} &\int_0^2 \frac{dx}{x^2 + 4x + 8} \\[0.4em] &\quad = \frac{1}{2}\int_1^2 \frac{du}{1 + u^2} \\[0.4em] &\quad = \frac{1}{2}\left(\arctan 2 - \frac{\pi}{4}\right) \\[0.4em] &\quad \approx 0.160875. \end{aligned}

Lessons Completing the square

Choosing a method

Look at the integrand and try the first row that fits.

You seeTryKey step
a function of g(x)g(x) times the factor g′(x)g'(x)substitutionu=g(x)u = g(x), du=g′(x) dxdu = g'(x)\,dx
a product such as xexx e^x or xsin⁡xx\sin x, or ln⁡x\ln x aloneintegration by partsuu is the factor that gets simpler
even powers of sin⁡x\sin x and cos⁡x\cos xhalf-angle identitiessin⁡2x=1−cos⁡2x2\sin^2 x = \frac{1 - \cos 2x}{2}
an odd power of sin⁡x\sin x or cos⁡x\cos xkeep one factoru=cos⁡xu = \cos x or u=sin⁡xu = \sin x
a2−x2\sqrt{a^2 - x^2}, a2+x2\sqrt{a^2 + x^2} or x2−a2\sqrt{x^2 - a^2}trigonometric substitutionx=asin⁡θx = a\sin\theta, atan⁡θa\tan\theta or asec⁡θa\sec\theta
a ratio of polynomials with a factorable denominatorpartial fractionsdivide first if it is not proper
an irreducible quadratic in a denominator or under a rootcompleting the squareu=x−hku = \frac{x - h}{k}

Some integrands, such as e−x2e^{-x^2}, sin⁡xx\frac{\sin x}{x} and 1ln⁡x\frac{1}{\ln x}, have no elementary antiderivative. Their definite integrals come from symmetry, numerical rules, series or the methods of the last unit.

Lessons Substitution, Integration by parts, The half-angle identity, Trigonometric substitution, Partial fractions, Completing the square

Structure, limits and approximation

Some integrals are settled by structure rather than by an antiderivative. Symmetry cancels or doubles whole terms, an improper integral is a limit of ordinary integrals, Simpson’s rule estimates a value with a guaranteed error, and a power series can be integrated one term at a time.

Symmetry

A function is even when f(−x)=f(x)f(-x) = f(x) and odd when f(−x)=−f(x)f(-x) = -f(x), for every xx in a domain symmetric about 00. For ff integrable on a symmetric interval [−a,a][-a, a],

∫−aaf(x) dx=0,∫−aaf(x) dx=2∫0af(x) dx,\begin{gathered} \int_{-a}^{a} f(x)\,dx = 0, \\[0.6em] \int_{-a}^{a} f(x)\,dx = 2\int_0^a f(x)\,dx, \end{gathered}

the first when ff is odd and the second when ff is even. A product of two even or two odd functions is even, and a product of an even and an odd function is odd. Any ff on [−a,a][-a, a] splits into an even part f(x)+f(−x)2\frac{f(x) + f(-x)}{2} and an odd part f(x)−f(−x)2\frac{f(x) - f(-x)}{2}, so read the parity of each term before calculating. Since x3x^3 is odd and cos⁡x\cos x is even,

∫−11(x3+cos⁡x) dx=0+2∫01cos⁡x dx=2sin⁡1≈1.682942.\begin{aligned} &\int_{-1}^{1} \big(x^3 + \cos x\big)\,dx \\[0.4em] &\quad = 0 + 2\int_0^1 \cos x\,dx \\[0.4em] &\quad = 2\sin 1 \approx 1.682942. \end{aligned}

For an improper integral the rule for odd functions applies only after the integral is known to converge. The integral of xx over (−∞,∞)(-\infty, \infty) diverges, even though ∫−RRx dx=0\int_{-R}^{R} x\,dx = 0 for every RR, and the integral of 1x\frac{1}{x} over [−1,1][-1, 1] diverges too.

Lessons Symmetry

Improper integrals

An integral over an unbounded interval, or of a function that is unbounded near a point, is defined as a limit. For ff continuous on [a,∞)[a, \infty), and for ff continuous on (a,b](a, b] and unbounded near aa,

∫a∞f(x) dx=lim⁡R→∞∫aRf(x) dx,∫abf(x) dx=lim⁡t→a+∫tbf(x) dx,\begin{aligned} &\int_a^\infty f(x)\,dx \\[0.4em] &\quad = \lim_{R\to\infty} \int_a^R f(x)\,dx, \\[0.8em] &\int_a^b f(x)\,dx \\[0.4em] &\quad = \lim_{t\to a^+} \int_t^b f(x)\,dx, \end{aligned}

and similarly at the other end. The integral converges when the limit exists and is finite, and diverges otherwise. Over (−∞,∞)(-\infty, \infty), or with a problem at both ends or at an interior point, split the integral so that each piece has one problem, and it converges only when every piece converges.

The p-test. For real pp,

IntegralConverges exactly whenValue
∫1∞dxxp\displaystyle\int_1^\infty \frac{dx}{x^p}p>1p > 11p−1\dfrac{1}{p - 1}
∫01dxxp\displaystyle\int_0^1 \frac{dx}{x^p}p<1p < 111−p\dfrac{1}{1 - p}

Both follow from ∫1Rx−p dx=R1−p−11−p\int_1^R x^{-p}\,dx = \frac{R^{1-p} - 1}{1 - p} for p≠1p \ne 1. In the boundary case p=1p = 1 both diverge, since ∫1Rdxx=ln⁡R→∞\int_1^R \frac{dx}{x} = \ln R \to \infty and ∫t1dxx=−ln⁡t→∞\int_t^1 \frac{dx}{x} = -\ln t \to \infty, so ∫0∞x−p dx\int_0^\infty x^{-p}\,dx diverges for every pp. For p>1p > 1 the remaining tail is

∫R∞dxxp=R1−pp−1.\int_R^\infty \frac{dx}{x^p} = \frac{R^{1-p}}{p - 1}.

For example p=2p = 2 gives ∫1Rx−2 dx=1−1R→1\int_1^R x^{-2}\,dx = 1 - \frac{1}{R} \to 1, and p=1.1p = 1.1 gives 1010, approached slowly.

  • Comparison test. If ff and gg are continuous with 0≤g≤f0 \le g \le f on [a,∞)[a, \infty), convergence of ∫a∞f dx\int_a^\infty f\,dx implies convergence of ∫a∞g dx\int_a^\infty g\,dx, and divergence of the second implies divergence of the first. Since x2≥xx^2 \ge x for x≥1x \ge 1, the integral ∫1∞e−x2 dx\int_1^\infty e^{-x^2}\,dx converges and is at most ∫1∞e−x dx=1e\int_1^\infty e^{-x}\,dx = \frac{1}{e}.
  • Absolute convergence. For ff continuous on [a,∞)[a, \infty), if ∫a∞∣f∣ dx\int_a^\infty \lvert f \rvert\,dx converges, so does ∫a∞f dx\int_a^\infty f\,dx.
  • Singularities inside. The fundamental theorem needs ff continuous on [a,b][a, b]. Applied blindly to 1x2\frac{1}{x^2} on [−1,1][-1, 1] it gives [−1x]−11=−2\big[-\frac{1}{x}\big]_{-1}^{1} = -2, a negative value for a positive integrand, while in fact ∫01dxx2\int_0^1 \frac{dx}{x^2} diverges and so does the whole integral. Any integral of 1x\frac{1}{x} over an interval containing 00 diverges in the same way.

Lessons When an infinite region has finite area

Simpson’s rule

Composite Simpson’s rule. Take an even number nn of subintervals of [a,b][a, b] with width h=b−anh = \frac{b - a}{n}, points xj=a+jhx_j = a + jh and values fj=f(xj)f_j = f(x_j). Then

Sn=h3[f0+4f1+2f2+4f3+2f4+⋯+2fn−2+4fn−1+fn].\begin{aligned} S_n = \frac{h}{3}\big[&f_0 + 4f_1 + 2f_2 \\[0.4em] &+ 4f_3 + 2f_4 + \cdots \\[0.4em] &+ 2f_{n-2} + 4f_{n-1} \\[0.4em] &+ f_n\big]. \end{aligned}

On each pair of subintervals the rule integrates the quadratic interpolant through the three points, which gives h3(f0+4f1+f2)\frac{h}{3}(f_0 + 4f_1 + f_2), and it is exact for every polynomial of degree at most 33.

Simpson’s error bound. If f(4)f^{(4)} is continuous on [a,b][a, b] and ∣f(4)(x)∣≤M\lvert f^{(4)}(x) \rvert \le M there, then

∣∫abf(x) dx−Sn∣≤(b−a) h4M180=M(b−a)5180 n4.\begin{aligned} &\left\lvert \int_a^b f(x)\,dx - S_n \right\rvert \\[0.4em] &\quad \le \frac{(b - a)\,h^4 M}{180} \\[0.4em] &\quad = \frac{M(b - a)^5}{180\,n^4}. \end{aligned}

Any valid bound MM can be used, and a smaller one needs fewer subintervals. The function e−x2e^{-x^2} has no elementary antiderivative. On [0,1][0, 1] its fourth derivative is (16x4−48x2+12) e−x2(16x^4 - 48x^2 + 12)\,e^{-x^2}, and since e−x2≤1e^{-x^2} \le 1 the crude bound M=16+48+12=76M = 16 + 48 + 12 = 76 holds. Two subintervals give S2=16(1+4e−1/4+e−1)≈0.747180S_2 = \frac{1}{6}\big(1 + 4e^{-1/4} + e^{-1}\big) \approx 0.747180. With n=100n = 100 the error is at most 76180⋅1004≈4.22×10−9\frac{76}{180 \cdot 100^4} \approx 4.22 \times 10^{-9} and S100≈0.7468241329S_{100} \approx 0.7468241329, so

∫01e−x2 dx=0.746824\int_0^1 e^{-x^2}\,dx = 0.746824

to six decimal places. The largest value of ∣f(4)∣\lvert f^{(4)} \rvert on [0,1][0, 1] is in fact 1212, at x=0x = 0.

Lessons Approximate an integral with a controlled error

Power series and term-by-term integration

A power series ∑n=0∞cn(x−a)n\sum_{n=0}^{\infty} c_n (x - a)^n has a radius of convergence RR with 0≤R≤∞0 \le R \le \infty. It converges absolutely for ∣x−a∣<R\lvert x - a \rvert < R and diverges for ∣x−a∣>R\lvert x - a \rvert > R.

  • Term-by-term integration. If R>0R > 0, then for ∣x−a∣<R\lvert x - a \rvert < R the integral from aa to xx of the series is ∑n=0∞cnn+1(x−a)n+1\sum_{n=0}^{\infty} \frac{c_n}{n + 1}(x - a)^{n+1}, a power series with the same radius RR. The series converges uniformly on every closed interval inside (a−R,a+R)(a - R, a + R), which justifies exchanging the sum and the integral there.
  • Weierstrass M-test. If each fnf_n is continuous on [α,β][\alpha, \beta] with ∣fn(x)∣≤Mn\lvert f_n(x) \rvert \le M_n there and ∑Mn\sum M_n converges, then ∑fn\sum f_n converges uniformly on [α,β][\alpha, \beta] and its integral is the sum of the integrals of the fnf_n.
  • Geometric tail bound. If an>0a_n > 0 and an+1≤r ana_{n+1} \le r\,a_n for all n≥N+1n \ge N + 1 with r<1r < 1, then ∑n>Nan≤aN+11−r\sum_{n > N} a_n \le \frac{a_{N+1}}{1 - r}.
SeriesValid for
ex=∑n=0∞xnn!e^x = \displaystyle\sum_{n=0}^{\infty} \frac{x^n}{n!}all xx
11−x=∑n=0∞xn\dfrac{1}{1 - x} = \displaystyle\sum_{n=0}^{\infty} x^n−1<x<1-1 < x < 1
ln⁡(1+x)=∑n=1∞(−1)n−1xnn\ln(1 + x) = \displaystyle\sum_{n=1}^{\infty} \frac{(-1)^{n-1}x^n}{n}−1<x≤1-1 < x \le 1
arctan⁡x=∑n=0∞(−1)nx2n+12n+1\arctan x = \displaystyle\sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{2n + 1}−1≤x≤1-1 \le x \le 1

The last two come from integrating the geometric series for 11+t\frac{1}{1 + t} and 11+t2\frac{1}{1 + t^2} term by term, and they also hold at the endpoints listed.

The function ex−1x\frac{e^x - 1}{x} has a removable singularity at 00. With the value 11 there it equals ∑n=1∞xn−1n!\sum_{n=1}^{\infty} \frac{x^{n-1}}{n!} for every xx, each term is at most 1n!\frac{1}{n!} on [0,1][0, 1], and ∑1n!\sum \frac{1}{n!} converges, so

∫01ex−1x dx=∑n=1∞1n n!=1+14+118+⋯ .\begin{aligned} &\int_0^1 \frac{e^x - 1}{x}\,dx = \sum_{n=1}^{\infty} \frac{1}{n\,n!} \\[0.4em] &\quad = 1 + \frac{1}{4} + \frac{1}{18} + \cdots. \end{aligned}

The terms an=1n n!a_n = \frac{1}{n\,n!} have ratio an+1an=n(n+1)2<1N+2\frac{a_{n+1}}{a_n} = \frac{n}{(n + 1)^2} < \frac{1}{N + 2} for n≥N+1n \ge N + 1. The partial sum S10≈1.3179021490S_{10} \approx 1.3179021490 therefore has a tail below 2.5×10−92.5 \times 10^{-9}, and the integral is 1.3179021.317902 to six decimal places.

Lessons Integrate an infinite expansion

Geometric construction

Each formula here comes from one recipe. Cut the quantity into small pieces whose size is known, such as a strip, a washer, a shell or a short piece of curve, write one piece in terms of the variable, and add the pieces with an integral. Simpler shapes inside and around the figure give quick bounds that check the answer.

Area between curves

For ff and gg continuous on [a,b][a, b], the region between their graphs for a≤x≤ba \le x \le b has area

A=∫ab∣f(x)−g(x)∣ dx.A = \int_a^b \lvert f(x) - g(x) \rvert\,dx.

Where f≥gf \ge g this is top minus bottom, ∫ab(f(x)−g(x)) dx\int_a^b \big(f(x) - g(x)\big)\,dx. Where the curves cross, split at the crossings and take top minus bottom on each piece. The limits are often the intersection points, the solutions of f(x)=g(x)f(x) = g(x). With horizontal strips, for xL(y)≤xR(y)x_L(y) \le x_R(y) on [c,d][c, d], the area is ∫cd(xR(y)−xL(y)) dy\int_c^d \big(x_R(y) - x_L(y)\big)\,dy.

The line y=2xy = 2x and the parabola y=x2y = x^2 meet where x(x−2)=0x(x - 2) = 0, and 2x≥x22x \ge x^2 between, so

A=∫02(2x−x2) dx=4−83=43.\begin{aligned} A &= \int_0^2 \big(2x - x^2\big)\,dx \\[0.4em] &= 4 - \frac{8}{3} = \frac{4}{3}. \end{aligned}

Horizontal strips give the same area, ∫04(y−y2) dy=43\int_0^4 \big(\sqrt{y} - \frac{y}{2}\big)\,dy = \frac{4}{3}. On [0,π][0, \pi] the difference sin⁡x−cos⁡x\sin x - \cos x changes sign at π4\frac{\pi}{4}, so the area between the two graphs is (2−1)+(2+1)=22(\sqrt{2} - 1) + (\sqrt{2} + 1) = 2\sqrt{2}, while the signed integral ∫0π(sin⁡x−cos⁡x) dx\int_0^\pi (\sin x - \cos x)\,dx is only 22.

Lessons Area between curves

Volumes by slicing and washers

If the cross-section of a solid perpendicular to the xx-axis at xx has area A(x)A(x), continuous on [a,b][a, b], the solid has volume V=∫abA(x) dxV = \int_a^b A(x)\,dx. Rotating the region r(x)≤y≤R(x)r(x) \le y \le R(x), a≤x≤ba \le x \le b, with rr and RR continuous and 0≤r≤R0 \le r \le R, about the xx-axis gives a solid of revolution whose slices are washers with outer radius RR and inner radius rr. Its volume is

V=π∫ab(R2−r2) dx,V = \pi\int_a^b \big(R^2 - r^2\big)\,dx,

with R=R(x)R = R(x) and r=r(x)r = r(x).

The disk method is the case r=0r = 0. The radii are distances from the axis of rotation, so about a horizontal line y=ky = k that does not cross the region, measure both radii from that line, and about the yy-axis use horizontal slices and integrate in yy. The washer area is π(R2−r2)\pi(R^2 - r^2), not π(R−r)2\pi(R - r)^2.

Rotating x2≤y≤xx^2 \le y \le x, 0≤x≤10 \le x \le 1, about the xx-axis gives R(x)=xR(x) = x and r(x)=x2r(x) = x^2, so

V=π∫01(x2−x4) dx=π3−π5=2π15.\begin{aligned} V &= \pi\int_0^1 \big(x^2 - x^4\big)\,dx \\[0.4em] &= \frac{\pi}{3} - \frac{\pi}{5} = \frac{2\pi}{15}. \end{aligned}

About the line y=−1y = -1 the radii become x+1x + 1 and x2+1x^2 + 1, and the volume is 7π15\frac{7\pi}{15}.

Lessons Washers

Volumes by cylindrical shells

Rotate the region g(x)≤y≤f(x)g(x) \le y \le f(x), a≤x≤ba \le x \le b, with 0≤a0 \le a and ff and gg continuous, about the yy-axis. A strip parallel to the axis at xx sweeps a cylindrical shell with radius xx, height h(x)=f(x)−g(x)h(x) = f(x) - g(x) and thickness dxdx, so

V=2π∫abx h(x) dx.V = 2\pi\int_a^b x\,h(x)\,dx.

About a vertical line x=cx = c with c≤ac \le a or c≥bc \ge b, the shell radius is ∣x−c∣\lvert x - c \rvert. Shells use strips parallel to the axis and washers use strips perpendicular to it, and both give the same volume. Rotating 0≤y≤x20 \le y \le x^2, 0≤x≤10 \le x \le 1, about the yy-axis gives 2π∫01x3 dx=π22\pi\int_0^1 x^3\,dx = \frac{\pi}{2} by shells and π∫01(1−y) dy=π2\pi\int_0^1 (1 - y)\,dy = \frac{\pi}{2} by washers.

Rotating 0≤y≤e−x20 \le y \le e^{-x^2}, 0≤x≤10 \le x \le 1, about the yy-axis and substituting u=−x2u = -x^2 gives

V=2π∫01xe−x2 dx=π(1−e−1)≈1.985865.\begin{aligned} V &= 2\pi\int_0^1 x e^{-x^2}\,dx \\[0.4em] &= \pi\big(1 - e^{-1}\big) \\[0.4em] &\approx 1.985865. \end{aligned}

The answer lies between the volumes πe\frac{\pi}{e} and π\pi of the cylinders of radius 11 with heights e−1e^{-1} and 11.

Lessons Cylindrical shells

Arc length

For ff continuously differentiable on [a,b][a, b], the graph of ff has length

L=∫ab1+f′(x)2 dx.L = \int_a^b \sqrt{1 + f'(x)^2}\,dx.

The arc length element ds=dx2+dy2=1+f′(x)2 dxds = \sqrt{dx^2 + dy^2} = \sqrt{1 + f'(x)^2}\,dx is the Pythagorean length of a short segment of the curve. A parametric curve (x(t),y(t))\big(x(t), y(t)\big), α≤t≤β\alpha \le t \le \beta, with continuous derivatives and traced once, has length ∫αβx′(t)2+y′(t)2 dt\int_\alpha^\beta \sqrt{x'(t)^2 + y'(t)^2}\,dt. The catenary y=cosh⁡xy = \cosh x on [0,1][0, 1] has length sinh⁡1≈1.175201\sinh 1 \approx 1.175201, because 1+sinh⁡2x=cosh⁡2x1 + \sinh^2 x = \cosh^2 x.

For y=x3/2y = x^{3/2} on [0,1][0, 1] the derivative y′=32xy' = \frac{3}{2}\sqrt{x} is continuous on [0,1][0, 1], with right derivative 00 at x=0x = 0. Substituting u=1+94xu = 1 + \frac{9}{4}x,

L=∫011+94x dx=1313−827≈1.439710.\begin{aligned} L &= \int_0^1 \sqrt{1 + \tfrac{9}{4}x}\,dx \\[0.4em] &= \frac{13\sqrt{13} - 8}{27} \\[0.4em] &\approx 1.439710. \end{aligned}

The length lies between the chord 2\sqrt{2} and the sum 22 of the two legs, as it should for a curve that rises steadily.

Lessons The length of a curve

Surface area of revolution

Rotating the graph of y=f(x)y = f(x), with f≥0f \ge 0 continuously differentiable on [a,b][a, b], about the xx-axis gives a lateral surface of area

S=2π∫abf 1+(f′)2 dx,S = 2\pi\int_a^b f\,\sqrt{1 + (f')^2}\,dx,

with ff and f′f' evaluated at xx.

A short piece of the curve sweeps a band close to a conical frustum, whose lateral surface has area π(r1+r2) ℓ\pi(r_1 + r_2)\,\ell, the circumference at the average radius times the slant length ℓ\ell. That gives the surface area element dS=2πy dsdS = 2\pi y\,ds. About the yy-axis, with 0≤a0 \le a, the radius is xx and the integrand is 2πx1+f′(x)22\pi x\sqrt{1 + f'(x)^2}. The formula counts only the curved lateral surface, not end disks. For the line y=rhxy = \frac{r}{h}x on [0,h][0, h] it gives the lateral area πrr2+h2\pi r\sqrt{r^2 + h^2} of a cone.

For y=xy = \sqrt{x} on [1,4][1, 4] the integrand simplifies, since x1+14x=124x+1\sqrt{x}\sqrt{1 + \frac{1}{4x}} = \frac{1}{2}\sqrt{4x + 1}, and

S=π∫144x+1 dx=π6(1717−55)≈30.846490.\begin{aligned} S &= \pi\int_1^4 \sqrt{4x + 1}\,dx \\[0.4em] &= \frac{\pi}{6}\big(17\sqrt{17} - 5\sqrt{5}\big) \\[0.4em] &\approx 30.846490. \end{aligned}

Lessons Surface area

Multiple integration and coordinates

A double integral adds f ΔAf\,\Delta A over small pieces of a plane region, and a triple integral adds f ΔVf\,\Delta V over small pieces of a solid. Fubini’s theorem computes them as iterated single integrals, and new coordinates fit a round or slanted region at the cost of a scale factor, the absolute value of a Jacobian determinant.

Double and iterated integrals

The double integral ∬Df dA\iint_D f\,dA is the limit of sums of f ΔAf\,\Delta A over small rectangles covering DD as their size tends to 00, and for f≥0f \ge 0 it is the volume under the surface z=f(x,y)z = f(x, y) over DD. The area of DD is ∬D1 dA\iint_D 1\,dA, and when it is positive the average value of ff is 1area(D)∬Df dA\frac{1}{\text{area}(D)}\iint_D f\,dA. For ff continuous on a rectangular region [a,b]×[c,d][a, b] \times [c, d],

∬[a,b]×[c,d]f dA=∫ab∫cdf(x,y) dy dx=∫cd∫abf(x,y) dx dy.\begin{aligned} &\iint_{[a,b]\times[c,d]} f\,dA \\[0.4em] &\quad = \int_a^b \int_c^d f(x, y)\,dy\,dx \\[0.4em] &\quad = \int_c^d \int_a^b f(x, y)\,dx\,dy. \end{aligned}

The inner integral holds the outer variable fixed. For a region a≤x≤ba \le x \le b, g1(x)≤y≤g2(x)g_1(x) \le y \le g_2(x) with g1≤g2g_1 \le g_2 continuous, and ff continuous on it,

∬Df dA=∫ab∫g1(x)g2(x)f(x,y) dy dx,\begin{aligned} &\iint_D f\,dA \\[0.4em] &= \int_a^b \int_{g_1(x)}^{g_2(x)} f(x, y)\,dy\,dx, \end{aligned}

and a region c≤y≤dc \le y \le d, h1(y)≤x≤h2(y)h_1(y) \le x \le h_2(y) works the same way with the roles exchanged. Over a rectangle a product g(x) k(y)g(x)\,k(y) of continuous functions integrates to the product of the two single integrals.

The surface z=x+yz = x + y over the unit square has inner integral ∫01(x+y) dy=x+12\int_0^1 (x + y)\,dy = x + \frac{1}{2}, so

∫01∫01(x+y) dy dx=∫01(x+12)dx=1.\begin{aligned} &\int_0^1 \int_0^1 (x + y)\,dy\,dx \\[0.4em] &\quad = \int_0^1 \left(x + \frac{1}{2}\right)dx = 1. \end{aligned}

Lessons Accumulate over two coordinates

Fubini’s theorem and reversing the order

Fubini’s theorem. If ff is continuous on a closed bounded region DD that can be written both as a≤x≤ba \le x \le b, g1(x)≤y≤g2(x)g_1(x) \le y \le g_2(x) and as c≤y≤dc \le y \le d, h1(y)≤x≤h2(y)h_1(y) \le x \le h_2(y), with continuous bounds, then both iterated integrals equal ∬Df dA\iint_D f\,dA. For an unbounded region or integrand the same holds when ∬D∣f∣ dA\iint_D \lvert f \rvert\,dA is finite, and for f≥0f \ge 0 the two orders always agree, finite or infinite, which is Tonelli’s theorem. Without absolute integrability the orders can disagree. For x2−y2(x2+y2)2\frac{x^2 - y^2}{(x^2 + y^2)^2} on 0<x,y≤10 < x, y \le 1, integrating first in yy gives π4\frac{\pi}{4} and integrating first in xx gives −π4-\frac{\pi}{4}.

To reverse the order, sketch the region, read the bounds in the new direction from the sketch, and check that both descriptions give the same region. The integral of ey2e^{y^2} in yy has no elementary antiderivative, but the triangle 0≤x≤y≤10 \le x \le y \le 1 can also be read as 0≤y≤10 \le y \le 1, 0≤x≤y0 \le x \le y. The new inner integral supplies the missing factor yy for the substitution u=y2u = y^2,

∫01∫x1ey2 dy dx=∫01∫0yey2 dx dy=∫01y ey2 dy=e−12≈0.859141.\begin{aligned} &\int_0^1 \int_x^1 e^{y^2}\,dy\,dx \\[0.4em] &\quad = \int_0^1 \int_0^y e^{y^2}\,dx\,dy \\[0.4em] &\quad = \int_0^1 y\,e^{y^2}\,dy \\[0.4em] &\quad = \frac{e - 1}{2} \approx 0.859141. \end{aligned}

Lessons Change the order

Polar coordinates

Polar coordinates are x=rcos⁡θx = r\cos\theta and y=rsin⁡θy = r\sin\theta with r≥0r \ge 0, so x2+y2=r2x^2 + y^2 = r^2. The annular sector between the radii rr and r+Δrr + \Delta r and the angles θ\theta and θ+Δθ\theta + \Delta\theta has area exactly (r+Δr2)Δr Δθ\left(r + \frac{\Delta r}{2}\right)\Delta r\,\Delta\theta, which gives the polar area element dA=r dr dθdA = r\,dr\,d\theta.

Double integral in polar coordinates. Let DD be the region α≤θ≤β\alpha \le \theta \le \beta, h1(θ)≤r≤h2(θ)h_1(\theta) \le r \le h_2(\theta), with 0<β−α≤2π0 < \beta - \alpha \le 2\pi and h1h_1 and h2h_2 continuous with 0≤h1≤h20 \le h_1 \le h_2. For ff continuous on DD, and with g(r,θ)=f(rcos⁡θ,rsin⁡θ)g(r, \theta) = f(r\cos\theta, r\sin\theta),

∬Df dA=∫αβ∫h1(θ)h2(θ)g(r,θ) r dr dθ.\begin{aligned} &\iint_D f\,dA \\[0.4em] &= \int_\alpha^\beta \int_{h_1(\theta)}^{h_2(\theta)} g(r, \theta)\,r\,dr\,d\theta. \end{aligned}

Polar coordinates suit radial symmetry, an integrand that depends on x2+y2x^2 + y^2 over a disk, an annulus or a sector. Over the annulus 1≤x2+y2≤41 \le x^2 + y^2 \le 4, for example, ∬dAx2+y2=2πln⁡2\iint \frac{dA}{x^2 + y^2} = 2\pi\ln 2. Over the unit disk,

∬x2+y2≤1(x2+y2) dA=∫02π∫01r2 r dr dθ=π2,\begin{aligned} &\iint_{x^2 + y^2 \le 1} \big(x^2 + y^2\big)\,dA \\[0.4em] &\quad = \int_0^{2\pi} \int_0^1 r^2\,r\,dr\,d\theta \\[0.4em] &\quad = \frac{\pi}{2}, \end{aligned}

and dividing by the area π\pi gives the average value 12\frac{1}{2}.

Lessons Let the circle choose the coordinates

Change of variables and the Jacobian

For a map (u,v)↦(x(u,v),y(u,v))(u, v) \mapsto \big(x(u, v), y(u, v)\big) with continuous partial derivatives, the Jacobian determinant is

∂(x,y)∂(u,v)=det⁡(xuxvyuyv)=xuyv−xvyu.\begin{aligned} \frac{\partial(x, y)}{\partial(u, v)} &= \det\begin{pmatrix} x_u & x_v \\ y_u & y_v \end{pmatrix} \\[0.4em] &= x_u y_v - x_v y_u. \end{aligned}

Its absolute value is the local area scale factor, and its sign only records orientation.

Change of variables theorem. Let SS be a closed bounded region of the uvuv-plane with piecewise smooth boundary, let the map TT be continuously differentiable on an open set containing SS and one-to-one with nonzero Jacobian on the interior of SS, and let R=T(S)R = T(S). For ff continuous on RR, and with J=∂(x,y)∂(u,v)J = \frac{\partial(x, y)}{\partial(u, v)},

∬Rf(x,y) dx dy=∬Sf(T(u,v)) ∣J∣ du dv.\begin{aligned} &\iint_R f(x, y)\,dx\,dy \\[0.4em] &= \iint_S f\big(T(u, v)\big)\,\lvert J \rvert\,du\,dv. \end{aligned}

When the new coordinates are given as functions u(x,y)u(x, y) and v(x,y)v(x, y), the inverse map has the reciprocal Jacobian, ∂(x,y)∂(u,v)=1/∂(u,v)∂(x,y)\frac{\partial(x, y)}{\partial(u, v)} = 1 \Big/ \frac{\partial(u, v)}{\partial(x, y)}.

CoordinatesAbsolute value of the Jacobian
polar, (rcos⁡θ,rsin⁡θ)(r\cos\theta, r\sin\theta)rr
linear, (x,y)=A (u,v)(x, y) = A\,(u, v) for a constant invertible matrix AA∣det⁡A∣\lvert \det A \rvert
spherical, as in the next topicr2sin⁡ϕr^2\sin\phi

The parallelogram DD with vertices (0,0)(0, 0), (12,12)\left(\frac{1}{2}, \frac{1}{2}\right), (1,0)(1, 0) and (12,−12)\left(\frac{1}{2}, -\frac{1}{2}\right) is 0≤x+y≤10 \le x + y \le 1, 0≤x−y≤10 \le x - y \le 1. With u=x+yu = x + y and v=x−yv = x - y it becomes the unit square, and ∂(u,v)∂(x,y)=−2\frac{\partial(u, v)}{\partial(x, y)} = -2, so the area scale is ∣∂(x,y)∂(u,v)∣=12\left\lvert \frac{\partial(x, y)}{\partial(u, v)} \right\rvert = \frac{1}{2} and the area of DD is 12\frac{1}{2}. Then

∬Dex+y dx dy=12∫01∫01eu dv du=e−12.\begin{aligned} &\iint_D e^{x+y}\,dx\,dy \\[0.4em] &\quad = \frac{1}{2}\int_0^1 \int_0^1 e^u\,dv\,du \\[0.4em] &\quad = \frac{e - 1}{2}. \end{aligned}

Lessons Build coordinates from the boundaries

Triple integrals and spherical coordinates

A triple integral ∭Ef dV\iiint_E f\,dV adds f ΔVf\,\Delta V over small boxes, and for ff continuous on a box it can be computed as an iterated integral in any of the six orders. A solid with continuous density δ≥0\delta \ge 0 has mass M=∭Eδ dVM = \iiint_E \delta\,dV. Spherical coordinates are

x=rsin⁡ϕcos⁡θ,y=rsin⁡ϕsin⁡θ,z=rcos⁡ϕ,\begin{aligned} x &= r\sin\phi\cos\theta, \\[0.4em] y &= r\sin\phi\sin\theta, \\[0.4em] z &= r\cos\phi, \end{aligned}

where r≥0r \ge 0 is the distance from the origin, the polar angle 0≤ϕ≤π0 \le \phi \le \pi is measured down from the positive zz-axis and the azimuth θ\theta runs over an interval of length at most 2π2\pi. Then x2+y2+z2=r2x^2 + y^2 + z^2 = r^2. A small cell has edges drdr, r dϕr\,d\phi and rsin⁡ϕ dθr\sin\phi\,d\theta, so the spherical volume element is

dV=r2sin⁡ϕ dr dϕ dθ.dV = r^2\sin\phi\,dr\,d\phi\,d\theta.

Many books write ρ\rho for this distance and keep rr for the distance from the zz-axis. The unit ball is 0≤r≤10 \le r \le 1, 0≤ϕ≤π0 \le \phi \le \pi, 0≤θ≤2π0 \le \theta \le 2\pi, with volume 4π3\frac{4\pi}{3}. With the radial density x2+y2+z2=r2x^2 + y^2 + z^2 = r^2, the integrand r2⋅r2sin⁡ϕr^2 \cdot r^2\sin\phi separates into one factor for each coordinate, and the mass is

M=∫01r4 dr⋅∫0πsin⁡ϕ dϕ⋅∫02πdθ=15⋅2⋅2π=4π5,\begin{aligned} M &= \int_0^1 r^4\,dr \cdot \int_0^\pi \sin\phi\,d\phi \\[0.4em] &\qquad \cdot \int_0^{2\pi} d\theta \\[0.4em] &= \frac{1}{5} \cdot 2 \cdot 2\pi = \frac{4\pi}{5}, \end{aligned}

and the average density is 35\frac{3}{5}.

Lessons A ball, a radial density and a volume cell

Fields, circulation and flux

A vector field is integrated along curves for work and circulation and across surfaces for flux. The fundamental theorem for line integrals, Green’s theorem and the divergence theorem each replace one of these integrals by a simpler one, and each needs its orientation and its conditions on smoothness and on the domain.

Line integrals

Let CC be an oriented piecewise smooth curve with a parametrization r(t)=(x(t),y(t))\mathbf{r}(t) = \big(x(t), y(t)\big), a≤t≤ba \le t \le b, that traces it once in the chosen direction, and let F=(P,Q)\mathbf{F} = (P, Q) be continuous on CC. The line integral is

∫CF⋅dr=∫abF(r(t))⋅r′(t) dt=∫CP dx+Q dy.\begin{aligned} &\int_C \mathbf{F} \cdot d\mathbf{r} \\[0.4em] &\quad = \int_a^b \mathbf{F}\big(\mathbf{r}(t)\big) \cdot \mathbf{r}'(t)\,dt \\[0.4em] &\quad = \int_C P\,dx + Q\,dy. \end{aligned}

For a force field it is the work along CC, and around a closed curve, written ∮C\oint_C, it is the circulation. The value does not depend on the parametrization as long as the direction is kept, reversing the direction changes its sign, and a curve made of pieces gets the sum of the pieces. The field (−y,x)(-y, x) gives 13\frac{1}{3} along y=x2y = x^2 from (0,0)(0, 0) to (1,1)(1, 1) and 00 along the straight segment, so its line integrals depend on the path.

Lessons When a path integral only needs its endpoints

Conservative fields and path independence

Fundamental theorem for line integrals. If ψ\psi is continuously differentiable on an open set DD and CC is a piecewise smooth curve in DD from AA to BB, then

∫C∇ψ⋅dr=ψ(B)−ψ(A).\int_C \nabla\psi \cdot d\mathbf{r} = \psi(B) - \psi(A).

A field is conservative on DD when F=∇ψ\mathbf{F} = \nabla\psi there for a potential function ψ\psi, and then P dx+Q dy=dψP\,dx + Q\,dy = d\psi is an exact differential. For a continuous field on a connected open set DD three statements are equivalent. The field is conservative, its line integrals are path independent, and its integral around every closed piecewise smooth curve in DD is 00.

If F=(P,Q)\mathbf{F} = (P, Q) is continuously differentiable and conservative on an open set, then Py=QxP_y = Q_x there. The converse holds on a simply connected open set, which in the plane means a connected open set with no holes, but not in general. The field (−yx2+y2,xx2+y2)\left(\frac{-y}{x^2 + y^2}, \frac{x}{x^2 + y^2}\right) has Py=QxP_y = Q_x on the plane without the origin, yet its circulation counterclockwise around the unit circle is 2π2\pi. A potential is found by integrating PP in xx and fixing the added function of yy with QQ. For example (2xy+1,x2)(2xy + 1, x^2) has the potential x2y+xx^2y + x, so every path from (0,0)(0, 0) to (1,2)(1, 2) gives 33.

The field (y,x)(y, x) is the gradient of ψ=xy\psi = xy on the whole plane, since y dx+x dy=d(xy)y\,dx + x\,dy = d(xy). Every piecewise smooth path from (0,0)(0, 0) to (1,1)(1, 1) therefore gives

∫C(y dx+x dy)=ψ(1,1)−ψ(0,0)=1,\begin{aligned} &\int_C (y\,dx + x\,dy) \\[0.4em] &\quad = \psi(1, 1) - \psi(0, 0) = 1, \end{aligned}

a closed path gives 00, and the reversed path gives −1-1.

Lessons When a path integral only needs its endpoints

Green’s theorem

Green’s theorem. Let CC be a positively oriented (counterclockwise, with the region on the left), piecewise smooth, simple closed curve, let DD be CC together with the region inside it, and let PP and QQ be continuously differentiable on an open set containing DD. Then

∮CP dx+Q dy=∬D(∂Q∂x−∂P∂y)dA.\begin{aligned} &\oint_C P\,dx + Q\,dy \\[0.4em] &\quad = \iint_D \left(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right) dA. \end{aligned}

The scalar curl Qx−PyQ_x - P_y measures circulation per unit area. Adding the circulations around small cells gives the theorem, because each shared interior edge is traversed once in each direction and cancels, leaving only the outer boundary. A clockwise boundary changes the sign. A single point of DD where PP or QQ is undefined breaks the hypothesis, as for the field with circulation 2π2\pi in the previous topic. Taking Qx−Py=1Q_x - P_y = 1 gives the area of DD as 12∮Cx dy−y dx\frac{1}{2}\oint_C x\,dy - y\,dx, which is πab\pi ab for the ellipse x=acos⁡tx = a\cos t, y=bsin⁡ty = b\sin t with a,b>0a, b > 0.

The field (−y,x)(-y, x) has scalar curl 1−(−1)=21 - (-1) = 2, so around the unit circle counterclockwise

∮C(−y dx+x dy)=∬D2 dA=2π.\begin{aligned} &\oint_C (-y\,dx + x\,dy) \\[0.4em] &\quad = \iint_D 2\,dA = 2\pi. \end{aligned}

Directly, the field has length 11 and points along the circle, so the work is the length 2π2\pi of the circle. Clockwise, the circulation is −2π-2\pi.

Lessons Add the turning inside

Flux and the divergence theorem

The flux of a field F\mathbf{F} through an oriented surface SS with unit normal n\mathbf{n} is the surface integral ∬SF⋅n dS\iint_S \mathbf{F} \cdot \mathbf{n}\,dS, where dSdS is the surface area element. For F=(F1,F2,F3)\mathbf{F} = (F_1, F_2, F_3) the divergence is

∇⋅F=∂F1∂x+∂F2∂y+∂F3∂z.\nabla \cdot \mathbf{F} = \frac{\partial F_1}{\partial x} + \frac{\partial F_2}{\partial y} + \frac{\partial F_3}{\partial z}.

Divergence theorem. Let EE be a bounded solid whose boundary SS is a closed piecewise smooth surface oriented by the outward unit normal n\mathbf{n}, and let F\mathbf{F} be continuously differentiable on an open set containing EE and SS. Then

∬SF⋅n dS=∭E∇⋅F dV.\begin{aligned} &\iint_S \mathbf{F} \cdot \mathbf{n}\,dS \\[0.4em] &\quad = \iiint_E \nabla \cdot \mathbf{F}\,dV. \end{aligned}

Adding the outward fluxes of small boxes gives the theorem, because the flux through each shared interior face is counted once in each direction and cancels. The inward normal changes the sign. The field must be continuously differentiable at every point of EE. The field r∣r∣3\frac{\mathbf{r}}{\lvert \mathbf{r} \rvert^3} has divergence 00 away from the origin, yet its outward flux through every sphere about the origin is 4π4\pi.

For F=(x,y,z)\mathbf{F} = (x, y, z) the divergence is 33, and through the unit sphere, the boundary of the unit ball BB, with the outward normal

∬SF⋅n dS=∭B3 dV=3⋅4π3=4π.\begin{aligned} \iint_S \mathbf{F} \cdot \mathbf{n}\,dS &= \iiint_B 3\,dV \\[0.4em] &= 3 \cdot \frac{4\pi}{3} = 4\pi. \end{aligned}

Directly, n=(x,y,z)\mathbf{n} = (x, y, z) on the unit sphere, so F⋅n=1\mathbf{F} \cdot \mathbf{n} = 1 and the flux is the area 4π4\pi. With the inward normal it is −4π-4\pi.

The theorems of this unit and the fundamental theorem of calculus all integrate a derivative over a region and get boundary values.

TheoremRegion integralBoundary term
Fundamental theorem of calculus∫abF′(x) dx\displaystyle\int_a^b F'(x)\,dxF(b)−F(a)F(b) - F(a)
Line integrals∫C∇ψ⋅dr\displaystyle\int_C \nabla\psi \cdot d\mathbf{r}ψ(B)−ψ(A)\psi(B) - \psi(A)
Green∬D(Qx−Py) dA\displaystyle\iint_D (Q_x - P_y)\,dA∮CP dx+Q dy\displaystyle\oint_C P\,dx + Q\,dy
Divergence∭E∇⋅F dV\displaystyle\iiint_E \nabla \cdot \mathbf{F}\,dV∬SF⋅n dS\displaystyle\iint_S \mathbf{F} \cdot \mathbf{n}\,dS

Lessons Count outward flow from the sources inside

Advanced transformations

Some single integrals are easiest in a larger setting. Squaring the Gaussian integral moves it into the plane, a parameter turns one integral into a family that can be differentiated, and the residue theorem evaluates real integrals with contours in the complex plane.

The Gaussian integral

∫−∞∞e−x2 dx=π.\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}.

Call the integral II. It is positive and finite, because e−x2≤e−∣x∣e^{-x^2} \le e^{-\lvert x \rvert} for ∣x∣≥1\lvert x \rvert \ge 1. Squaring it as a product of independent integrals in xx and yy gives an improper double integral over the whole plane, and polar coordinates evaluate it,

I2=∬R2e−(x2+y2) dA=∫02π∫0∞e−r2 r dr dθ=2π⋅12=π.\begin{aligned} I^2 &= \iint_{\mathbb{R}^2} e^{-(x^2 + y^2)}\,dA \\[0.4em] &= \int_0^{2\pi} \int_0^\infty e^{-r^2}\,r\,dr\,d\theta \\[0.4em] &= 2\pi \cdot \frac{1}{2} = \pi. \end{aligned}

The integrand is positive and continuous, so Tonelli’s theorem allows the product to be written as a double integral and computed in Cartesian or in polar coordinates. The positive root gives I=πI = \sqrt{\pi}. The square I2I^2 is not ∫−∞∞e−2x2 dx\int_{-\infty}^{\infty} e^{-2x^2}\,dx, which equals π/2\sqrt{\pi/2}.

IntegralValueConditions
∫0∞e−x2 dx\displaystyle\int_0^\infty e^{-x^2}\,dxπ2\dfrac{\sqrt{\pi}}{2}
∫−∞∞e−ax2 dx\displaystyle\int_{-\infty}^{\infty} e^{-ax^2}\,dxπa\sqrt{\dfrac{\pi}{a}}a>0a > 0
∫−∞∞x2e−x2 dx\displaystyle\int_{-\infty}^{\infty} x^2 e^{-x^2}\,dxπ2\dfrac{\sqrt{\pi}}{2}
∫−∞∞1σ2π e−(x−μ)2/(2σ2) dx\displaystyle\int_{-\infty}^{\infty} \frac{1}{\sigma\sqrt{2\pi}}\,e^{-(x - \mu)^2/(2\sigma^2)}\,dx11σ>0\sigma > 0
Γ(12)=∫0∞t−1/2e−t dt\displaystyle\Gamma\left(\tfrac{1}{2}\right) = \int_0^\infty t^{-1/2}e^{-t}\,dtπ\sqrt{\pi}

Lessons Square the question

Differentiation under the integral sign

Leibniz integral rule. If f(x,t)f(x, t) and ∂f∂t(x,t)\frac{\partial f}{\partial t}(x, t) are continuous on [a,b]×[c,d][a, b] \times [c, d], then F(t)=∫abf(x,t) dxF(t) = \int_a^b f(x, t)\,dx is differentiable on [c,d][c, d], with one-sided derivatives at cc and dd, and

F′(t)=∫ab∂f∂t(x,t) dx.F'(t) = \int_a^b \frac{\partial f}{\partial t}(x, t)\,dx.

Dominated version. For an improper integral, or an integrand that is not continuous at the ends, suppose that f(⋅,t)f(\cdot, t) is integrable over the interval of integration for each tt in an interval JJ, that ∂f∂t\frac{\partial f}{\partial t} exists and is continuous for xx inside the interval and tt in JJ, and that ∣∂f∂t(x,t)∣≤g(x)\left\lvert \frac{\partial f}{\partial t}(x, t) \right\rvert \le g(x) for all such xx and tt, with ∫g dx\int g\,dx finite. Then FF is differentiable on JJ with the same formula, one-sided at an endpoint of JJ that belongs to JJ.

With variable limits a(t)a(t) and b(t)b(t), differentiable with values in [x0,x1][x_0, x_1], and ff and ∂f∂t\frac{\partial f}{\partial t} continuous on [x0,x1]×[c,d][x_0, x_1] \times [c, d],

ddt∫a(t)b(t)f(x,t) dx=f(b(t),t) b′(t)−f(a(t),t) a′(t)+∫a(t)b(t)∂f∂t(x,t) dx.\begin{aligned} &\frac{d}{dt}\int_{a(t)}^{b(t)} f(x, t)\,dx \\[0.4em] &\quad = f\big(b(t), t\big)\,b'(t) \\[0.4em] &\qquad - f\big(a(t), t\big)\,a'(t) \\[0.4em] &\qquad + \int_{a(t)}^{b(t)} \frac{\partial f}{\partial t}(x, t)\,dx. \end{aligned}

To evaluate ∫01x−1ln⁡x dx\int_0^1 \frac{x - 1}{\ln x}\,dx, embed it in the parameter family

F(a)=∫01xa−1ln⁡x dxF(a) = \int_0^1 \frac{x^a - 1}{\ln x}\,dx

for a≥0a \ge 0. The integrand has removable endpoint limits, 00 as x→0+x \to 0^+ and aa as x→1−x \to 1^-. Its derivative in aa is xax^a, and 0≤xa≤10 \le x^a \le 1 on (0,1)(0, 1) is an integrable bound. So F′(a)=∫01xa dx=1a+1F'(a) = \int_0^1 x^a\,dx = \frac{1}{a + 1}, with a right derivative at 00. The known value F(0)=0F(0) = 0 fixes the constant in F(a)=ln⁡(a+1)F(a) = \ln(a + 1), and

∫01x−1ln⁡x dx=F(1)=ln⁡2≈0.693147.\begin{aligned} &\int_0^1 \frac{x - 1}{\ln x}\,dx = F(1) \\[0.4em] &\quad = \ln 2 \approx 0.693147. \end{aligned}

In the same way, differentiating ∫0∞e−tx dx=1t\int_0^\infty e^{-tx}\,dx = \frac{1}{t} for t>0t > 0 repeatedly, with the bound xne−t0xx^n e^{-t_0x} for t≥t0>0t \ge t_0 > 0, gives ∫0∞xne−tx dx=n!tn+1\int_0^\infty x^n e^{-tx}\,dx = \frac{n!}{t^{n+1}} for every whole number nn.

Lessons Make a family, then differentiate the family

Contour integration and residues

A complex function ff is analytic on an open set when it has a complex derivative at every point of the set. For a piecewise smooth curve z(t)z(t), a≤t≤ba \le t \le b, the contour integral is ∫Cf(z) dz=∫abf(z(t)) z′(t) dt\int_C f(z)\,dz = \int_a^b f\big(z(t)\big)\,z'(t)\,dt. Counterclockwise around a circle of radius ρ>0\rho > 0 about pp, the basic integral is ∮dzz−p=2πi\oint \frac{dz}{z - p} = 2\pi i.

  • Cauchy’s theorem. If ff is analytic on an open set containing a simple closed piecewise smooth curve CC and its inside, then ∮Cf(z) dz=0\oint_C f(z)\,dz = 0.

  • Residue theorem. Let CC be a positively oriented simple closed piecewise smooth curve, and let ff be analytic on an open set containing CC and its inside, except at finitely many isolated singularities z1,…,znz_1, \ldots, z_n inside CC and none on CC. Then

    ∮Cf(z) dz=2πi∑k=1nRes⁡(f,zk).\begin{aligned} &\oint_C f(z)\,dz \\[0.4em] &\quad = 2\pi i\sum_{k=1}^{n} \operatorname{Res}(f, z_k). \end{aligned}

The residue Res⁡(f,p)\operatorname{Res}(f, p) is the coefficient of 1z−p\frac{1}{z - p} in the Laurent expansion of ff about pp.

Singularity at ppResidueConditions
simple polelim⁡z→p (z−p) f(z)\displaystyle\lim_{z\to p}\,(z - p)\,f(z)
simple zero of a denominatorg(p)h′(p)\dfrac{g(p)}{h'(p)} for f=ghf = \dfrac{g}{h}gg and hh analytic at pp, h(p)=0h(p) = 0, h′(p)≠0h'(p) \ne 0
pole of order mm1(m−1)!lim⁡z→pdm−1dzm−1[(z−p)mf(z)]\displaystyle\frac{1}{(m - 1)!}\lim_{z\to p}\frac{d^{m-1}}{dz^{m-1}}\big[(z - p)^m f(z)\big]

To evaluate a real integral over (−∞,∞)(-\infty, \infty), close the segment [−R,R][-R, R] with the semicircle ΓR\Gamma_R of radius RR in the upper half-plane, and show that the integral over the arc tends to 00 as R→∞R \to \infty.

  • ML estimate. If ∣f∣≤M\lvert f \rvert \le M on a curve Γ\Gamma of length LL, then ∣∫Γf(z) dz∣≤ML\left\lvert \int_\Gamma f(z)\,dz \right\rvert \le ML.
  • Jordan’s lemma. If a>0a > 0 and MRM_R is the largest value of ∣g∣\lvert g \rvert on ΓR\Gamma_R, then ∣∫ΓRg(z) eiaz dz∣≤πMRa\left\lvert \int_{\Gamma_R} g(z)\,e^{iaz}\,dz \right\rvert \le \frac{\pi M_R}{a}, which tends to 00 when MR→0M_R \to 0.
  • Rational integrands. If PP and QQ are polynomials, the degree of QQ is at least the degree of PP plus 22 and QQ has no real zeros, the arc integral tends to 00 by the ML estimate and ∫−∞∞P(x)Q(x) dx\int_{-\infty}^{\infty} \frac{P(x)}{Q(x)}\,dx is 2πi2\pi i times the sum of the residues of PQ\frac{P}{Q} in the upper half-plane.

For example eiz1+z2\frac{e^{iz}}{1 + z^2} has residue e−12i\frac{e^{-1}}{2i} at ii. Since ∣eiz∣≤1\lvert e^{iz} \rvert \le 1 in the upper half-plane, the ML estimate makes the arc vanish, and taking real parts gives ∫−∞∞cos⁡x1+x2 dx=πe\int_{-\infty}^{\infty} \frac{\cos x}{1 + x^2}\,dx = \frac{\pi}{e}.

The function f(z)=11+z4f(z) = \frac{1}{1 + z^4} has simple poles at the four fourth roots of −1-1. For R>1R > 1 the upper semicircle encloses p1=eiπ/4p_1 = e^{i\pi/4} and p2=e3iπ/4p_2 = e^{3i\pi/4}, where Res⁡(f,p)=14p3=−p4\operatorname{Res}(f, p) = \frac{1}{4p^3} = -\frac{p}{4} because p4=−1p^4 = -1. The two residues add up to −i22-\frac{i}{2\sqrt{2}}, and on the arc ∣1+z4∣≥R4−1\lvert 1 + z^4 \rvert \ge R^4 - 1, so the arc integral is at most πRR4−1\frac{\pi R}{R^4 - 1} in absolute value, which tends to 00. The integral over the real line is 2πi⋅(−i22)=π22\pi i \cdot \left(-\frac{i}{2\sqrt{2}}\right) = \frac{\pi}{\sqrt{2}}, and since the integrand is even,

∫0∞dx1+x4=π22≈1.110721.\begin{aligned} \int_0^\infty \frac{dx}{1 + x^4} &= \frac{\pi}{2\sqrt{2}} \\[0.4em] &\approx 1.110721. \end{aligned}

Lessons Leave the real line, collect the poles, return