MATH-283 (Calculus III) Curriculum

Session 14 - Double Integrals

Handout Worksheet
Fall Semester 2026
Session Objectives & Overview

By the end of this session, you will be able to:

  • Formulate, interpret, and evaluate iterated double integrals over rectangular domains using Fubini's Theorem, apply the product separability theorem ∬Rg(x)h(y) dA=(∫abg(x) dx)(∫cdh(y) dy)\iint_R g(x)h(y)\,dA = \left(\int_a^b g(x)\,dx\right)\left(\int_c^d h(y)\,dy\right)∬R​g(x)h(y)dA=(∫ab​g(x)dx)(∫cd​h(y)dy), and compute the spatial average value favg=1Area(R)∬Rf(x,y) dAf_{\text{avg}} = \frac{1}{\text{Area}(R)}\iint_R f(x, y)\,dAfavg​=Area(R)1β€‹βˆ¬R​f(x,y)dA across planar design envelopes.
  • Classify, sketch, and evaluate double integrals over general planar regions by distinguishing between Type I (vertical slices: g1(x)≀y≀g2(x)g_1(x) \le y \le g_2(x)g1​(x)≀y≀g2​(x)) and Type II (horizontal slices: h1(y)≀x≀h2(y)h_1(y) \le x \le h_2(y)h1​(y)≀x≀h2​(y)) geometries, determining curve intersections to establish outer constant bounds.
  • Execute the 4-step order-switching protocol (Extract β†’\toβ†’ Sketch β†’\toβ†’ Invert β†’\toβ†’ Re-slice) to convert between Type I and Type II iterated integrals, partition compound multi-curve boundaries into disjoint sub-regions, and unlock closed-form solutions for non-elementary integrands such as ex2e^{x^2}ex2, sin⁑yy\frac{\sin y}{y}ysiny​, and 1+y3\sqrt{1+y^3}1+y3​.

Micro-Lecture

1. Spatial Accumulation & The Volume Under a Surface

  • Mathematical Foundation: In single-variable calculus, definite integrals accumulate quantities along a 1-dimensional interval ∫abf(x) dx\int_a^b f(x)\,dx∫ab​f(x)dx. When physical phenomena vary across a 2-dimensional surface, integration generalizes to the double integral ∬Df(x,y) dA\iint_D f(x, y)\,dA∬D​f(x,y)dA, defined in Multiple Integrals as the infinite limit of a double Riemann sum. Geometrically, if f(x,y)β‰₯0f(x, y) \ge 0f(x,y)β‰₯0, the double integral calculates the exact volume of the cylindrical solid lying between the planar region DDD and the surface z=f(x,y)z = f(x, y)z=f(x,y). If f(x,y)=1f(x, y) = 1f(x,y)=1, the integral yields the planar surface area Area(D)=∬D1 dA\text{Area}(D) = \iint_D 1\,dAArea(D)=∬D​1dA.
  • Engineering Context: Civil and structural engineers model concrete bridge foundations subject to non-uniform soil reaction pressure; integrating the contact stress distribution Οƒ(x,y)\sigma(x, y)Οƒ(x,y) over the footing footprint calculates the total structural bearing load supported by the ground.

2. Constant Domains & Fubini's Iterated Slicing

  • Mathematical Foundation: On a closed rectangular domain R=[a,b]Γ—[c,d]R = [a, b] \times [c, d]R=[a,b]Γ—[c,d], all boundary edges are constant lines. As proven in Double Integrals over a Rectangular Region, Fubini's Theorem establishes that double integrals equal iterated single integrals evaluated in either coordinate order: ∬Rf dA=∫cd∫abf dx dy=∫ab∫cdf dy dx\iint_R f\,dA = \int_c^d \int_a^b f\,dx\,dy = \int_a^b \int_c^d f\,dy\,dx∬R​fdA=∫cdβ€‹βˆ«ab​fdxdy=∫abβ€‹βˆ«cd​fdydx. Choosing the optimal slicing direction eliminates tedious integration by parts when one variable acts as a chain-rule multiplier. Furthermore, dividing the double integral by the rectangular area produces the spatial average value favg=1Area(R)∬Rf dAf_{\text{avg}} = \frac{1}{\text{Area}(R)}\iint_R f\,dAfavg​=Area(R)1β€‹βˆ¬R​fdA.
  • Engineering Context: Entertainment and electrical engineers evaluate radiant light emissions from flat-panel LED displays, using spatial averages to calculate overall luminaire efficiency and mitigate optical hotspot banding across modular video walls.

3. General Boundaries & Slicing Inversion

  • Mathematical Foundation: When a domain is bounded by curved or angled profiles, it is classified as a general region (formulated in Double Integrals over a General Region). In Type I domains (vertical slices), outer limits are constants x∈[a,b]x \in [a, b]x∈[a,b] and inner limits are functions y∈[g1(x),g2(x)]y \in [g_1(x), g_2(x)]y∈[g1​(x),g2​(x)]. In Type II domains (horizontal slices), outer limits are constants y∈[c,d]y \in [c, d]y∈[c,d] and inner limits are functions x∈[h1(y),h2(y)]x \in [h_1(y), h_2(y)]x∈[h1​(y),h2​(y)]. When an integrand lacks an elementary antiderivative (e.g., ex2e^{x^2}ex2), engineers employ Switching the Order of Integration via the 4-step re-slicing protocol to invert boundary equations and evaluate the system.
  • Engineering Context: Ocean and aerospace engineers evaluate moments of inertia and hydrostatic waterplane stability curves on ship hulls and tapered composite wings, switching integration orders to eliminate algebraic radical singularities in automated finite element routines.

Skill Block 1

Worked Example A: Polynomial Iterated Double Integral over a Rectangle

Evaluate the double integral:

∬R(4x3+6xy2) dAwhereΒ R=[0,2]Γ—[1,3]\iint_R (4x^3 + 6xy^2) \, dA \quad \text{where } R = [0, 2] \times [1, 3]∬R​(4x3+6xy2)dAwhereΒ R=[0,2]Γ—[1,3]

Worked Solution:
  1. Formulate the Iterated Integral:
    State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):

    ∬Rf(x,y) dA=∫ab∫cdf(x,y) dy dx\iint_R f(x, y) \, dA = \int_{a}^{b} \int_{c}^{d} f(x, y) \, dy \, dx∬R​f(x,y)dA=∫abβ€‹βˆ«cd​f(x,y)dydx
    Substitute the domain limits x∈[0,2]x \in [0, 2]x∈[0,2] and y∈[1,3]y \in [1, 3]y∈[1,3]:
    I=∫02∫13(4x3+6xy2) dy dxI = \int_{0}^{2} \int_{1}^{3} (4x^3 + 6xy^2) \, dy \, dxI=∫02β€‹βˆ«13​(4x3+6xy2)dydx

  2. Evaluate the Inner Integral (dydydy):
    Treat xxx as a constant and integrate with respect to yyy:

    ∫13(4x3+6xy2) dy=[4x3y+2xy3]y=1y=3\int_{1}^{3} (4x^3 + 6xy^2) \, dy = \left[ 4x^3 y + 2xy^3 \right]_{y=1}^{y=3}∫13​(4x3+6xy2)dy=[4x3y+2xy3]y=1y=3​
    Substitute the boundary limits:
    =(4x3(3)+2x(3)3)βˆ’(4x3(1)+2x(1)3)= \Big(4x^3(3) + 2x(3)^3\Big) - \Big(4x^3(1) + 2x(1)^3\Big)=(4x3(3)+2x(3)3)βˆ’(4x3(1)+2x(1)3)
    =(12x3+54x)βˆ’(4x3+2x)=8x3+52x= (12x^3 + 54x) - (4x^3 + 2x) = 8x^3 + 52x=(12x3+54x)βˆ’(4x3+2x)=8x3+52x

  3. Evaluate the Outer Integral (dxdxdx):
    Integrate the resulting single-variable function with respect to xxx:

    I=∫02(8x3+52x) dx=[2x4+26x2]02I = \int_{0}^{2} (8x^3 + 52x) \, dx = \left[ 2x^4 + 26x^2 \right]_{0}^{2}I=∫02​(8x3+52x)dx=[2x4+26x2]02​
    Evaluate at the limits:
    I=(2(2)4+26(2)2)βˆ’0=2(16)+26(4)=32+104=136I = \Big(2(2)^4 + 26(2)^2\Big) - 0 = 2(16) + 26(4) = 32 + 104 = 136I=(2(2)4+26(2)2)βˆ’0=2(16)+26(4)=32+104=136

