MATH-283 (Calculus III) Curriculum

Session 16 - Triple Integrals

Handout Worksheet
Fall Semester 2026
Session Objectives & Overview

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

  • Formulate and evaluate triple iterated integrals over 3D rectangular boxes B=[a,b]×[c,d]×[r,s]B = [a, b] \times [c, d] \times [r, s]B=[a,b]×[c,d]×[r,s] in Cartesian coordinates, applying Fubini's Theorem across all six differential orderings and leveraging algebraic product separability ∭Bg(x)h(y)k(z) dV=(∫g)(∫h)(∫k)\iiint_B g(x)h(y)k(z) \, dV = \left(\int g\right)\left(\int h\right)\left(\int k\right)∭B​g(x)h(y)k(z)dV=(∫g)(∫h)(∫k) to compute enclosed volumes, total mass, and 3D spatial average values.
  • Construct rigorous integration limits for general non-rectangular 3D regions E⊂R3E \subset \mathbb{R}^3E⊂R3 by executing the 4-step spatial projection protocol: identifying upper and lower bounding surfaces along the inner axis (z=u1(x,y)z = u_1(x, y)z=u1​(x,y) to z=u2(x,y)z = u_2(x, y)z=u2​(x,y)), projecting the 3D boundary onto the 2D coordinate shadow plane DxyD_{xy}Dxy​, establishing boundary curves, and bounding the outer variable with constants.
  • Strategically reorder the differentials of a triple integral (dz dy dx→dx dz dydz \, dy \, dx \to dx \, dz \, dydzdydx→dxdzdy) to bypass non-elementary antiderivatives (such as ∫ex2 dx\int e^{x^2} \, dx∫ex2dx), inverting surface functions and analyzing trace intersection curves to determine whether the domain remains unified or partitions into multiple piecewise iterated integrals.

Micro-Lecture

1. Spatial Volume Accumulation & The 3D Riemann Sum

  • Mathematical Foundation: In single-variable calculus, definite integrals accumulate quantities along a 1D line segment; in multivariable calculus, double integrals accumulate over a 2D planar patch. As established in Multiple Integrals, when physical properties vary continuously across 3D space, we partition a bounded volume E⊂R3E \subset \mathbb{R}^3E⊂R3 into an array of infinitesimal rectangular cuboids with side dimensions Δx,Δy,Δz\Delta x, \Delta y, \Delta zΔx,Δy,Δz. The volume of each cuboid is ΔV=Δx Δy Δz\Delta V = \Delta x \, \Delta y \, \Delta zΔV=ΔxΔyΔz. The triple integral is defined as the infinite refinement limit of the triple Riemann sum:
    ∭Ef(x,y,z) dV=lim⁡l→∞m→∞n→∞∑i=1l∑j=1m∑k=1nf(xijk∗,yijk∗,zijk∗) Δxi Δyj Δzk\iiint_{E} f(x, y, z) \, dV = \lim_{\substack{l \to \infty \\ m \to \infty \\ n \to \infty}} \sum_{i=1}^{l} \sum_{j=1}^{m} \sum_{k=1}^{n} f(x_{ijk}^*, y_{ijk}^*, z_{ijk}^*) \, \Delta x_i \, \Delta y_j \, \Delta z_k∭E​f(x,y,z)dV=l→∞m→∞n→∞​lim​i=1∑l​j=1∑m​k=1∑n​f(xijk∗​,yijk∗​,zijk∗​)Δxi​Δyj​Δzk​
    If f(x,y,z)=1f(x, y, z) = 1f(x,y,z)=1, the integral yields the spatial volume V(E)=∭E1 dVV(E) = \iiint_E 1 \, dVV(E)=∭E​1dV. If f(x,y,z)=ρ(x,y,z)f(x, y, z) = \rho(x, y, z)f(x,y,z)=ρ(x,y,z) represents a mass density (in kg/m3\text{kg/m}^3kg/m3) or charge density (in C/m3\text{C/m}^3C/m3), the integral calculates the total accumulated mass MMM or total enclosed electric charge QQQ.
  • Engineering Context: Mechanical and aerospace engineers utilize triple integrals to compute mass, centers of gravity, and mass moments of inertia for irregular structural components, ensuring aerodynamic stability and preventing rotational flutter in turbine blades.

2. Fubini's Theorem in Three Dimensions

  • Mathematical Foundation: When the spatial domain is a rectangular box B=[a,b]×[c,d]×[r,s]B = [a, b] \times [c, d] \times [r, s]B=[a,b]×[c,d]×[r,s], the limits of integration for all three coordinate variables are fixed, independent constants. As established in Triple Integrals over a Rectangular Region, Fubini's Theorem guarantees that if f(x,y,z)f(x, y, z)f(x,y,z) is continuous across BBB, the triple integral evaluates as an iterated integral in any of the 3!=63! = 63!=6 possible differential permutations (dz dy dxdz \, dy \, dxdzdydx, dz dx dydz \, dx \, dydzdxdy, dy dz dxdy \, dz \, dxdydzdx, etc.). Furthermore, if the integrand factors into single-variable components f(x,y,z)=g(x)h(y)k(z)f(x, y, z) = g(x)h(y)k(z)f(x,y,z)=g(x)h(y)k(z), the triple integral decouples into the product of three independent single integrals:
    ∭Bg(x)h(y)k(z) dV=(∫abg(x) dx)(∫cdh(y) dy)(∫rsk(z) dz)\iiint_{B} g(x)h(y)k(z) \, dV = \left(\int_{a}^{b} g(x) \, dx\right) \left(\int_{c}^{d} h(y) \, dy\right) \left(\int_{r}^{s} k(z) \, dz\right)∭B​g(x)h(y)k(z)dV=(∫ab​g(x)dx)(∫cd​h(y)dy)(∫rs​k(z)dz)
  • Engineering Context: Electrical and computer engineers evaluate total thermal dissipation and electromagnetic flux in cubic microchip packaging and semiconductor silicon substrates where material properties vary along orthogonal fabrication axes.

3. General Spatial Regions & The Shadow Projection Method

  • Mathematical Foundation: For general solids with curved faces or inclined planar walls, boundaries depend functionally on neighboring variables. As detailed in Triple Integrals over a General Region, setting up a triple integral requires a systematic "inside-out" perspective: identify the entering floor surface z=u1(x,y)z = u_1(x, y)z=u1​(x,y) and exiting ceiling surface z=u2(x,y)z = u_2(x, y)z=u2​(x,y) along the inner axis, then project the boundary intersection onto the coordinate shadow plane DxyD_{xy}Dxy​. The 2D shadow DxyD_{xy}Dxy​ is then treated as a standard double integral region with curve-to-curve middle limits and constant outer limits. When an integrand contains non-elementary single-variable functions, reorienting the spatial sweep order (as documented in Changing the Order of Integration in Triple Integrals) generates algebraic multipliers that make closed-form evaluation possible.
  • Engineering Context: Civil and environmental engineers use general triple integrals to calculate bulk earthwork excavation volumes and hydraulic reservoir storage capacity, while chemical engineers model pollutant destruction yields across packed catalyst beds.

Skill Block 1

Worked Example A: Evaluating a Monomial Triple Integral over a Rectangular Box

Evaluate the triple integral:

∭B12xyz2 dV\iiint_{B} 12xyz^2 \, dV∭B​12xyz2dV
where BBB is the rectangular box defined by 0≤x≤10 \le x \le 10≤x≤1, 1≤y≤21 \le y \le 21≤y≤2, and 0≤z≤20 \le z \le 20≤z≤2.

Worked Solution:
  1. Set Up the Iterated Integral:
    State Fubini's Theorem for Triple Integrals (defined in Triple Integrals over a Rectangular Region):

    ∭Bf(x,y,z) dV=∫x=ab∫y=cd∫z=rsf(x,y,z) dz dy dx\iiint_{B} f(x, y, z) \, dV = \int_{x=a}^{b} \int_{y=c}^{d} \int_{z=r}^{s} f(x, y, z) \, dz \, dy \, dx∭B​f(x,y,z)dV=∫x=ab​∫y=cd​∫z=rs​f(x,y,z)dzdydx
    Substitute the given boundaries:
    I=∫01∫12∫0212xyz2 dz dy dxI = \int_{0}^{1} \int_{1}^{2} \int_{0}^{2} 12xyz^2 \, dz \, dy \, dxI=∫01​∫12​∫02​12xyz2dzdydx

  2. Evaluate the Innermost Integral (dzdzdz):
    Treat xxx and yyy as constants during partial integration:

    ∫0212xyz2 dz=[12xyz33]02=[4xyz3]02=4xy(2)3−0=32xy\int_{0}^{2} 12xyz^2 \, dz = \left[ 12xy \frac{z^3}{3} \right]_{0}^{2} = \left[ 4xyz^3 \right]_{0}^{2} = 4xy(2)^3 - 0 = 32xy∫02​12xyz2dz=[12xy3z3​]02​=[4xyz3]02​=4xy(2)3−0=32xy

  3. Evaluate the Middle Integral (dydydy):
    Substitute the inner result into the yyy-layer, holding xxx constant:

    ∫1232xy dy=[16xy2]12=16x(2)2−16x(1)2=64x−16x=48x\int_{1}^{2} 32xy \, dy = \left[ 16xy^2 \right]_{1}^{2} = 16x(2)^2 - 16x(1)^2 = 64x - 16x = 48x∫12​32xydy=[16xy2]12​=16x(2)2−16x(1)2=64x−16x=48x

  4. Evaluate the Outermost Integral (dxdxdx):
    Substitute the middle result into the xxx-layer:

    ∫0148x dx=[24x2]01=24(1)2−0=24\int_{0}^{1} 48x \, dx = \left[ 24x^2 \right]_{0}^{1} = 24(1)^2 - 0 = 24∫01​48xdx=[24x2]01​=24(1)2−0=24
    Thus, ∭B12xyz2 dV=24\iiint_B 12xyz^2 \, dV = 24∭B​12xyz2dV=24.

