Strang Calculus Multiple Integral Manual - ch14
Strang Calculus Multiple Integral Manual - ch14
The double integral f (x, y)dA gives the volume between R and t h e surface z = f (x,y) The base is first .
cut into small s q u a r e s of area AA. The volume above the ith piece is approximately f (x*,yi) A A . The limit of
the sum C f ( x i , y i ) A A is the volume integral. Three properties of double integrals are h ( f g ) d A = $$ f d A +
+ $$ g d A and $$ cfdA = c $$ f d A and $/a
f dA= f dA +
f d A if R splits i n t o S and T.
If R is the rectangle 0 5 x 5 4,4 5 y 5 6, the integral $$ x dA can be computed two ways. One is x dy dx,
when the inner integral is xy]: = 2x. The outer integral gives x2]t = 16. When the x integral comes first it
equals $ x dx = ix'l! = 8. Then the y integral equals 8y]q = 16. This is the volume between t h e base
rectangle and t h e p l a n e x = x.
The area R is $$ l d y dx. When R is the triangle between x = 0, y = 22, and y = 1, the inner limits on y
1 1
are hc and 1. This is the length of a t h i n vertical strip. The (outer) limits on x are 0 and 2.The area is ;I.
1 Now the strip is horizontal and the outer integral
In the opposite order, the (inner) limits on x are 0 and 2y.
is :$ f 1 When the density is p(x, y) , the total mass in the region R is
dy = 3. $sp dx dy. The moments are
My = $$ px dx dy and Mx = $$ py dx dy. The centroid has Z = M,/M.
25 so
1 bx
$0
2
dy dx + $1 so
b(2-X)
dy d~ = :$ $y;i(ylb)
dx dy = b 25 f (a, b) - f (a, 0) - f (0, b) + f (O,O)
27 1,'Jt(22 - 3y + 1)dx dy = 4 29 1
: f (x)dz = 1: JJ(~'
ldy dx 31 5 0 , 0 0 0 ~
3 S $ ~ $ ~ x 2 d x d y = ~3 5 2 J o ' / a$Ow
Idx dy = a4
i-A *-L
5 7 ~ C ~ = l Cf b(-, l
n n
a)
is exact for f = 1,x, y, xy 39 Volume 8.5 4 1 Volumes in 2,2 ln(1 + 4)
43 $: $: xydx dy = $: &dy = 1n2;J01 Jo1 xzYdydx = e d x = In2 5,'
45 With long rectangles C yiAA = C AA = 1but s$y dA = i
2 $: 2xy dx = x 2 y ] = ; 2e (e2 -
~ (e2 - I ) ~J2 dy = (e2 - l)$]ie= (e2 - 1)(2e2- 2) = 2(e2 - 112;
$: =71.2 e =
11 11' 2
.L.sZe*
y
= ln2e-ln2 = I n $ = 1.
4 $:ye~~dx= e z ~ ] =
f e2y - e ~ ; $ d ( e ~ y - e ~ )=d [ie2~-e']:
~ = i e 2 - e + -,
21 ' s 3 ,/A=d
2 4 T % T y ] x=
24- - 2 4 m ; the x integral is [$(6 +
- $(3 2 ~ ) ~ ' ~=] $g3l2 +
!-~ - - $4312 + $-
+ +
Note! 3 22 y is not zero in the region of integration.
1
6 The region is above y = x3 and below y = x (from 0 to 1). Area = $ i ( z - x3)dx = 1% - %]A 4
1
= z.
8 The region is below the parabola y = 1 - x2 and above its mirror image y = x2 - 1.
8
1
+
Area = $-,(I - x2 - x2 1)dx = 122- $x3]?1 = 9.
10 The area is all below the axis y = 0, where horirontal strips crws from x = y to x = lyl (which is -y). Note
0
that the y integral stops at y =O. Area = $ ! l $ ' Y d x d y = - 2 y d y = [-y2]El = 1.
12 The strips in Problem 6 from y = x3 up to x are changed to strips from x = y across to z = y113.The outer
integral on y is by chance also from O to 1. Area = J,'(y1/3 - y)dy = [gy4/3
- fy2]Q= f.
14 Between the upper parabola y = 1 - x2 in Problem 8 and the x axis,the strips now cross from the
1 ,p=i
left side x = -4- to the right side x = + d G . This half of the area is So $-- - - dx dy =
$', 24I-y' dy = - $ ( I - y)3/211= $.The other half has strips from left side to right side of y = x2 - 1
or x = fd-. This area is ~ _ O ~ $ - dy ~ d x $).
(also
16 The triangle in Problem 10 had sides x = y, x = -y, and y = -1. Now the strips are vertical. They go
from y = -1 up to y = x on the left side: area = $_O1$f d dx = J : ~ ( X +
1)dx = ? ( x l)2]?1 +
= f . The
strips go from -1 up to y = - z on the right side: area = So: $ dy dx = # ( - x
yl
1)dx = f . +
Check: f f = 1. +
18 The triangle has corners at (0,O) and (-1,0) and (-1, -1). Its area is /_Ol$;' dy dx = J~$_;Y dx dy(= 4).
20 The triangle has corners at (0,O) and (2,4) and (4,4). Horinontal strips go from x = to x = y : :
area = $'jY dx dy = 4. Vertical strips are of two kinds: from y = x up to y = 22 or to y = 4.
O2
~rea=$,$' ye dydx+$,'$~dydx=2+2=4.
22 (Hard Problem) The boundary lines are y = $ x from (-2, -1) to (0,0), and y = -22 from (0,O) to
(1, -2), and y = or x = -3y - 5 from (-2, -1) to (1, -2). (This is the hardest one: note first the
slope - 4.)
Vertical strips go from the third line up to the first or second: area = $:2fi-6,
.. dy dz+
1: 4::- dy dx = + f = 8.
Horinontal strips cross from the k t or third Iines to the second:
q + q = f'*
5 )3
,
area = $ _ ; ' J ' ~dx' ~dy~ + J_O~$;'/~ dx dy =
24 The top of the triangle is (a, b). &om x = 0 to a the vertical strips lead to $, lblc
a bx/a
dy dx =
45 28 3
~ , " ~ ~ d ~ = ~ . ~ ahn d e9 = n& =~2 . = ~ =
56 The area of the quarter-circle 4. The moment is iero around the axis y = 0 (by symmetry): jZ = 0.
The other moment, with a factor 2 that accounts for symmetry of left and right, is
lo&/2 $Gy d y d ~ = 2 $ ~ ( ~ - $ ) d x = 2 [ f - $]Pl2=
$. Then g =
38 The integral $
: x2dx dy has the usual midpoint error -%
for the integral of z2 (see Section 5.8).
The y integral So dy = 1 is done exactly. So the error is
1
-&
(and the same for $$ dy). The
integral of zy is computed exactly. Errors decrease with exponent p = 2, the order of accuracy.
40 The exact integral is I,'$,' = 2 J';l4 j"* = 2 $;I4 sec 8 dB = 2 ln(sec 8 + tan 8 )]:I4 =
2 h(,h 1)- +
42 The exact integral is I,' I,' ez s h r y dz dy = I,1( e - 1)sin dy = (- COB K Y ) ] =
~ ? (e - 1)-
+
We change variables to improve the limits of integration. The disk z2 y2 5 9 becomes the rectangle
0 I. r I 8,O I.8 I. 2*. The inner limits of $ ' d y dz are y = f a .In polar coordinates this area integral
becomes /I r dr d8 = 9%
A polar rectangle has sides dr and r dB. Two sides are not straight but the angles are still 90'. The area
bet-n the circles r = 1 and r = 3 and the rays 0 = 0 and 8 = r / 4 is B(S2 - 12)= 1. The integral $$ z dy dx
changes to $$ r2cos 0 dr dB. This is the moment around the y axis. Then I is the ratio My/M. This is the z
coordinate of the centroid, and it is the average value of x.
In a rotation through a , the point that reaches (u,v ) starts at x = u cos a - v sin a , y = u sin a v cos a. +
A rectangle in the uv plane comes from a in zy. The areas are equal so the stretching factor is
J = 1. This is the determinant of the matrix
$1(u eos a - v sin a)2du dv.
cosa -sina
sina cosa I
. The moment of inertia I$ z2dz dy changes to
For single integrals dz changes to (dx/du)du. For double integrals dx d y changes to J du du with J =
a ( x , y ) / a ( u ; v ) .The stretching factor J is the determinant.of the 2 by 2 matrix ax/au ax/& , ~h~
functions z(u, v ) and y(u, v ) connect an zy region R to a uv region S , and $JR dx dy
ay/au ay/&
=-IIQ
J d u d v = a& of I€.
I
For polar coordinates z = u cos v and y = u sin v (or r sin 8). For z = u,y = u+ 4v the 2 by 2 determinant is
J = 4. A square in the uv plane comes from a parallelogram in xy. In the opposite direction the change has
u = z and u = i(y 1 This J is constant because this change of variables is linear.
- z ) and a new J = 3.
I$,"'"$~rdrdO=f SS=qnarter-~irclewithu~Oandv~0;$,'$~~dudv
6 R is symmetric across the y axis; I,' loQ u du dv = $ divided by area gives (ti, ir) = (4/3r,4 / 3 4
dz;dzyy region R* becomes R in the x*y* plane; dz dy = dz'dy* when region moves
72 ~ , ' l ~ $ ; : ~ ~
J=l az/ar*
ay/ar*
&/a*
ay/ae* I=/ cos 8*
sine*
-r* sin 8*
r*cosB*
1
I, r*dr*d8*
=* JI4
3~14 1
1 I
23 J = 2 1 = 3, area lo so3du dv = 3; J =
1 1
2e2u+v e2u+v
eu+2v - 3e3u+3', J: $
: 3e3u+3vdUdv =
J: (eS+~'- e3')dv = i(c6 - 2e3 + 1)
25 comers (2, Y) = (o,o), (1, o), (1,f (I)), (0, f (0));( t . 1) gives 2 = Y = f (+I; J = I 1
vf.(u)
0
f (u)
27 ~2 = s:l4 $011 '
e-'a dr dB = s:I4 (e- 11.ina 0 - 1)do
29 r = r2dr dB/ r dr dB =:j :a3 sin3 B dB/ra2 = 81 9
2r 1 2
I, r r dr dB = 5 so
33 Along the right side; along the bottom; at the bottom right corner
35 /J zy dz dy = I,' +
(u cos a - v sin a)(usin a v cos a)du dv = (cos2a - sin2 a )
3 7 ~ ~ ~ ~ ~ r ~ t r d r d B = ~~ (Q 5z =~c o
- s4a -~s i) n a , y = s i n a + c o s a g o e s t o u = 1 , v = 1
Is dx dy are lP2J'Yy dx dy for the lower triangle plus IT212 1-,/s dx dy for the circular top.
. -
4 (See Problem 36 of Section 14.1) Jrl4
3r/4 so(r sin B)r dr dB
1
= [fI;[- COB B]$? = 9;divide by area f to reach
a
--
g=.*=.&. 2
3 2
24 Problem 18 has J= I 1 1
I
= 1. So the area of R is l x area of unit square = 1. Problem 20 has
uv
J = I -2u 2v
I
= 2(u2 d), +
and integration over the square gives area of R =
I,' :$ + So
4 . Check in z, y coordinates: area of R = 2 1 (1- z2)dz = $.
2(u2 v2)du dv = 5
34 (a) False (forgot the stretching factor J) (b) False (z can be larger than x2) (c) False (forgot to divide
by the area) (d) Tkue (odd function integrated over symmetric interval) (e) False (the straight-sided
region is a traperoid: angle from 0 to 4 and radius from r l to r2 yields area i ( r z - r f ) sin 6 cos 4).
54 44
36 /$ pdA = So2u $'5 r 2 ( r dr dB) = as+. This is the polar m o m e n t of inertia lowith density p = 1.
38 $1 f dA = f (P)$$ dA is the M e a n Value T h e o r e m for double integrals (compare Property 7, Section
5.6). If f = z or f = y, choose P = centroid (Z,y).
Six important solid shapes are a box, prism, cone, cylinder, t e t r a h e d r o n , a n d sphere. The integral
$$$ dx dy dz adds the volume dx d y d z of small boxes. For computation it becomes t h r e e single integrals.
I
The inner integral dz is the length of a line through the solid. The variables y and z are held constant. The
double integral $$ dz dy is the a r e a of a slice, with z held constant. Then the z integral adds up the volumes of
slices.
If the solid region V is bounded by the planes z = 0, y = 0, z = 0, and z + 2y + 32 = 1, the limits on the
inner z integral are 0 and 1- 2y - 3s. The limits on y are 0 and z(1-
1 1 In
3s). The limits on z are 0 and 3.
the new variables u = x, v = 2y, w = 32, the equation of the outer boundary is u v + +
w = 1. The volume
1 mom dx = du and dy = dv/2 and dz = dw/J, the volume of an xyz box
of the tetrahedron in uuw space is 8.
is dx dy dz = 81d u dv dw. So the volume of V is g1
g.
1$: $:dz dy dz =
30 I y I z 5 z 5 1 and all other orders zzy, yzz,zzy,zyz; all six contain (O,O,O); to contain (1,0,1)
14.3 DipIe h tegrah (page 540)
23 Horizontal slices are circles of area lrr2 = r ( 4 - I); volume = lolr(4 - z)da = 8 r ; centroid
4
has z = 0 , =~0 , r = ~ t Z l ( -
4 z)dz/8r =
25 I = f gives zeros; = fl
f dy dz, = l lo.
: f dx dz, B'I = f dx I,'
27 J!
', $!l(y2 + z2)dx dy dz = y;JJJ x2dV = 3 JJJ(x - T ) 2 d V =
I!, i;
29 1: dx dy dz = 6 S 1 aapesoidal rule is second-order; correct for 1,x, y, z, xy, xz, yz, xyz
1 1 1
2 ~ h e a r e a o f 0 < x sy < z 5 1 i s $ o 1x Y dxdydx.Thefourfacesarex=O,y=x,z= y , a = l .
4 I,' x dx dy dz = J', f dy dz = I,' f dz = &. Divide by the volume a to find f = a;
1
1 S
J,' E :J 1 2 1 '
y dx dy dz = loJo d d y dz = Jo %dz = and j? = It; by symmetry E = 8 .
6 Volume of half-cube = J', l!,
J', dz dy dz = 4.
8 $', /_II dx dy dz = / +
: 2(z 1)dz = [(z l)l]; = 3. +
10 $
,! dz dy dz = I!,Ji1(y 1)dy dz = + w d r= /!, [w]~,
= (tetrahedron).
+ +
#
12 The plane faces are x = 0, y = 0, z = 0, and 2x y r = 4 (which goes through 3 points). The volume
2 4-2s
(4 - 22 - y)dy dx = v d z = [-w]i
2 4-2s 4-2s-y
is $0 /o 0
I dz dy dx = lo/, 12 = -ST-*
= A? 16
Check: Multiply standard volume ) by (4)(4)(2) =
(4-y-2)/2
y.
Check: Double the volume in Problem 11.
14 Put dz last and stop at z = 1: lolo lo
1 4-z 1 4-2
dx dy dz = Jo I, Y d y dz =
$0'
&.- = [ (4;:)'~; = '-3' - 37 12 m*
16 (Still tetrahedron of Problem 12: volume still ). Limits of integration: the top vertex
falh from (0,0,4) onto the y axis at (0, -4,O). The corner (2,0,0) stays on the x axis.
The corner (O,4,0) swings up to (0,0,4). The volume integral is J
: J!~ dz dydz = 7. 16
1 dl-%'
18 The plane z = x cuts the circular base in half, leaving z 2 0. Volume = lo /: dz dy dx =
/ , ' 2 ~ 1 / ~ d x = [ - $ ( l - x ) l o = 3.
2 312 1 2
20 Lying along the x axis the cylinder goes-from x = 0 to x = 6. Its slices are circular disks y2 (z - 1)2= 1 +
resting on the x axis. Volume = lo1- A-
'+-dz
Jl-yl dy dx = still 6 r .
.
22 Change variables to X = f , Y = f , Z = f; then dXdYdZ =
-
Volume = "2".
abc dXdYdZ = l/l
1 Centroid (z,9,g) = ( a x , by,cZ) = (f,q,
f ). (Recall volume ) and centroid ( i , !, 4) of standard
tetrahedron: this is Example 2.)
24 (a) Change variables to X = j,Y = f ,Z = $. Then the solid is Xa Y2 Z2 = 1, a unit sphere of volume + +
,
4%
. Therefore the original volume is (4)(2)( $) = 9 9.
(b) The hypervolume in 4 dimensions is 1
following the pattern of 1 for interval, f for triangle, f for tetrahedron.
m,
26 Average of f =
1/Az
I,
jVf (x, y, a)dV/ J dV = integral of f(x?y,y,.)divided by the volume.
28 Volume of unit cube = Ez1 A
C : I ~ ~ ( A X=) ~1.
SO In one variable, the midpoint rule is correct for the functions 1 and x. In three variables it is correct for
1, X? Y? a? xY?m ?YZ?W E *
52 Simpson's Rule has coefficients ), $ , ) over a unit interval. In three dimensions the 8 corners of the cube will
have coefficients (a)3 = A.
The center will have ( i ) 3 = E.
The centers of the 12 edges will have
= A.
The centers of the 6 faces have = &. (Check: 8(1) 64 12(4) 6(16) = 216.) + + +
When N3 cubes an stacked together, with N small cubes each way, there are only 2N 1 meshpoints +
+
along each direction. This makes (2N 1)3points or about 8 per cube. (Visualire the 8 new points
of the cube as having x , y, z equal to rero or f .)
< <
The solid region 1 p 5 2,O 5 # r ,0 S 8 <_ 2 r is a hollow sphere. Its volume is 4r(a3 - l3)/s. &om
the 4 and t9 integrals the area of a spherical shell at radius p equals 47rp2. Newton discovered that the outside
gravitational attraction of a sphere is the same as for an equal mass located at t h e center.
~ ( I , B , z ) = ( D , $ , o ) ; ( ~ , ~ , B ) = ( D 4, ;( ,r ~
, e), ~ ) = ( 5 , ~ 0 ~ - ~ $ , 5 ) ; ( p , # , e ) = ( 5 J Z , ; , ~ 0 ~ - ~ 9 )
6 ( xy , ) = ( , 1 ) ( 8 z) = ( 1 8 z = r on the positive x axis (z2 0, y = 0(= 8), r = 0)
10 z = cos t, y = $sin t, z = g s n t. The unit sphere intersects the plane y = z.
i
+ + +
12 The surface z = 1 r2 = 1 z2 y2 is a paraboloid (parabola rotated around the z axis). The region is
above the halfdisk 0 5 r 5 1 , O 5 8 5 r. The volume is tr.
<
1 4 This is the volume of a half-cylinder (because of 0 B 5 r) : height 7r, radius r, volume $lr4.
16 The upper surface p = 2 is the top of a sphere. The lower surface p = see 4 is the plane z = p cos 4 = 1.
(The angle # = % is the meeting of sphere and plane, where sec 4 = 2.) The volume is
2~$:'~(-) sin 4 d 4 = 2aI-) cos 4 - &]:I3 = 27rI-t - -+ ) + )I -
- 3'
&
18 The region 1 5 p < 3 is a hollow sphere (spherical shell). The limits 0 < ) 5 f keep the part that lies
above a 45' cone. The volume is y(1-q).
20 Fkom the unit ball p 5 1 keep the part above the cone ) = 1 radian and inside the wedge 0 8 5 1 radian. <
!
Volume = $', sin )d) = z1( l - coo I).
+
23 The curve p = 1- cos is a cardioid in the zz plane (like r = 1- cos 9 in the z y plane). So we have a
cardioid of revolution. Its volume is $
as in Problem 9.3.35.
24 Mass = Jo
2n
So J0
n R
psin)(p+ 1)dp d) dB = JcR3 +2xR2.
1 1
26 Newton's ochievcment The cosine law (see hint) gives cos a = -P . Then integrate :
J$/(~% + &)dV.
uses the same u = D2 - 2pD cos
The second integral is
++
$J J f 5= s.
So
d u 2 D -
%-=
29D 9
The first integral over q5
1
p2 = q2 as in the text: $+- ,,
[*I:I: = &(- - 1 D+p )= D(Df-p,, ..The 9 integral gives 2x and then the p integral is
:$ i*ra D'-p'
2D p2 dp = w. The two integrals give
as Newton hoped and expected.
28 The s m h movement produces a right triangle with hypotenuse AD and almost the same angle a. So the
new small side Aq is A D cos a.
30 f $ q d A = 47rp2D + g.
Divide by 4xp2 to find = D
2
+
for the shell. Then the integral over p gives
Resource: Calculus
Gilbert Strang
The following may not correspond to a particular course on MIT OpenCourseWare, but has been
provided by the author as an individual learning resource.
For information about citing these materials or our Terms of Use, visit: https://ocw.mit.edu/terms.