Double Integrals
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 β¬Rβg(x)h(y)dA=(β«abβg(x)dx)(β«cdβh(y)dy), and compute the spatial average value favgβ=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)) and Type II (horizontal slices: h1β(y)β€xβ€h2β(y)) geometries, determining curve intersections to establish outer constant bounds.
- Execute the 4-step order-switching protocol (Extract β Sketch β Invert β 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 ex2, ysinyβ, and 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 β«abβf(x)dx. When physical phenomena vary across a 2-dimensional surface, integration generalizes to the double integral β¬Dβf(x,y)dA, defined in Multiple Integrals as the infinite limit of a double Riemann sum. Geometrically, if f(x,y)β₯0, the double integral calculates the exact volume of the cylindrical solid lying between the planar region D and the surface z=f(x,y). If f(x,y)=1, the integral yields the planar surface area Area(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) 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], 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: β¬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β=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] and inner limits are functions yβ[g1β(x),g2β(x)]. In Type II domains (horizontal slices), outer limits are constants yβ[c,d] and inner limits are functions xβ[h1β(y),h2β(y)]. When an integrand lacks an elementary antiderivative (e.g., 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
Double Integrals over Rectangular Regions & Spatial Averages
General Overview
Skill Block 1 covers the evaluation of double integrals over rectangular domains using Fubini's Theorem. Students learn to translate double integrals into iterated single integrals, strategically select whether to integrate dx or dy first to simplify calculus mechanics, exploit algebraic product separability, and calculate the spatial average value of a two-variable function across a planar area.
Conceptual Framework
- Double Integral Definition over a Rectangle:
For a continuous function f(x,y) on R=[a,b]Γ[c,d]={(x,y)βR2β£aβ€xβ€b,cβ€yβ€d}, the double integral is:β¬Rβf(x,y)dA=Ξxβ0Ξyβ0βlimβi=1βmβj=1βnβf(xijββ,yijββ)ΞxΞy(Defined in Double Integrals over a Rectangular Region.) - Fubini's Theorem for Rectangular Regions:β¬Rβf(x,y)dA=β«cdβ(β«abβf(x,y)dx)dy=β«abβ(β«cdβf(x,y)dy)dx(Defined in Double Integrals over a Rectangular Region.)
- Product Separability Property:
If f(x,y)=g(x)β h(y), then:β¬Rβg(x)h(y)dA=(β«abβg(x)dx)(β«cdβh(y)dy) - Spatial Average Value Formula:
The spatial average value of f(x,y) over a region R is:favgβ=Area(R)1ββ¬Rβf(x,y)dAwhere for a rectangle, Area(R)=(bβa)(dβc). (Defined in Double Integrals over a Rectangular Region.)
Key Equations & Structural Rules
- Outer Limits Are Always Constants: In any valid iterated integral, the outer limits of integration must strictly be numeric constants.
- Inner Variable Substitution Integrity: When evaluating the inner integral with respect to y, the variable x is treated strictly as an algebraic constant. When evaluating with respect to x, y is treated as a constant.
- Fubini Order Selection Strategy: Before calculating, inspect the integrand for cross-variable couplings like xsin(xy) or yexy. Always choose the inner variable whose differentiation matches the exterior factor, converting integration by parts into simple u-substitution.
Worked Example A: Polynomial Iterated Double Integral over a Rectangle
Evaluate the double integral:
π‘ Reveal Worked Solution
Formulate the Iterated Integral:
State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):β¬Rβf(x,y)dA=β«abββ«cdβf(x,y)dydxSubstitute the domain limits xβ[0,2] and yβ[1,3]:I=β«02ββ«13β(4x3+6xy2)dydxEvaluate the Inner Integral (dy):
Treat x as a constant and integrate with respect to y:β«13β(4x3+6xy2)dy=[4x3y+2xy3]y=1y=3βSubstitute the boundary limits:=(4x3(3)+2x(3)3)β(4x3(1)+2x(1)3)=(12x3+54x)β(4x3+2x)=8x3+52xEvaluate the Outer Integral (dx):
Integrate the resulting single-variable function with respect to x: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=136
Worked Example B: Strategic Order Selection via Fubini's Shortcut
Evaluate the double integral:
π‘ Reveal Worked Solution
Analyze Slicing Strategy:
State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):β¬Rβf(x,y)dA=β«abββ«cdβf(x,y)dydx=β«cdββ«abβf(x,y)dxdy- If we integrate with respect to x first, β«xsin(xy)dx requires integration by parts because x appears both outside and inside the sine function.
- If we integrate with respect to y first, x acts as a constant multiplier. Since βyββ(xy)=x, the factor x is the exact derivative of the inner argument, yielding an immediate single-step antiderivative!
Formulate and Evaluate Inner Integral (dy):
Set up the iterated integral with dy on the inside:I=β«01ββ«0Οβxsin(xy)dydxIntegrate with respect to y, using β«xsin(xy)dy=βcos(xy):β«0Οβxsin(xy)dy=[βcos(xy)]y=0y=Οβ=βcos(Οx)β(βcos(0))=1βcos(Οx)Evaluate Outer Integral (dx):
Substitute into the outer integral:I=β«01β(1βcos(Οx))dx=[xβΟsin(Οx)β]01βEvaluate at the limits:I=(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 dydx to dxdy, 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 y first eliminates the need for integration by parts.
- (c) Referencing Optoelectronic LED Display Panel Average Irradiance, explain physically why the spatial average irradiance Eavgβ divides the total double integral by the rectangular surface area Area(R)=(bβa)(dβc).
Review the Conceptual Framework above and consult the definitions in Double Integrals over a Rectangular Region. For part (c), examine the Spatial Average Value Formula and the discussion in the linked optoelectronics engineering application.
Answer: * (a) Fubini's Theorem holds whenever f(x,y) is continuous on R=[a,b]Γ[c,d]. The limits are constants because rectangular boundaries represent geometric lines parallel to the coordinate axes (x=a,x=b,y=c,y=d).
- (b) The variable x is held constant during y-integration, acting as the exact derivative multiplier of the argument xy (βyββ(xy)=x), enabling direct antiderivation β«xsin(xy)dy=βcos(xy).
- (c) The double integral β¬RβEdA evaluates total emitted radiant flux (watts); dividing by the panel area (square meters) calculates the uniform continuous energy density (watts per square meter) distributed across the display face.
Step-by-Step Execution
- (a) State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):
β¬Rβf(x,y)dA=β«cdββ«abβf(x,y)dxdy=β«abββ«cdβf(x,y)dydxThis identity is valid for all functions f(x,y) that are continuous on the closed rectangle R. The limits are fixed constants because a rectangle in Cartesian space is explicitly defined by intervals xβ[a,b] and yβ[c,d].
- (b) In the expression β«xsin(xy)dy, x is independent of y. By the chain rule, dydβ[βcos(xy)]=xsin(xy). If integrating with respect to x first, x is the integration variable and would require integration by parts with u=x and dv=sin(xy)dx.
- (c) State the Spatial Average Value Formula:
Eavgβ=Area(R)1ββ¬RβE(x,y)dAIn optoelectronic display engineering, light emission varies across the panel due to heating and circuit resistance. The double integral accumulates total radiant power (W). Dividing by the geometric surface area Area(R)=(0.5)(0.4)=0.2Β m2 normalizes the total power to express the mean flux per unit area (W/m2), which is used to calibrate LED driver brightness and prevent visual hotspot banding.
Practice Problem 2
Evaluate the double integral:
Apply Fubini's Theorem by setting up an iterated integral with constant limits. Integrate with respect to x first (or y first), evaluate the antiderivative at the limits, and then divide by the rectangular area using the Spatial Average Value Formula (defined in Double Integrals over a Rectangular Region).
Answer:
Step-by-Step Execution
Formulate the Iterated Integral:
State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):β¬Rβf(x,y)dA=β«13ββ«02β(3x2β2y)dxdyEvaluate the Inner Integral (dx):
Treat y as a constant:β«02β(3x2β2y)dx=[x3β2yx]x=0x=2β=((2)3β2y(2))β0=8β4yEvaluate the Outer Integral (dy):
β«13β(8β4y)dy=[8yβ2y2]13β=(8(3)β2(3)2)β(8(1)β2(1)2)=(24β18)β(8β2)=6β6=0Compute the Spatial Average Value:
State the Spatial Average Value Formula:favgβ=Area(R)1ββ¬Rβf(x,y)dAThe area of the rectangular domain is:Area(R)=(2β0)(3β1)=2Γ2=4favgβ=41β(0)=0
Practice Problem 3
Evaluate the double integral by selecting the strategically simpler order of integration:
Compare the inner integrals: integrating with respect to y first requires integration by parts (β«yexydy), whereas integrating with respect to x first treats y as a constant multiplier (β«yexydx=exy). State Fubini's Theorem for Rectangular Regions and integrate with respect to x first.
Answer:
Step-by-Step Execution
Formulate the Strategically Optimal Iterated Integral:
State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):β¬RβyexydA=β«02ββ«01βyexydxdyEvaluate the Inner Integral (dx):
Holding y constant, y acts as the exact derivative multiplier of xy (βxββ(xy)=y):β«01βyexydx=[exy]x=0x=1β=ey(1)βey(0)=eyβ1Evaluate the Outer Integral (dy):
β«02β(eyβ1)dy=[eyβy]02β=(e2β2)β(e0β0)=e2β2β1=e2β3
Practice Problem 4
Evaluate the double integral:
Use Fubini's Theorem to set up iterated integrals. Notice that the term 2xcosy is product-separable (g(x)h(y)), allowing decoupled integration. State Fubini's Theorem for Rectangular Regions and evaluate.
Answer:
Step-by-Step Execution
Formulate the Iterated Integral:
State Fubini's Theorem for Rectangular Regions (defined in Double Integrals over a Rectangular Region):β¬Rβ(2xcosyβ3y2)dA=β«0Ο/2ββ«03β(2xcosyβ3y2)dxdyEvaluate the Inner Integral (dx):
Treat y as a constant:β«03β(2xcosyβ3y2)dx=[x2cosyβ3y2x]x=0x=3β=(9cosyβ9y2)β0=9cosyβ9y2Evaluate the Outer Integral (dy):
β«0Ο/2β(9cosyβ9y2)dy=[9sinyβ3y3]0Ο/2βEvaluate at the limits:=(9sin(2Οβ)β3(2Οβ)3)β0=9(1)β3(8Ο3β)=9β83Ο3β(Note: Multiplying over the interval xβ[0,3], 3Γ3y2=9y2, whose antiderivative is 3y3=83Ο3β.)
π Regroup 1
- Order Selection Heuristic: How can you determine within 5 seconds whether dydx or dxdy 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
Double Integrals over General Regions (Type I vs. Type II)
General Overview
Skill Block 2 expands double integration to non-rectangular, general planar domains. Students learn to classify regions as Type I (vertical slices bounded by functions of x) or Type II (horizontal slices bounded by functions of y), compute intersection points to define the outer constant interval, and set up iterated integrals with functional inner limits.
Conceptual Framework
- Type I General Region (Vertical Slices):
A planar domain bounded by two vertical lines and two continuous functions of x:D={(x,y)βR2β£aβ€xβ€b,g1β(x)β€yβ€g2β(x)}β¬Dβf(x,y)dA=β«abβ(β«g1β(x)g2β(x)βf(x,y)dy)dx(Defined in Double Integrals over a General Region.) - Type II General Region (Horizontal Slices):
A planar domain bounded by two horizontal lines and two continuous functions of y:D={(x,y)βR2β£cβ€yβ€d,h1β(y)β€xβ€h2β(y)}β¬Dβf(x,y)dA=β«cdβ(β«h1β(y)h2β(y)βf(x,y)dx)dy(Defined in Double Integrals over a General Region.) - Planar Area via Double Integration:
Setting f(x,y)=1 calculates the geometric area of the domain:Area(D)=β¬Dβ1dA
Key Equations & Structural Rules
- Golden Rule of Multiple Integrals: The outer limits must ALWAYS be pure constants; the inner limits may be functions of the outer variable, but NEVER the inner variable itself.
- Intersection Equating Protocol: To find the outer constant limits of a region enclosed by two curves, equate the equations (g1β(x)=g2β(x) or h1β(y)=h2β(y)) and solve for the real roots.
- Top/Bottom Identification: Always test an interior test point between the intersection roots to verify which curve lies strictly above (or to the right of) the other.
Worked Example A: Type I General Region (Parabola and Slanted Line)
Evaluate the double integral:
π‘ Reveal Worked Solution
Find Intersection Points (Outer Bounds):
Set the two boundary curves equal to locate where they intersect:x2=2xβΉx2β2x=0βΉx(xβ2)=0The outer constant limits are x=0 and x=2.Establish the Functional Inner Limits:
Test an intermediate value x=1: y=2(1)=2 and y=(1)2=1. Since 2>1, the line y=2x is the upper boundary and the parabola y=x2 is the lower boundary. Thus, D is a Type I Region:D={(x,y)βR2β£0β€xβ€2,x2β€yβ€2x}Formulate and Evaluate Inner Integral (dy):
State the Type I Double Integral Formula (defined in Double Integrals over a General Region):β¬Dβf(x,y)dA=β«abββ«g1β(x)g2β(x)βf(x,y)dydxI=β«02ββ«x22xβ(x+y)dydxEvaluate with respect to y, treating x as constant:β«x22xβ(x+y)dy=[xy+21βy2]y=x2y=2xβ=(x(2x)+21β(2x)2)β(x(x2)+21β(x2)2)=(2x2+2x2)β(x3+21βx4)=4x2βx3β21βx4Evaluate Outer Integral (dx):
β«02β(4x2βx3β21βx4)dx=[34βx3β41βx4β101βx5]02βSubstitute x=2:=(34β(8)β41β(16)β101β(32))β0=332ββ4β516βFind a common denominator (15):=15160ββ1560ββ1548β=1552β
Worked Example B: Type II General Region (Horizontal Slices)
Evaluate the double integral:
π‘ Reveal Worked Solution
Identify the Type II Boundary Structure:
Because the boundary equations explicitly express x as functions of y, this region is naturally formulated as a Type II Region (defined in Double Integrals over a General Region):- Outer constant limits for y: 0β€yβ€2
- Inner limits for x: Left boundary is x=y2, right boundary is x=y
Formulate the Iterated Integral:
State the Type II Double Integral Formula:β¬Dβf(x,y)dA=β«cdββ«h1β(y)h2β(y)βf(x,y)dxdyI=β«02ββ«y2yβ(x+y)dxdyEvaluate Inner Integral (dx):
Treat y as a constant:β«y2yβ(x+y)dx=[21βx2+yx]x=y2x=yβ=(21βy2+y(y))β(21β(y2)2+y(y2))=23βy2βy3β21βy4Evaluate Outer Integral (dy):
I=β«02β(23βy2βy3β21βy4)dy=[21βy3β41βy4β101βy5]02βEvaluate at y=2:=(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=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.
Consult the Conceptual Framework above and the definitions in Double Integrals over a General Region. For part (c), examine the boundary geometry of the submerged parabolic floodgate.
Answer: * (a) Type I has constant outer bounds for xβ[a,b] and functional inner bounds for yβ[g1β(x),g2β(x)]. Type II has constant outer bounds for yβ[c,d] and functional inner bounds for xβ[h1β(y),h2β(y)].
- (b) The roots x=0,2 define the outer interval endpoints. At x=1, y=2(1)=2>12=1, confirming y=2x is the upper boundary.
- (c) Hydrostatic pressure P(y)=Οg(2βy) is constant along horizontal lines of equal depth, making horizontal strip integration (dxdy) clean and uncoupled.
Step-by-Step Execution
- (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)}, where x has constant limits and y has functional limits.
- Type II (Horizontal Slices): D={(x,y)β£cβ€yβ€d,h1β(y)β€xβ€h2β(y)}, where y has constant limits and x has functional limits.
- (b) Setting g1β(x)=g2β(x)βΉx2=2xβΉx(xβ2)=0 finds the geometric intersection points (0,0) and (2,4), which bound the domain horizontally. Testing x=1β(0,2) yields g2β(1)=2 and g1β(1)=1; since 2>1, y=2x is strictly above y=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=y2 and on the right by x=4, a Type II horizontal slice lies at a uniform water depth (2βy), so P(y)=Οg(2βy) factors out of the inner dx integral as a pure constant.
Practice Problem 6
Evaluate the double integral over the Type I region:
State the Type I Double Integral Formula (defined in Double Integrals over a General Region). Set up the iterated integral with y going from x2 to x and x going from 0 to 1. Remember that β«y1/2dy=32βy3/2.
Answer:
Step-by-Step Execution
Formulate the Type I Iterated Integral:
State the Type I Double Integral Formula:β¬Dβf(x,y)dA=β«abββ«g1β(x)g2β(x)βf(x,y)dydxI=β«01ββ«x2xβy1/2dydxEvaluate the Inner Integral (dy):
Treat x as a constant:β«x2xβy1/2dy=[32βy3/2]y=x2y=xβ=32β(x3/2β(x2)3/2)=32β(x3/2βx3)Evaluate the Outer Integral (dx):
I=32ββ«01β(x3/2βx3)dx=32β[52βx5/2β41βx4]01βEvaluate at the limits:=32β(52ββ41β)=32β(208β5β)=32β(203β)=202β=101β
Practice Problem 7
Evaluate the double integral over the Type II region:
State the Type II Double Integral Formula (defined in Double Integrals over a General Region). Set up horizontal slices where x goes from y to yβ and y goes from 0 to 1.
Answer:
Step-by-Step Execution
Formulate the Type II Iterated Integral:
State the Type II Double Integral Formula:β¬Dβf(x,y)dA=β«cdββ«h1β(y)h2β(y)βf(x,y)dxdyI=β«01ββ«yyββyβdxdyEvaluate the Inner Integral (dx):
Since yβ contains no x, it is treated as a constant:β«yyββyβdx=yβ[x]yyββ=yβ(yββy)=yβy3/2Evaluate the Outer Integral (dy):
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 xy-plane, confirming that Type I and Type II evaluations yield identical numeric results).
Practice Problem 8
An optoelectronic display module features a wedge-shaped active zone R bounded by y=0, y=x, and x=2 (spatial units in meters). Emitted irradiance is given by:
- Compute the total radiant power β¬RβEdA as a Type I double integral.
- Find the spatial average irradiance Eavgβ across the wedge in terms of E0β.
The wedge is bounded horizontally by xβ[0,2] and vertically by yβ[0,x]. Set up the Type I integral, evaluate the inner integral in y, then integrate in x. Use the formula for the area of a triangle (Area=21βbaseΓheight) to find Eavgβ=Area1ββ¬RβEdA.
Answer:
Step-by-Step Execution
Evaluate the Total Irradiance Integral:
State the Type I Double Integral Formula (defined in Double Integrals over a General Region):β¬RβE(x,y)dA=β«02ββ«0xβE0β(1β0.2x)dydxSince the integrand contains no y, it acts as a constant during inner integration:β«0xβE0β(1β0.2x)dy=E0β(1β0.2x)[y]0xβ=E0β(xβ0.2x2)Now integrate with respect to x from 0 to 2:β¬RβEdA=E0ββ«02β(xβ51βx2)dx=E0β[21βx2β151βx3]02βEvaluate at x=2:=E0β(24ββ158β)=E0β(2β158β)=E0β(1530β8β)=1522βE0βCompute Spatial Average Irradiance (Eavgβ):
State the Spatial Average Value Formula:Eavgβ=Area(R)1ββ¬RβEdAThe region R is a right triangle with base 2Β m and height 2Β m:Area(R)=21β(2)(2)=2Β m2Eavgβ=21β(1522βE0β)=1511βE0ββ0.733E0βΒ W/m2
π Regroup 2
- Common Slicing Trap: Why does setting limits like β«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 y introduces piecewise square roots or vertical step discontinuities, whereas solving for x yields single continuous functions).
Skill Block 3
Switching the Order of Integration & Non-Elementary Antiderivatives
General Overview
Skill Block 3 masters the transformation of double integrals between Type I and Type II representations. Students execute the 4-step re-slicing protocol, algebraically invert boundary functions, decompose multi-segment curves into unions of sub-regions, and solve previously "impossible" integrals involving non-elementary antiderivatives (ex2, ysinyβ, 1+y3β).
Conceptual Framework
- The Core Slicing Inversion Principle:
Because inner limits depend on the outer variable, switching the order of integration requires inverting the boundary equations and redrawing the geometric slicing frame:β«abββ«g1β(x)g2β(x)βf(x,y)dydx=β«cdββ«h1β(y)h2β(y)βf(x,y)dxdy(Defined in Switching the Order of Integration.) - The 4-Step Re-slicing Protocol:
- Extract: Read the given bounds to identify current orientation (aβ€xβ€b,g1β(x)β€yβ€g2β(x)).
- Sketch: Plot the boundary equations and identify all vertices on the coordinate plane.
- Invert: Solve the inner equations for the opposite variable (x=gβ1(y)).
- Re-slice: Determine the new outer constant span and inner entrance/exit boundaries.
- Decomposition of Non-Monotonic Boundaries:
If the left or right boundary changes character at an intermediate height y=k, partition the domain:β¬DβfdA=β¬D1ββfdA+β¬D2ββfdA
Key Equations & Structural Rules
- Never Just Swap Limits: Swapping β«01ββ«y1βfdxdyββ«y1ββ«01βfdydx is mathematically illegal because the outer integral would contain a variable y.
- Non-Elementary Antiderivative Detection: When you see integrands like ex2, sin(y2), cos(x2), or ysinyβ whose inner antiderivative cannot be expressed in elementary terms, immediately plan to switch the order of integration.
Worked Example A: Dual-Type Geometric Slicing on a Triangular Domain
Let R be the triangular region with vertices (0,0),(3,0), and (3,2). Evaluate β¬Rβ(2x+3y)dA by setting it up first as a Type I integral, then verifying the result as a Type II integral.
π‘ Reveal Worked Solution
Identify the Boundary Lines:
The triangle has three boundary segments:- Bottom edge: y=0 along the x-axis from x=0 to x=3.
- Right edge: x=3 from y=0 to y=2.
- Hypotenuse edge connecting (0,0) to (3,2): slope m=3β02β0β=32ββΉy=32βxβΊx=23βy.
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β€32βxI=β«03ββ«032βxβ(2x+3y)dydxEvaluate inner integral (dy):β«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=2x2Evaluate outer integral (dx):I=β«03β2x2dx=[32βx3]03β=32β(27)=18Method 2: Type II Setup (Horizontal Slices):
State the Type II Double Integral Formula (defined in Double Integrals over a General Region):- Outer bounds for y: yβ[0,2]
- Horizontal slices enter on the left at the line x=23βy and exit on the right at the vertical line x=3:
I=β«02ββ«23βy3β(2x+3y)dxdyEvaluate inner integral (dx):β«23βy3β(2x+3y)dx=[x2+3yx]x=23βyx=3β=(9+9y)β((23βy)2+3y(23βy))=9+9yβ(49βy2+29βy2)=9+9yβ427βy2Evaluate outer integral (dy):I=β«02β(9+9yβ427βy2)dy=[9y+29βy2β49βy3]02β=9(2)+29β(4)β49β(8)=18+18β18=18Both methods produce the exact same value of 18.
Worked Example B: Overcoming a Non-Elementary Antiderivative
Evaluate the double integral by switching the order of integration:
π‘ Reveal Worked Solution
Execute Step 1 (Extract Original Bounds):
From the integral limits:TypeΒ II:0β€yβ€1,yβ€xβ€1Notice that β«ex2dx has no elementary antiderivative.Execute Step 2 & 3 (Sketch Domain & Invert Curves):
The boundary equations are y=0, y=1, x=1, and x=y.
Inverting x=y gives y=x. The region is a triangle with vertices at (0,0),(1,0), and (1,1).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 x: 0β€xβ€1
- Vertical slices enter at the bottom line y=0 and exit at the slanted line y=x:
0β€yβ€xFormulate the new integral:I=β«01ββ«0xβex2dydx
Evaluate Switched Integral:
The inner integrand ex2 contains no y, so it is treated as a constant:β«0xβex2dy=ex2[y]0xβ=xex2Now evaluate the outer integral using u-substitution: let u=x2βΉdu=2xdxβΉxdx=2duβ: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)ββ«abβfdxdy)?
- (b) In Worked Example B, describe how sketching the triangular region 0β€yβ€1,yβ€xβ€1 allows you to reformulate the limits as 0β€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?
Review the Conceptual Framework above and consult the 4-Step Re-slicing Protocol in Switching the Order of Integration. For part (c), look at how inverting y=xββx=y2 eliminates radical exponents.
Answer: * (a) Outer limits must always be fixed numeric constants so the integral evaluates to a physical scalar. Putting variables on the outer integral produces an invalid variable expression.
- (b) The domain vertices are (0,0),(1,0),(1,1). Slicing vertically means x spans the constant range [0,1] while vertical arrows enter at y=0 and exit at y=x.
- (c) Inverting y=xβ to x=y2 eliminates radical square roots from the inner antiderivative, replacing fractional powers with integer polynomial powers (x2β3x3β).
Step-by-Step Execution
- (a) State the Golden Rule of Iterated Integrals (defined in Switching the Order of Integration):
β¬DβfdAβRThe outer integral must always have constant limits [a,b] or [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β€1 describe horizontal arrows entering at x=y and terminating at the vertical line x=1, starting from y=0 up to y=1. When sketched, the domain is the lower right triangle bounded by y=0,x=1, and y=x. Slicing vertically, x ranges continuously from 0 to 1, and for each x, the vertical ray runs from the baseline y=0 up to the line y=x.
- (c) In aerospace structural analysis, the boundary curve is y=xβ. When integrating in Type I chords, intermediate steps involve fractional powers β«x2xβdx=β«x5/2dx. Re-slicing horizontally as Type II inverts the boundary to the clean polynomial x=y2, so the inner integral integrates pure integer powers β«y24βx2dx=[3x3β]y24β=364βy6β, eliminating radical fractions.
Practice Problem 10
Let R be bounded by y=x3, y=4, and the vertical lines x=β2 and x=2.
- Convert the bounds to a Type II region, splitting the domain into two sub-regions R1β and R2β as necessary.
- Evaluate the total geometric area β¬Rβ1dA.
The cubic curve y=x3 passes through (β2,β8), (0,0), (34β,4), and (2,8). Because the horizontal line y=4 intersects the cubic at x=34β, the left/right boundaries change at y=4. Split the region horizontally into R1β (from y=β8 to y=4) and R2β (from y=4 to y=8). Invert the curve to x=y1/3.
Answer:
Step-by-Step Execution
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/3The domain spans vertically from y=(β2)3=β8 to y=(2)3=8. At y=4, the curve crosses x=41/3β1.587:- Lower Sub-region R1β (from y=β8 to y=4):
Left boundary is the vertical line x=β2; right boundary is the curve x=y1/3.R1β={(x,y)βR2β£β8β€yβ€4,β2β€xβ€y1/3} - Upper Sub-region R2β (from y=4 to y=8):
Left boundary is the curve x=y1/3; right boundary is the vertical line x=2.R2β={(x,y)βR2β£4β€yβ€8,y1/3β€xβ€2}
- Lower Sub-region R1β (from y=β8 to y=4):
Evaluate the Area Integral as the Sum of Two Type II Integrals:
Area(R)=β«β84ββ«β2y1/3β1dxdy+β«48ββ«y1/32β1dxdyEvaluate R1β Integral:
β«β84β[x]β2y1/3βdy=β«β84β(y1/3+2)dy=[43βy4/3+2y]β84β=(43β(4)4/3+2(4))β(43β(β8)4/3+2(β8))Since 44/3=4Γ41/3 and (β8)4/3=((β8)1/3)4=(β2)4=16:=(3(41/3)+8)β(43β(16)β16)=3(41/3)+8β(12β16)=3(41/3)+12Evaluate R2β Integral:
β«48β[x]y1/32βdy=β«48β(2βy1/3)dy=[2yβ43βy4/3]48β=(2(8)β43β(8)4/3)β(2(4)β43β(4)4/3)Since 84/3=(81/3)4=24=16:=(16β43β(16))β(8β3(41/3))=(16β12)β8+3(41/3)=3(41/3)β4Sum the Components:
TotalΒ Area=(3(41/3)+12)+(3(41/3)β4)=6(41/3)+8=634β+8β17.524
Practice Problem 11
Evaluate the double integral by switching the order of integration:
The inner integral β«ysinyβdy is non-elementary. Extract the Type I limits: 0β€xβ€Ο,xβ€yβ€Ο. Sketch the triangular region and re-slice as Type II: y ranges from 0 to Ο, and for each y, x ranges from 0 to y. Apply the 4-Step Re-slicing Protocol (defined in Switching the Order of Integration).
Answer:
Step-by-Step Execution
Extract Original Boundaries & Recognize Non-Elementary Integrand:
The original Type I bounds are:0β€xβ€Ο,xβ€yβ€ΟBecause β«ysinyβdy cannot be solved in elementary terms, we must switch to Type II.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,Ο), and (Ο,Ο).
- The vertical extent is 0β€yβ€Ο.
- A horizontal slice enters at the left vertical line x=0 and exits at the line x=y:
0β€xβ€yFormulate the Type II integral:I=β«0Οββ«0yβysinyβdxdy
Evaluate the Switched Iterated Integral:
Evaluate inner integral (dx):β«0yβysinyβdx=ysinyβ[x]0yβ=ysinyβ(yβ0)=sinyEvaluate outer integral (dy):I=β«0Οβsinydy=[βcosy]0Οβ=βcos(Ο)β(βcos(0))=β(β1)β(β1)=1+1=2
π Regroup 3
- Order Reversal as an Unlocking Key: Notice how in Problems 9 and 11, functions that appeared mathematically "unintegrable" (ex2 and ysinyβ) instantly dissolved into high-school u-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:
- Double integrals extend accumulation into 2D space, calculating volumes, areas, total physical charges, and spatial averages.
- Over rectangular regions, Fubini's Theorem permits free order swapping; choose the order that creates chain-rule derivatives and avoids integration by parts.
- 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.
- 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.