◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Switching the Order of Integration

Formulas & Definitions
Core Definition

Switching the order of integration is an analytical technique used to transform a double integral between Type I (dy dxdy \, dxdydx) and Type II (dx dydx \, dydxdy) representations across a general planar region DDD. Because the inner limits are functions of the outer variable, one cannot simply exchange the integral signs; instead, one must algebraically invert the boundary curves and re-slice the geometric domain:

∫ab∫g1(x)g2(x)f(x,y) dy dx=∫cd∫h1(y)h2(y)f(x,y) dx dy\int_{a}^{b} \int_{g_1(x)}^{g_2(x)} f(x, y) \, dy \, dx = \int_{c}^{d} \int_{h_1(y)}^{h_2(y)} f(x, y) \, dx \, dy∫ab​∫g1​(x)g2​(x)​f(x,y)dydx=∫cd​∫h1​(y)h2​(y)​f(x,y)dxdy

Concept Explanation

In multivariable calculus, many iterated integrals cannot be evaluated in their original formulation because the inner single-variable integrand possesses no elementary antiderivative (e.g., ∫ex2 dx\int e^{x^2} \, dx∫ex2dx, ∫sin⁡(y2) dy\int \sin(y^2) \, dy∫sin(y2)dy, or ∫sin⁡yy dy\int \frac{\sin y}{y} \, dy∫ysiny​dy). By reversing the slicing direction across the planar region, the opposite variable is integrated first, frequently producing an algebraic multiplier that unlocks standard uuu-substitution or power-rule integration.

The 4-Step Re-slicing Protocol

To safely and rigorously reverse the order of integration, always execute the following four steps:

  1. Step 1: Extract the Original Boundary Inequalities
    Read the given iterated integral from outer constants to inner functions:

    Type I:a≤x≤b,g1(x)≤y≤g2(x)\text{Type I:} \quad a \le x \le b, \quad g_1(x) \le y \le g_2(x)Type I:a≤x≤b,g1​(x)≤y≤g2​(x)
    Type II:c≤y≤d,h1(y)≤x≤h2(y)\text{Type II:} \quad c \le y \le d, \quad h_1(y) \le x \le h_2(y)Type II:c≤y≤d,h1​(y)≤x≤h2​(y)

  2. Step 2: Accurately Sketch the Region DDD
    Plot all boundary lines and curves (x=a,x=b,y=g1(x),y=g2(x)x = a, x = b, y = g_1(x), y = g_2(x)x=a,x=b,y=g1​(x),y=g2​(x)) on the Cartesian coordinate plane. Identify and label all vertex intersection points. Shade the enclosed region DDD. Skipping this sketch is the leading cause of student boundary calculation errors.

  3. Step 3: Algebraically Invert the Boundary Functions
    Solve each inner boundary equation for the opposite variable:

    • If originally given y=g(x)y = g(x)y=g(x), solve for xxx in terms of yyy: x=g−1(y)=h(y)x = g^{-1}(y) = h(y)x=g−1(y)=h(y).
    • Pay careful attention to algebraic domain restrictions (e.g., y=x2  ⟹  x=+yy = x^2 \implies x = +\sqrt{y}y=x2⟹x=+y​ for the right branch, x=−yx = -\sqrt{y}x=−y​ for the left branch).
  4. Step 4: Re-slice the Domain and Formulate New Limits

    • If transforming from Type I to Type II (Vertical →\to→ Horizontal Slices): Find the minimum and maximum yyy-values across the entire shaded region to establish the new outer constant limits c≤y≤dc \le y \le dc≤y≤d. Then cast a horizontal arrow from left to right through the region to identify the left entrance curve x=h1(y)x = h_1(y)x=h1​(y) and the right exit curve x=h2(y)x = h_2(y)x=h2​(y).
    • If transforming from Type II to Type I (Horizontal →\to→ Vertical Slices): Find the minimum and maximum xxx-values across the shaded region for the outer constant limits a≤x≤ba \le x \le ba≤x≤b. Cast a vertical arrow upward from bottom to top to identify the bottom entrance curve y=g1(x)y = g_1(x)y=g1​(x) and top exit curve y=g2(x)y = g_2(x)y=g2​(x).
    • If the entrance or exit curves change algebraic definition over the domain, split the integral into the sum of two or more disjoint iterated integrals.