Worked Example B: Strategic Order Selection via Fubini's Shortcut

Evaluate the double integral:

∬Rxsin⁑(xy) dAwhereΒ R=[0,1]Γ—[0,Ο€]\iint_R x \sin(xy) \, dA \quad \text{where } R = [0, 1] \times [0, \pi]∬R​xsin(xy)dAwhereΒ R=[0,1]Γ—[0,Ο€]
demonstrating why integrating with respect to yyy first is superior to integrating with respect to xxx first.

Worked Solution:
  1. Analyze Slicing Strategy:
    State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):

    ∬Rf(x,y) dA=∫ab∫cdf(x,y) dy dx=∫cd∫abf(x,y) dx dy\iint_R f(x, y) \, dA = \int_{a}^{b} \int_{c}^{d} f(x, y) \, dy \, dx = \int_{c}^{d} \int_{a}^{b} f(x, y) \, dx \, dy∬R​f(x,y)dA=∫abβ€‹βˆ«cd​f(x,y)dydx=∫cdβ€‹βˆ«ab​f(x,y)dxdy

    • If we integrate with respect to xxx first, ∫xsin⁑(xy) dx\int x \sin(xy) \, dx∫xsin(xy)dx requires integration by parts because xxx appears both outside and inside the sine function.
    • If we integrate with respect to yyy first, xxx acts as a constant multiplier. Since βˆ‚βˆ‚y(xy)=x\frac{\partial}{\partial y}(xy) = xβˆ‚yβˆ‚β€‹(xy)=x, the factor xxx is the exact derivative of the inner argument, yielding an immediate single-step antiderivative!
  2. Formulate and Evaluate Inner Integral (dydydy):
    Set up the iterated integral with dydydy on the inside:

    I=∫01∫0Ο€xsin⁑(xy) dy dxI = \int_{0}^{1} \int_{0}^{\pi} x \sin(xy) \, dy \, dxI=∫01β€‹βˆ«0π​xsin(xy)dydx
    Integrate with respect to yyy, using ∫xsin⁑(xy) dy=βˆ’cos⁑(xy)\int x \sin(xy) \, dy = -\cos(xy)∫xsin(xy)dy=βˆ’cos(xy):
    ∫0Ο€xsin⁑(xy) dy=[βˆ’cos⁑(xy)]y=0y=Ο€=βˆ’cos⁑(Ο€x)βˆ’(βˆ’cos⁑(0))=1βˆ’cos⁑(Ο€x)\int_{0}^{\pi} x \sin(xy) \, dy = \left[ -\cos(xy) \right]_{y=0}^{y=\pi} = -\cos(\pi x) - (-\cos(0)) = 1 - \cos(\pi x)∫0π​xsin(xy)dy=[βˆ’cos(xy)]y=0y=π​=βˆ’cos(Ο€x)βˆ’(βˆ’cos(0))=1βˆ’cos(Ο€x)

  3. Evaluate Outer Integral (dxdxdx):
    Substitute into the outer integral:

    I=∫01(1βˆ’cos⁑(Ο€x)) dx=[xβˆ’sin⁑(Ο€x)Ο€]01I = \int_{0}^{1} \Big(1 - \cos(\pi x)\Big) \, dx = \left[ x - \frac{\sin(\pi x)}{\pi} \right]_{0}^{1}I=∫01​(1βˆ’cos(Ο€x))dx=[xβˆ’Ο€sin(Ο€x)​]01​
    Evaluate at the limits:
    I=(1βˆ’sin⁑(Ο€)Ο€)βˆ’(0βˆ’sin⁑(0)Ο€)=(1βˆ’0)βˆ’(0βˆ’0)=1I = \left( 1 - \frac{\sin(\pi)}{\pi} \right) - \left( 0 - \frac{\sin(0)}{\pi} \right) = (1 - 0) - (0 - 0) = 1I=(1βˆ’Ο€sin(Ο€)​)βˆ’(0βˆ’Ο€sin(0)​)=(1βˆ’0)βˆ’(0βˆ’0)=1

Practice Problem 1: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) Under what mathematical conditions does Fubini's Theorem permit swapping the integration order from dy dxdy \, dxdydx to dx dydx \, dydxdy, and why must all four limits of integration remain constants for a rectangular region?
  • (b) In Worked Example B, explain why integrating with respect to yyy first eliminates the need for integration by parts.
  • (c) Referencing Optoelectronic LED Display Panel Average Irradiance, explain physically why the spatial average irradiance EavgE_{\text{avg}}Eavg​ divides the total double integral by the rectangular surface area Area(R)=(bβˆ’a)(dβˆ’c)\text{Area}(R) = (b-a)(d-c)Area(R)=(bβˆ’a)(dβˆ’c).
Solution Key:
  • (a) State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):
    ∬Rf(x,y) dA=∫cd∫abf(x,y) dx dy=∫ab∫cdf(x,y) dy dx\iint_R f(x, y) \, dA = \int_c^d \int_a^b f(x, y) \, dx \, dy = \int_a^b \int_c^d f(x, y) \, dy \, dx∬R​f(x,y)dA=∫cdβ€‹βˆ«ab​f(x,y)dxdy=∫abβ€‹βˆ«cd​f(x,y)dydx
    This identity is valid for all functions f(x,y)f(x, y)f(x,y) that are continuous on the closed rectangle RRR. The limits are fixed constants because a rectangle in Cartesian space is explicitly defined by intervals x∈[a,b]x \in [a, b]x∈[a,b] and y∈[c,d]y \in [c, d]y∈[c,d].
  • (b) In the expression ∫xsin⁑(xy) dy\int x \sin(xy)\,dy∫xsin(xy)dy, xxx is independent of yyy. By the chain rule, ddy[βˆ’cos⁑(xy)]=xsin⁑(xy)\frac{d}{dy}[-\cos(xy)] = x\sin(xy)dyd​[βˆ’cos(xy)]=xsin(xy). If integrating with respect to xxx first, xxx is the integration variable and would require integration by parts with u=xu = xu=x and dv=sin⁑(xy) dxdv = \sin(xy)\,dxdv=sin(xy)dx.
  • (c) State the Spatial Average Value Formula:
    Eavg=1Area(R)∬RE(x,y) dAE_{\text{avg}} = \frac{1}{\text{Area}(R)} \iint_R E(x, y) \, dAEavg​=Area(R)1β€‹βˆ¬R​E(x,y)dA
    In optoelectronic display engineering, light emission varies across the panel due to heating and circuit resistance. The double integral accumulates total radiant power (W\text{W}W). Dividing by the geometric surface area Area(R)=(0.5)(0.4)=0.2Β m2\text{Area}(R) = (0.5)(0.4) = 0.2\text{ m}^2Area(R)=(0.5)(0.4)=0.2Β m2 normalizes the total power to express the mean flux per unit area (W/m2\text{W/m}^2W/m2), which is used to calibrate LED driver brightness and prevent visual hotspot banding.
Show your steps and calculations here.

Practice Problem 2

Evaluate the double integral:

∬R(3x2βˆ’2y) dAwhereΒ R=[0,2]Γ—[1,3]\iint_R (3x^2 - 2y) \, dA \quad \text{where } R = [0, 2] \times [1, 3]∬R​(3x2βˆ’2y)dAwhereΒ R=[0,2]Γ—[1,3]
and determine the spatial average value favgf_{\text{avg}}favg​ of the function across the region RRR.