Worked Example B: Average Value of a 3D Linear Potential Field

Find the average value f‾\overline{f}f​ of the scalar field f(x,y,z)=x+2y−zf(x, y, z) = x + 2y - zf(x,y,z)=x+2y−z over the rectangular box B=[0,2]×[0,3]×[1,2]B = [0, 2] \times [0, 3] \times [1, 2]B=[0,2]×[0,3]×[1,2].

Worked Solution:
  1. Calculate Geometric Box Volume:
    State the Volume of a Rectangular Box Formula:

    V(B)=(b−a)(d−c)(s−r)=(2−0)(3−0)(2−1)=2×3×1=6V(B) = (b - a)(d - c)(s - r) = (2 - 0)(3 - 0)(2 - 1) = 2 \times 3 \times 1 = 6V(B)=(b−a)(d−c)(s−r)=(2−0)(3−0)(2−1)=2×3×1=6

  2. Evaluate the Triple Integral over the Box:

    I=∫02∫03∫12(x+2y−z) dz dy dxI = \int_{0}^{2} \int_{0}^{3} \int_{1}^{2} (x + 2y - z) \, dz \, dy \, dxI=∫02​∫03​∫12​(x+2y−z)dzdydx

    Inner zzz-integral:

    ∫12(x+2y−z) dz=[(x+2y)z−z22]12=(2(x+2y)−2)−(1(x+2y)−12)=x+2y−32\int_{1}^{2} (x + 2y - z) \, dz = \left[ (x + 2y)z - \frac{z^2}{2} \right]_{1}^{2} = \left(2(x + 2y) - 2\right) - \left(1(x + 2y) - \frac{1}{2}\right) = x + 2y - \frac{3}{2}∫12​(x+2y−z)dz=[(x+2y)z−2z2​]12​=(2(x+2y)−2)−(1(x+2y)−21​)=x+2y−23​

    Middle yyy-integral:

    ∫03(x+2y−32) dy=[xy+y2−32y]03=3x+9−92=3x+92\int_{0}^{3} \left(x + 2y - \frac{3}{2}\right) \, dy = \left[ xy + y^2 - \frac{3}{2}y \right]_{0}^{3} = 3x + 9 - \frac{9}{2} = 3x + \frac{9}{2}∫03​(x+2y−23​)dy=[xy+y2−23​y]03​=3x+9−29​=3x+29​

    Outer xxx-integral:

    ∫02(3x+92) dx=[32x2+92x]02=32(4)+92(2)=6+9=15\int_{0}^{2} \left(3x + \frac{9}{2}\right) \, dx = \left[ \frac{3}{2}x^2 + \frac{9}{2}x \right]_{0}^{2} = \frac{3}{2}(4) + \frac{9}{2}(2) = 6 + 9 = 15∫02​(3x+29​)dx=[23​x2+29​x]02​=23​(4)+29​(2)=6+9=15

  3. Compute the Spatial Average Value:
    State the Average Value Formula for a 3D Function:

    f‾=1V(B)∭Bf(x,y,z) dV=156=2.5\overline{f} = \frac{1}{V(B)} \iiint_{B} f(x, y, z) \, dV = \frac{15}{6} = 2.5f​=V(B)1​∭B​f(x,y,z)dV=615​=2.5

Practice Problem 1: Concept & Reading Check

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

  • (a) State the formal relationship between the differential volume element dVdVdV and Cartesian increments dx,dy,dzdx, dy, dzdx,dy,dz, and explain why Fubini's Theorem permits all six integration orders to yield identical values on a rectangular box.
  • (b) In Worked Example A, explain why the variables xxx and yyy can be treated as constants during the innermost integration with respect to zzz.
  • (c) In the linked engineering application Microprocessor Heat Sink Volumetric Thermal Dissipation (Computer & Mechanical Engineering), how does integrating a 3D volumetric thermal generation density q(x,y,z)q(x, y, z)q(x,y,z) differ conceptually from finding the geometric volume of the heat-spreader block?
Solution Key:
  • (a) The differential volume element dV=dz dy dxdV = dz \, dy \, dxdV=dzdydx represents the infinitesimal volume of a 3D coordinate prism. Under Fubini's Theorem for continuous functions on rectangular boxes, the boundaries of each variable are independent constants, meaning the order of accumulation does not affect the continuous summation limit of the Riemann sum across the domain.
  • (b) When performing partial integration with respect to zzz, the coordinates xxx and yyy represent independent orthogonal directions held stationary during the vertical line-slicing pass. Consequently, any expression containing only xxx and yyy acts as a constant scaling multiplier with respect to zzz.
  • (c) Evaluating ∭B1 dV\iiint_B 1 \, dV∭B​1dV sums unit volume increments, yielding the spatial capacity V(B)V(B)V(B) in cubic centimeters (cm3\text{cm}^3cm3). In contrast, integrating the non-uniform power density ∭Bq(x,y,z) dV\iiint_B q(x, y, z) \, dV∭B​q(x,y,z)dV weights each spatial differential element dVdVdV by its local thermal dissipation rate in W/cm3\text{W/cm}^3W/cm3, accumulating the total thermal heat dissipation rate PtotalP_{\text{total}}Ptotal​ in watts (W\text{W}W).
Show your steps and calculations here.

Practice Problem 2

Evaluate the triple integral over the rectangular box:

I=∫−15∫24∫01(x+yz2) dz dy dxI = \int_{-1}^{5} \int_{2}^{4} \int_{0}^{1} (x + yz^2) \, dz \, dy \, dxI=∫−15​∫24​∫01​(x+yz2)dzdydx

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

    I=∫−15∫24[∫01(x+yz2) dz]dy dxI = \int_{-1}^{5} \int_{2}^{4} \left[ \int_{0}^{1} (x + yz^2) \, dz \right] dy \, dxI=∫−15​∫24​[∫01​(x+yz2)dz]dydx

  2. Evaluate the Innermost Integral (dzdzdz):
    Treat xxx and yyy as constants:

    ∫01(x+yz2) dz=[xz+yz33]01=(x(1)+13y(1)3)−0=x+13y\int_{0}^{1} (x + yz^2) \, dz = \left[ xz + y \frac{z^3}{3} \right]_{0}^{1} = \left( x(1) + \frac{1}{3}y(1)^3 \right) - 0 = x + \frac{1}{3}y∫01​(x+yz2)dz=[xz+y3z3​]01​=(x(1)+31​y(1)3)−0=x+31​y

  3. Evaluate the Middle Integral (dydydy):
    Substitute the result into the yyy-integral over [2,4][2, 4][2,4]:

    ∫24(x+13y) dy=[xy+16y2]24=(4x+166)−(2x+46)=2x+126=2x+2\int_{2}^{4} \left( x + \frac{1}{3}y \right) \, dy = \left[ xy + \frac{1}{6}y^2 \right]_{2}^{4} = \left( 4x + \frac{16}{6} \right) - \left( 2x + \frac{4}{6} \right) = 2x + \frac{12}{6} = 2x + 2∫24​(x+31​y)dy=[xy+61​y2]24​=(4x+616​)−(2x+64​)=2x+612​=2x+2

  4. Evaluate the Outermost Integral (dxdxdx):
    Substitute the linear expression into the xxx-integral over [−1,5][-1, 5][−1,5]:

    ∫−15(2x+2) dx=[x2+2x]−15=(52+2(5))−((−1)2+2(−1))=(25+10)−(1−2)=35−(−1)=36\int_{-1}^{5} (2x + 2) \, dx = \left[ x^2 + 2x \right]_{-1}^{5} = \left( 5^2 + 2(5) \right) - \left( (-1)^2 + 2(-1) \right) = (25 + 10) - (1 - 2) = 35 - (-1) = 36∫−15​(2x+2)dx=[x2+2x]−15​=(52+2(5))−((−1)2+2(−1))=(25+10)−(1−2)=35−(−1)=36

Show your steps and calculations here.

Practice Problem 3

Evaluate the triple integral over the unit cube B=[0,1]×[0,1]×[0,1]B = [0, 1] \times [0, 1] \times [0, 1]B=[0,1]×[0,1]×[0,1]:

I=∫01∫01∫01ex+2y+3z dz dy dxI = \int_{0}^{1} \int_{0}^{1} \int_{0}^{1} e^{x + 2y + 3z} \, dz \, dy \, dxI=∫01​∫01​∫01​ex+2y+3zdzdydx

