Chapter 6: Applications of Double Integration
Section 6.3: Surface Area
Example 6.3.12
Obtain the surface integral of gx,y=x y over the part of the top half of the ellipsoid 3 x2+5 y2+7 z2=1 that sits above the first-quadrant region bounded by the ellipse 3 x2+5 y2=1/2. See Example 6.2.9.
Solution
Mathematical Solution
The surface is defined by Fx,y=7−21 x2−35⁢y2/7 , the top half of the ellipsoid 3 x2+5 y2+7 z2=1, so the surface-area element is
dσ=1+Fx2+Fy2dA=17⁢12⁢x2+10⁢y2−73⁢x2+5⁢y2−1dA
Iterating in the order dy dx results in the integral
∫016⁢∫0110⁢10−60 x2x⁢y 17⁢12⁢x2+10⁢y2−73⁢x2+5⁢y2−1 ⅆyⅆx ≐ 0.002381642202
where 3 x2+5 y2=1/2, the equation of the ellipse, has been solved for y=yx, and xR=1/6.
As a function of y, the integrand is sufficiently complicated that Maple cannot find an antiderivative for it. Hence, the iterated integral is evaluated numerically.
Maple Solution - Interactive
Solve 3 x2+5 y2+7 z2=1 for z=zx,y
Context Panel: Assign to a Name≻q
3 x2+5 y2+7 z2=1→assign to a nameq
Type the name q and press the Enter key.
Context Panel: Solve≻Obtain Solutions for≻z
Context Panel: Assign to a Name≻Z
q
3⁢x2+5⁢y2+7⁢z2=1
→solutions for z
17⁢−21⁢x2−35⁢y2+7,−17⁢−21⁢x2−35⁢y2+7
→assign to a name
Z
The simplest approach is to employ the task template in Table 6.3.12(a), after noting that the plane region R is an ellipse.
Tools≻Tasks≻Browse:
Calculus - Vector≻Integration≻Surface Integration≻Surface Defined over an Ellipse
Surface Integral on a Surface Defined over an Ellipse
Integrand
fx,y,z=
Surface
zx,y=
Equation of Ellipse:
From θ= to θ=
Table 6.3.12(a) Task template for surface integration over a sector of an ellipse
A solution from first principles is given in Table 6.3.12(b).
Solve the equation of the ellipse for y=yx
Write the equation of the ellipse.
Context Panel: Solve≻Obtain Solutions for≻y
Context Panel: Assign to a Name≻Y
3 x2+5 y2=1/2→solutions for y110⁢−60⁢x2+10,−110⁢−60⁢x2+10→assign to a nameY
Obtain x=xL and x=xR for the ellipse
Write the equation yT=0 and press the Enter key.
Context Panel: Solve≻Obtain Solutions for≻x
Y1=0
110⁢−60⁢x2+10=0
→solutions for x
−16⁢6,16⁢6
Obtain λ=1+fx2+fy2
Calculus palette: Partial-derivative template
Context Panel: Assign Name
λ=1+∂∂ x Z12+∂∂ y Z12→assign
Write an appropriate iterated integral and evaluate numerically
Calculus palette: Iterated double-integral template
Context Panel: 2-D Math≻Convert To≻Inert Form
Context Panel: Approximate≻10 (digits)
∫01/6∫0Y1x y λ ⅆy ⅆx
∫016⁢6∫0110⁢−60⁢x2+10x⁢y⁢1+9⁢x2−21⁢x2−35⁢y2+7+25⁢y2−21⁢x2−35⁢y2+7ⅆyⅆx
→at 10 digits
0.002381642203
Table 6.3.12(b) Solution from first principles
Maple Solution - Coded
Initialize
Install the Student MultivariateCalculus package.
withStudent:-MultivariateCalculus:
Define the bounding surfaces zx,y=Z± and the bounds of the region R as yx=Y±, and X=±1/6. Note the use of the solve and eval commands.
q≔3 x2+5 y2=1/2:Z≔solve3 x2+5 y2+7 z2=1,z:Y≔solveq,y:X≔solveevalq,y=0, x:
Use the diff command to obtain the partial derivatives with respect to x and y.
Apply the simplify command.
λ≔simplify1+diffZ1,x2+diffZ1,y2
17⁢7⁢12⁢x2+10⁢y2−73⁢x2+5⁢y2−1
Form the integral via the MultiInt command with a pre-defined domain option
Evaluate the integral with the evalf command
S≔MultiIntx y λ,x,y=SectorEllipseq,r,θ,0,π/2,output=integral
∫012⁢π∫0160⁢6⁢10⁢tan⁡θ2+1110+16⁢tan⁡θ217⁢r3⁢cos⁡θ⁢sin⁡θ⁢7⁢12⁢r2⁢cos⁡θ2+10⁢r2⁢sin⁡θ2−73⁢r2⁢cos⁡θ2+5⁢r2⁢sin⁡θ2−1ⅆrⅆθ
evalfS = 0.002381642202
Use the SurfaceInt command from the Student VectorCalculus package
Q≔Student:-VectorCalculus:-SurfaceIntx y,x,y,z=Surfacex,y,Z1,x,y=SectorEllipseq,r,θ,0,π/2,output=integral
∫012⁢π∫0160⁢6⁢10⁢tan⁡θ2+1110+16⁢tan⁡θ2r3⁢cos⁡θ⁢sin⁡θ⁢1+9⁢r2⁢cos⁡θ2−21⁢r2⁢cos⁡θ2−35⁢r2⁢sin⁡θ2+7+25⁢r2⁢sin⁡θ2−21⁢r2⁢cos⁡θ2−35⁢r2⁢sin⁡θ2+7ⅆrⅆθ
evalfQ = 0.002381642202
Because the region of integration is an ellipse, Maple has automatically switched to polar coordinates. Without the list of coordinates r,θ included in the definition of the circle, Maple would use x for r and y for θ, writing a deceptive iterated integral. The integrand of this surface integral is so complicated that any attempt to find a closed-form for the inner integral has been abandoned in favor of a numeric evaluation of the double integral.
A solution from first principles that uses the top-level Int command is given below.
Intx y λ,y=0..Y1,x=0..X1=evalfIntx y λ,y=0..Y1,x=0..X1
∫016⁢6∫0110⁢−60⁢x2+1017⁢x⁢y⁢7⁢12⁢x2+10⁢y2−73⁢x2+5⁢y2−1ⅆyⅆx=0.002381642203
<< Previous Example Section 6.3 Next Example >>
© 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