Solution Key:
  1. Formulate the Iterated Integral:
    State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):

    ∬Rf(x,y) dA=∫13∫02(3x2βˆ’2y) dx dy\iint_R f(x, y) \, dA = \int_{1}^{3} \int_{0}^{2} (3x^2 - 2y) \, dx \, dy∬R​f(x,y)dA=∫13β€‹βˆ«02​(3x2βˆ’2y)dxdy

  2. Evaluate the Inner Integral (dxdxdx):
    Treat yyy as a constant:

    ∫02(3x2βˆ’2y) dx=[x3βˆ’2yx]x=0x=2=((2)3βˆ’2y(2))βˆ’0=8βˆ’4y\int_{0}^{2} (3x^2 - 2y) \, dx = \left[ x^3 - 2yx \right]_{x=0}^{x=2} = \Big((2)^3 - 2y(2)\Big) - 0 = 8 - 4y∫02​(3x2βˆ’2y)dx=[x3βˆ’2yx]x=0x=2​=((2)3βˆ’2y(2))βˆ’0=8βˆ’4y

  3. Evaluate the Outer Integral (dydydy):

    ∫13(8βˆ’4y) dy=[8yβˆ’2y2]13=(8(3)βˆ’2(3)2)βˆ’(8(1)βˆ’2(1)2)\int_{1}^{3} (8 - 4y) \, dy = \left[ 8y - 2y^2 \right]_{1}^{3} = \Big(8(3) - 2(3)^2\Big) - \Big(8(1) - 2(1)^2\Big)∫13​(8βˆ’4y)dy=[8yβˆ’2y2]13​=(8(3)βˆ’2(3)2)βˆ’(8(1)βˆ’2(1)2)
    =(24βˆ’18)βˆ’(8βˆ’2)=6βˆ’6=0= (24 - 18) - (8 - 2) = 6 - 6 = 0=(24βˆ’18)βˆ’(8βˆ’2)=6βˆ’6=0

  4. Compute the Spatial Average Value:
    State the Spatial Average Value Formula:

    favg=1Area(R)∬Rf(x,y) dAf_{\text{avg}} = \frac{1}{\text{Area}(R)} \iint_R f(x, y) \, dAfavg​=Area(R)1β€‹βˆ¬R​f(x,y)dA
    The area of the rectangular domain is:
    Area(R)=(2βˆ’0)(3βˆ’1)=2Γ—2=4\text{Area}(R) = (2 - 0)(3 - 1) = 2 \times 2 = 4Area(R)=(2βˆ’0)(3βˆ’1)=2Γ—2=4
    favg=14(0)=0f_{\text{avg}} = \frac{1}{4}(0) = 0favg​=41​(0)=0

Show your steps and calculations here.

Practice Problem 3

Evaluate the double integral by selecting the strategically simpler order of integration:

∬Ryexy dAwhereΒ R=[0,1]Γ—[0,2]\iint_R y e^{xy} \, dA \quad \text{where } R = [0, 1] \times [0, 2]∬R​yexydAwhereΒ R=[0,1]Γ—[0,2]

Solution Key:
  1. Formulate the Strategically Optimal Iterated Integral:
    State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):

    ∬Ryexy dA=∫02∫01yexy dx dy\iint_R y e^{xy} \, dA = \int_{0}^{2} \int_{0}^{1} y e^{xy} \, dx \, dy∬R​yexydA=∫02β€‹βˆ«01​yexydxdy

  2. Evaluate the Inner Integral (dxdxdx):
    Holding yyy constant, yyy acts as the exact derivative multiplier of xyxyxy (βˆ‚βˆ‚x(xy)=y\frac{\partial}{\partial x}(xy) = yβˆ‚xβˆ‚β€‹(xy)=y):

    ∫01yexy dx=[exy]x=0x=1=ey(1)βˆ’ey(0)=eyβˆ’1\int_{0}^{1} y e^{xy} \, dx = \left[ e^{xy} \right]_{x=0}^{x=1} = e^{y(1)} - e^{y(0)} = e^y - 1∫01​yexydx=[exy]x=0x=1​=ey(1)βˆ’ey(0)=eyβˆ’1

  3. Evaluate the Outer Integral (dydydy):

    ∫02(eyβˆ’1) dy=[eyβˆ’y]02=(e2βˆ’2)βˆ’(e0βˆ’0)=e2βˆ’2βˆ’1=e2βˆ’3\int_{0}^{2} (e^y - 1) \, dy = \left[ e^y - y \right]_{0}^{2} = \Big(e^2 - 2\Big) - \Big(e^0 - 0\Big) = e^2 - 2 - 1 = e^2 - 3∫02​(eyβˆ’1)dy=[eyβˆ’y]02​=(e2βˆ’2)βˆ’(e0βˆ’0)=e2βˆ’2βˆ’1=e2βˆ’3

Show your steps and calculations here.

Practice Problem 4

Evaluate the double integral:

∬R(2xcos⁑yβˆ’3y2) dAwhereΒ R=[0,3]Γ—[0,Ο€2]\iint_R (2x \cos y - 3y^2) \, dA \quad \text{where } R = [0, 3] \times \left[0, \frac{\pi}{2}\right]∬R​(2xcosyβˆ’3y2)dAwhereΒ R=[0,3]Γ—[0,2π​]

Solution Key:
  1. Formulate the Iterated Integral:
    State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):

    ∬R(2xcos⁑yβˆ’3y2) dA=∫0Ο€/2∫03(2xcos⁑yβˆ’3y2) dx dy\iint_R (2x \cos y - 3y^2) \, dA = \int_{0}^{\pi/2} \int_{0}^{3} (2x \cos y - 3y^2) \, dx \, dy∬R​(2xcosyβˆ’3y2)dA=∫0Ο€/2β€‹βˆ«03​(2xcosyβˆ’3y2)dxdy

  2. Evaluate the Inner Integral (dxdxdx):
    Treat yyy as a constant:

    ∫03(2xcos⁑yβˆ’3y2) dx=[x2cos⁑yβˆ’3y2x]x=0x=3=(9cos⁑yβˆ’9y2)βˆ’0=9cos⁑yβˆ’9y2\int_{0}^{3} (2x \cos y - 3y^2) \, dx = \left[ x^2 \cos y - 3y^2 x \right]_{x=0}^{x=3} = \Big(9 \cos y - 9y^2\Big) - 0 = 9 \cos y - 9y^2∫03​(2xcosyβˆ’3y2)dx=[x2cosyβˆ’3y2x]x=0x=3​=(9cosyβˆ’9y2)βˆ’0=9cosyβˆ’9y2

  3. Evaluate the Outer Integral (dydydy):

    ∫0Ο€/2(9cos⁑yβˆ’9y2) dy=[9sin⁑yβˆ’3y3]0Ο€/2\int_{0}^{\pi/2} (9 \cos y - 9y^2) \, dy = \left[ 9 \sin y - 3y^3 \right]_{0}^{\pi/2}∫0Ο€/2​(9cosyβˆ’9y2)dy=[9sinyβˆ’3y3]0Ο€/2​
    Evaluate at the limits:
    =(9sin⁑(Ο€2)βˆ’3(Ο€2)3)βˆ’0=9(1)βˆ’3(Ο€38)=9βˆ’3Ο€38= \left( 9 \sin\left(\frac{\pi}{2}\right) - 3\left(\frac{\pi}{2}\right)^3 \right) - 0 = 9(1) - 3\left(\frac{\pi^3}{8}\right) = 9 - \frac{3\pi^3}{8}=(9sin(2π​)βˆ’3(2π​)3)βˆ’0=9(1)βˆ’3(8Ο€3​)=9βˆ’83Ο€3​
    (Note: Multiplying over the interval x∈[0,3]x \in [0, 3]x∈[0,3], 3Γ—3y2=9y23 \times 3y^2 = 9y^23Γ—3y2=9y2, whose antiderivative is 3y3=3Ο€383y^3 = \frac{3\pi^3}{8}3y3=83Ο€3​.)


Show your steps and calculations here.

Regroup 1

  • Order Selection Heuristic: How can you determine within 5 seconds whether dy dxdy \, dxdydx or dx dydx \, dydxdy is the faster calculation on a rectangle? (Look for cross-variable products where one variable can act as a chain-rule multiplier).
  • Constant Bounds Integrity: Why is Fubini's Theorem so straightforward over rectangles compared to curved domains? (Because the region is defined entirely by constant endpoints, so switching differentials requires zero change to boundary equations).

Skill Block 2

Worked Example A: Type I General Region (Parabola and Slanted Line)

Evaluate the double integral:

∬D(x+y) dA\iint_D (x + y) \, dA∬D​(x+y)dA
where DDD is the region bounded by the parabola y=x2y = x^2y=x2 and the line y=2xy = 2xy=2x.