Solution Key:
  1. Apply the Product Separability Property:
    State the Algebraic Separability Theorem (defined in Triple Integrals over a Rectangular Region):

    ∭Bg(x)h(y)k(z) dV=(∫abg(x) dx)(∫cdh(y) dy)(∫rsk(z) dz)\iiint_{B} g(x)h(y)k(z) \, dV = \left( \int_{a}^{b} g(x) \, dx \right) \left( \int_{c}^{d} h(y) \, dy \right) \left( \int_{r}^{s} k(z) \, dz \right)∭B​g(x)h(y)k(z)dV=(∫ab​g(x)dx)(∫cd​h(y)dy)(∫rs​k(z)dz)
    Decompose the exponential integrand:
    ex+2y+3z=ex⋅e2y⋅e3ze^{x + 2y + 3z} = e^x \cdot e^{2y} \cdot e^{3z}ex+2y+3z=ex⋅e2y⋅e3z
    Thus:
    I=(∫01ex dx)(∫01e2y dy)(∫01e3z dz)I = \left( \int_{0}^{1} e^x \, dx \right) \left( \int_{0}^{1} e^{2y} \, dy \right) \left( \int_{0}^{1} e^{3z} \, dz \right)I=(∫01​exdx)(∫01​e2ydy)(∫01​e3zdz)

  2. Evaluate Each Single-Variable Factor:

    • xxx-factor: ∫01ex dx=[ex]01=e1−e0=e−1\int_{0}^{1} e^x \, dx = [e^x]_{0}^{1} = e^1 - e^0 = e - 1∫01​exdx=[ex]01​=e1−e0=e−1
    • yyy-factor: ∫01e2y dy=[12e2y]01=12(e2−e0)=e2−12\int_{0}^{1} e^{2y} \, dy = \left[ \frac{1}{2} e^{2y} \right]_{0}^{1} = \frac{1}{2}(e^2 - e^0) = \frac{e^2 - 1}{2}∫01​e2ydy=[21​e2y]01​=21​(e2−e0)=2e2−1​
    • zzz-factor: ∫01e3z dz=[13e3z]01=13(e3−e0)=e3−13\int_{0}^{1} e^{3z} \, dz = \left[ \frac{1}{3} e^{3z} \right]_{0}^{1} = \frac{1}{3}(e^3 - e^0) = \frac{e^3 - 1}{3}∫01​e3zdz=[31​e3z]01​=31​(e3−e0)=3e3−1​
  3. Multiply the Evaluated Terms:

    I=(e−1)⋅(e2−12)⋅(e3−13)=(e−1)(e2−1)(e3−1)6I = (e - 1) \cdot \left(\frac{e^2 - 1}{2}\right) \cdot \left(\frac{e^3 - 1}{3}\right) = \frac{(e - 1)(e^2 - 1)(e^3 - 1)}{6}I=(e−1)⋅(2e2−1​)⋅(3e3−1​)=6(e−1)(e2−1)(e3−1)​

Show your steps and calculations here.

Practice Problem 4

Evaluate the triple integral over the rectangular region EEE:

∭E(3−x+2yz) dVwhere E=[0,2]×[−1,1]×[1,2]\iiint_{E} (3 - x + 2yz) \, dV \quad \text{where } E = [0, 2] \times [-1, 1] \times [1, 2]∭E​(3−x+2yz)dVwhere E=[0,2]×[−1,1]×[1,2]

Solution Key:
  1. Formulate the Iterated Triple Integral:
    Apply Fubini's Theorem for Triple Integrals:

    I=∫02∫−11∫12(3−x+2yz) dz dy dxI = \int_{0}^{2} \int_{-1}^{1} \int_{1}^{2} (3 - x + 2yz) \, dz \, dy \, dxI=∫02​∫−11​∫12​(3−x+2yz)dzdydx

  2. Evaluate the Innermost Integral (dzdzdz):
    Treat xxx and yyy as constants:

    ∫12(3−x+2yz) dz=[(3−x)z+yz2]12=((3−x)(2)+y(4))−((3−x)(1)+y(1))\int_{1}^{2} (3 - x + 2yz) \, dz = \left[ (3 - x)z + yz^2 \right]_{1}^{2} = \Big((3 - x)(2) + y(4)\Big) - \Big((3 - x)(1) + y(1)\Big)∫12​(3−x+2yz)dz=[(3−x)z+yz2]12​=((3−x)(2)+y(4))−((3−x)(1)+y(1))
    =(6−2x+4y)−(3−x+y)=3−x+3y= (6 - 2x + 4y) - (3 - x + y) = 3 - x + 3y=(6−2x+4y)−(3−x+y)=3−x+3y

  3. Evaluate the Middle Integral (dydydy):
    Substitute the expression into the yyy-layer over [−1,1][-1, 1][−1,1]:

    ∫−11(3−x+3y) dy=[(3−x)y+32y2]−11\int_{-1}^{1} (3 - x + 3y) \, dy = \left[ (3 - x)y + \frac{3}{2}y^2 \right]_{-1}^{1}∫−11​(3−x+3y)dy=[(3−x)y+23​y2]−11​
    Notice that 3y3y3y is an odd function over the symmetric interval [−1,1][-1, 1][−1,1], so its integral is zero:
    =((3−x)(1)+32)−((3−x)(−1)+32)=(3−x)−(−(3−x))=2(3−x)=6−2x= \left( (3 - x)(1) + \frac{3}{2} \right) - \left( (3 - x)(-1) + \frac{3}{2} \right) = (3 - x) - (-(3 - x)) = 2(3 - x) = 6 - 2x=((3−x)(1)+23​)−((3−x)(−1)+23​)=(3−x)−(−(3−x))=2(3−x)=6−2x

  4. Evaluate the Outermost Integral (dxdxdx):

    ∫02(6−2x) dx=[6x−x2]02=(6(2)−(2)2)−0=12−4=8\int_{0}^{2} (6 - 2x) \, dx = \left[ 6x - x^2 \right]_{0}^{2} = \left( 6(2) - (2)^2 \right) - 0 = 12 - 4 = 8∫02​(6−2x)dx=[6x−x2]02​=(6(2)−(2)2)−0=12−4=8


Show your steps and calculations here.

Regroup 1

  • Clarification on Iterated Order: Fubini's Theorem guarantees that on a box with constant boundaries, any of the six differential orders (dz dy dx,dx dz dydz\,dy\,dx, dx\,dz\,dydzdydx,dxdzdy, etc.) yields identical results. Choose orders that simplify algebra early (e.g., integrating odd terms over symmetric intervals like [−1,1][-1, 1][−1,1]).
  • Separability Requirement: Separability applies only when the integrand factors into g(x)h(y)k(z)g(x)h(y)k(z)g(x)h(y)k(z) and all six integration limits are constants. Never attempt to separate when limits contain variables.

Skill Block 2

Worked Example A: First-Octant Tetrahedron with Scaled Boundary Plane

Evaluate the triple integral:

∭Ez dV\iiint_{E} z \, dV∭E​zdV
where EEE is the solid tetrahedron in the first octant bounded by the planes x=0x = 0x=0, y=0y = 0y=0, z=0z = 0z=0, and 2x+y+z=22x + y + z = 22x+y+z=2.

Worked Solution:
  1. Determine Inner zzz-Limits:
    State the Type I Region Definition (defined in Triple Integrals over a General Region):

    u1(x,y)≤z≤u2(x,y)u_1(x, y) \le z \le u_2(x, y)u1​(x,y)≤z≤u2​(x,y)
    In the first octant, the floor is the xyxyxy-plane (z=0z = 0z=0). The top boundary is the slanted plane 2x+y+z=2  ⟹  z=2−2x−y2x + y + z = 2 \implies z = 2 - 2x - y2x+y+z=2⟹z=2−2x−y.
    Thus, the inner limits are:
    0≤z≤2−2x−y0 \le z \le 2 - 2x - y0≤z≤2−2x−y

  2. Project onto the xyxyxy-Plane to Form 2D Shadow DDD:
    Set z=0z = 0z=0 in the plane equation:

    2x+y=2  ⟹  y=2−2x2x + y = 2 \implies y = 2 - 2x2x+y=2⟹y=2−2x
    In the first octant (x≥0,y≥0x \ge 0, y \ge 0x≥0,y≥0), the region DDD is a triangle bounded by y=0y = 0y=0 and y=2−2xy = 2 - 2xy=2−2x.
    Middle limits: 0≤y≤2−2x0 \le y \le 2 - 2x0≤y≤2−2x.
    Outer limits: Set y=0  ⟹  2x=2  ⟹  0≤x≤1y = 0 \implies 2x = 2 \implies 0 \le x \le 1y=0⟹2x=2⟹0≤x≤1.

  3. Set Up and Evaluate the Iterated Integral:

    I=∫01∫02−2x∫02−2x−yz dz dy dxI = \int_{0}^{1} \int_{0}^{2-2x} \int_{0}^{2-2x-y} z \, dz \, dy \, dxI=∫01​∫02−2x​∫02−2x−y​zdzdydx

    Inner zzz-integral:

    ∫02−2x−yz dz=[z22]02−2x−y=12(2−2x−y)2\int_{0}^{2-2x-y} z \, dz = \left[ \frac{z^2}{2} \right]_{0}^{2-2x-y} = \frac{1}{2}(2 - 2x - y)^2∫02−2x−y​zdz=[2z2​]02−2x−y​=21​(2−2x−y)2

    Middle yyy-integral: Use uuu-substitution: let u=2−2x−y  ⟹  dy=−duu = 2 - 2x - y \implies dy = -duu=2−2x−y⟹dy=−du.
    When y=0,u=2−2xy = 0, u = 2 - 2xy=0,u=2−2x; when y=2−2x,u=0y = 2 - 2x, u = 0y=2−2x,u=0.

    ∫2−2x012u2(−du)=∫02−2x12u2 du=[u36]02−2x=(2−2x)36=8(1−x)36=43(1−x)3\int_{2-2x}^{0} \frac{1}{2}u^2 (-du) = \int_{0}^{2-2x} \frac{1}{2}u^2 \, du = \left[ \frac{u^3}{6} \right]_{0}^{2-2x} = \frac{(2 - 2x)^3}{6} = \frac{8(1 - x)^3}{6} = \frac{4}{3}(1 - x)^3∫2−2x0​21​u2(−du)=∫02−2x​21​u2du=[6u3​]02−2x​=6(2−2x)3​=68(1−x)3​=34​(1−x)3

    Outer xxx-integral:

    ∫0143(1−x)3 dx=[−49(1−x)4]01=0−(−49(1)4)=49\int_{0}^{1} \frac{4}{3}(1 - x)^3 \, dx = \left[ -\frac{4}{9}(1 - x)^4 \right]_{0}^{1} = 0 - \left(-\frac{4}{9}(1)^4\right) = \frac{4}{9}∫01​34​(1−x)3dx=[−94​(1−x)4]01​=0−(−94​(1)4)=94​
    Thus, ∭Ez dV=49\iiint_E z \, dV = \frac{4}{9}∭E​zdV=94​.