Engineering Applications

  • Aerospace & Mechanical Engineering: Computing structural torsional polar moments of inertia Jz=∬D(x2+y2) dAJ_z = \iint_D (x^2 + y^2) \, dAJz​=∬D​(x2+y2)dA and mass distributions for tapered aerodynamic airfoils and turbine blades where chordwise integration is algebraically intractable.
  • Ocean & Naval Engineering: Evaluating waterplane second moments of area Ixx=∬Dy2 dAI_{xx} = \iint_D y^2 \, dAIxx​=∬D​y2dA to calculate the transverse metacentric height GM‾\overline{GM}GM and righting arm buoyancy stability for ship hulls with parabolic waterlines.
  • Electrical & Signal Engineering: Evaluating 2D spatial convolution integrals and aperture field distributions where non-separable kernel boundaries become elementary after coordinate axis reversal.

Example Problems

Example 1: Resolving a Non-Elementary Antiderivative (ex2e^{x^2}ex2)

Evaluate the double integral by switching the order of integration:

I=∫01∫y1ex2 dx dyI = \int_{0}^{1} \int_{y}^{1} e^{x^2} \, dx \, dyI=∫01​∫y1​ex2dxdy

Solution

Step 1: Extract Original Boundaries (Type II)
Outer limits: 0≤y≤10 \le y \le 10≤y≤1.
Inner limits: y≤x≤1y \le x \le 1y≤x≤1.
Notice that ∫ex2 dx\int e^{x^2} \, dx∫ex2dx has no elementary antiderivative; switching order is mandatory.

Step 2: Sketch the Planar Region
The region DDD is bounded on the left by the slanted line x=yx = yx=y (or y=xy = xy=x), on the right by the vertical line x=1x = 1x=1, and below by the horizontal line y=0y = 0y=0. The vertices are (0,0),(1,0),(0, 0), (1, 0),(0,0),(1,0), and (1,1)(1, 1)(1,1), forming a right triangle in the first quadrant.

Step 3: Re-slice Vertically (Type I Formulation)

  • Outer constant limits for xxx: The triangle extends from x=0x = 0x=0 to x=1x = 1x=1.
  • Inner functional limits for yyy: For any fixed x∈[0,1]x \in [0, 1]x∈[0,1], a vertical ray enters the triangle at the bottom line y=0y = 0y=0 and exits at the top line y=xy = xy=x. Thus:
    D={(x,y)∈R2∣0≤x≤1, 0≤y≤x}D = \{(x, y) \in \mathbb{R}^2 \mid 0 \le x \le 1, \, 0 \le y \le x\}D={(x,y)∈R2∣0≤x≤1,0≤y≤x}

Step 4: Set up and Evaluate the Type I Integral

I=∫01∫0xex2 dy dxI = \int_{0}^{1} \int_{0}^{x} e^{x^2} \, dy \, dxI=∫01​∫0x​ex2dydx

Evaluate Inner Integral (dydydy):
Since ex2e^{x^2}ex2 does not depend on yyy, it acts as a constant:

∫0xex2 dy=ex2[y]0x=xex2\int_{0}^{x} e^{x^2} \, dy = e^{x^2} [y]_{0}^{x} = x e^{x^2}∫0x​ex2dy=ex2[y]0x​=xex2

Evaluate Outer Integral (dxdxdx):
Now use uuu-substitution: let u=x2  ⟹  du=2x dx  ⟹  x dx=du2u = x^2 \implies du = 2x \, dx \implies x \, dx = \frac{du}{2}u=x2⟹du=2xdx⟹xdx=2du​.
When x=0,u=0x = 0, u = 0x=0,u=0; when x=1,u=1x = 1, u = 1x=1,u=1.