Worked Solution:
  1. Find Intersection Points (Outer Bounds):
    Set the two boundary curves equal to locate where they intersect:

    x2=2xβ€…β€ŠβŸΉβ€…β€Šx2βˆ’2x=0β€…β€ŠβŸΉβ€…β€Šx(xβˆ’2)=0x^2 = 2x \implies x^2 - 2x = 0 \implies x(x - 2) = 0x2=2x⟹x2βˆ’2x=0⟹x(xβˆ’2)=0
    The outer constant limits are x=0x = 0x=0 and x=2x = 2x=2.

  2. Establish the Functional Inner Limits:
    Test an intermediate value x=1x = 1x=1: y=2(1)=2y = 2(1) = 2y=2(1)=2 and y=(1)2=1y = (1)^2 = 1y=(1)2=1. Since 2>12 > 12>1, the line y=2xy = 2xy=2x is the upper boundary and the parabola y=x2y = x^2y=x2 is the lower boundary. Thus, DDD is a Type I Region:

    D={(x,y)∈R2∣0≀x≀2, x2≀y≀2x}D = \{(x, y) \in \mathbb{R}^2 \mid 0 \le x \le 2, \, x^2 \le y \le 2x\}D={(x,y)∈R2∣0≀x≀2,x2≀y≀2x}

  3. Formulate and Evaluate Inner Integral (dydydy):
    State the Type I Double Integral Formula (defined in Double Integrals over a General Region):

    ∬Df(x,y) dA=∫ab∫g1(x)g2(x)f(x,y) dy dx\iint_D f(x, y) \, dA = \int_{a}^{b} \int_{g_1(x)}^{g_2(x)} f(x, y) \, dy \, dx∬D​f(x,y)dA=∫abβ€‹βˆ«g1​(x)g2​(x)​f(x,y)dydx
    I=∫02∫x22x(x+y) dy dxI = \int_{0}^{2} \int_{x^2}^{2x} (x + y) \, dy \, dxI=∫02β€‹βˆ«x22x​(x+y)dydx
    Evaluate with respect to yyy, treating xxx as constant:
    ∫x22x(x+y) dy=[xy+12y2]y=x2y=2x=(x(2x)+12(2x)2)βˆ’(x(x2)+12(x2)2)\int_{x^2}^{2x} (x + y) \, dy = \left[ xy + \frac{1}{2}y^2 \right]_{y=x^2}^{y=2x} = \Big(x(2x) + \frac{1}{2}(2x)^2\Big) - \Big(x(x^2) + \frac{1}{2}(x^2)^2\Big)∫x22x​(x+y)dy=[xy+21​y2]y=x2y=2x​=(x(2x)+21​(2x)2)βˆ’(x(x2)+21​(x2)2)
    =(2x2+2x2)βˆ’(x3+12x4)=4x2βˆ’x3βˆ’12x4= (2x^2 + 2x^2) - (x^3 + \frac{1}{2}x^4) = 4x^2 - x^3 - \frac{1}{2}x^4=(2x2+2x2)βˆ’(x3+21​x4)=4x2βˆ’x3βˆ’21​x4

  4. Evaluate Outer Integral (dxdxdx):

    ∫02(4x2βˆ’x3βˆ’12x4) dx=[43x3βˆ’14x4βˆ’110x5]02\int_{0}^{2} \left( 4x^2 - x^3 - \frac{1}{2}x^4 \right) \, dx = \left[ \frac{4}{3}x^3 - \frac{1}{4}x^4 - \frac{1}{10}x^5 \right]_{0}^{2}∫02​(4x2βˆ’x3βˆ’21​x4)dx=[34​x3βˆ’41​x4βˆ’101​x5]02​
    Substitute x=2x = 2x=2:
    =(43(8)βˆ’14(16)βˆ’110(32))βˆ’0=323βˆ’4βˆ’165= \left( \frac{4}{3}(8) - \frac{1}{4}(16) - \frac{1}{10}(32) \right) - 0 = \frac{32}{3} - 4 - \frac{16}{5}=(34​(8)βˆ’41​(16)βˆ’101​(32))βˆ’0=332β€‹βˆ’4βˆ’516​
    Find a common denominator (151515):
    =16015βˆ’6015βˆ’4815=5215= \frac{160}{15} - \frac{60}{15} - \frac{48}{15} = \frac{52}{15}=15160β€‹βˆ’1560β€‹βˆ’1548​=1552​

Worked Example B: Type II General Region (Horizontal Slices)

Evaluate the double integral:

∬D(x+y) dA\iint_D (x + y) \, dA∬D​(x+y)dA
where DDD is the region bounded by x=y2x = y^2x=y2 and x=yx = yx=y from y=0y = 0y=0 to y=2y = 2y=2.

Worked Solution:
  1. Identify the Type II Boundary Structure:
    Because the boundary equations explicitly express xxx as functions of yyy, this region is naturally formulated as a Type II Region (defined in Double Integrals over a General Region):

    • Outer constant limits for yyy: 0≀y≀20 \le y \le 20≀y≀2
    • Inner limits for xxx: Left boundary is x=y2x = y^2x=y2, right boundary is x=yx = yx=y
  2. Formulate the Iterated Integral:
    State the Type II Double Integral Formula:

    ∬Df(x,y) dA=∫cd∫h1(y)h2(y)f(x,y) dx dy\iint_D f(x, y) \, dA = \int_{c}^{d} \int_{h_1(y)}^{h_2(y)} f(x, y) \, dx \, dy∬D​f(x,y)dA=∫cdβ€‹βˆ«h1​(y)h2​(y)​f(x,y)dxdy
    I=∫02∫y2y(x+y) dx dyI = \int_{0}^{2} \int_{y^2}^{y} (x + y) \, dx \, dyI=∫02β€‹βˆ«y2y​(x+y)dxdy

  3. Evaluate Inner Integral (dxdxdx):
    Treat yyy as a constant:

    ∫y2y(x+y) dx=[12x2+yx]x=y2x=y=(12y2+y(y))βˆ’(12(y2)2+y(y2))\int_{y^2}^{y} (x + y) \, dx = \left[ \frac{1}{2}x^2 + yx \right]_{x=y^2}^{x=y} = \left( \frac{1}{2}y^2 + y(y) \right) - \left( \frac{1}{2}(y^2)^2 + y(y^2) \right)∫y2y​(x+y)dx=[21​x2+yx]x=y2x=y​=(21​y2+y(y))βˆ’(21​(y2)2+y(y2))
    =32y2βˆ’y3βˆ’12y4= \frac{3}{2}y^2 - y^3 - \frac{1}{2}y^4=23​y2βˆ’y3βˆ’21​y4

  4. Evaluate Outer Integral (dydydy):

    I=∫02(32y2βˆ’y3βˆ’12y4) dy=[12y3βˆ’14y4βˆ’110y5]02I = \int_{0}^{2} \left( \frac{3}{2}y^2 - y^3 - \frac{1}{2}y^4 \right) \, dy = \left[ \frac{1}{2}y^3 - \frac{1}{4}y^4 - \frac{1}{10}y^5 \right]_{0}^{2}I=∫02​(23​y2βˆ’y3βˆ’21​y4)dy=[21​y3βˆ’41​y4βˆ’101​y5]02​
    Evaluate at y=2y = 2y=2:
    =(12(8)βˆ’14(16)βˆ’110(32))βˆ’0=4βˆ’4βˆ’3.2=βˆ’3.2=βˆ’165= \left( \frac{1}{2}(8) - \frac{1}{4}(16) - \frac{1}{10}(32) \right) - 0 = 4 - 4 - 3.2 = -3.2 = -\frac{16}{5}=(21​(8)βˆ’41​(16)βˆ’101​(32))βˆ’0=4βˆ’4βˆ’3.2=βˆ’3.2=βˆ’516​

Practice Problem 5: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) What is the structural distinction between a Type I region and a Type II region in terms of which variable has constant outer limits versus functional inner limits?
  • (b) In Worked Example A, explain how setting x2=2xx^2 = 2xx2=2x establishes the outer constant boundaries, and how testing an intermediate point confirms which curve is on top.
  • (c) Referencing Hydrostatic Pressure Load on a Parabolic Dam Floodgate, explain why fluid pressure depth is integrated along horizontal slices (Type II) rather than vertical slices.