Worked Example B: Spatial Volume Enclosed by a Parabolic Cylinder

Find the volume of the solid region EEE bounded by the parabolic cylinder y=x2y = x^2y=x2 and the planes z=0z = 0z=0 and y+z=4y + z = 4y+z=4.

Worked Solution:
  1. Determine Inner zzz-Limits:
    The solid extends from the floor z=0z = 0z=0 up to the ceiling plane y+z=4  ⟹  z=4−yy + z = 4 \implies z = 4 - yy+z=4⟹z=4−y.
    Thus, 0≤z≤4−y0 \le z \le 4 - y0≤z≤4−y.

  2. Determine 2D Projection Shadow DDD in the xyxyxy-Plane:
    Setting z=0z = 0z=0 in the ceiling plane gives the line y=4y = 4y=4. The lateral wall is the parabolic cylinder y=x2y = x^2y=x2.
    Middle limits (yyy): for any fixed xxx, yyy goes from the parabola up to the line: x2≤y≤4x^2 \le y \le 4x2≤y≤4.
    Outer limits (xxx): the intersection points are x2=4  ⟹  −2≤x≤2x^2 = 4 \implies -2 \le x \le 2x2=4⟹−2≤x≤2.

  3. Formulate and Evaluate the Volume Integral:
    State the Volume via Triple Integration Formula:

    V=∭E1 dV=∫−22∫x24∫04−y1 dz dy dxV = \iiint_{E} 1 \, dV = \int_{-2}^{2} \int_{x^2}^{4} \int_{0}^{4-y} 1 \, dz \, dy \, dxV=∭E​1dV=∫−22​∫x24​∫04−y​1dzdydx

    Inner zzz-integral:

    ∫04−y1 dz=4−y\int_{0}^{4-y} 1 \, dz = 4 - y∫04−y​1dz=4−y

    Middle yyy-integral:

    ∫x24(4−y) dy=[4y−y22]x24=(16−8)−(4x2−x42)=8−4x2+x42\int_{x^2}^{4} (4 - y) \, dy = \left[ 4y - \frac{y^2}{2} \right]_{x^2}^{4} = \left( 16 - 8 \right) - \left( 4x^2 - \frac{x^4}{2} \right) = 8 - 4x^2 + \frac{x^4}{2}∫x24​(4−y)dy=[4y−2y2​]x24​=(16−8)−(4x2−2x4​)=8−4x2+2x4​

    Outer xxx-integral: Exploiting the symmetry of the even integrand:

    V=2∫02(8−4x2+x42) dx=2[8x−43x3+x510]02V = 2 \int_{0}^{2} \left( 8 - 4x^2 + \frac{x^4}{2} \right) \, dx = 2 \left[ 8x - \frac{4}{3}x^3 + \frac{x^5}{10} \right]_{0}^{2}V=2∫02​(8−4x2+2x4​)dx=2[8x−34​x3+10x5​]02​
    V=2(16−323+165)=2(240−160+4815)=2(12815)=25615V = 2 \left( 16 - \frac{32}{3} + \frac{16}{5} \right) = 2 \left( \frac{240 - 160 + 48}{15} \right) = 2\left(\frac{128}{15}\right) = \frac{256}{15}V=2(16−332​+516​)=2(15240−160+48​)=2(15128​)=15256​

Practice Problem 5: Concept & Reading Check

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

  • (a) Define a Type I 3D region, and describe the geometric meaning of the 2D projection shadow DxyD_{xy}Dxy​ in the xyxyxy-plane.
  • (b) In Worked Example A, explain how setting the ceiling surface equation equal to z=0z = 0z=0 reveals the boundary line of the 2D shadow region.
  • (c) In the linked engineering application Aerospace Fuel Tank Center of Mass and Slosh Baffle (Aerospace & Mechanical Engineering), why does determining the vertical center of mass coordinate zˉ\bar{z}zˉ require evaluating two separate triple integrals (MxyM_{xy}Mxy​ and MMM)?
Solution Key:
  • (a) A Type I region in R3\mathbb{R}^3R3 is defined as E={(x,y,z)∣(x,y)∈Dxy, u1(x,y)≤z≤u2(x,y)}E = \{(x, y, z) \mid (x, y) \in D_{xy}, \, u_1(x, y) \le z \le u_2(x, y)\}E={(x,y,z)∣(x,y)∈Dxy​,u1​(x,y)≤z≤u2​(x,y)}. The domain DxyD_{xy}Dxy​ is the orthogonal geometric projection (or shadow) of the 3D solid onto the xyxyxy-coordinate plane, which serves as the 2D domain over which the middle and outer single integrals are bounded.
  • (b) In the first octant, the bottom boundary of the tetrahedron is the horizontal plane z=0z = 0z=0. The boundary where the ceiling surface terminates into the floor occurs along the trace curve where z=0z = 0z=0. Substituting z=0z = 0z=0 into 2x+y+z=22x + y + z = 22x+y+z=2 gives 2x+y=2  ⟹  y=2−2x2x + y = 2 \implies y = 2 - 2x2x+y=2⟹y=2−2x, which forms the hypotenuse boundary of the triangular shadow in the xyxyxy-plane.
  • (c) Center of mass represents the weighted average spatial location of mass throughout the solid. To calculate the vertical coordinate zˉ=MxyM\bar{z} = \frac{M_{xy}}{M}zˉ=MMxy​​, engineers must evaluate the first moment about the xyxyxy-plane Mxy=∭Ezρ(x,y,z) dVM_{xy} = \iiint_E z \rho(x, y, z) \, dVMxy​=∭E​zρ(x,y,z)dV (which weights each parcel of mass by its vertical lever arm zzz) and normalize it by the total accumulated mass M=∭Eρ(x,y,z) dVM = \iiint_E \rho(x, y, z) \, dVM=∭E​ρ(x,y,z)dV.
Show your steps and calculations here.

Practice Problem 6

Evaluate the triple integral over the general region EEE:

∭E(5x−3y) dV\iiint_{E} (5x - 3y) \, dV∭E​(5x−3y)dV
where EEE is the solid tetrahedron in the first octant bounded by x=0x = 0x=0, y=0y = 0y=0, z=0z = 0z=0, and x+y+z=1x + y + z = 1x+y+z=1, using the integration order dz dy dxdz \, dy \, dxdzdydx.

