Week 5
Week 5
Week 5
1 X i j k
f( , , ) .
n3 n n n
(i/n,j/n,k/n)∈E
As in two dimensions, triple integrals can be evaluated by iterated 1D integral computations. Here
is a simple example:
1 If E is the box {x ∈ [1, 2], y ∈ [0, 1], z ∈ [0, 1]} and f (x, y, z) = 24x2 y 3 z.
R1 R1 R1
0 0 0
24x2 y 3 z dz dy dx .
R1 R1
We start from the core 0
24x2 y 3 z dz = 12x3 y 3 , then integrate the middle layer, 0 12x3 y 3 dy
R2
= 3x2 and finally and finally handle the outer layer: 1 3x2 dx =7.
For the most inner integral, x = x0 and y = y0 are fixed. The integral is integrating up the
function z → f (x0 , y0 , z) along the part intersecting the body. After completing the middle
integral, we have computed the integral on the plane z = const intersected with R. The
most outer integral sums up all these 2-dimensional sections.
In calculus, two important reductions are used to compute triple integrals. In single variable
calculus, one reduces the problem directly to a one dimensional integral by slicing the body
along an axes. The slices are 2-dimensional. In multi-variable calculus, we usually reduce the
problem to an integration problem in two dimensions. This is more flexible.
The single variable method slices the solid along a line. If g(z) is the double
Rb
integral along the two dimensional slice, then a g(z) dz. The multi-variable
method sees the solid sandwiched between the graphs of two functions g(x, y) and
h(x, y) over a common two dimensional region R. The integral reduces to a double
R R R h(x,y)
integral [
R g(x,y)
f (x, y, z) dz]dA.
below the graph RofRa function f (x, y) and above a region R, considered part of the xy-plane.
It is the integral R
f (x, y) dA. We actually have expressed this now as a triple integral. It
is more natural to think of volume as a triple integral also when considering physical units
as we would measure volume using cubic meters for example.
1
3 Find the volume of the unit sphere. Solution: The sphere is sandwiched between the graphs
of two functions obtained by solving for z. Let R be the unit disc in the xy plane. If we use
the sandwich method, we get
Z Z Z √1−x2 −y2
V = [ √ 1dz]dA .
R − 1−x2 −y 2
RR p
which gives a double integral 2 1 − x2 − y 2 dA which is of course best solved in polar
R 2π R 1 √ R
coordinates. We have 0 0 1 − r2 r drdθ = 4π/3.
With the washer method which is √ in this case also called disc method, we slice along
2 2
the z axes and get a disc of radius 1 − z with area R 1 π(1 − z 2). This is a method suitable
for single variable calculus because we get directly −1 π(1 − z ) dz = 4π/3.
4
RRR
The mass of a body with density ρ(x, y, z) is defined as R
ρ(x, y, z) dV . For bodies with
constant density ρ, the mass is ρV , where V is the volume. Compute the mass of a body
which is bounded by the parabolic cylinder z = 4 − x2 , and the planes x = 0, y = 0, y =
6, z = 0 if the density of the body is z. Solution:
Z 2 Z 6 Z 4−x2 Z 2Z 6
z dz dy dx = (4 − x2 )2 /2 dydx
0
Z 20 0 5
0 0
3
x 8x
= 6 (4 − x2 )2 /2 dx = 6( − + 16x)|20 = 2 · 512/5
0 5 3
2
The problem of computing volumes has been tackled early in mathematics:
Assume, we hit 5277 of n=10000 times. The volume so measured is 0.5277. The actual volume
of 1/8’th of the sphere is π/6 = 0.524. For n → ∞ the Monte Carlo computation gives the
actual volume. The Monte-Carlo integral is stronger than the Riemann integral. The law of large
3
numbers in probability theory proves it to be equivalent to the Lebesgue integral and allows to
measure much more sets than solids with piecewise smooth boundaries.
2) Is there an ”integral’ which can measure every solid in space and which has the property
that the volume of a rotated or translated body remains the same? No! Many sets turn out to
be ”crazy” in the sense that one can not measure their volume. An example is the paradox of
Banach and Tarski which tells that one can slice up the unit ball x2 + y 2 + z 2 ≤ 1 into 5 pieces
A, B, C, D, E, rotate and translate them in space so that the pieces A, B, C fit together to be a
unit ball and D, E form a second unit ball. Since the volume has obviously doubled and volume
should be additive, some of the sets A, B, C, D, E do not have a well defined volume.
Homework
1 Evaluate the triple integral
Z 5Z zZ 2y
2
ze−y dxdydz .
0 0 0
E = {x2 + y 2 ≤ z 2 0 ≤ z ≤ 5 } ,
which has the z-axis as its center of symmetry.
4 Integrate f (x, y, z) = x2 +y 2 −z over the tetrahedron with vertices
(0, 0, 0), (4, 4, 0), (0, 4, 0), (0, 0, 12).
5 This is a classic by Archimedes: What is the volume of the body
obtained by intersecting the solid cylinders x2 + z 2 ≤ 1 and y 2 +
z 2 ≤ 1?
4
Math S21a: Multivariable calculus Oliver Knill, Summer 2017
ZZ ZZ
f (x, y, z) dxdydz = g(r, θ, z) r drdθdz
T (R) R
Spherical coordinates use ρ, the distance to the origin as well as two Euler angles: θ the polar
angle and φ, the angle between the vector and the z axis. The coordinate change is
The integration factor measures the volume of a spherical wedge which is dρ, ρ sin(φ) dθ, ρdφ =
ρ2 sin(φ)dθdφdρ.
ZZ ZZ
f (x, y, z) dxdydz = g(ρ, θ, z) ρ2 sin(φ) dρdθdφ
T (R) R
1
The momentRofR Rinertia of a body G with respect to an axis L is defined as the
triple integral G
r(x, y, z)2 dzdydx, where r(x, y, z) = ρ sin(φ) is the distance
from the axis L.
4 R5 8πR5
= (− cos(φ) + cos(φ)3 /3)|π0 (L5 /5)(2π) = · · 2π = .
3 5 15
If the sphere rotates with angular velocity ω, then Iω 2 /2 is the kinetic energy of that sphere.
Example: the moment of inertia of the earth is 8·1037 kgm2 . The angular velocity is ω = 2π/day =
2π/(86400s). The rotational energy is 8 · 1037 kgm2 /(7464960000s2 ) ∼ 1029 J ∼ 2.51024 kcal.
3 Find the volume and the center of mass of a diamond, √ the intersection of the unit sphere
with the cone given in cylindrical coordinates as z = 3r.
Solution: we use spherical coordinates to find the center of mass
Z 1 Z 2π Z π/6
1
x = =0 ρ3 sin2 (φ) cos(θ) dφdθdρ
0 0 0 V
Z 1 Z 2π Z π/6
1
y = ρ3 sin2 (φ) sin(θ) dφdθdρ = 0
0 0 0 V
Z 1 Z 2π Z π/6
1 2π
z = ρ3 cos(φ) sin(φ) dφdθdρ =
0 0 0 V 32V
2
RRR 2
4 Find R
z dV for the solid obtained by intersecting {1 ≤ x2 + y 2 + z 2 ≤ 4 } with the
double cone {z 2 ≥ x2 + y 2 }.
Solution: since the result for the double cone is twice the result for the single cone, we
work with the diamond shaped region R in {z > 0} and multiply the result at the end with
2. In spherical coordinates, the solid R is given by 1 ≤ ρ ≤ 2 and 0 ≤ φ ≤ π/4. With
z = ρ cos(φ), we have
Z 2 Z 2π Z π/4
ρ4 cos2 (φ) sin(φ) dφdθdρ
1 0 0
5 5
2 1 − cos3 (φ)) π/4 31
=( − )2π( |0 = 2π (1 − 2−3/2 ) .
5 5 3 5
√ 3
The result for the double cone is 4π(31/5)(1 − 1/ 2 ) .
Remarks: There are other coordinate systems besides Euclidean, cylindrical and spherical. One
of them are torus coordinates, where T (r, φ, θ) = (1+r cos(φ)) cos(θ), (1+r cos(φ) sin(θ), r sin(φ)),
a coordinate system which works inside the solid torus r ≤ 1. Are there spherical coordinates
in higher dimensions? Yes, there are. They are called hyperspherical coordinates. In four
dimensions (the space of quaternions) for example we would have a third angle ψ and get
(x, y, z, w) = (ρ sin(ψ) sin(φ) sin(θ), ρ sin(ψ) sin(φ) cos(θ), ρ sin(ψ) cos(φ), ρ cos(ψ)) .
The four dimensional case is especially interesting because one can write the sphere S 3 in four
dimensions as the set of pairs of complex numbers z, w with |z|2 + |w|2 = 1. The 3 sphere is
special because it is the group SU (2) of all unitary 2 × 2 matrices of determinant 1. It is also the
set of all quaternions of length 1. The quaternions are historically interesting for multivariable
calculus because they predated vector calculus we teach today and incorporate both the dot and
cross product.
3
Homework
1 Assume the density of a solid E = x2 + y 2 − z 2 < 1, −1 <
z < 1 is given by the forth power of the distance to the z-axes:
σ(x, y, z) = (x2 + y 2)2. Find its mass
Z Z Z
M= (x2 + y 2)2 dxdydz .
E
RRR 2 2
2 Find the moment of inertia E (x + y ) dV of the body E
whose volume is given by the integral
Z π/4 Z π/2 Z 3
Vol(E) = ρ2 sin(φ) dρdθdφ .
0 0 0
4
Math S21a: Multivariable calculus Oliver Knill, Summer 2017
An example is the dipol field F~ (x, y) = hx − 1, yi/((x − 1)2 + y 2 )3/2 − hx + 1, yi/((x + 1)2 + y 2 )3/2
generated by a positive and negative point charge. Here is a picture
When is a vector field a gradient field? F~ (x, y) = hP (x, y), Q(x, y)i = ∇f (x, y) implies Qx (x, y) =
Py (x, y). If this does not hold at some point, F~ is no gradient field.
Clairaut test: If Qx (x, y) − Py (x, y) is not zero at some point, then F~ (x, y) =
hP (x, y), Q(x, y)i is not a gradient field.
We will see that curl(F~ ) = Qx − Py = 0 is also sufficient for F~ to be a gradient field if F~ is defined
everywhere. How do we get f the function with F~ = ∇f ? We will look at examples in class.
2 Is the vector field F~ (x, y) = hxy, 2xy 2 i a gradient field? Solution: No: Qx − Py = 2y 2 − x
is not zero.
1
Vector fields in weather forecast On weather maps, one can see isoterms, curves of constant
temperature or isobars, curves p(x, y) = c of constant pressure. These are level curves. The wind
velocity F~ (x, y) is close but not always exactly perpendicular to the isobars, the lines of equal
pressure p. In reality, the scalar pressure field p and the velocity field F~ also depend on time. The
equations which describe the weather dynamics are called the Navier Stokes equations
d ~ ~
F + F · ∇F~ = ν∆F~ − ∇p + f, divF~ = 0
dt
(where ∆ and div are defined later. This is an other example of a partial differential equation.
It is one of the millenium problems to prove that these equations have smooth solutions in space.
Vector fields are important in differential equations. We look at some examples in population
dynamics and mechanics. You can skip this motivational part:
3 Let x(t) denote the population of a ”prey species” like tuna fish and y(t) is the population
size of a ”predator” like sharks. We have x′ (t) = ax(t) + bx(t)y(t) with positive a, b because
both more predators and more prey species will lead to prey consumption. The rate of
change of y(t) is y ′ (t) = −cy(t) + dxy, where c, d are positive. This can be written using a
vector field ~r′ = F~ (~r(t)). We have a negative sign in the first part because predators would
die out without food. The second term is explained because both more predators as well as
more prey leads to a growth of predators through reproduction. A concrete example is the
Volterra-Lodka system
ẋ = 0.4x − 0.4xy
ẏ = −0.1y + 0.2xy ,
4 A class vector fields important in mechanics are Hamiltonian fields: If H(x, y) is a function
of two variables, then hHy (x, y), −Hx (x, y)i is called a Hamiltonian vector field. An ex-
ample is the harmonic oscillator H(x, y) = (x2 +y 2 )/2. Its vector field (Hy (x, y), −Hx (x, y)) =
(y, −x). The flow lines of a Hamiltonian vector fields are located on the level curves of H
(as you have shown in th homework with the chain rule).
5 Newton’s law m~r′′ = F relates the acceleration ~r′′ of a body with the force F acting at
the point. For example, if x(t) is the position of a mass point in [−1, 1] attached at two
springs and the mass is m = 2, then the point experiences a force (−x + (−x)) = −2x so
2
that mx′′ = 2x or x′′ (t) = −x(t). If we introduce y(t) = x′ (t) of t, then x′ (t) = y(t) and
y ′ (t) = −x(t). Of course y is the velocity of the mass point, so a pair (x, y), thought of as
an initial condition, describes the system so that nature knows what the future evolution of
the system has to be given that data.
6 We don’t yet know yet the curve t 7→ ~r(t) = hx(t), y(t)i, but we know the tangents ~r ′ (t) =
hx′ (t), y ′ (t)i = hy(t), −x(t)i. In other words, we know a direction at each point. The
equation (x′ = y, y ′ = −x) is called a system of ordinary differential equations (ODE’s)
More generally, the problem when studying ODE’s is to find solutions x(t), y(t) of equations
x′ (t) = f (x(t), y(t)), y ′ (t) = g(x(t), y(t)). Here we look for curves x(t), y(t) so that at any
given point (x, y), the tangent vector (x′ (t), y ′ (t)) is (y, −x). You can check by differentiation
that the circles (x(t), y(t)) = (r sin(t), r cos(t)) are solutions. They form a family of curves.
7 If x(t) is the angle of a pendulum, then the gravity acting on it produces a force G(x) =
−gm sin(x), where m is the mass of the pendulum and where g is a constant. For example,
if x = 0 (pendulum at bottom) or x = π (pendulum at the top), then the force is zero. The
Newton equation ”mass times acceleration = force” gives
ẍ(t) = −g sin(x(t)) .
8 The equation of motion for the pendulum ẍ(t) = −g sin(x(t)) can be written with y = ẋ
also as
d
(x(t), y(t)) = (y(t), −g sin(x(t)) .
dt
Each possible motion of the pendulum x(t) is described by a curve ~r(t) = (x(t), y(t)). Writ-
ing down explicit formulas for (x(t), y(t)) is in this case not possible with known functions
like sin, cos, exp, log etc. However, one still can understand the curves.
Curves on the top of the picture represent situations, where the velocity y is large. They describe
the pendulum spinning around fast in the clockwise direction. Curves starting near the point
(0, 0), where the pendulum is at a stable rest, describe small oscillations of the pendulum.
3
Homework
1 a) Draw the gradient vector field of f (x, y) = sin(x2 − y 2).
b) Draw the gradient vector field of f (x, y) = (x − 1)2 + (y − 2)2.
Hint: In both cases, draw a contour map of f and use gradients
to draw the vector field F (x, y) = ∇f .
2 The vector field F~ (x, y) = hx/(x2 + y 2)(3/2), y/(x2 + y 2)(3/2)i
appears in electrostatics. Find a function f (x, y) such that F~ =
∇f .
~ (x, y) = hP (x, y), Q(x, y)i = hxy, x2i a
3 a) Is the vector field F
gradient field?
b) Is the vector field F~ (x, y) = hP (x, y), Q(x, y)i = hsin(x) +
y, cos(y) + xi a gradient field?
In both cases, find f (x, y) satisfying ∇f (x, y) = F~ (x, y) or give
a reason, why it does not exist.
4 Which of the following vector fields F~ = hP, Qi can be written as
F~ = hP, Qi = hfx, fy i? Make use of Clairaut’s identity Qx = Py ,
to see whether f exists. If yes, find f .
a) F~ (x, y) = hx11, y 9i.
b) F~ (x, y) = hy 9, x7i.
c) F~ (x, y) = h10y + 10x, 10x + 10yi.
d) F~ (x, y) = h9 − y 2 + 4x3y 3, −2xy + 3x4y 2i.
4
Math S21a: Multivariable calculus Oliver Knill, Summer 2017
1 Let C : t 7→ ~r(t) = hcos(t), sin(t)i be a circle parametrized by t ∈ [0, 2π] and let F~ (x, y) =
R
~
h−y, xi. Calculate the line integral I = C F~ (~r) · dr.
R 2π R
Solution: We have I = 0 F~ (~r(t)) · ~r ′ (t) dt = 0 (− sin(t), cos(t)) · (− sin(t), cos(t)) dt =
2π
R 2π 2
0
sin (t) + cos2 (t) dt = 2π
2 Let ~r(t) be a curve given in polar coordinates as ~r(t) = cos(t), φ(t) = t defined on [0, π]. Let
R
F~ be the vector field F~ (x, y) = (−xy, 0). Calculate the line integral C F~ · dr. ~ Solution: In
2
Cartesian coordinates, the curve is r(t) = (cos (t), cos(t) sin(t)). The velocity vector is then
~r ′ (t) = h−2 sin(t) cos(t), − sin2 (t) + cos2 (t)) = (x(t), y(t)i. The line integral is
Z π Z π
~ ′
F (~r(t)) · ~r (t) dt = (cos3 (t) sin(t), 0) · (−2 sin(t) cos(t), − sin2 (t) + cos2 (t)) dt
0
Z π 0
The first generalization of the fundamental theorem of calculus to higher dimensions is the fun-
damental theorem of line integrals.
In other words, the line integral is the potential difference between the end points ~r(b) and ~r(a),
if F~ is a gradient field.
1
3 Let f (x, y, z) be the temperature distribution in a room and let ~r(t) the path of a fly in the
room, then f (~r(t)) is the temperature, the fly experiences at the point ~r(t) at time t. The
change of temperature for the fly is dtd f (~r(t)). The line-integral of the temperature gradient
∇f along the path of the fly coincides with the temperature difference between the end point
and initial point.
4 If ~r(t) is parallel to the level curve of f , then d/dtf (~r(t)) = 0 and ~r ′ (t) orthogonal to
∇f (~r(t)).
5 If ~r(t) is orthogonal to the level curve, then |d/dtf (~r(t))| = |∇f ||~r ′ (t)| and ~r ′ (t) is parallel
to ∇f (~r(t)).
The proof of the fundamental theorem uses the chain rule in the second equality and the funda-
mental theorem of calculus in the third equality of the following identities:
Z b Z b Z b
d
F~ (~r(t)) · ~r (t) dt =
′ ′
∇f (~r(t)) · ~r (t) dt = f (~r(t)) dt = f (~r(b)) − f (~r(a)) .
a a a dt
For a gradient field, the line-integral along any closed curve is zero.
When is a vector field a gradient field? F~ (x, y) = ∇f (x, y) implies Py (x, y) = Qx (x, y). If this
does not hold at some point, F~ = hP, Qi is no gradient field. This is called the component test
or Clairot test. We will see later that the condition curl(F~ ) = Qx − Py = 0 implies that the field
is conservative, if the region satisfies a certain property.
2
Remarks: The fundamental theorem of line integrals works in any dimension. You can formulate
and check it yourself. The reason is that curves, vector fields, chain rule and integration along
curves are easy to generalize to any dimensions. We will see next week that if R is a region
“without holes” then F~ is a gradient field if and only if curl(F~ ) = 0 everywhere in R. A region
R is called simply connected, if every curve in R can be contracted to a point in a continuous
way and every two points can be connected by a path. A disc is an example of a simply connected
region, an annular region is an example which is not. Any region with a hole is not simply
connected. For simply connected regions, the existence of a gradient field is equivalent to the field
having curl zero everywhere.
A device which implements a non gradient force field
is called a perpetual motion machine. It realizes a
force field for which the energy gain is positive along
some closed loop. The first law of thermodynamics
forbids the existence of such a machine. It is infor-
mative to contemplate some of the ideas people have
come up and to see why they don’t work. Here is an
example: consider a O-shaped pipe which is filled only
on the right side with water. A wooden ball falls on
the right hand side in the air and moves up in the water.
Why are there no ”perpetual motion machines”. Ben-
jamin Peirce refers in his book ”A system of analytic me-
chanics” of 1855 to the ”antropic principle”: ”Such a
series of motions would receive the technical name of a
’perpetual motion’ by which is to be understood, that of
a system which would constantly return to the same po-
sition, with an increase of power, unless a portion of the
power were drawn off in some way and appropriated, if it
were desired, to some species of work. A constitution of
the fixed forces, such as that here supposed and in which
a perpetual motion would possible, may not, perhaps, be
incompatible with the unbounded power of the Creator;
but, if it had been introduced into nature, it would have
proved destructive to human belief, in the spiritual ori-
gin of force, and the necessity of a First Cause superior
to matter, and would have subjected the grand plans of
Divine benevolence to the will and caprice of man”.
Here is a futile attempt featured in Youtube.
3
Non-conservative
fields can also be
generated by optical
illusion as M.C. Es-
cher did. The illusion
suggests the existence
of a force field which
is not conservative.
Can you figure out
how Escher’s pictures
”work”?
Homework
1 Let C be the space curve ~r(t) = hcos(t), sin(t), ti for t ∈ [0, π/4]
R
~
and let F (x, y, z) = hy, x, 15i. Calculate the line integral C F~ ·
~
dr.
2 What is the work done by moving in the force field F ~ (x, y) =
h2x3 + 1, 2y 4i along the parabola y = x2 from (−1, 1) to (1, 1)?
a) compute directly b) use the theorem.
3 Let F~ be the vector field F~ (x, y) = h−y, xi/2. Compute the line
integral of F along an ellipse ~r(t) = ha cos(t), b sin(t)i with width
2a and height 2b. The result should depend on a and b.
4 After summer school, you relax in a Jacuzzi and move along curve
C given by part of the curve x20 + y 20 = 1 in the first quadrant,
oriented counter clockwise. The hot water in the tub has the
velocity F~ (x, y) = hx + y, y 4 + xi. Calculate the line integral
R
~ ~
C F · dr, the energy you gain from the fluid force when dislocating
from (1, 0) to (0, 1). Be lazy.
5 Find a closed curve C : ~r(t) for which the vector field