Solution Key:
  • (a) State the Classification of General Planar Regions (defined in Double Integrals over a General Region):
    • Type I (Vertical Slices): D={(x,y)∣a≀x≀b, g1(x)≀y≀g2(x)}D = \{(x, y) \mid a \le x \le b, \, g_1(x) \le y \le g_2(x)\}D={(x,y)∣a≀x≀b,g1​(x)≀y≀g2​(x)}, where xxx has constant limits and yyy has functional limits.
    • Type II (Horizontal Slices): D={(x,y)∣c≀y≀d, h1(y)≀x≀h2(y)}D = \{(x, y) \mid c \le y \le d, \, h_1(y) \le x \le h_2(y)\}D={(x,y)∣c≀y≀d,h1​(y)≀x≀h2​(y)}, where yyy has constant limits and xxx has functional limits.
  • (b) Setting g1(x)=g2(x)β€…β€ŠβŸΉβ€…β€Šx2=2xβ€…β€ŠβŸΉβ€…β€Šx(xβˆ’2)=0g_1(x) = g_2(x) \implies x^2 = 2x \implies x(x - 2) = 0g1​(x)=g2​(x)⟹x2=2x⟹x(xβˆ’2)=0 finds the geometric intersection points (0,0)(0, 0)(0,0) and (2,4)(2, 4)(2,4), which bound the domain horizontally. Testing x=1∈(0,2)x = 1 \in (0, 2)x=1∈(0,2) yields g2(1)=2g_2(1) = 2g2​(1)=2 and g1(1)=1g_1(1) = 1g1​(1)=1; since 2>12 > 12>1, y=2xy = 2xy=2x is strictly above y=x2y = x^2y=x2 over the entire interior.
  • (c) In hydraulic engineering, hydrostatic pressure is governed by fluid depth below the surface. Because the floodgate is bounded on the left by x=y2x = y^2x=y2 and on the right by x=4x = 4x=4, a Type II horizontal slice lies at a uniform water depth (2βˆ’y)(2 - y)(2βˆ’y), so P(y)=ρg(2βˆ’y)P(y) = \rho g(2 - y)P(y)=ρg(2βˆ’y) factors out of the inner dxdxdx integral as a pure constant.
Show your steps and calculations here.

Practice Problem 6

Evaluate the double integral over the Type I region:

∬Dy dAwhereΒ D={(x,y)∈R2∣0≀x≀1, x2≀y≀x}\iint_D \sqrt{y} \, dA \quad \text{where } D = \{(x, y) \in \mathbb{R}^2 \mid 0 \le x \le 1, \, x^2 \le y \le x\}∬D​y​dAwhereΒ D={(x,y)∈R2∣0≀x≀1,x2≀y≀x}

Solution Key:
  1. Formulate the Type I Iterated Integral:
    State the Type I Double Integral Formula:

    ∬Df(x,y) dA=∫ab∫g1(x)g2(x)f(x,y) dy dx\iint_D f(x, y) \, dA = \int_{a}^{b} \int_{g_1(x)}^{g_2(x)} f(x, y) \, dy \, dx∬D​f(x,y)dA=∫abβ€‹βˆ«g1​(x)g2​(x)​f(x,y)dydx
    I=∫01∫x2xy1/2 dy dxI = \int_{0}^{1} \int_{x^2}^{x} y^{1/2} \, dy \, dxI=∫01β€‹βˆ«x2x​y1/2dydx

  2. Evaluate the Inner Integral (dydydy):
    Treat xxx as a constant:

    ∫x2xy1/2 dy=[23y3/2]y=x2y=x=23(x3/2βˆ’(x2)3/2)=23(x3/2βˆ’x3)\int_{x^2}^{x} y^{1/2} \, dy = \left[ \frac{2}{3} y^{3/2} \right]_{y=x^2}^{y=x} = \frac{2}{3} \left( x^{3/2} - (x^2)^{3/2} \right) = \frac{2}{3} \left( x^{3/2} - x^3 \right)∫x2x​y1/2dy=[32​y3/2]y=x2y=x​=32​(x3/2βˆ’(x2)3/2)=32​(x3/2βˆ’x3)

  3. Evaluate the Outer Integral (dxdxdx):

    I=23∫01(x3/2βˆ’x3) dx=23[25x5/2βˆ’14x4]01I = \frac{2}{3} \int_{0}^{1} (x^{3/2} - x^3) \, dx = \frac{2}{3} \left[ \frac{2}{5}x^{5/2} - \frac{1}{4}x^4 \right]_{0}^{1}I=32β€‹βˆ«01​(x3/2βˆ’x3)dx=32​[52​x5/2βˆ’41​x4]01​
    Evaluate at the limits:
    =23(25βˆ’14)=23(8βˆ’520)=23(320)=220=110= \frac{2}{3} \left( \frac{2}{5} - \frac{1}{4} \right) = \frac{2}{3} \left( \frac{8 - 5}{20} \right) = \frac{2}{3} \left( \frac{3}{20} \right) = \frac{2}{20} = \frac{1}{10}=32​(52β€‹βˆ’41​)=32​(208βˆ’5​)=32​(203​)=202​=101​

Show your steps and calculations here.

Practice Problem 7

Evaluate the double integral over the Type II region:

∬Dy dAwhereΒ D={(x,y)∈R2∣0≀y≀1, y≀x≀y}\iint_D \sqrt{y} \, dA \quad \text{where } D = \{(x, y) \in \mathbb{R}^2 \mid 0 \le y \le 1, \, y \le x \le \sqrt{y}\}∬D​y​dAwhereΒ D={(x,y)∈R2∣0≀y≀1,y≀x≀y​}

Solution Key:
  1. Formulate the Type II Iterated Integral:
    State the Type II Double Integral Formula:

    ∬Df(x,y) dA=∫cd∫h1(y)h2(y)f(x,y) dx dy\iint_D f(x, y) \, dA = \int_{c}^{d} \int_{h_1(y)}^{h_2(y)} f(x, y) \, dx \, dy∬D​f(x,y)dA=∫cdβ€‹βˆ«h1​(y)h2​(y)​f(x,y)dxdy
    I=∫01∫yyy dx dyI = \int_{0}^{1} \int_{y}^{\sqrt{y}} \sqrt{y} \, dx \, dyI=∫01β€‹βˆ«yy​​y​dxdy

  2. Evaluate the Inner Integral (dxdxdx):
    Since y\sqrt{y}y​ contains no xxx, it is treated as a constant:

    ∫yyy dx=y[x]yy=y(yβˆ’y)=yβˆ’y3/2\int_{y}^{\sqrt{y}} \sqrt{y} \, dx = \sqrt{y} [x]_{y}^{\sqrt{y}} = \sqrt{y} (\sqrt{y} - y) = y - y^{3/2}∫yy​​y​dx=y​[x]yy​​=y​(yβ€‹βˆ’y)=yβˆ’y3/2

  3. Evaluate the Outer Integral (dydydy):

    I=∫01(yβˆ’y3/2) dy=[12y2βˆ’25y5/2]01=12βˆ’25=5βˆ’410=110I = \int_{0}^{1} (y - y^{3/2}) \, dy = \left[ \frac{1}{2}y^2 - \frac{2}{5}y^{5/2} \right]_{0}^{1} = \frac{1}{2} - \frac{2}{5} = \frac{5 - 4}{10} = \frac{1}{10}I=∫01​(yβˆ’y3/2)dy=[21​y2βˆ’52​y5/2]01​=21β€‹βˆ’52​=105βˆ’4​=101​
    (Notice that Practice Problems 6 and 7 represent the exact same geometric region in the xyxyxy-plane, confirming that Type I and Type II evaluations yield identical numeric results).

Show your steps and calculations here.

Practice Problem 8

An optoelectronic display module features a wedge-shaped active zone RRR bounded by y=0y = 0y=0, y=xy = xy=x, and x=2x = 2x=2 (spatial units in meters). Emitted irradiance is given by:

E(x,y)=E0(1βˆ’0.2x)(W/m2)E(x, y) = E_0 (1 - 0.2x) \quad (\text{W/m}^2)E(x,y)=E0​(1βˆ’0.2x)(W/m2)

  1. Compute the total radiant power ∬RE dA\iint_R E \, dA∬R​EdA as a Type I double integral.
  2. Find the spatial average irradiance EavgE_{\text{avg}}Eavg​ across the wedge in terms of E0E_0E0​.