Solution Key:
  1. Set Up the Type I Boundaries:
    State the Type I Region Definition (defined in Triple Integrals over a General Region):

    0≤z≤1−x−y,0≤y≤1−x,0≤x≤10 \le z \le 1 - x - y, \quad 0 \le y \le 1 - x, \quad 0 \le x \le 10≤z≤1−x−y,0≤y≤1−x,0≤x≤1
    Formulate the ordered iterated integral:
    I=∫01∫01−x∫01−x−y(5x−3y) dz dy dxI = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} (5x - 3y) \, dz \, dy \, dxI=∫01​∫01−x​∫01−x−y​(5x−3y)dzdydx

  2. Evaluate the Innermost Integral (dzdzdz):
    Treat xxx and yyy as constants:

    ∫01−x−y(5x−3y) dz=[(5x−3y)z]01−x−y=(5x−3y)(1−x−y)\int_{0}^{1-x-y} (5x - 3y) \, dz = \Big[ (5x - 3y)z \Big]_{0}^{1-x-y} = (5x - 3y)(1 - x - y)∫01−x−y​(5x−3y)dz=[(5x−3y)z]01−x−y​=(5x−3y)(1−x−y)
    Expand the product:
    =5x(1−x)−5xy−3y(1−x)+3y2=(5x−5x2)−(5x+3−3x)y+3y2=(5x−5x2)−(2x+3)y+3y2= 5x(1 - x) - 5xy - 3y(1 - x) + 3y^2 = (5x - 5x^2) - (5x + 3 - 3x)y + 3y^2 = (5x - 5x^2) - (2x + 3)y + 3y^2=5x(1−x)−5xy−3y(1−x)+3y2=(5x−5x2)−(5x+3−3x)y+3y2=(5x−5x2)−(2x+3)y+3y2

  3. Evaluate the Middle Integral (dydydy):

    ∫01−x((5x−5x2)−(2x+3)y+3y2) dy=[(5x−5x2)y−2x+32y2+y3]01−x\int_{0}^{1-x} \left( (5x - 5x^2) - (2x + 3)y + 3y^2 \right) \, dy = \left[ (5x - 5x^2)y - \frac{2x+3}{2}y^2 + y^3 \right]_{0}^{1-x}∫01−x​((5x−5x2)−(2x+3)y+3y2)dy=[(5x−5x2)y−22x+3​y2+y3]01−x​
    Substitute y=1−xy = 1 - xy=1−x and factor (1−x)2(1 - x)^2(1−x)2:
    =5x(1−x)2−2x+32(1−x)2+(1−x)3=(1−x)2[5x−2x+32+(1−x)]= 5x(1 - x)^2 - \frac{2x+3}{2}(1 - x)^2 + (1 - x)^3 = (1 - x)^2 \left[ 5x - \frac{2x+3}{2} + (1 - x) \right]=5x(1−x)2−22x+3​(1−x)2+(1−x)3=(1−x)2[5x−22x+3​+(1−x)]
    =(1−x)2[4x+1−x−32]=(1−x)2(3x−12)= (1 - x)^2 \left[ 4x + 1 - x - \frac{3}{2} \right] = (1 - x)^2 \left( 3x - \frac{1}{2} \right)=(1−x)2[4x+1−x−23​]=(1−x)2(3x−21​)
    Expand the polynomial:
    =(1−2x+x2)(3x−12)=3x3−132x2+4x−12= (1 - 2x + x^2)\left(3x - \frac{1}{2}\right) = 3x^3 - \frac{13}{2}x^2 + 4x - \frac{1}{2}=(1−2x+x2)(3x−21​)=3x3−213​x2+4x−21​

  4. Evaluate the Outermost Integral (dxdxdx):

    ∫01(3x3−132x2+4x−12) dx=[34x4−136x3+2x2−12x]01\int_{0}^{1} \left( 3x^3 - \frac{13}{2}x^2 + 4x - \frac{1}{2} \right) \, dx = \left[ \frac{3}{4}x^4 - \frac{13}{6}x^3 + 2x^2 - \frac{1}{2}x \right]_{0}^{1}∫01​(3x3−213​x2+4x−21​)dx=[43​x4−613​x3+2x2−21​x]01​
    =34−136+2−12=9−26+24−612=112= \frac{3}{4} - \frac{13}{6} + 2 - \frac{1}{2} = \frac{9 - 26 + 24 - 6}{12} = \frac{1}{12}=43​−613​+2−21​=129−26+24−6​=121​

Show your steps and calculations here.

Practice Problem 7

Find the volume of the first-octant tetrahedron bounded by the intercept plane:

x2+y1+z3=1\frac{x}{2} + \frac{y}{1} + \frac{z}{3} = 12x​+1y​+3z​=1
with x≥0,y≥0,z≥0x \ge 0, y \ge 0, z \ge 0x≥0,y≥0,z≥0, setting up and evaluating a triple integral in rectangular coordinates.

Solution Key:
  1. Determine Integration Limits:
    State the Volume via Triple Integration Formula:

    V=∭E1 dVV = \iiint_{E} 1 \, dVV=∭E​1dV

    • Inner zzz-bounds: Solve for z  ⟹  z=3−32x−3yz \implies z = 3 - \frac{3}{2}x - 3yz⟹z=3−23​x−3y. Floor is z=0z = 0z=0.
    • Middle yyy-bounds: Set z=0  ⟹  x2+y=1  ⟹  0≤y≤1−x2z = 0 \implies \frac{x}{2} + y = 1 \implies 0 \le y \le 1 - \frac{x}{2}z=0⟹2x​+y=1⟹0≤y≤1−2x​.
    • Outer xxx-bounds: Set y=0  ⟹  x2=1  ⟹  0≤x≤2y = 0 \implies \frac{x}{2} = 1 \implies 0 \le x \le 2y=0⟹2x​=1⟹0≤x≤2.
      V=∫02∫01−x/2∫03−3x/2−3y1 dz dy dxV = \int_{0}^{2} \int_{0}^{1 - x/2} \int_{0}^{3 - 3x/2 - 3y} 1 \, dz \, dy \, dxV=∫02​∫01−x/2​∫03−3x/2−3y​1dzdydx
  2. Evaluate Inner zzz-Integral:

    ∫03−3x/2−3y1 dz=3−32x−3y=3(1−x2)−3y\int_{0}^{3 - 3x/2 - 3y} 1 \, dz = 3 - \frac{3}{2}x - 3y = 3\left(1 - \frac{x}{2}\right) - 3y∫03−3x/2−3y​1dz=3−23​x−3y=3(1−2x​)−3y

  3. Evaluate Middle yyy-Integral:

    ∫01−x/2(3(1−x2)−3y) dy=[3(1−x2)y−32y2]01−x/2\int_{0}^{1 - x/2} \left( 3\left(1 - \frac{x}{2}\right) - 3y \right) \, dy = \left[ 3\left(1 - \frac{x}{2}\right)y - \frac{3}{2}y^2 \right]_{0}^{1 - x/2}∫01−x/2​(3(1−2x​)−3y)dy=[3(1−2x​)y−23​y2]01−x/2​
    Substitute y=1−x2y = 1 - \frac{x}{2}y=1−2x​:
    =3(1−x2)2−32(1−x2)2=32(1−x2)2=32(1−x+x24)=32−32x+38x2= 3\left(1 - \frac{x}{2}\right)^2 - \frac{3}{2}\left(1 - \frac{x}{2}\right)^2 = \frac{3}{2}\left(1 - \frac{x}{2}\right)^2 = \frac{3}{2}\left(1 - x + \frac{x^2}{4}\right) = \frac{3}{2} - \frac{3}{2}x + \frac{3}{8}x^2=3(1−2x​)2−23​(1−2x​)2=23​(1−2x​)2=23​(1−x+4x2​)=23​−23​x+83​x2

  4. Evaluate Outermost xxx-Integral:

    V=∫02(32−32x+38x2) dx=[32x−34x2+18x3]02V = \int_{0}^{2} \left( \frac{3}{2} - \frac{3}{2}x + \frac{3}{8}x^2 \right) \, dx = \left[ \frac{3}{2}x - \frac{3}{4}x^2 + \frac{1}{8}x^3 \right]_{0}^{2}V=∫02​(23​−23​x+83​x2)dx=[23​x−43​x2+81​x3]02​
    =(32(2)−34(4)+18(8))−0=3−3+1=1= \left( \frac{3}{2}(2) - \frac{3}{4}(4) + \frac{1}{8}(8) \right) - 0 = 3 - 3 + 1 = 1=(23​(2)−43​(4)+81​(8))−0=3−3+1=1
    (Note: This matches the geometric pyramid volume formula V=16⋅a⋅b⋅c=16(2)(1)(3)=1V = \frac{1}{6} \cdot a \cdot b \cdot c = \frac{1}{6}(2)(1)(3) = 1V=61​⋅a⋅b⋅c=61​(2)(1)(3)=1.)

Show your steps and calculations here.

Practice Problem 8

Let E={(x,y,z)∈R3∣x≥0,y≥0,z≥0,x+y+z≤1}E = \{(x, y, z) \in \mathbb{R}^3 \mid x \ge 0, y \ge 0, z \ge 0, x + y + z \le 1\}E={(x,y,z)∈R3∣x≥0,y≥0,z≥0,x+y+z≤1} be a solid tetrahedron with variable mass density ρ(x,y,z)=x+y+z\rho(x, y, z) = x + y + zρ(x,y,z)=x+y+z. Compute the total mass:

M=∭Eρ(x,y,z) dVM = \iiint_{E} \rho(x, y, z) \, dVM=∭E​ρ(x,y,z)dV

