Calculus III: Iterated Integrals and Fubini's Theorem

In engineering fields, from calculating heat transfer in a plate to determining the stress over a beam's cross-section, you frequently encounter problems requiring integration over two-dimensional regions. Iterated integrals and Fubini's Theorem are the core tools that make these calculations tractable, allowing you to break down a complex double integral into a sequence of simpler single-variable integrals. Mastering these concepts is essential for modeling continuous systems in mechanical, civil, and electrical engineering, as well as for success in advanced calculus courses and certification exams.

Foundations of Double and Iterated Integrals

A double integral, denoted ${\iint }_{R}f(x,y)dA$, represents the signed volume under the surface $z=f(x,y)$ over a planar region $R$. When $f(x,y)=1$, the double integral yields the area of $R$. The direct evaluation of a double integral from its Riemann sum definition is often impractical. Instead, we use an iterated integral, which computes the volume by integrating sequentially with respect to one variable at a time while treating the other as a constant. For a rectangular region $R=[a,b]\times [c,d]$, this looks like integrating first in $y$ then in $x$: ${\int }_a^b\left({\int }_c^{d}f(x,y)dy\right)dx$. The key question is: when does this iterative process equal the original double integral? This is precisely what Fubini's Theorem addresses.

Fubini's Theorem: The Bridge to Practical Computation

Fubini's Theorem provides the conditions under which a double integral can be evaluated as an iterated integral. Its standard statement for rectangular regions is: If $f(x,y)$ is continuous on the rectangle $R=[a,b]\times [c,d]$, then the double integral equals either iterated integral order. Mathematically, this is expressed as:

$${\iint }_{R}f(x,y)dA={\int }_a^b{\int }_c^{d}f(x,y)dydx={\int }_c^d{\int }_a^{b}f(x,y)dxdy.$$

The theorem extends to non-rectangular regions, but the limits of integration become functions that describe the region's boundary. The critical condition is that $f$ must be integrable over $R$, which is guaranteed if $f$ is continuous or has only a finite number of discontinuities on a "well-behaved" region. For engineering applications, you will almost always work with continuous functions over closed, bounded regions, so Fubini's Theorem applies. This theorem is your license to compute by choosing the most convenient order of integration.

Evaluating Iterated Integrals: A Step-by-Step Methodology

Evaluating an iterated integral involves performing integration from the inside out, treating the outer variable as a constant during the inner integration. Consider the region $R$ bounded by $y=x^2$ and $y=1$ for $0\le x\le 1$. To compute ${\iint }_R(x+2y)dA$ with integration order $dydx$, you first describe $R$ as: $0\le x\le 1$ and $x^2\le y\le 1$.

The iterated integral is: $${\int }_0^1{\int }_{x^2}^1(x+2y)dydx.$$

Step 1: Integrate with respect to $y$. Treat $x$ as constant: $${\int }_{x^2}^1(x+2y)dy={\left[xy+y^2\right]}_{y=x^2}^{y=1}=(x\cdot 1+1^2)-(x\cdot x^2+(x^2{)}^2)=(x+1)-(x^3+x^4).$$

Step 2: Integrate the result with respect to $x$: $${\int }_0^1\left(x+1-x^3-x^4\right)dx={\left[\frac{x^2}{2}+x-\frac{x^4}{4}-\frac{x^5}{5}\right]}_0^1=\frac{1}{2}+1-\frac{1}{4}-\frac{1}{5}=\frac{10+20-5-4}{20}=\frac{21}{20}.$$

This step-by-step approach, where you carefully handle the variable limits, is fundamental. In exam settings, a common strategy is to sketch the region first to avoid limit errors.

Changing the Order of Integration

Sometimes, the given order of integration leads to a difficult or impossible inner integral. Changing the order of integration means swapping $dxdy$ to $dydx$ or vice versa, which can simplify the antiderivative. This requires re-determining the limits of integration based on the geometry of the region $R$.

The process is methodical:

1. Identify the current limits from the iterated integral. For ${\int }_0^1{\int }_{x^2}^{1}f(x,y)dydx$, the limits describe $R$: $0\le x\le 1$ and $x^2\le y\le 1$.
2. Sketch the region. Plot the curves $y=x^2$ and $y=1$ from $x=0$ to $x=1$. The region is bounded on the left by $x=0$, on the right by $x=1$, below by $y=x^2$, and above by $y=1$.
3. Redescribe $R$ for the opposite order. To integrate $dxdy$, you need $y$ as the outer variable. The total range in $y$ is from the minimum to maximum $y$ in $R$: from $y=0$ to $y=1$. For a fixed $y$, $x$ runs from the left boundary to the right boundary. The left boundary is $x=0$, and the right boundary is given by $y=x^2$, so $x=\sqrt{y}$. Thus, $R$ is described as $0\le y\le 1$ and $0\le x\le \sqrt{y}$.
4. Write the new iterated integral: ${\int }_0^1{\int }_0^{\sqrt{y}}f(x,y)dxdy$.

This geometric reinterpretation is a powerful skill. On engineering exams, problems often test your ability to change order by providing an integral that is hard to evaluate in its original form.

Strategic Order Selection for Simplification

Changing the order of integration can transform an intractable problem into a straightforward one. It simplifies computation significantly in these common scenarios:

• When the inner integral has no elementary antiderivative in the given order. For example, ${\int }_0^1{\int }_y^1{e}^{x^2}dxdy$ is impossible to integrate as is because $e^{x^2}$ has no elementary antiderivative with respect to $x$. Changing the order by sketching the region $0\le y\le 1$, $y\le x\le 1$ (a triangle) and rewriting as ${\int }_0^1{\int }_0^x{e}^{x^2}dydx={\int }_0^{1}x{e}^{x^2}dx$ makes it solvable via substitution.
• When the integrand separates more cleanly. If $f(x,y)=g(x)h(y)$, choosing limits that are constants for the inner integral can allow you to factor the integrals, as in ${\int }_a^{b}g(x)dx\cdot {\int }_c^{d}h(y)dy$.
• In engineering applications like calculating moments of inertia, where the region geometry might be simpler to describe in one order. For a region bounded by curves like $x=y^2$ and $x=2y$, integrating in $dxdy$ order often leads to constant limits for the inner integral, reducing algebraic complexity.

Always assess the integrand and the region before committing to an order. This strategic thinking saves time and reduces errors in both homework and timed assessments.

Common Pitfalls

1. Incorrect Limit Setup from Misreading the Region: The most frequent error is deriving wrong new limits when changing order without sketching the region. Correction: Always draw the region defined by the original limits. Label the bounding curves and find their intersections. Describe the region in the new order by stating the outer variable's constant range and the inner variable's functional range.
2. Forgetting the Conditions of Fubini's Theorem: Applying iterated integration to a function with a discontinuity curve inside the region without adjusting the limits can lead to incorrect results. Correction: Verify the function is continuous over the closed region. If discontinuities exist, split the region into subregions where the function is continuous.
3. Algebra Errors in Variable Manipulation: When changing order, solving for $x$ in terms of $y$ (or vice versa) can lead to mistakes, such as taking the wrong branch of a square root. Correction: From your sketch, identify which function is the left boundary and which is the right (or top and bottom). For the region $0\le y\le 1$, $-\sqrt{y}\le x\le \sqrt{y}$, if the physical context requires a positive $x$, you would use $0\le x\le \sqrt{y}$.
4. Integrating in the Wrong Sequence: Performing the outer integration before the inner one is a procedural error. Correction: Remember the golden rule: integrate from the inside out. The inner differential tells you which variable to integrate first while treating all others as constants.

Summary

• Fubini's Theorem is the workhorse for evaluating double integrals, stating that under continuity conditions, ${\iint }_{R}f(x,y)dA$ can be computed as an iterated integral in either order.
• Evaluating an iterated integral requires integrating from the inside out, carefully applying the fundamental theorem of calculus to limits that may be functions.
• Changing the order of integration is a strategic technique that involves re-describing the region of integration geometrically to swap the sequence $dxdy$ to $dydx$ or vice versa.
• Determining new limits mandates sketching the region defined by the original limits and then expressing the boundaries in terms of the new outer variable.
• Order simplification is most beneficial when the inner integral is overly complex or impossible in the original order, often turning a daunting problem into a routine calculation.