Solution Key:
  1. Evaluate the Total Irradiance Integral:
    State the Type I Double Integral Formula (defined in Double Integrals over a General Region):

    ∬RE(x,y) dA=∫02∫0xE0(1βˆ’0.2x) dy dx\iint_R E(x, y) \, dA = \int_{0}^{2} \int_{0}^{x} E_0 (1 - 0.2x) \, dy \, dx∬R​E(x,y)dA=∫02β€‹βˆ«0x​E0​(1βˆ’0.2x)dydx
    Since the integrand contains no yyy, it acts as a constant during inner integration:
    ∫0xE0(1βˆ’0.2x) dy=E0(1βˆ’0.2x)[y]0x=E0(xβˆ’0.2x2)\int_{0}^{x} E_0 (1 - 0.2x) \, dy = E_0 (1 - 0.2x) [y]_0^x = E_0 (x - 0.2x^2)∫0x​E0​(1βˆ’0.2x)dy=E0​(1βˆ’0.2x)[y]0x​=E0​(xβˆ’0.2x2)
    Now integrate with respect to xxx from 000 to 222:
    ∬RE dA=E0∫02(xβˆ’15x2) dx=E0[12x2βˆ’115x3]02\iint_R E \, dA = E_0 \int_{0}^{2} \left( x - \frac{1}{5}x^2 \right) \, dx = E_0 \left[ \frac{1}{2}x^2 - \frac{1}{15}x^3 \right]_{0}^{2}∬R​EdA=E0β€‹βˆ«02​(xβˆ’51​x2)dx=E0​[21​x2βˆ’151​x3]02​
    Evaluate at x=2x = 2x=2:
    =E0(42βˆ’815)=E0(2βˆ’815)=E0(30βˆ’815)=2215E0= E_0 \left( \frac{4}{2} - \frac{8}{15} \right) = E_0 \left( 2 - \frac{8}{15} \right) = E_0 \left( \frac{30 - 8}{15} \right) = \frac{22}{15}E_0=E0​(24β€‹βˆ’158​)=E0​(2βˆ’158​)=E0​(1530βˆ’8​)=1522​E0​

  2. Compute Spatial Average Irradiance (EavgE_{\text{avg}}Eavg​):
    State the Spatial Average Value Formula:

    Eavg=1Area(R)∬RE dAE_{\text{avg}} = \frac{1}{\text{Area}(R)} \iint_R E \, dAEavg​=Area(R)1β€‹βˆ¬R​EdA
    The region RRR is a right triangle with base 2Β m2\text{ m}2Β m and height 2Β m2\text{ m}2Β m:
    Area(R)=12(2)(2)=2Β m2\text{Area}(R) = \frac{1}{2}(2)(2) = 2\text{ m}^2Area(R)=21​(2)(2)=2Β m2
    Eavg=12(2215E0)=1115E0β‰ˆ0.733E0Β W/m2E_{\text{avg}} = \frac{1}{2} \left( \frac{22}{15}E_0 \right) = \frac{11}{15}E_0 \approx 0.733E_0\text{ W/m}^2Eavg​=21​(1522​E0​)=1511​E0β€‹β‰ˆ0.733E0​ W/m2


Show your steps and calculations here.

Regroup 2

  • Common Slicing Trap: Why does setting limits like ∫x22x∫02f dx dy\int_{x^2}^{2x} \int_0^2 f\,dx\,dy∫x22xβ€‹βˆ«02​fdxdy cause an automatic failure in engineering software? (Because variables appear in the outer limits, leaving a non-constant function instead of a physical scalar number).
  • Choosing Type I vs. Type II: When is Type II strictly required over Type I? (When solving boundary equations for yyy introduces piecewise square roots or vertical step discontinuities, whereas solving for xxx yields single continuous functions).

Skill Block 3

Worked Example A: Dual-Type Geometric Slicing on a Triangular Domain

Let RRR be the triangular region with vertices (0,0),(3,0),(0, 0), (3, 0),(0,0),(3,0), and (3,2)(3, 2)(3,2). Evaluate ∬R(2x+3y) dA\iint_R (2x + 3y) \, dA∬R​(2x+3y)dA by setting it up first as a Type I integral, then verifying the result as a Type II integral.

Worked Solution:
  1. Identify the Boundary Lines:
    The triangle has three boundary segments:

    • Bottom edge: y=0y = 0y=0 along the xxx-axis from x=0x = 0x=0 to x=3x = 3x=3.
    • Right edge: x=3x = 3x=3 from y=0y = 0y=0 to y=2y = 2y=2.
    • Hypotenuse edge connecting (0,0)(0, 0)(0,0) to (3,2)(3, 2)(3,2): slope m=2βˆ’03βˆ’0=23β€…β€ŠβŸΉβ€…β€Šy=23xβ€…β€ŠβŸΊβ€…β€Šx=32ym = \frac{2 - 0}{3 - 0} = \frac{2}{3} \implies y = \frac{2}{3}x \iff x = \frac{3}{2}ym=3βˆ’02βˆ’0​=32β€‹βŸΉy=32​x⟺x=23​y.
  2. Method 1: Type I Setup (Vertical Slices):
    State the Type I Double Integral Formula (defined in Double Integrals over a General Region):

    x∈[0,3],0≀y≀23xx \in [0, 3], \quad 0 \le y \le \frac{2}{3}xx∈[0,3],0≀y≀32​x
    I=∫03∫023x(2x+3y) dy dxI = \int_{0}^{3} \int_{0}^{\frac{2}{3}x} (2x + 3y) \, dy \, dxI=∫03β€‹βˆ«032​x​(2x+3y)dydx
    Evaluate inner integral (dydydy):
    ∫023x(2x+3y) dy=[2xy+32y2]023x=2x(23x)+32(23x)2=43x2+32(49x2)=43x2+23x2=2x2\int_{0}^{\frac{2}{3}x} (2x + 3y) \, dy = \left[ 2xy + \frac{3}{2}y^2 \right]_{0}^{\frac{2}{3}x} = 2x\left(\frac{2}{3}x\right) + \frac{3}{2}\left(\frac{2}{3}x\right)^2 = \frac{4}{3}x^2 + \frac{3}{2}\left(\frac{4}{9}x^2\right) = \frac{4}{3}x^2 + \frac{2}{3}x^2 = 2x^2∫032​x​(2x+3y)dy=[2xy+23​y2]032​x​=2x(32​x)+23​(32​x)2=34​x2+23​(94​x2)=34​x2+32​x2=2x2
    Evaluate outer integral (dxdxdx):
    I=∫032x2 dx=[23x3]03=23(27)=18I = \int_{0}^{3} 2x^2 \, dx = \left[ \frac{2}{3}x^3 \right]_{0}^{3} = \frac{2}{3}(27) = 18I=∫03​2x2dx=[32​x3]03​=32​(27)=18

  3. Method 2: Type II Setup (Horizontal Slices):
    State the Type II Double Integral Formula (defined in Double Integrals over a General Region):

    • Outer bounds for yyy: y∈[0,2]y \in [0, 2]y∈[0,2]
    • Horizontal slices enter on the left at the line x=32yx = \frac{3}{2}yx=23​y and exit on the right at the vertical line x=3x = 3x=3:
      I=∫02∫32y3(2x+3y) dx dyI = \int_{0}^{2} \int_{\frac{3}{2}y}^{3} (2x + 3y) \, dx \, dyI=∫02β€‹βˆ«23​y3​(2x+3y)dxdy
      Evaluate inner integral (dxdxdx):
      ∫32y3(2x+3y) dx=[x2+3yx]x=32yx=3=(9+9y)βˆ’((32y)2+3y(32y))\int_{\frac{3}{2}y}^{3} (2x + 3y) \, dx = \left[ x^2 + 3yx \right]_{x=\frac{3}{2}y}^{x=3} = \Big(9 + 9y\Big) - \left( \left(\frac{3}{2}y\right)^2 + 3y\left(\frac{3}{2}y\right) \right)∫23​y3​(2x+3y)dx=[x2+3yx]x=23​yx=3​=(9+9y)βˆ’((23​y)2+3y(23​y))
      =9+9yβˆ’(94y2+92y2)=9+9yβˆ’274y2= 9 + 9y - \left( \frac{9}{4}y^2 + \frac{9}{2}y^2 \right) = 9 + 9y - \frac{27}{4}y^2=9+9yβˆ’(49​y2+29​y2)=9+9yβˆ’427​y2
      Evaluate outer integral (dydydy):
      I=∫02(9+9yβˆ’274y2) dy=[9y+92y2βˆ’94y3]02I = \int_{0}^{2} \left( 9 + 9y - \frac{27}{4}y^2 \right) \, dy = \left[ 9y + \frac{9}{2}y^2 - \frac{9}{4}y^3 \right]_{0}^{2}I=∫02​(9+9yβˆ’427​y2)dy=[9y+29​y2βˆ’49​y3]02​
      =9(2)+92(4)βˆ’94(8)=18+18βˆ’18=18= 9(2) + \frac{9}{2}(4) - \frac{9}{4}(8) = 18 + 18 - 18 = 18=9(2)+29​(4)βˆ’49​(8)=18+18βˆ’18=18
      Both methods produce the exact same value of 181818.