Solution Key:
  1. Formulate the Mass Integral:
    State the Total Mass Formula with Variable Density (defined in Triple Integrals over a General Region):

    M=∭Eρ(x,y,z) dV=∫01∫01−x∫01−x−y(x+y+z) dz dy dxM = \iiint_{E} \rho(x, y, z) \, dV = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} (x + y + z) \, dz \, dy \, dxM=∭E​ρ(x,y,z)dV=∫01​∫01−x​∫01−x−y​(x+y+z)dzdydx

  2. Evaluate Inner zzz-Integral:
    Treat xxx and yyy as constants:

    ∫01−x−y(x+y+z) dz=[(x+y)z+12z2]01−x−y=(x+y)(1−x−y)+12(1−x−y)2\int_{0}^{1-x-y} (x + y + z) \, dz = \left[ (x + y)z + \frac{1}{2}z^2 \right]_{0}^{1-x-y} = (x + y)(1 - x - y) + \frac{1}{2}(1 - x - y)^2∫01−x−y​(x+y+z)dz=[(x+y)z+21​z2]01−x−y​=(x+y)(1−x−y)+21​(1−x−y)2
    Express (x+y)(x + y)(x+y) as 1−(1−x−y)1 - (1 - x - y)1−(1−x−y):
    =(1−(1−x−y))(1−x−y)+12(1−x−y)2=(1−x−y)−12(1−x−y)2= \Big(1 - (1 - x - y)\Big)(1 - x - y) + \frac{1}{2}(1 - x - y)^2 = (1 - x - y) - \frac{1}{2}(1 - x - y)^2=(1−(1−x−y))(1−x−y)+21​(1−x−y)2=(1−x−y)−21​(1−x−y)2

  3. Evaluate Middle yyy-Integral:
    Use uuu-substitution: let u=1−x−y  ⟹  dy=−duu = 1 - x - y \implies dy = -duu=1−x−y⟹dy=−du.
    When y=0,u=1−xy = 0, u = 1 - xy=0,u=1−x; when y=1−x,u=0y = 1 - x, u = 0y=1−x,u=0.

    ∫1−x0(u−12u2)(−du)=∫01−x(u−12u2) du=[12u2−16u3]01−x=12(1−x)2−16(1−x)3\int_{1-x}^{0} \left( u - \frac{1}{2}u^2 \right) (-du) = \int_{0}^{1-x} \left( u - \frac{1}{2}u^2 \right) \, du = \left[ \frac{1}{2}u^2 - \frac{1}{6}u^3 \right]_{0}^{1-x} = \frac{1}{2}(1 - x)^2 - \frac{1}{6}(1 - x)^3∫1−x0​(u−21​u2)(−du)=∫01−x​(u−21​u2)du=[21​u2−61​u3]01−x​=21​(1−x)2−61​(1−x)3

  4. Evaluate Outermost xxx-Integral:
    Use www-substitution: let w=1−x  ⟹  dx=−dww = 1 - x \implies dx = -dww=1−x⟹dx=−dw.
    When x=0,w=1x = 0, w = 1x=0,w=1; when x=1,w=0x = 1, w = 0x=1,w=0.

    M=∫10(12w2−16w3)(−dw)=∫01(12w2−16w3) dw=[16w3−124w4]01=16−124=4−124=324=18M = \int_{1}^{0} \left( \frac{1}{2}w^2 - \frac{1}{6}w^3 \right) (-dw) = \int_{0}^{1} \left( \frac{1}{2}w^2 - \frac{1}{6}w^3 \right) \, dw = \left[ \frac{1}{6}w^3 - \frac{1}{24}w^4 \right]_{0}^{1} = \frac{1}{6} - \frac{1}{24} = \frac{4 - 1}{24} = \frac{3}{24} = \frac{1}{8}M=∫10​(21​w2−61​w3)(−dw)=∫01​(21​w2−61​w3)dw=[61​w3−241​w4]01​=61​−241​=244−1​=243​=81​


Show your steps and calculations here.

Regroup 2

  • Common Failure Mode: Outer Variable Dependency: Never allow inner variables (yyy or zzz) to appear in the outer limits (dxdxdx). The outermost integral must always collapse to pure numerical constants, guaranteeing a single scalar real number output.
  • Visualizing Shadows: If students struggle with the 2D projection DDD, instruct them to shine a mental flashlight straight down the zzz-axis onto the xyxyxy-floor. The perimeter where light is blocked defines the curve equations for the middle integral.

Skill Block 3

Worked Example A: Resolving a Non-Elementary Antiderivative via Axis Reordering

Consider the iterated triple integral:

I=∫14∫01∫2y2ex2z dx dy dzI = \int_{1}^{4} \int_{0}^{1} \int_{2y}^{2} \frac{e^{x^2}}{z} \, dx \, dy \, dzI=∫14​∫01​∫2y2​zex2​dxdydz
Rewrite the integral in the order dy dx dzdy \, dx \, dzdydxdz and evaluate its exact value.

Worked Solution:
  1. Analyze Original Bounds:
    Extract the inequalities:

    1≤z≤4,0≤y≤1,2y≤x≤21 \le z \le 4, \quad 0 \le y \le 1, \quad 2y \le x \le 21≤z≤4,0≤y≤1,2y≤x≤2
    The innermost integral ∫ex2 dx\int e^{x^2} \, dx∫ex2dx has no elementary antiderivative. Changing the order of the inner two variables to dy dxdy \, dxdydx is required.

  2. Re-slice the 2D Domain in the xyxyxy-Plane:
    The outer variable zzz remains on [1,4][1, 4][1,4]. In the xyxyxy-plane, the region is bounded by x=2yx = 2yx=2y, x=2x = 2x=2, and y=0y = 0y=0.
    Solve for yyy in terms of xxx:

    x=2y  ⟹  y=x2x = 2y \implies y = \frac{x}{2}x=2y⟹y=2x​
    For the new order dy dxdy \, dxdydx:

    • For a fixed xxx, yyy goes from y=0y = 0y=0 up to y=x2y = \frac{x}{2}y=2x​.
    • xxx spans the constant interval 0≤x≤20 \le x \le 20≤x≤2.
  3. Set Up and Evaluate the Rewritten Integral:

    I=∫14∫02∫0x/2ex2z dy dx dzI = \int_{1}^{4} \int_{0}^{2} \int_{0}^{x/2} \frac{e^{x^2}}{z} \, dy \, dx \, dzI=∫14​∫02​∫0x/2​zex2​dydxdz

    Inner yyy-integral:

    ∫0x/2ex2z dy=[ex2zy]0x/2=xex22z\int_{0}^{x/2} \frac{e^{x^2}}{z} \, dy = \left[ \frac{e^{x^2}}{z} y \right]_{0}^{x/2} = \frac{x e^{x^2}}{2z}∫0x/2​zex2​dy=[zex2​y]0x/2​=2zxex2​

    Middle xxx-integral: 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=2,u=4x = 2, u = 4x=2,u=4.

    ∫02xex22z dx=12z∫0412eu du=14z[eu]04=e4−14z\int_{0}^{2} \frac{x e^{x^2}}{2z} \, dx = \frac{1}{2z} \int_{0}^{4} \frac{1}{2} e^u \, du = \frac{1}{4z} [e^u]_{0}^{4} = \frac{e^4 - 1}{4z}∫02​2zxex2​dx=2z1​∫04​21​eudu=4z1​[eu]04​=4ze4−1​

    Outer zzz-integral:

    ∫14e4−14z dz=e4−14[ln⁡∣z∣]14=e4−14(ln⁡4−ln⁡1)=(e4−1)ln⁡44\int_{1}^{4} \frac{e^4 - 1}{4z} \, dz = \frac{e^4 - 1}{4} \Big[ \ln |z| \Big]_{1}^{4} = \frac{e^4 - 1}{4} (\ln 4 - \ln 1) = \frac{(e^4 - 1)\ln 4}{4}∫14​4ze4−1​dz=4e4−1​[ln∣z∣]14​=4e4−1​(ln4−ln1)=4(e4−1)ln4​

Worked Example B: Reordering across Intersecting Cylindrical and Planar Boundaries

Rewrite the triple integral ∫01∫01−x2∫01−xf(x,y,z) dz dy dx\int_{0}^{1} \int_{0}^{1-x^2} \int_{0}^{1-x} f(x, y, z) \, dz \, dy \, dx∫01​∫01−x2​∫01−x​f(x,y,z)dzdydx into the order dx dy dzdx \, dy \, dzdxdydz.

