Chapter 5: Double Integration
Section 5.3: Regions with Curved Boundaries
Essentials
A double integral of fx,y over a plane region that can be described by the inequalities a≤x≤b,yBx≤y≤yTx can be iterated as
∫ab∫y=yBxy=yTxfx,y dy dx
Figures 5.3.(1, 2) show some of the regions for which this holds.
use plots,Student:-VectorCalculus in module() local p1,p2,p3,p4,p5,V,q,YT,YB; YT:=3+(2/3)*sqrt(-x^2+8*x-7): YB:=3-(2/3)*sqrt(-x^2+8*x-7): V:=RootedVector(root=[4,1],<0,4>): p1:=shadebetween(YB,YT,x=1..7,scaling=constrained,view=[0..7,0..5],color=red): p2:=plot([YT,YB],x=1..7,color=[black,green],thickness=3): p3:=textplot({[1.7,4.9,typeset(y[T])],[1.7,1.2,typeset(y[B])],[1.5,.3,typeset(x=a)],[6.5,.3,typeset(x=b)]}): p4:=plot([[[1,0],[1,3]],[[7,3],[7,0]]],linestyle=dot,color=blue): p5:=PlotVector(V,color=black): q:=display(p5,p1,p2,p3,p4,labels=[x,y],tickmarks=[0,0]); print(q); end module: end use:
Figure 5.3.1
use plots, Student:-VectorCalculus in module() local p1,p2,p3,p4,p5,p6,q,V,YB,YT; V:=RootedVector(root=[4,1],<0,7>): YT:=6+(2/3)*sqrt(-x^2+8*x-7): YB:=3-(2/3)*sqrt(-x^2+8*x-7): p1:=shadebetween(YB,YT,x=1..7,scaling=constrained,view=[0..7,0..8],color=red): p2:=plot([YT,YB],x=1..7,color=[black,green],thickness=3): p3:=textplot({[1.7,7.9,typeset(y[T])],[1.7,1.2,typeset(y[B])],[1.5,.3,typeset(x=a)],[6.5,.3,typeset(x=b)]}): p4:=plot([[[1,3],[1,6]],[[7,3],[7,6]]],linestyle=dot,thickness=5,color=gold): p5:=plot([[[1,0],[1,3]],[[7,3],[7,0]]],linestyle=dot,color=blue): p6:=PlotVector(V,color=black): q:=display(p6,p1,p2,p3,p4,p5,labels=[x,y],tickmarks=[0,0]); print(q); end module: end use:
Figure 5.3.2
A double integral of fx,y over a plane region that can be described by the inequalities xLy≤x≤xRy,c≤y≤d can be iterated as
∫cd∫x=xLyx=xRyfx,y dx dy
Figures 5.3.(3,4) show some of the regions for which this holds.
use plots, Student:-VectorCalculus in module() local XR,XL,V,q,p1,p2,p3,p4,p5,p6; V:=RootedVector(root=[1,4],<4,0>): XR,XL:=3+(2/3)*sqrt(-y^2+8*y-7), 3-(2/3)*sqrt(-y^2+8*y-7): p1:=implicitplot(9*(x-3)^2+4*(y-4)^2 <= 36,x=0..6,y=0..8,filledregions=true,scaling=constrained,view=[0..6,0..8]): p2:=implicitplot(x=XL,x=1..3,y=1..7,color=green,thickness=3,gridrefine=3): p3:=implicitplot(x=XR,x=3..5,y=1..7,color=black,thickness=3,gridrefine=3): p4:=plot([[[0,7],[3,7]],[[3,1],[0,1]]],linestyle=dot,color=black): p5:=PlotVector(V,color=black): p6:=textplot({[1,7.5,typeset(y=d)],[1,.5,typeset(y=c)],[1,6.,typeset(x[L])],[5,2,typeset(x[R])]}): q:=display(p5,p1,p2,p3,p4,p6,labels=[x,y],tickmarks=[0,0]): print(q); end module: end use:
Figure 5.3.3
use plots, Student:-VectorCalculus, plottools in module() local p1,p2,p3,p4,p5,p6,q,V,YB,YT,p124,pp; V:=RootedVector(root=[1,4],<7,0>): YT:=6+(2/3)*sqrt(-x^2+8*x-7): YB:=3-(2/3)*sqrt(-x^2+8*x-7): p1:=shadebetween(YB,YT,x=1..7,scaling=constrained,color=red): p2:=plot([YT,YB],x=1..7,color=[black,green],thickness=3): p4:=plot([[[1,3],[1,6]],[[7,3],[7,6]]],linestyle=dot,thickness=5,color=gold): p124:=display(p1,p2,p4): pp:=rotate(p124,-Pi/2,[4,4]): p3:=textplot({[1.,6,typeset(x[L])],[7.7,1.2,typeset(x[R])],[1,.5,typeset(y=c)],[1,7.5,typeset(y=d)]}): p5:=plot([[[0,1],[3,1]],[[3,7],[0,7]]],linestyle=dot,color=blue): p6:=PlotVector(V,color=black): q:=display(p6,pp,p3,p5,labels=[x,y],tickmarks=[0,0],view=[0..9,0..8]); print(q); end module: end use:
Figure 5.3.4
Iterated Integrals in Maple
Table 5.3.1 lists two ways to iterate a double integral of fx,y at top level. It uses either the int or Int commands, the first leading to immediate evaluation; the second, to the unevaluated (or inert) integral.
∫x=ax=b∫y=yBxy=yTxfx,y dy dx
Intfx,y,y=yBx..yTx,x=a..b=intfx,y,y=yBx..yTx,x=a..b
IntIntfx,y,y=yBx..yTx,x=a..b=intintfx,y,y=yBx..yTx,x=a..b
∫y=cy=d∫x=xLyx=xRyfx,y dx dy
Intfx,y,x=xLy..xRy,y=c..d=intfx,y,x=xLy..xRy,y=c..d
IntIntfx,y,x=xLy..xRy,y=c..d=intintfx,y,x=xLy..xRy,y=c..d
Table 5.3.1 Iterating a double integral with top-level commands
The first syntax wherein a single call to either Int or int is made with a list as second argument is the preference throughout this Study Guide. Unfortunately, the iterated double-integral template in the Calculus palette is still a place-holder for the iteration of the int command.
Table 5.3.2 details the MultiInt command in the Student MultivariateCalculus package.
The MultiInt command in the Student MultivariateCalculus package behaves much like the preferred syntax for the Int command, namely, that the arguments after the integrand need only be equations connecting a variable name with a range. These equations are not placed in a list, but the order in which they are to be given is the same as for the Int command.
After the equations specifying the region of integration, the MultiInt command can take the additional parameters output = integral or output = steps. The first such parameter causes the command to return the unevaluated integral; the second, the detailed steps by which the integral is evaluated. The option output = value is the default, and it is not needed to obtain the immediate evaluation of the integral.
The MultiInt command itself has been updated to recognize the same pre-defined domains that the modified int command in the VectorCalculus packages have recognized since their inception. The two-dimensional domains are Circle, Ellipse, Rectangle, Sector, and Triangle; the three-dimensional, are Parallelepiped, Sphere, and Tetrahedron. The Region option is used to define an arbitrary region in either two or three dimensions. The syntax for implementing integration over one of these domains is that of the modified int command in the VectorCalculus packages: [vars] = Domain. See Table 5.3.4 for details.
There is access to the MultiInt command via the Context Panel once the Student MultivariateCalculus package has been loaded. For example, bring up the Context Panel on fx,y, an expression in x and y, to serve as the integrand, and follow the path
Student Multivariate Calculus≻Integrate≻Iterated
The pane shown in Figure 5.3.1 appears and is used to fix the coordinate system and the names of the variables in that system. Clicking the OK button launches the pane shown in Figure 5.3.2 where the order of iteration is established, and the limits of integration imposed. The drop-down box to the right of "Output" provides three options, namely, the value of the integral, the inert integral (shown to the right of Figure 5.3.2), and a stepwise evaluation of the iterated integral.
Unfortunately, Context Panel access to the MultiInt command does not include access to the updated recognition of pre-defined domains.
Table 5.3.2 The MultiInt command in the Student MultivariateCalculus package
Figure 5.3.1 Coordinates
Figure 5.3.2 Limits and output
∫ab∫uxvxfx,y dy dx
Table 5.3.3 contains the visualization task template that not only evaluates an iterated double integral, but also displays the planar region over which the integral takes place. In Table 5.3.3, the double integral of fx,y=x y is integrated over the region bounded by y=x2 and y=x. The left-hand graph shows the region, and the arrow guides the eye in visualizing that the first (inner) integration is in the y-direction. The right-hand graph shows the integrand over the region of integration. If fx,y≥0, then this graph shows the solid whose volume is computed by the iterated integral.
The display in the table shows the end-state of the task template. To re-initialize it so that, for example, the buttons under "Value of Integral" work, re-select the Area Element. (If one of these buttons is pressed first, the Area Element will reset to "Select dA.")
Tools≻Tasks≻Browse:
Calculus - Multivariate≻Integration≻Visualizing Regions of Integration≻Cartesian 2-D
Evaluate ∬RΨx,y dA and Graph R
Area Element dA
Select dAdy dxdx dy
, Ψ=
Value of Integral
G=
b=
g=
a=
Bounding Curves
"Volume"
Table 5.3.3 Cartesian 2-D visualization task-template
Table 5.3.4 lists examples of the correct syntax for the two-dimensional domains recognized by the MultiInt command. In each case, the optional parameter output = integral causes the command to return the unevaluated integral.
Region
Syntax
Comments
Circle
MultiIntfx,y,x,y=Circlea,b,r
Center: a,b; radius: r
Ellipse - equation
MultiIntfx,y,x,y=Ellipsea x2+b y2=c
a,b,c positive
Ellipse - parameters
MultiIntfx,y,x,y=Ellipseu,v,a,b,θ
Center: u,v; lengths of semimajor and semiminor axes: a,b; Angle of rotation: θ (optional)
Sector -
MultiIntfx,y,x,y=SectorCirclea,b,r,α,β
Start angle: α; End angle: β
Ellipse
MultiIntfx,y,x,y=SectorEllipseu,v,a,b,α,β
Rectangle
MultiIntfx,y,x,y=Rectanglea..b,c..d
a≤x≤b,c≤y≤d
Triangle
MultiIntfx,y,x,y=Trianglev1,v2,v3
v1,v2,v3are the vertices of the triangle, given as vectors.
(dy dx)
MultiIntfx,y,x,y=Regiona..b,yBx..yTx
∫ab∫yBxyTxfx,y dy dx
(dx dy)
MultiIntfx,y,y,x=Regionc..d,xLy..xRy
∫cd∫xLyxRyfx,y dx dy
Table 5.3.4 Regions understood by the MultiInt command in the Student MultivariateCalculus package
Note the behavior of the Region option for integrating over a general region. The order of the integration is controlled by the list of variable names that follows the integrand. The order of the ranges is precisely the order of the ranges as they appear on the iterated integrals. This is the opposite of the behavior for Int and MultiInt where the range for the inner integral is always given first!
Tables 5.3.(5 - 9) contain examples of task-templates that are based on the commands in Table 5.3.4.
Calculus - Vector≻Integration≻Multiple integration≻2D≻
The center of the disk is given in Cartesian coordinates.
Integrate f over a Disk
fx,yorfr,θ} =
Center h,k: h=, k= Radius =
From θ= to θ=
dθ dr
Table 5.3.5 Task template for integration over a disk
Integrate fx,y over an Ellipse
fx,y=
Equation of ellipse:
Table 5.3.6 Task template for integration over an ellipse
Integrate fx,y over a Rectangle
Inner Integral: ≤x≤ Outer Integral: ≤y≤
dx dy
Table 5.3.7 Task template for integration over a rectangle
Integrate fx,y over a Triangle
Vertices
x1=
x2=
x3=
y1=
y2=
y3=
dy dx
Table 5.3.8 Task template for integration over a triangle
Integrate fx,y over a General Region
∫x=ax=b∫y=uxy=vxf ⅆy ⅆx
vx=
ux=
Table 5.3.9 Task template for integration over a general 2-D region
Examples
In Examples 5.3.(1 - 10), integrate fx,y=1+2 x2+3 y2 over the given region R.
Example 5.3.1
R={x,y | 0≤x≤1,x2≤y≤x}
Example 5.3.2
R={x,y | −2≤x≤3,0≤y≤6+x−x2}
Example 5.3.3
R is the interior and boundary of the triangle whose vertices are 0,0,3,1,2,5.
Example 5.3.4
R={x,y | 0≤x≤1,x3/2≤y≤2−2−x2}
Example 5.3.5
R={x,y | 1−y2≤x≤4−1−y2,−1≤y≤1}
Example 5.3.6
R={x,y | −3≤x≤2,x2−6≤y≤−x}
Example 5.3.7
R is the interior and boundary of the ellipse 4 x2+9 y2=36.
Example 5.3.8
R is the region bounded by the ellipse 4 x2+9 y2=36 and the line y=2−2 x/3.
Example 5.3.9
R is the region bounded by the ellipse 4 x2+9 y2=36 and the line y=2−1+5/3x.
Example 5.3.10
R={x,y | 0≤x≤π,cosx≤y≤1+sinx}
<< Previous Section Table of Contents Next Section >>
© Maplesoft, a division of Waterloo Maple Inc., 2024. All rights reserved. This product is protected by copyright and distributed under licenses restricting its use, copying, distribution, and decompilation.
For more information on Maplesoft products and services, visit www.maplesoft.com
Download Help Document