Worked Example B: Overcoming a Non-Elementary Antiderivative

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

Worked Solution:
  1. Execute Step 1 (Extract Original Bounds):
    From the integral limits:

    TypeΒ II:0≀y≀1,y≀x≀1\text{Type II:} \quad 0 \le y \le 1, \quad y \le x \le 1TypeΒ II:0≀y≀1,y≀x≀1
    Notice that ∫ex2 dx\int e^{x^2} \, dx∫ex2dx has no elementary antiderivative.

  2. Execute Step 2 & 3 (Sketch Domain & Invert Curves):
    The boundary equations are y=0y = 0y=0, y=1y = 1y=1, x=1x = 1x=1, and x=yx = yx=y.
    Inverting x=yx = yx=y gives y=xy = xy=x. The region is a triangle with vertices at (0,0),(1,0),(0, 0), (1, 0),(0,0),(1,0), and (1,1)(1, 1)(1,1).

  3. Execute Step 4 (Re-slice as Type I):
    Apply the 4-Step Re-slicing Protocol (defined in Switching the Order of Integration):

    • Outer constant limits for xxx: 0≀x≀10 \le x \le 10≀x≀1
    • Vertical slices enter at the bottom line y=0y = 0y=0 and exit at the slanted line y=xy = xy=x:
      0≀y≀x0 \le y \le x0≀y≀x
      Formulate the new integral:
      I=∫01∫0xex2 dy dxI = \int_{0}^{1} \int_{0}^{x} e^{x^2} \, dy \, dxI=∫01β€‹βˆ«0x​ex2dydx
  4. Evaluate Switched Integral:
    The inner integrand ex2e^{x^2}ex2 contains no yyy, so it is treated 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
    Now evaluate the outer integral using 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​:
    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)

Practice Problem 9: Concept & Reading Check

Open this Skill Block on the Memex to locate the information needed to answer the following:

  • (a) Why is it strictly invalid to switch the order of integration over a general region by merely swapping the functional integral signs (e.g., writing ∫g1(x)g2(x)∫abf dx dy\int_{g_1(x)}^{g_2(x)} \int_a^b f\,dx\,dy∫g1​(x)g2​(x)β€‹βˆ«ab​fdxdy)?
  • (b) In Worked Example B, describe how sketching the triangular region 0≀y≀1,y≀x≀10 \le y \le 1, y \le x \le 10≀y≀1,y≀x≀1 allows you to reformulate the limits as 0≀x≀1,0≀y≀x0 \le x \le 1, 0 \le y \le x0≀x≀1,0≀y≀x.
  • (c) In Aircraft Wing Torsional Polar Moment of Inertia, why does re-slicing a tapered airfoil cross-section from vertical chords to horizontal spanwise strips simplify the calculation?
Solution Key:
  • (a) State the Golden Rule of Iterated Integrals (defined in Switching the Order of Integration):
    ∬Df dA∈R\iint_D f \, dA \in \mathbb{R}∬D​fdA∈R
    The outer integral must always have constant limits [a,b][a, b][a,b] or [c,d][c, d][c,d]. Merely moving the functional limits to the outer integral leaves variables in the final evaluated result, violating the fundamental theorem of calculus.
  • (b) The original limits y≀x≀1y \le x \le 1y≀x≀1 describe horizontal arrows entering at x=yx = yx=y and terminating at the vertical line x=1x = 1x=1, starting from y=0y = 0y=0 up to y=1y = 1y=1. When sketched, the domain is the lower right triangle bounded by y=0,x=1,y = 0, x = 1,y=0,x=1, and y=xy = xy=x. Slicing vertically, xxx ranges continuously from 000 to 111, and for each xxx, the vertical ray runs from the baseline y=0y = 0y=0 up to the line y=xy = xy=x.
  • (c) In aerospace structural analysis, the boundary curve is y=xy = \sqrt{x}y=x​. When integrating in Type I chords, intermediate steps involve fractional powers ∫x2x dx=∫x5/2 dx\int x^2 \sqrt{x}\,dx = \int x^{5/2}\,dx∫x2x​dx=∫x5/2dx. Re-slicing horizontally as Type II inverts the boundary to the clean polynomial x=y2x = y^2x=y2, so the inner integral integrates pure integer powers ∫y24x2 dx=[x33]y24=64βˆ’y63\int_{y^2}^4 x^2\,dx = [\frac{x^3}{3}]_{y^2}^4 = \frac{64 - y^6}{3}∫y24​x2dx=[3x3​]y24​=364βˆ’y6​, eliminating radical fractions.
Show your steps and calculations here.

Practice Problem 10

Let RRR be bounded by y=x3y = x^3y=x3, y=4y = 4y=4, and the vertical lines x=βˆ’2x = -2x=βˆ’2 and x=2x = 2x=2.

  1. Convert the bounds to a Type II region, splitting the domain into two sub-regions R1R_1R1​ and R2R_2R2​ as necessary.
  2. Evaluate the total geometric area ∬R1 dA\iint_R 1 \, dA∬R​1dA.