Worked Solution:
  1. Extract Original Inequalities:

    0≤x≤1,0≤y≤1−x2,0≤z≤1−x0 \le x \le 1, \quad 0 \le y \le 1 - x^2, \quad 0 \le z \le 1 - x0≤x≤1,0≤y≤1−x2,0≤z≤1−x
    Bounding surfaces: x=0,y=0,z=0x = 0, y = 0, z = 0x=0,y=0,z=0, parabolic cylinder y=1−x2  ⟹  x=1−yy = 1 - x^2 \implies x = \sqrt{1 - y}y=1−x2⟹x=1−y​, and plane z=1−x  ⟹  x=1−zz = 1 - x \implies x = 1 - zz=1−x⟹x=1−z.

  2. Determine New Inner xxx-Limits:
    Looking along the xxx-axis from x=0x = 0x=0, the forward boundary is constrained by two distinct surfaces:

    0≤x≤min⁡(1−z,1−y)0 \le x \le \min(1 - z, \sqrt{1 - y})0≤x≤min(1−z,1−y​)

  3. Find the Trace Intersection in the yzyzyz-Plane:
    Equate the two surface equations to determine where the boundary switches:

    1−z=1−y  ⟹  (1−z)2=1−y  ⟹  y=1−(1−z)2=2z−z21 - z = \sqrt{1 - y} \implies (1 - z)^2 = 1 - y \implies y = 1 - (1 - z)^2 = 2z - z^21−z=1−y​⟹(1−z)2=1−y⟹y=1−(1−z)2=2z−z2
    In the yzyzyz-shadow (0≤z≤1,0≤y≤10 \le z \le 1, 0 \le y \le 10≤z≤1,0≤y≤1), the curve y=2z−z2y = 2z - z^2y=2z−z2 partitions the domain into two regions:

    • Region 1 (where 0≤y≤2z−z20 \le y \le 2z - z^20≤y≤2z−z2): The plane is the active bound: 0≤x≤1−z0 \le x \le 1 - z0≤x≤1−z.
    • Region 2 (where 2z−z2≤y≤12z - z^2 \le y \le 12z−z2≤y≤1): The cylinder is the active bound: 0≤x≤1−y0 \le x \le \sqrt{1 - y}0≤x≤1−y​.
  4. Formulate the Split Iterated Integrals:

    ∫01∫02z−z2∫01−zf(x,y,z) dx dy dz+∫01∫2z−z21∫01−yf(x,y,z) dx dy dz\int_{0}^{1} \int_{0}^{2z-z^2} \int_{0}^{1-z} f(x, y, z) \, dx \, dy \, dz + \int_{0}^{1} \int_{2z-z^2}^{1} \int_{0}^{\sqrt{1-y}} f(x, y, z) \, dx \, dy \, dz∫01​∫02z−z2​∫01−z​f(x,y,z)dxdydz+∫01​∫2z−z21​∫01−y​​f(x,y,z)dxdydz

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 changing the order of integration mandatory when facing an integrand such as ex2z\frac{e^{x^2}}{z}zex2​ with inner differential dxdxdx?
  • (b) In Worked Example B, explain how the intersection trace y=2z−z2y = 2z - z^2y=2z−z2 indicates that the reordered integral must be split into two pieces.
  • (c) In the linked engineering application Chemical Reactor Continuous Yield Optimization (Chemical & Environmental Engineering), how does changing the integration order from dx dy dzdx \, dy \, dzdxdydz to dy dx dzdy \, dx \, dzdydxdz enable chemical engineers to avoid expensive numerical simulations?
Solution Key:
  • (a) The function ex2e^{x^2}ex2 does not possess an elementary antiderivative. In the order dx dy dzdx \, dy \, dzdxdydz, the inner integration cannot be completed using standard functions. Reordering to dy dx dzdy \, dx \, dzdydxdz integrates yyy first across the interval [0,x/2][0, x/2][0,x/2], generating the variable factor xxx, which immediately allows uuu-substitution with u=x2,du=2x dxu = x^2, du = 2x \, dxu=x2,du=2xdx.
  • (b) The forward extent of the solid along the xxx-axis is bounded by two competing geometric surfaces: the inclined plane x=1−zx = 1 - zx=1−z and the parabolic cylinder x=1−yx = \sqrt{1 - y}x=1−y​. Setting them equal yields the trace curve y=2z−z2y = 2z - z^2y=2z−z2. Because one surface cuts in earlier for y≤2z−z2y \le 2z - z^2y≤2z−z2 and the other for y≥2z−z2y \ge 2z - z^2y≥2z−z2, the 2D shadow must be partitioned into two distinct sub-domains.
  • (c) Analytical solutions eliminate the need to execute 3D finite-difference or finite-volume Navier-Stokes grid calculations for every reactor geometry. Having an explicit closed-form formula allows engineers to optimize catalyst bed length, width, and height in milliseconds within automated design loops.
Show your steps and calculations here.

Practice Problem 10

Consider the solid tetrahedron EEE in the first octant bounded by x=0,y=0,z=0x = 0, y = 0, z = 0x=0,y=0,z=0, and x+y+z=1x + y + z = 1x+y+z=1. Starting from the iterated integral in the order dz dy dxdz \, dy \, dxdzdydx:

I=∫01∫01−x∫01−x−yxyz dz dy dxI = \int_{0}^{1} \int_{0}^{1-x} \int_{0}^{1-x-y} xyz \, dz \, dy \, dxI=∫01​∫01−x​∫01−x−y​xyzdzdydx
rewrite the integral in the order dx dz dydx \, dz \, dydxdzdy and evaluate its exact value.

Solution Key:
  1. Determine New Bounds for Order dx dz dydx \, dz \, dydxdzdy:
    State the Re-slicing Protocol in the Projection Plane (defined in Changing the Order of Integration in Triple Integrals):

    • Inner xxx-bounds: Isolate xxx from x+y+z=1  ⟹  0≤x≤1−y−zx + y + z = 1 \implies 0 \le x \le 1 - y - zx+y+z=1⟹0≤x≤1−y−z.
    • Shadow in yzyzyz-plane: Set x=0  ⟹  y+z=1x = 0 \implies y + z = 1x=0⟹y+z=1.
    • Middle zzz-bounds: Solving for zzz gives 0≤z≤1−y0 \le z \le 1 - y0≤z≤1−y.
    • Outer yyy-bounds: Intercept is at y=1  ⟹  0≤y≤1y = 1 \implies 0 \le y \le 1y=1⟹0≤y≤1.
      I=∫01∫01−y∫01−y−zxyz dx dz dyI = \int_{0}^{1} \int_{0}^{1-y} \int_{0}^{1-y-z} xyz \, dx \, dz \, dyI=∫01​∫01−y​∫01−y−z​xyzdxdzdy
  2. Evaluate Inner xxx-Integral:
    Treat yyy and zzz as constants:

    ∫01−y−zxyz dx=[12x2yz]01−y−z=12yz(1−y−z)2\int_{0}^{1-y-z} xyz \, dx = \left[ \frac{1}{2}x^2 yz \right]_{0}^{1-y-z} = \frac{1}{2}yz(1 - y - z)^2∫01−y−z​xyzdx=[21​x2yz]01−y−z​=21​yz(1−y−z)2

  3. Evaluate Middle zzz-Integral:
    Expand (1−y−z)2=(1−y)2−2(1−y)z+z2(1 - y - z)^2 = (1 - y)^2 - 2(1 - y)z + z^2(1−y−z)2=(1−y)2−2(1−y)z+z2:

    ∫01−y(12y(1−y)2z−y(1−y)z2+12yz3)dz\int_{0}^{1-y} \left( \frac{1}{2}y(1-y)^2 z - y(1-y)z^2 + \frac{1}{2}yz^3 \right) dz∫01−y​(21​y(1−y)2z−y(1−y)z2+21​yz3)dz
    =[14y(1−y)2z2−13y(1−y)z3+18yz4]01−y= \left[ \frac{1}{4}y(1-y)^2 z^2 - \frac{1}{3}y(1-y)z^3 + \frac{1}{8}yz^4 \right]_{0}^{1-y}=[41​y(1−y)2z2−31​y(1−y)z3+81​yz4]01−y​
    Substitute z=1−yz = 1 - yz=1−y:
    =(14−13+18)y(1−y)4=(6−8+324)y(1−y)4=124y(1−y)4= \left( \frac{1}{4} - \frac{1}{3} + \frac{1}{8} \right) y(1 - y)^4 = \left( \frac{6 - 8 + 3}{24} \right) y(1 - y)^4 = \frac{1}{24}y(1 - y)^4=(41​−31​+81​)y(1−y)4=(246−8+3​)y(1−y)4=241​y(1−y)4

  4. Evaluate Outermost yyy-Integral:
    Use uuu-substitution: let u=1−y  ⟹  y=1−uu = 1 - y \implies y = 1 - uu=1−y⟹y=1−u and dy=−dudy = -dudy=−du.
    When y=0,u=1y = 0, u = 1y=0,u=1; when y=1,u=0y = 1, u = 0y=1,u=0.

    I=∫10124(1−u)u4(−du)=∫01124(u4−u5) du=124[u55−u66]01=124(130)=1720I = \int_{1}^{0} \frac{1}{24}(1 - u)u^4 (-du) = \int_{0}^{1} \frac{1}{24}(u^4 - u^5) \, du = \frac{1}{24} \left[ \frac{u^5}{5} - \frac{u^6}{6} \right]_{0}^{1} = \frac{1}{24}\left(\frac{1}{30}\right) = \frac{1}{720}I=∫10​241​(1−u)u4(−du)=∫01​241​(u4−u5)du=241​[5u5​−6u6​]01​=241​(301​)=7201​

Show your steps and calculations here.

Practice Problem 11

Rewrite the volume integral:

V=∫02∫0x∫0x−y1 dz dy dxV = \int_{0}^{2} \int_{0}^{x} \int_{0}^{x-y} 1 \, dz \, dy \, dxV=∫02​∫0x​∫0x−y​1dzdydx
into the order dx dz dydx \, dz \, dydxdzdy and evaluate the resulting integral to determine the volume.