I=∫0112eu du=[12eu]01=12(e1−e0)=12(e−1)I = \int_{0}^{1} \frac{1}{2} e^u \, du = \left[ \frac{1}{2} e^u \right]_{0}^{1} = \frac{1}{2}(e^1 - e^0) = \frac{1}{2}(e - 1)I=∫01​21​eudu=[21​eu]01​=21​(e1−e0)=21​(e−1)

Example 2: Resolving a Rational Transcendental Integrand (sin⁡yy\frac{\sin y}{y}ysiny​)

Evaluate the double integral by switching the order of integration:

I=∫0π∫xπsin⁡yy dy dxI = \int_{0}^{\pi} \int_{x}^{\pi} \frac{\sin y}{y} \, dy \, dxI=∫0π​∫xπ​ysiny​dydx

Solution

Step 1: Extract Original Boundaries (Type I)
Outer limits: 0≤x≤π0 \le x \le \pi0≤x≤π.
Inner limits: x≤y≤πx \le y \le \pix≤y≤π.
As written, the inner integral requires finding ∫sin⁡yy dy\int \frac{\sin y}{y} \, dy∫ysiny​dy, which is the non-elementary sine integral function Si(y)\text{Si}(y)Si(y).

Step 2: Sketch the Planar Region
The lower boundary is the line y=xy = xy=x, the upper boundary is the horizontal line y=πy = \piy=π, and the left boundary is the vertical line x=0x = 0x=0. The region is an upper triangle with vertices (0,0),(0,π),(0, 0), (0, \pi),(0,0),(0,π), and (π,π)(\pi, \pi)(π,π).

Step 3: Re-slice Horizontally (Type II Formulation)

  • Outer constant limits for yyy: The triangle extends vertically from y=0y = 0y=0 to y=πy = \piy=π.
  • Inner functional limits for xxx: For any fixed y∈[0,π]y \in [0, \pi]y∈[0,π], a horizontal ray enters at the vertical line x=0x = 0x=0 and exits at the line x=yx = yx=y. Thus:
    D={(x,y)∈R2∣0≤y≤π, 0≤x≤y}D = \{(x, y) \in \mathbb{R}^2 \mid 0 \le y \le \pi, \, 0 \le x \le y\}D={(x,y)∈R2∣0≤y≤π,0≤x≤y}

Step 4: Set up and Evaluate the Type II Integral

I=∫0π∫0ysin⁡yy dx dyI = \int_{0}^{\pi} \int_{0}^{y} \frac{\sin y}{y} \, dx \, dyI=∫0π​∫0y​ysiny​dxdy

Evaluate Inner Integral (dxdxdx):

∫0ysin⁡yy dx=sin⁡yy[x]0y=sin⁡yy(y−0)=sin⁡y\int_{0}^{y} \frac{\sin y}{y} \, dx = \frac{\sin y}{y} [x]_{0}^{y} = \frac{\sin y}{y}(y - 0) = \sin y∫0y​ysiny​dx=ysiny​[x]0y​=ysiny​(y−0)=siny

Evaluate Outer Integral (dydydy):

∫0πsin⁡y dy=[−cos⁡y]0π=−cos⁡(π)−(−cos⁡(0))=−(−1)+1=2\int_{0}^{\pi} \sin y \, dy = [-\cos y]_{0}^{\pi} = -\cos(\pi) - (-\cos(0)) = -(-1) + 1 = 2∫0π​sinydy=[−cosy]0π​=−cos(π)−(−cos(0))=−(−1)+1=2


Connections & References

  • Parent Concepts: Double Integrals over a General Region, Multiple Integrals
  • Sub-concepts & Extensions:
    • Double Integrals in Polar Coordinates
    • Changing the Order of Integration in Triple Integrals
    • Change of Variables in Multiple Integrals
  • Course Links:
    • Calc-3 Session 14 - Double Integrals
    • Calc-3 Session 15 - Double Integrals in Polar Coordinates
    • Calc-3 Session 16 - Triple Integrals