Reading an Integral One Variable at a Time
The product measure from the previous tutorial lets us integrate functions on \(X\times Y\). Iterated integration is a way to compute such an integral by first fixing one coordinate, integrating over the other, and then integrating the result over the fixed coordinate. The order matters for how the calculation is set up, even when the theorems of this course guarantee that either order gives the same answer.
Throughout, take \((X,\mathcal{F},\mu)\) and \((Y,\mathcal{G},\nu)\) to be sigma-finite measure spaces, and take functions on \(X\times Y\) to be measurable with respect to \(\mathcal{F}\otimes\mathcal{G}\). The sigma-finiteness hypothesis is essential here: it is among the hypotheses of the Tonelli and Fubini theorems established earlier in this course. We will use those theorems by name rather than re-prove them.
The sections here are sections of a function: for fixed \(x\), the function of \(y\) is \(y\mapsto f(x,y)\), and for fixed \(y\), the function of \(x\) is \(x\mapsto f(x,y)\). The earlier section results for product-measurable sets, together with Tonelli’s Theorem, ensure that the inner integral of a nonnegative measurable function is measurable as a function of the remaining coordinate. In particular, the iterated integrals are well-defined as extended nonnegative integrals.
For a signed function, care is needed. Fubini’s Theorem applies when \(f\) is integrable on the product space, meaning \(\int_{X\times Y}|f|\,d(\mu\otimes\nu)<\infty\). It then gives section integrability almost everywhere and equality of both iterated integrals with the product-space integral. If the exceptional set of outer coordinates where a section is not integrable is nonempty, the inner integral can be assigned any value there without changing the outer integral. For a nonnegative function, Tonelli’s Theorem instead permits infinite values and does not require the product-space integral to be finite.
A Product of Separate Functions
A useful first case is a function that separates into one factor depending only on \(x\) and another depending only on \(y\). The result below turns a two-variable integral into a product of one-variable integrals. Its hypotheses ensure that the factors and their integrals are finite; this avoids ambiguity about multiplying zero and infinity.
Proof. The function \((x,y)\mapsto u(x)v(y)\) is product-measurable. For each fixed \(x\), nonnegative homogeneity of the integral gives $$ \int_Yu(x)v(y)\,d\nu(y)=u(x)\int_Yv(y)\,d\nu(y). $$ Integrating this equality over \(X\), and using nonnegative homogeneity once more, yields $$ \int_X\left(\int_Yu(x)v(y)\,d\nu(y)\right)d\mu(x) =\left(\int_Yv\,d\nu\right)\left(\int_Xu\,d\mu\right). $$ Tonelli’s Theorem applies because the spaces are sigma-finite and the product function is nonnegative and measurable. It identifies this iterated integral with the product-space integral and also identifies the iterated integral in the other order with that same value. This proves the result.
Worked Example: Integrating a Sum on a Rectangle
Let \(R=[0,2]\times[1,3]\), with two-dimensional Lebesgue measure, and set \(f(x,y)=x+2y\). For each fixed \(x\), integrate first in \(y\): $$ \int_1^3(x+2y)\,dy =\left[xy+y^2\right]_1^3 =(3x+9)-(x+1) =2x+8. $$ Therefore, $$ \int_0^2(2x+8)\,dx =\left[x^2+8x\right]_0^2 =4+16 =20. $$ In the reverse order, for fixed \(y\), $$ \int_0^2(x+2y)\,dx =\left[\frac{x^2}{2}+2yx\right]_0^2 =2+4y. $$ Thus $$ \int_1^3(2+4y)\,dy =\left[2y+2y^2\right]_1^3 =(6+18)-(2+2) =20. $$ The rectangle has constant coordinate bounds, so reversing the order does not change the region being integrated over.
Changing the Order by Slicing a Region
For a nonrectangular region, the bounds in an iterated integral describe its sections. Fixing \(x\) gives a vertical slice; fixing \(y\) gives a horizontal slice. The following result formalizes this procedure by extending a function on a region to be zero outside it. It is a direct application of Fubini’s Theorem, not a new version of that theorem.
Proof. Extend \(f\) to \(F:X\times Y\to\mathbb{R}\) by setting \(F(x,y)=f(x,y)\) on \(E\) and \(F(x,y)=0\) outside \(E\). This extension is measurable: for any Borel set \(C\subseteq\mathbb{R}\), its preimage under \(F\) is assembled from the measurable part of \(E\) on which \(f\) lies in \(C\), and, if \(0\in C\), the measurable set \(E^c\). Also, $$ \int_{X\times Y}|F|\,d(\mu\otimes\nu)=\int_E|f|\,d(\mu\otimes\nu)<\infty. $$ Fubini’s Theorem therefore applies to \(F\). For each \(x\), \(F(x,y)\) equals \(f(x,y)\) on \(E_x\) and zero outside \(E_x\), so its section integral is the integral of \(f(x,\cdot)\) over \(E_x\), wherever that section is integrable. The corresponding statement holds for each horizontal section \(E^y\). Finally, the product-space integral of \(F\) is the integral of \(f\) over \(E\). Substituting these three identifications into Fubini’s equalities proves the result.
When the integrand is nonnegative and its integral may be infinite, use Tonelli’s Theorem in the same extension-by-zero method. The sigma-finiteness hypotheses on both factors still apply. For geometric calculations, it is especially important to describe the region itself before writing bounds: the bounds must include exactly the points of the region, up to boundary sets that do not affect the integral in the examples below.
Worked Example: Integrating Over a Triangle in Both Orders
Let \(E=\{(x,y):0\leq y\leq x\leq1\}\), and integrate \(f(x,y)=x+y\). With \(x\) fixed, \(y\) ranges from \(0\) to \(x\), while \(x\) ranges from \(0\) to \(1\). Hence $$ \int_0^1\int_0^x(x+y)\,dy\,dx =\int_0^1\left[xy+\frac{y^2}{2}\right]_{y=0}^{y=x}dx =\int_0^1\frac{3x^2}{2}\,dx =\frac12. $$ To reverse the order, fix \(y\). The inequalities \(y\leq x\leq1\) show that \(x\) ranges from \(y\) to \(1\), and \(y\) ranges from \(0\) to \(1\). Therefore $$ \int_0^1\int_y^1(x+y)\,dx\,dy =\int_0^1\left[\frac{x^2}{2}+yx\right]_{x=y}^{x=1}dy =\int_0^1\left(\frac12+y-\frac{3y^2}{2}\right)dy =\frac12+\frac12-\frac12 =\frac12. $$ The two descriptions of the triangle are \(0\leq y\leq x\leq1\) and \(0\leq y\leq1,\ y\leq x\leq1\). Deriving the second from the defining inequalities is the key step in reversing the order.
Worked Example: An Integral over the Positive Quadrant
Consider \(f(x,y)=e^{-x}e^{-y}\) on \(E=[0,\infty)\times[0,\infty)\), with Lebesgue measure in both coordinates. Each factor is nonnegative, measurable, and has integral $$ \int_0^\infty e^{-t}\,dt =\lim_{b\to\infty}\left[-e^{-t}\right]_0^b =1. $$ The factorization theorem gives $$ \int_E e^{-x}e^{-y}\,d(\lambda\otimes\lambda) =\left(\int_0^\infty e^{-x}\,dx\right) \left(\int_0^\infty e^{-y}\,dy\right) =1\cdot1 =1. $$ The same value can be seen directly in either order. For example, $$ \int_0^\infty\left(\int_0^\infty e^{-x}e^{-y}\,dy\right)dx =\int_0^\infty e^{-x}\cdot1\,dx =1. $$ The positive quadrant has infinite area, but that does not prevent the nonnegative function from having a finite integral. Tonelli’s Theorem permits the iterated calculation, and the factorization computes its value.
Choosing an Order and Checking the Hypotheses
For a region described by inequalities, the practical task is to express each slice accurately. A reliable process is:
These specify which pairs \((x,y)\) are included before any integration bounds are chosen.
Solve the inequalities for the range of the other coordinate, and identify the outer coordinate’s range.
Confirm that every point in the region is included and that points outside it are excluded. If the slice changes form, split the region into pieces.
Use Tonelli for nonnegative measurable functions, including when the value may be infinite; use Fubini for integrable signed functions.
Changing the order is a choice of description, not permission to ignore the theorem’s hypotheses. In particular, for a signed function, separate iterated integrals may fail to exist or may not be interchangeable if absolute integrability is absent. For nonnegative functions on sigma-finite spaces, Tonelli allows either order, but one or both answers may be infinite. Always check measurability and the relevant integrability or nonnegativity condition before treating the two orders as equal.
Check Your Understanding
Use the definitions and results in this tutorial to answer the following questions.
- What does the inner integral in an iterated integral hold fixed?
- Which hypotheses allow Tonelli’s Theorem to justify iterated integration when a nonnegative integral might be infinite?
- Why is absolute integrability needed when applying Fubini’s Theorem to a signed function?
- For the region \(0\leq y\leq x\leq1\), what are the bounds when integrating in \(x\) first?
- How does extending a function by zero outside a measurable region turn integration over that region into a product-space integral?