Solution Key:
  1. Extract Original Inequalities and Invert for xxx:
    From the original integral:

    0≤z≤x−y,0≤y≤x,0≤x≤20 \le z \le x - y, \quad 0 \le y \le x, \quad 0 \le x \le 20≤z≤x−y,0≤y≤x,0≤x≤2
    The lower bound for xxx is the slanted plane z=x−y  ⟹  x=y+zz = x - y \implies x = y + zz=x−y⟹x=y+z.
    The upper bound for xxx is the vertical cut-off plane x=2x = 2x=2.
    Thus, the inner limits are:
    y+z≤x≤2y + z \le x \le 2y+z≤x≤2

  2. Determine 2D Projection onto the yzyzyz-Plane:
    Set the inner limits equal to find the boundary curve:

    y+z=2  ⟹  z=2−yy + z = 2 \implies z = 2 - yy+z=2⟹z=2−y
    In the first quadrant of the yzyzyz-plane (y≥0,z≥0y \ge 0, z \ge 0y≥0,z≥0):

    • Middle zzz-bounds: 0≤z≤2−y0 \le z \le 2 - y0≤z≤2−y.
    • Outer yyy-bounds: 0≤y≤20 \le y \le 20≤y≤2.
  3. Formulate and Evaluate the Rewritten Integral:

    V=∫02∫02−y∫y+z21 dx dz dyV = \int_{0}^{2} \int_{0}^{2-y} \int_{y+z}^{2} 1 \, dx \, dz \, dyV=∫02​∫02−y​∫y+z2​1dxdzdy

    Inner xxx-integral:

    ∫y+z21 dx=[x]y+z2=2−(y+z)=(2−y)−z\int_{y+z}^{2} 1 \, dx = [x]_{y+z}^{2} = 2 - (y + z) = (2 - y) - z∫y+z2​1dx=[x]y+z2​=2−(y+z)=(2−y)−z

    Middle zzz-integral:

    ∫02−y((2−y)−z) dz=[(2−y)z−z22]02−y=(2−y)2−12(2−y)2=12(2−y)2\int_{0}^{2-y} ((2 - y) - z) \, dz = \left[ (2 - y)z - \frac{z^2}{2} \right]_{0}^{2-y} = (2 - y)^2 - \frac{1}{2}(2 - y)^2 = \frac{1}{2}(2 - y)^2∫02−y​((2−y)−z)dz=[(2−y)z−2z2​]02−y​=(2−y)2−21​(2−y)2=21​(2−y)2

    Outer yyy-integral: Use uuu-substitution: let u=2−y  ⟹  dy=−duu = 2 - y \implies dy = -duu=2−y⟹dy=−du.
    When y=0,u=2y = 0, u = 2y=0,u=2; when y=2,u=0y = 2, u = 0y=2,u=0.

    V=∫2012u2(−du)=∫0212u2 du=[u36]02=86=43V = \int_{2}^{0} \frac{1}{2}u^2 (-du) = \int_{0}^{2} \frac{1}{2}u^2 \, du = \left[ \frac{u^3}{6} \right]_{0}^{2} = \frac{8}{6} = \frac{4}{3}V=∫20​21​u2(−du)=∫02​21​u2du=[6u3​]02​=68​=34​

Show your steps and calculations here.

Practice Problem 12

Given the triple integral representing the volume of solid EEE:

V=∫01/2∫01/2∫01−2x+2y1 dz dy dxV = \int_{0}^{1/2} \int_{0}^{1/2} \int_{0}^{1-2x+2y} 1 \, dz \, dy \, dxV=∫01/2​∫01/2​∫01−2x+2y​1dzdydx
rewrite the integral in the order dx dz dydx \, dz \, dydxdzdy (splitting the region into two iterated integrals if necessary) and evaluate the volume.

Solution Key:
  1. Analyze Original Boundaries:
    Original limits:

    0≤x≤12,0≤y≤12,0≤z≤1−2x+2y0 \le x \le \frac{1}{2}, \quad 0 \le y \le \frac{1}{2}, \quad 0 \le z \le 1 - 2x + 2y0≤x≤21​,0≤y≤21​,0≤z≤1−2x+2y
    Solving the ceiling surface for xxx:
    2x=1+2y−z  ⟹  x=1+2y−z22x = 1 + 2y - z \implies x = \frac{1 + 2y - z}{2}2x=1+2y−z⟹x=21+2y−z​
    However, xxx is also bounded by the constant plane x=1/2x = 1/2x=1/2. Thus, the upper bound is:
    x=min⁡(12, 1+2y−z2)x = \min\left(\frac{1}{2}, \, \frac{1 + 2y - z}{2}\right)x=min(21​,21+2y−z​)

  2. Locate the Intersection Trace Curve:
    Set the two constraints equal:

    1+2y−z2=12  ⟹  1+2y−z=1  ⟹  z=2y\frac{1 + 2y - z}{2} = \frac{1}{2} \implies 1 + 2y - z = 1 \implies z = 2y21+2y−z​=21​⟹1+2y−z=1⟹z=2y
    In the yzyzyz-plane where 0≤y≤1/20 \le y \le 1/20≤y≤1/2:

    • For 0≤z≤2y0 \le z \le 2y0≤z≤2y: 1+2y−z2≥12\frac{1 + 2y - z}{2} \ge \frac{1}{2}21+2y−z​≥21​, so the constant wall is the active upper bound: 0≤x≤120 \le x \le \frac{1}{2}0≤x≤21​.
    • For 2y≤z≤1+2y2y \le z \le 1 + 2y2y≤z≤1+2y: the slanted plane is the active upper bound: 0≤x≤1+2y−z20 \le x \le \frac{1 + 2y - z}{2}0≤x≤21+2y−z​.
  3. Formulate the Split Iterated Integrals:

    V=∫01/2∫02y∫01/21 dx dz dy+∫01/2∫2y1+2y∫01+2y−z21 dx dz dyV = \int_{0}^{1/2} \int_{0}^{2y} \int_{0}^{1/2} 1 \, dx \, dz \, dy + \int_{0}^{1/2} \int_{2y}^{1+2y} \int_{0}^{\frac{1+2y-z}{2}} 1 \, dx \, dz \, dyV=∫01/2​∫02y​∫01/2​1dxdzdy+∫01/2​∫2y1+2y​∫021+2y−z​​1dxdzdy

  4. Evaluate the First Integral (V1V_1V1​):

    ∫01/21 dx=12\int_{0}^{1/2} 1 \, dx = \frac{1}{2}∫01/2​1dx=21​
    ∫02y12 dz=[12z]02y=y\int_{0}^{2y} \frac{1}{2} \, dz = \left[ \frac{1}{2}z \right]_{0}^{2y} = y∫02y​21​dz=[21​z]02y​=y
    ∫01/2y dy=[y22]01/2=12(14)=18\int_{0}^{1/2} y \, dy = \left[ \frac{y^2}{2} \right]_{0}^{1/2} = \frac{1}{2}\left(\frac{1}{4}\right) = \frac{1}{8}∫01/2​ydy=[2y2​]01/2​=21​(41​)=81​

  5. Evaluate the Second Integral (V2V_2V2​):

    ∫01+2y−z21 dx=1+2y−z2\int_{0}^{\frac{1+2y-z}{2}} 1 \, dx = \frac{1 + 2y - z}{2}∫021+2y−z​​1dx=21+2y−z​
    ∫2y1+2y1+2y−z2 dz\int_{2y}^{1+2y} \frac{1 + 2y - z}{2} \, dz∫2y1+2y​21+2y−z​dz
    Use uuu-substitution: let u=1+2y−z  ⟹  dz=−duu = 1 + 2y - z \implies dz = -duu=1+2y−z⟹dz=−du.
    When z=2y,u=1z = 2y, u = 1z=2y,u=1; when z=1+2y,u=0z = 1 + 2y, u = 0z=1+2y,u=0.
    ∫10u2(−du)=∫01u2 du=[u24]01=14\int_{1}^{0} \frac{u}{2} (-du) = \int_{0}^{1} \frac{u}{2} \, du = \left[ \frac{u^2}{4} \right]_{0}^{1} = \frac{1}{4}∫10​2u​(−du)=∫01​2u​du=[4u2​]01​=41​
    Now perform the outermost integration with respect to yyy:
    ∫01/214 dy=[14y]01/2=18\int_{0}^{1/2} \frac{1}{4} \, dy = \left[ \frac{1}{4}y \right]_{0}^{1/2} = \frac{1}{8}∫01/2​41​dy=[41​y]01/2​=81​

  6. Combine the Components:

    V=V1+V2=18+18=14V = V_1 + V_2 = \frac{1}{8} + \frac{1}{8} = \frac{1}{4}V=V1​+V2​=81​+81​=41​


Show your steps and calculations here.

Regroup 3

  • Trace Intersection Verification: Always check boundary values at the intersection trace (e.g., when z=2yz = 2yz=2y, both formulas for xxx yield 1/21/21/2). This ensures continuity across the partitioned domains.
  • Why Split? Splitting occurs because geometry does not care about our coordinate axes; when multiple surfaces clip a volume simultaneously, piecewise integration is mathematically required.

Synthesis Wrap-up

  • Core Takeaways:
    1. Triple integrals over rectangular boxes are computed layer-by-layer; Fubini's Theorem guarantees all 6 order permutations yield identical results.
    2. For general regions, follow the 4-step projection protocol: surface ceiling/floor →\to→ 2D projection shadow DDD →\to→ boundary curves →\to→ outer constants.
    3. When an inner antiderivative is non-elementary, reordering the differentials can introduce variable multipliers that unlock standard uuu-substitution.
  • Knewton Alta Homework Preview:
    • Triple Integrals over a Rectangular Region and a General Region
    • Changing the Order of Integration of Triple Integrals