Solution Key:
  1. Convert to Type II by Inverting Boundaries & Splitting Domain:
    State the Boundary Inversion Protocol (defined in Switching the Order of Integration):

    y=x3β€…β€ŠβŸΉβ€…β€Šx=y1/3y = x^3 \implies x = y^{1/3}y=x3⟹x=y1/3
    The domain spans vertically from y=(βˆ’2)3=βˆ’8y = (-2)^3 = -8y=(βˆ’2)3=βˆ’8 to y=(2)3=8y = (2)^3 = 8y=(2)3=8. At y=4y = 4y=4, the curve crosses x=41/3β‰ˆ1.587x = 4^{1/3} \approx 1.587x=41/3β‰ˆ1.587:

    • Lower Sub-region R1R_1R1​ (from y=βˆ’8y = -8y=βˆ’8 to y=4y = 4y=4):
      Left boundary is the vertical line x=βˆ’2x = -2x=βˆ’2; right boundary is the curve x=y1/3x = y^{1/3}x=y1/3.
      R1={(x,y)∈R2βˆ£βˆ’8≀y≀4,β€‰βˆ’2≀x≀y1/3}R_1 = \{(x, y) \in \mathbb{R}^2 \mid -8 \le y \le 4, \, -2 \le x \le y^{1/3}\}R1​={(x,y)∈R2βˆ£βˆ’8≀y≀4,βˆ’2≀x≀y1/3}
    • Upper Sub-region R2R_2R2​ (from y=4y = 4y=4 to y=8y = 8y=8):
      Left boundary is the curve x=y1/3x = y^{1/3}x=y1/3; right boundary is the vertical line x=2x = 2x=2.
      R2={(x,y)∈R2∣4≀y≀8, y1/3≀x≀2}R_2 = \{(x, y) \in \mathbb{R}^2 \mid 4 \le y \le 8, \, y^{1/3} \le x \le 2\}R2​={(x,y)∈R2∣4≀y≀8,y1/3≀x≀2}
  2. Evaluate the Area Integral as the Sum of Two Type II Integrals:

    Area(R)=βˆ«βˆ’84βˆ«βˆ’2y1/31 dx dy+∫48∫y1/321 dx dy\text{Area}(R) = \int_{-8}^{4} \int_{-2}^{y^{1/3}} 1 \, dx \, dy + \int_{4}^{8} \int_{y^{1/3}}^{2} 1 \, dx \, dyArea(R)=βˆ«βˆ’84β€‹βˆ«βˆ’2y1/3​1dxdy+∫48β€‹βˆ«y1/32​1dxdy

    Evaluate R1R_1R1​ Integral:

    βˆ«βˆ’84[x]βˆ’2y1/3 dy=βˆ«βˆ’84(y1/3+2) dy=[34y4/3+2y]βˆ’84\int_{-8}^{4} [x]_{-2}^{y^{1/3}} \, dy = \int_{-8}^{4} (y^{1/3} + 2) \, dy = \left[ \frac{3}{4}y^{4/3} + 2y \right]_{-8}^{4}βˆ«βˆ’84​[x]βˆ’2y1/3​dy=βˆ«βˆ’84​(y1/3+2)dy=[43​y4/3+2y]βˆ’84​
    =(34(4)4/3+2(4))βˆ’(34(βˆ’8)4/3+2(βˆ’8))= \left( \frac{3}{4}(4)^{4/3} + 2(4) \right) - \left( \frac{3}{4}(-8)^{4/3} + 2(-8) \right)=(43​(4)4/3+2(4))βˆ’(43​(βˆ’8)4/3+2(βˆ’8))
    Since 44/3=4Γ—41/34^{4/3} = 4 \times 4^{1/3}44/3=4Γ—41/3 and (βˆ’8)4/3=((βˆ’8)1/3)4=(βˆ’2)4=16(-8)^{4/3} = ((-8)^{1/3})^4 = (-2)^4 = 16(βˆ’8)4/3=((βˆ’8)1/3)4=(βˆ’2)4=16:
    =(3(41/3)+8)βˆ’(34(16)βˆ’16)=3(41/3)+8βˆ’(12βˆ’16)=3(41/3)+12= \left( 3(4^{1/3}) + 8 \right) - \left( \frac{3}{4}(16) - 16 \right) = 3(4^{1/3}) + 8 - (12 - 16) = 3(4^{1/3}) + 12=(3(41/3)+8)βˆ’(43​(16)βˆ’16)=3(41/3)+8βˆ’(12βˆ’16)=3(41/3)+12

    Evaluate R2R_2R2​ Integral:

    ∫48[x]y1/32 dy=∫48(2βˆ’y1/3) dy=[2yβˆ’34y4/3]48\int_{4}^{8} [x]_{y^{1/3}}^{2} \, dy = \int_{4}^{8} (2 - y^{1/3}) \, dy = \left[ 2y - \frac{3}{4}y^{4/3} \right]_{4}^{8}∫48​[x]y1/32​dy=∫48​(2βˆ’y1/3)dy=[2yβˆ’43​y4/3]48​
    =(2(8)βˆ’34(8)4/3)βˆ’(2(4)βˆ’34(4)4/3)= \left( 2(8) - \frac{3}{4}(8)^{4/3} \right) - \left( 2(4) - \frac{3}{4}(4)^{4/3} \right)=(2(8)βˆ’43​(8)4/3)βˆ’(2(4)βˆ’43​(4)4/3)
    Since 84/3=(81/3)4=24=168^{4/3} = (8^{1/3})^4 = 2^4 = 1684/3=(81/3)4=24=16:
    =(16βˆ’34(16))βˆ’(8βˆ’3(41/3))=(16βˆ’12)βˆ’8+3(41/3)=3(41/3)βˆ’4= \left( 16 - \frac{3}{4}(16) \right) - \left( 8 - 3(4^{1/3}) \right) = (16 - 12) - 8 + 3(4^{1/3}) = 3(4^{1/3}) - 4=(16βˆ’43​(16))βˆ’(8βˆ’3(41/3))=(16βˆ’12)βˆ’8+3(41/3)=3(41/3)βˆ’4

    Sum the Components:

    TotalΒ Area=(3(41/3)+12)+(3(41/3)βˆ’4)=6(41/3)+8=643+8β‰ˆ17.524\text{Total Area} = \Big(3(4^{1/3}) + 12\Big) + \Big(3(4^{1/3}) - 4\Big) = 6(4^{1/3}) + 8 = 6\sqrt[3]{4} + 8 \approx 17.524TotalΒ Area=(3(41/3)+12)+(3(41/3)βˆ’4)=6(41/3)+8=634​+8β‰ˆ17.524

Show your steps and calculations here.

Practice Problem 11

Evaluate the double integral by switching the order of integration:

∫0Ο€βˆ«xΟ€sin⁑yy dy dx\int_{0}^{\pi} \int_{x}^{\pi} \frac{\sin y}{y} \, dy \, dx∫0Ο€β€‹βˆ«xπ​ysiny​dydx

Solution Key:
  1. Extract Original Boundaries & Recognize Non-Elementary Integrand:
    The original Type I bounds are:

    0≀x≀π,x≀y≀π0 \le x \le \pi, \quad x \le y \le \pi0≀x≀π,x≀y≀π
    Because ∫sin⁑yy dy\int \frac{\sin y}{y} \, dy∫ysiny​dy cannot be solved in elementary terms, we must switch to Type II.

  2. Re-slice as a Type II General Region:
    Apply the 4-Step Re-slicing Protocol (defined in Switching the Order of Integration):

    • The triangular domain has vertices at (0,0),(0,Ο€),(0, 0), (0, \pi),(0,0),(0,Ο€), and (Ο€,Ο€)(\pi, \pi)(Ο€,Ο€).
    • The vertical extent is 0≀y≀π0 \le y \le \pi0≀y≀π.
    • A horizontal slice enters at the left vertical line x=0x = 0x=0 and exits at the line x=yx = yx=y:
      0≀x≀y0 \le x \le y0≀x≀y
      Formulate 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
  3. Evaluate the Switched Iterated Integral:
    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):
    I=∫0Ο€sin⁑y dy=[βˆ’cos⁑y]0Ο€=βˆ’cos⁑(Ο€)βˆ’(βˆ’cos⁑(0))=βˆ’(βˆ’1)βˆ’(βˆ’1)=1+1=2I = \int_{0}^{\pi} \sin y \, dy = \left[ -\cos y \right]_{0}^{\pi} = -\cos(\pi) - (-\cos(0)) = -(-1) - (-1) = 1 + 1 = 2I=∫0π​sinydy=[βˆ’cosy]0π​=βˆ’cos(Ο€)βˆ’(βˆ’cos(0))=βˆ’(βˆ’1)βˆ’(βˆ’1)=1+1=2


Show your steps and calculations here.

Regroup 3

  • Order Reversal as an Unlocking Key: Notice how in Problems 9 and 11, functions that appeared mathematically "unintegrable" (ex2e^{x^2}ex2 and sin⁑yy\frac{\sin y}{y}ysiny​) instantly dissolved into high-school uuu-substitutions the moment the slicing direction was switched.
  • Region Decomposition Diagnostic: How do you know when a region must be split into multiple integrals during order switching? (Whenever a horizontal or vertical scanning line transitions between different bounding formulas, such as switching from a curve to a flat boundary wall).

Synthesis Wrap-up

  • Core Takeaways from Session 14:
    1. Double integrals extend accumulation into 2D space, calculating volumes, areas, total physical charges, and spatial averages.
    2. Over rectangular regions, Fubini's Theorem permits free order swapping; choose the order that creates chain-rule derivatives and avoids integration by parts.
    3. Over general regions, outer limits must always be constants; inner limits must be functions. Slicing direction (Type I vertical vs. Type II horizontal) dictates boundary formulation.
    4. Switching the order of integration is not just a mathematical curiosityβ€”it is a vital tool for overcoming non-elementary antiderivatives and simplifying multi-component engineering boundary models.
  • Knewton Alta Homework Alignment:
    • Double Integrals over Rectangular Regions Part 1: Basic iterated integration, Fubini's Theorem, and product separability.
    • Double Integrals over Rectangular Regions Part 2: Spatial average values and applications to physical mass/charge distributions.
    • Double Integrals over General Regions: Setting up and evaluating Type I and Type II regions, and executing the 4-step order-switching protocol.