Introduction to Vector Calculus

Similar documents
Math Vector Calculus II

Section 4.1. the k-th coordinate function of f. We call the system of equations . =.

1.1 Single Variable Calculus versus Multivariable Calculus Rectangular Coordinate Systems... 4

Math 461 Homework 8. Paul Hacking. November 27, 2018

Archive of Calculus IV Questions Noel Brady Department of Mathematics University of Oklahoma

Practice Midterm Exam 1. Instructions. You have 60 minutes. No calculators allowed. Show all your work in order to receive full credit.

Math 461 Homework 8 Paul Hacking November 27, 2018

Practice Final Solutions

Vector Calculus handout

Multivariable Calculus : Change of Variables in Multiple Integrals

SOLUTIONS TO SECOND PRACTICE EXAM Math 21a, Spring 2003

Math 234 Exam 3 Review Sheet

1 + f 2 x + f 2 y dy dx, where f(x, y) = 2 + 3x + 4y, is

MATH 52 FINAL EXAM SOLUTIONS

MATH 12 CLASS 23 NOTES, NOV Contents 1. Change of variables: the Jacobian 1

3. Interpret the graph of x = 1 in the contexts of (a) a number line (b) 2-space (c) 3-space

Review Sheet for the Final

Geometry and Motion Selected answers to Sections A and C Dwight Barkley 2016

14.1. Multiple Integration. Iterated Integrals and Area in the Plane. Iterated Integrals. Iterated Integrals. MAC2313 Calculus III - Chapter 14

Vector Calculus. Dr. D. Sukumar

Name: SOLUTIONS Date: 11/9/2017. M20550 Calculus III Tutorial Worksheet 8

7a3 2. (c) πa 3 (d) πa 3 (e) πa3

Arc Length. Philippe B. Laval. Today KSU. Philippe B. Laval (KSU) Arc Length Today 1 / 12

Substitutions in Multiple Integrals

Notes on multivariable calculus

Exercises for Multivariable Differential Calculus XM521

Math 1272 Solutions for Fall 2005 Final Exam

1 The Differential Geometry of Surfaces

18.02 Multivariable Calculus Fall 2007

1. (a) The volume of a piece of cake, with radius r, height h and angle θ, is given by the formula: [Yes! It s a piece of cake.]

PRACTICE PROBLEMS. Please let me know if you find any mistakes in the text so that i can fix them. 1. Mixed partial derivatives.

Vector Calculus, Maths II

Math 107. Rumbos Fall Solutions to Review Problems for Exam 3

In general, the formula is S f ds = D f(φ(u, v)) Φ u Φ v da. To compute surface area, we choose f = 1. We compute

Math 212-Lecture Integration in cylindrical and spherical coordinates

WORKSHEET #13 MATH 1260 FALL 2014

Change of Variables, Parametrizations, Surface Integrals

Distance Formula in 3-D Given any two points P 1 (x 1, y 1, z 1 ) and P 2 (x 2, y 2, z 2 ) the distance between them is ( ) ( ) ( )

Practice Problems for Exam 3 (Solutions) 1. Let F(x, y) = xyi+(y 3x)j, and let C be the curve r(t) = ti+(3t t 2 )j for 0 t 2. Compute F dr.

(You may need to make a sin / cos-type trigonometric substitution.) Solution.

AP Calculus (BC) Chapter 10 Test No Calculator Section. Name: Date: Period:

Contents. MATH 32B-2 (18W) (L) G. Liu / (TA) A. Zhou Calculus of Several Variables. 1 Multiple Integrals 3. 2 Vector Fields 9

SOME PROBLEMS YOU SHOULD BE ABLE TO DO

MATHEMATICS: SPECIALIST UNITS 3C AND 3D FORMULA SHEET 2015

THE DIFFERENTIAL GEOMETRY OF PARAMETRIC PRIMITIVES

Sections minutes. 5 to 10 problems, similar to homework problems. No calculators, no notes, no books, no phones. No green book needed.

n=0 ( 1)n /(n + 1) converges, but not

Multivariable Calculus Lecture #11 Notes

3. On the grid below, sketch and label graphs of the following functions: y = sin x, y = cos x, and y = sin(x π/2). π/2 π 3π/2 2π 5π/2

Exercise. Exercise 1.1. MA112 Section : Prepared by Dr.Archara Pacheenburawana 1

Math 23b Practice Final Summer 2011

Multiple Choice. Compute the Jacobian, (u, v), of the coordinate transformation x = u2 v 4, y = uv. (a) 2u 2 + 4v 4 (b) xu yv (c) 3u 2 + 7v 6

Study guide for Exam 1. by William H. Meeks III October 26, 2012

Section 14.1 Vector Functions and Space Curves

DO NOT BEGIN THIS TEST UNTIL INSTRUCTED TO START

Math Exam IV - Fall 2011

Math 20C Homework 2 Partial Solutions

Look out for typos! Homework 1: Review of Calc 1 and 2. Problem 1. Sketch the graphs of the following functions:

The Jacobian. v y u y v (see exercise 46). The total derivative is also known as the Jacobian Matrix of the transformation T (u; v) :

Topic 5.6: Surfaces and Surface Elements

MTHE 227 Problem Set 10 Solutions. (1 y2 +z 2., 0, 0), y 2 + z 2 < 4 0, Otherwise.

Overview of vector calculus. Coordinate systems in space. Distance formula. (Sec. 12.1)

REVIEW 2, MATH 3020 AND MATH 3030

G G. G. x = u cos v, y = f(u), z = u sin v. H. x = u + v, y = v, z = u v. 1 + g 2 x + g 2 y du dv

A different parametric curve ( t, t 2 ) traces the same curve, but this time the par-

Name. FINAL- Math 13

Ideas from Vector Calculus Kurt Bryan

SOLUTIONS TO THE FINAL EXAM. December 14, 2010, 9:00am-12:00 (3 hours)

Derivatives and Integrals

vand v 3. Find the area of a parallelogram that has the given vectors as adjacent sides.

(a) The points (3, 1, 2) and ( 1, 3, 4) are the endpoints of a diameter of a sphere.

The Change-of-Variables Formula for Double Integrals

Solution. This is a routine application of the chain rule.

EXAM 2 ANSWERS AND SOLUTIONS, MATH 233 WEDNESDAY, OCTOBER 18, 2000

10.1 Review of Parametric Equations

Name: ID: Math 233 Exam 1. Page 1

Practice problems. 1. Given a = 3i 2j and b = 2i + j. Write c = i + j in terms of a and b.

MATH Final Review

MATH 255 Applied Honors Calculus III Winter Midterm 1 Review Solutions

CALCULUS 3 February 6, st TEST

MATH 332: Vector Analysis Summer 2005 Homework

Created by T. Madas SURFACE INTEGRALS. Created by T. Madas

Week 3: Differential Geometry of Curves

Name: ID: Math 233 Exam 1. Page 1

Exam 1 Review SOLUTIONS

Chapter 14: Vector Calculus

Map Projections. Jim Emery. Revision 7/21/ The Earth Ellipsoid 2. 2 The Fundamental Quadratic Forms of a Surface 6

ACCESS TO SCIENCE, ENGINEERING AND AGRICULTURE: MATHEMATICS 1 MATH00030 SEMESTER / Functions

Workbook. MAT 397: Calculus III

The Calculus of Vec- tors

MIDTERM EXAMINATION. Spring MTH301- Calculus II (Session - 3)

Math 5378, Differential Geometry Solutions to practice questions for Test 2

Faculty of Engineering, Mathematics and Science School of Mathematics

MATH 32A: MIDTERM 2 REVIEW. sin 2 u du z(t) = sin 2 t + cos 2 2

Solutions for the Practice Final - Math 23B, 2016

A1. Let r > 0 be constant. In this problem you will evaluate the following integral in two different ways: r r 2 x 2 dx

10.1 Curves Defined by Parametric Equation

Fundamental Theorems of Vector

Definition 1.1 Let a and b be numbers, a smaller than b. Then the set of all numbers between a and b :

f(p i )Area(T i ) F ( r(u, w) ) (r u r w ) da

Transcription:

Introduction to Vector Calculus Phillip Zenor Department of Mathematics Auburn University Edward E. Slaminka Department of Mathematics Auburn University Donald Thaxton Department of Physics Auburn University October 15, 2008 Copyright c 1992

4.4 Arc Length for Curves in Other Coordinate Systems 39 4.4 Arc Length for Curves in Other Coordinate Systems Arc Length for Polar Coordinates Suppose that u(t) (r(t), θ(t)), a t b, is a parametrization of a path in the plane given by polar coordinates. Then (x(t), y(t)) s(t) P ( u(t)) is a parametrization for the path in rectangular coordinates. The derivative is given by d d s(t) (x(t), y(t)) dt dt cos(θ(t)) sin(θ(t)) r(t) sin(θ(t)) r (t). r(t) cos(θ(t)) θ (t) Thus x (t) y (t) s (t) r (t) cos(θ(t)) θ (t)r(t) sin(θ(t)) r (t) sin(θ(t)) + θ (t)r(t) cos(θ(t)). We now have s (t) (x (t)) 2 + (y (t)) 2 ( ((r (t) cos(θ(t)) θ (t)r(t) sin(θ(t))) 2 + (r (t) sin(θ(t)) + θ (t)r(t) cos(θ(t))) 2) 1/2. Suppressing t and expanding, we obtain s ( r 2 cos 2 θ + θ 2 r 2 sin 2 θ 2r θ r cos θ sin θ + r 2 sin 2 θ + θ 2 r 2 cos 2 θ + 2r θ r sin θ cos θ ) 1/2 ( r 2 (cos 2 θ + sin 2 θ) + θ 2 r 2 (sin 2 θ + cos 2 θ) ) 1/2 (r (t)) 2 + (θ (t)r(t)) 2. Thus the length of the curve, where u(a) A, u(b) B, and the image of u is C, is given by L B A C ds b a (r (t)) 2 + (θ (t)r(t)) 2 dt. (1) EXAMPLE 1: Let C be the graph with polar equation r k, for 0 θ π. Thus C is a vertical line in rθ space. The graph P (C) is the top half of the circle with radius k. (See Figure 1.) Figure 1. The graph P (C) is the top half of the circle with radius k.

40 A parametrization for C in polar coordinates is u(t) (k, t), 0 t π. The arc length is given by B AC ds, where A P ( u(0)) (k, 0) and B P ( u(π)) ( k, 0). r (t) 0 and θ (t) 1. Using Equation (1), we obtain 1.2 0.8 0.4 1.2 1.4 Figure 2. The graph of the curve in Example 2. L ( k,0) (k,0) C ds π 0 π 0 k 2 dt πk. (r (t)) 2 + (θ (t)r(t)) 2 dt EXAMPLE 2: Let C be the graph with polar equation r e θ, for 0 θ ln 2. A parametrization for C in polar coordinates is u(t) (e t, t), 0 t ln 2. Thus r (t) e t and θ (t) 1. The arc length is given by B AC ds, where A P (u(0)) (1, 0) and B P (u(ln 2)) (2 cos(ln 2), 2 sin(ln 2)). Using Equation (1), we obtain L (2 cos(ln 2),2 sin(ln 2)) (1,0) d s ln 2 0 (r (t)) 2 + (θ (t)r(t)) 2 dt ln 2 0 2 e 2t + e 2t dt ln 2 0 e t dt 2 2 2 2. Arc Length for Cylindrical Coordinates In a similar fashion, we can compute arc length in cylindrical coordinates in R 3. If u(t) (r(t), θ(t), z(t)) is a parametrization of a path in R 3, then (x(t), y(t), z(t)) s(t) C( u(t)) is a parametrization of the path in rectangular coordinates. The derivative is given by d d s(t) (x(t), y(t), z(t)) dt dt cos(θ(t)) r(t) sin(θ(t)) 0 sin(θ(t)) r(t) cos(θ(t)) 0 0 0 1 r (t) θ (t) z (t).

4.4 Arc Length for Curves in Other Coordinate Systems 41 Thus x (t) y (t) z (t) s (t) Suppressing the variable t, s x 2 + y 2 + z 2 r (t) cos(θ(t)) θ (t)r(t) sin(θ(t)) r (t) sin(θ(t)) + θ (t)r(t) cos(θ(t)) z (t) (r cos θ θ r sin θ) 2 + (r sin θ + θ r cos θ) 2 + z 2 r 2 2 cos 2 θ + θ 2 r 2 sin 2 θ 2r θ r cos θ sin θ+. r 2 sin 2 θ + θ 2 r 2 cos 2 θ + 2r θ r sin θ cos θ + z 2 r 2 (cos 2 θ + sin 2 θ) + θ 2 r 2 (cos 2 θ + sin 2 θ) + z 2 r 2 + r 2 θ 2 + z 2. Thus the length of the curve is given by L B A C dr b a r 2 (t) + r 2 (t)θ 2 + z 2 (t) dt. (2) EXAMPLE 3: Let u(t) (1, 2πt, t), for 0 t 3, be a parametrization of a helix in cylindrical coordinates. Find the length of the helix. See Figure 3. Solution: s(t) C( u(t)). By Equation (2), s (t) r 2 (t) + r 2 (t)θ 2 (t) + z 2 (t) Figure 3. The helix from Example 3. Thus the arc length is given by L (2π) 2 + 1 4π 2 + 1 3 0 4π 2 + 1 dt 3 4π 2 + 1. Arc Length for Spherical Coordinates For spherical coordinates we follow a similar procedure. Let u(t) (ρ(t), φ(t), θ(t)) be a parametrization of a path in R 3 in

42 spherical coordinates. Then (x(t), y(t), z(t)) s(t) S( u(t)) will be a parametrization of the path in rectangular coordinates. The derivative is given by d d s(t) (x(t), y(t), z(t)) dt dt cos θ(t) sin φ(t) ρ(t) cos θ(t) cos φ(t) ρ(t) sin θ(t) sin φ(t) sin θ(t) sin φ(t) ρ(t) sin θ(t) cos φ(t) ρ(t) cos θ(t) sin φ(t) cos φ(t) ρ(t) sin φ(t) 0 ρ (t) φ (t) θ (t). Suppressing the variable t, x ρ cos θ sin φ + φ ρ cos θ cos φ θ ρ sin θ sin φ y z s ρ sin θ sin φ + φ ρ sin θ cos φ + θ ρ cos θ sin φ ρ cos φ φ ρ sin φ. We now have s x 2 + y 2 + z 2 ( (ρ cos θ sin φ + φ ρ cos θ cos φ θ ρ sin θ sin φ) 2 +(ρ sin θ sin φ + φ ρ sin θ cos φ + θ ρ cos φ sin φ) 2 +(ρ cos φ θ ρ sin φ) 2) ρ 2 + ρ 2 φ 2 + ρ 2 θ 2 sin 2 φ. EXERCISES 4.4 By expanding and simplifying, we obtain L B A C dr In Exercises 1 5, find the length, in xy space, of the graphs with the given polar equation. 1. r e 2θ, 0 θ ln 3. 2. r 2 cos(θ), 0 θ 2π. 3. r 3 sec(θ), 0 θ π 4. b a ρ 2 (t) + ρ 2 (t)φ 2 (t) + ρ 2 (t)θ 2 (t) sin 2 φ(t) dt. 4. r cos 2 (θ/2), 0 θ 2π. 5. r sin 2 (θ/2), 0 θ 2π. 6. Let h(t) (ρ(t), φ(t), θ(t)) describe, in spherical coordinates, the location of a particle at time t. Show that the magnitude of the particle s velocity in xyz space is given by ρ 2 (t) + ρ 2 (t)φ 2 (t) + ρ 2 (t)θ 2 (t) sin 2 φ(t).

4.5 Change of Area with Linear Transformations 43 7. Find the arc length of the helix given (in xyz coordinates) by r(t) (cos(2πt), sin(2πt), t), 0 t 4, using both the xyz coordinate version of the arc length integral and the cylindrical coordinate version. 8. The path parametrized by cos(2πt) sin ( ) tπ 4 r(t) sin(2πt) sin ( ) tπ 4 cos ( ), 0 t 4. tπ 4 9. Let r(t) (1, 2πt, t), 0 t 4, represent the helix of Exercise 7 in cylindrical coordinates. a. Show that the helix can be represented as ( t2 ( ) ) 1 u(t) + 1, Arcsin, 2πt, t2 + 1 0 t 4, in spherical ρφθ coordinates. b. Compute the arc length in spherical coordinates, and compare it to Exercise 7. in xyz coordinates describes a spiral along the unit sphere centered at the origin, starting at the north pole and descending to the south pole. Compute the arc length of this path using both xyz coordinates and spherical coordinates. 4.5 Change of Area with Linear Transformations In this section, we introduce the idea of the rate that a linear transformation changes area or volume. First, it is helpful to recall some facts about areas of parallelograms and volumes of parallelepipeds. Suppose that the vectors A and B are drawn emanating from a common point in R 3 forming adjacent edges of a parallelogram P. Then A B is the area of P. If the vectors A (a 1, b 1 ) and B (b 1, b 2 ) are drawn emanating from the origin in R 2 forming adjacent edges of a parallelogram P, then (a 1, a 2, 0) (b 1, b 2, 0) a 1 b 2 a 2 b 1 is the area of P. ( ) det a1 a 2 a b 1 b 1 b 2 a 2 b 1 (which is the area of P). 2 If the vectors A (a 1, a 2, a 3 ), B (b 1, b 2, b 3 ) and C (c 1, c 2, c 3 ) are drawn emanating from the origin and they do not lie in a common plane (they are not co-planer), then they form the adjacent edges of a parallelepiped, P. The magnitude of

44 their triple product is the volume of P. Computationally, the volume of P is A ( B C) a 1 a 2 a 3 det b 1 b 2 b 3 c 1 c 2 c 3. The following theorem is be useful in this and subsequent sections. v Theorem 1 ( a1 a det 2 b 1 b 2 ) ( a1 b det 1 a 2 b 2 ) and det a 1 a 2 a 3 b 1 b 2 b 3 c 1 c 2 c 3 det a 1 b 1 c 1 a 2 b 2 c 2 a 3 b 3 c 3. y T(u, v) ua + vb u While this theorem is easily proven with direct computation, it seems rather remarkable that the parallelepiped determined by the vectors (a 1, a 2, a 3 ), (b 1, b 2, b 3 ), and (c 1, c 2, c 3 ) has the same volume as the parallelepiped determined by (a 1, b 1, c 1 ), (a 2, b 2, c 2 ), and (a 3, b 3, c 3 ). Let T be defined by A B Figure 1. The image of the unit square 0 u 1, 0 v 1 is the parallelogram with adjacent edges A and B drawn emanating from the origin. x T (u, v) u A + v B, where A and B are vectors in R 3. If A and B are drawn emanating from the origin, they are adjacent sides of a parallelogram P with area A B. Let R be the unit square in uv space with adjacent sides the vectors (1, 0) and (0, 1). The area of R is one square unit, and T (R) P, which has an area of A B. (See Figure 1.) Thus T will take one square unit of area onto a parallelogram having area A B. It turns out that A B can properly be thought of as the rate that T changes area. If C is a set in uv space, then the area of T (C) is Area(C) A B. If A (a 1, a 2 ) and B (b 1, b 2 ) are in R 2, then (a 1, a 2, 0) (b 1, b 2, 0) a 1 b 2 a 2 b 1 detdt gives the area of the parallelogram determined by A and B. Thus, if T is a linear transformation from R 2 into R 2, then detd( T ) is the rate that T changes area.

4.5 Change of Area with Linear Transformations 45 In the same fashion, if T is a linear transformation from R 3 into R 3, then det(dt ) is the rate that T changes volume. To recapitulate: If T (r) r A is a linear transformation from R into R n, then A is the rate that T changes length. If T ( r) A T r is a linear transformation from R 2 into R 2, then deta T is the rate that T changes area. v If T ( r) A T r is a linear transformation from R 2 into R 3, then T (1, 0) T (0, 1) is the rate that T changes area. u If T ( r) A T r is a linear transformation from R 3 into R 3, then deta T is the rate that T changes volume. T(u, v) ua + vb EXAMPLE 1: Let A (a, 0) and B (0, b), where a and b are positive numbers. Let T be defined by T (u, v) u A + v B (au, bv). y (0, b) (a, 0) Let C be the unit circle u 2 + v 2 1. Then T (C) is the ellipse with equation ( x 2 ( a) + y ) 2 b 1. See Figure 2. The area bounded by C is π and the area bounded by T (C) is A B π abπ. Thus the area bounded by the ellipse ( x 2 ( a) + y ) 2 b 1 is abπ. This agrees with the value obtained by integration. Figure 2. The image of the unit disc u 2 + v 2 1 is the ellipse ` x 2 + ` y 2 1 a b x EXAMPLE 2: Let C be the unit circle in uv space. 2u + 3v (a) Let T (u, v) u u + v. Then T (1, 0) (2, 1, 1) and T (0, 1) (3, 0, 1). The rate that T changes area is (2, 1, 1) (3, 0, 1) 11. The area of T (C) is π 11. (b) Let T (u, v) ( 2u + 3v u ). The rate that T changes area is ( det DT ) ( ) 2 3 det 3. 1 0 Since the area of the unit circle is π, the area of T (C) is 3π.

46 u w T(u,v,w)(u, 2v,3w) v EXAMPLE 3: Let S denote the surface bounded between the graphs of v u 2 and v u. Let T (u, v) (2u + v, v u). Find the area of the image of S. Solution: The area of S is 1 0 u u2 du 1 6. The rate that T ( ) changes area is 2 1 det 3. Thus the area of T 1 1 (S) (3) ( ) 1 6 1 2. EXAMPLE 4: Let T be the transformation defined by T (u, v, w) (au, bv, cw). z Then T takes the unit sphere u 2 + v 2 + w 2 1 onto the ellipsoid E with equation ( x 2 ( a) + y ) 2 ( b + z ) 2 c 1. See Figure 3 for the case that a 1, b 2, and c 3. DT a 0 0 0 b 0 0 0 c. x Figure 3. T (u, v, w) (u, 2v, 3w) takes the unit sphere onto the ellipse x 2 + y2 4 + z2 9 1. y Thus T changes volume at a rate of abc, since the unit cube goes onto a box with sides of length a, b, and c. Since the volume bounded by a unit sphere is ( 4 3) π, the volume bounded by the ellipsoid E is 4 3 abcπ. If T is a linear transformation from R 2 into R, then T takes an object with area onto an object with no area. Therefore, the rate that T changes area is zero. This is consistent with the fact that T (1, 0) and T (0, 1) must point in the same or opposite directions (they are either a positive number, a negative number, or zero.) Indeed, if the domain of T is R 2, then: If the range of T is R 1, then the rate that T changes area is zero. If the range of T is R 2 or R 3, then T (0, 1) T (1, 0) 0 if and only if the image of T is a line or the origin. Similarly, if the domain of T is R 3, then: If the range of T is either R or R 2, then the rate that T changes volume is zero. If the range is R 3 then deta T 0 if and only if the image of T is a plane, a line, or the origin.

4.5 Change of Area with Linear Transformations 47 Summary If T is a linear transformation from R 2 into R 2, then the rate that T changes area is det(a T ) det(dt ). If T is a linear transformation from R 3 into R 3, then the rate that T changes volume is det(a T ) det(dt ). If T (u, v) u A + v B is a linear transformation from R 2 into R 3, then the rate that T changes area is A B. If T is a linear transformation from R 3 into R 2 or R, then the rate that T changes volume is 0. If T is a linear transformation from R 2 into R, then the rate that T changes volume is 0. EXERCISES 4.5 In Exercises 1 5, determine the rate that each transformation changes area. 1. T (u, v) (3u, 2v). 2. T (u, v) (u + 3v, v u). 3. T (u, v) (u v, 2v, u + v). 4. T (s, t) (s t, 3t + s, t s). 5. T (u, v) u + v. In Exercises 6 10, determine the rate that each transformation changes volume. 6. T (u, v, w) (2u v, u + w, u + v + w). 7. T (r, s, t) (2r + s t, r s 3t, r + s + t). 8. T (u, v, w) (2u w + v, u + v 22w). 9. T (u, v, w) (u + v + w, u + v + w). 10. T (u, v, w) u + v + w. 11. Find a linear transformation from the uv plane into the xy plane that takes the circle u 2 +v 2 1 onto the ellipse E with equation 4x 2 + 9y 2 36. Find the area bounded by E using the techniques of this section. 12. Find the area bounded by the ellipse 2x 2 +3y 2 5. 13. Find a linear transformation from uvw space onto xyz space that takes the unit sphere centered at the origin onto the ellipsoid E with equation 3x 2 + 4y 2 + 2z 2 1. Find the volume bounded by E. 14. Let R be the rectangle bounded between the lines u 2, u 2, v 0, and v 4 and let T (u, v) ( u + v, u 2v, 2u + 4v). Find the area of T (R). 15. Let R be the rectangle bounded between the lines u 2, u 2, v 2, and v 6 and let T (u, v) ( u + v, 2u 2v). Find the area of T (R).

48 16. Let R be the rectangle bounded between the lines u 2, u 2, v 2, and v 6 and let T (u, v) ( u + v, u 2v). Find the area of T (R). 17. Let R be the region bounded between the graphs of v u 2 and v u 3 and let T (u, v) (u + v, u 2v, 2u + v). Find the area of T (R). 18. Let R be the region bounded between the graphs of u v 2 and u v 3 and let T (u, v) (3u + v, u 2v). Find the area of T (R). 19. Let C be the box bounded between the planes u 2, u 2, v 3, v 4, w 0, and w 10. Let T (u, v, w) ( u + v + 2w, u 2v, 2u + 4v + w). Find the volume of T (C). 20. Let E be the ellipsoid u2 4 + y2 9 + z2 1, and let T (u, v) (u + v + 2w, u 2v, 2u + v + w). Find the volume of T (E). 21. Let V be the solid obtained by rotating the region bounded between the u axis and the graph of v u 2 +1, 1 u 1 about the u axis. Let T (u, v, w) (u + v + 2w, u 2v, 2u + v + w). Find the volume of T (E). f(t) f(t+h) f(t+h)-f(t) 4.6 The Jacobian When a function from R n into R m is not linear, then the rate that the function changes area or volume becomes a local property. In this section, we learn how a function changes area or volume at a point. This idea is not really new. Let f be a function from R into R n. Recall the geometric interpretation of the derivative, in which f (t) represents the rate that f changes arc length at t. To see that this is a reasonable interpretation, let [t, t + h] be an interval in the domain of f. (See Figure 1.) If h is small, then f(t + h) f(t) is an approximation of the length of f([t, t + h]), and Figure 1. f(t + h) f(t) approximates the arc length from f(t) to f(t + h). f (t) lim h0 length of f([t, t + h]) length of [t, t + h] lim h0 f(t + h) f(t) h. The rate that f changes length at t is called the Jacobian of f at t. We denote the Jacobian of f at t by J f(t). Notice that if g is a function from a subset of R n into R and if r is a parametrization for a curve C with endpoints A and B and domain [a, b], then B A C g d r b a g( r(t))j r(t) dt. The Jacobian can also be defined for functions with domain in R n

4.6 The Jacobian 49 for any positive integer n. However, we define it only for functions with domain in R 2 or R 3. Let D be a subset of R 2 (uv space), and let f be a function from D into R 2 or R 3. We want the Jacobian of f at (u, v) to be the rate that f changes area at (u, v). Let S h be the square having sides of length h with vertices (u, v), (u + h, v), (u, v + h) and (u + h, v + h). Assume that S h is a subset of D. We approximate the surface f(s h ) with the parallelogram P h with the adjacent sides the vectors f(u + h, v) f(u, v) and f(u, v + h) f(u, v), drawn emanating from f(u, v). See Figure 2. f(u+h,v) (u, v+h) f (u,v) (u+h,v) f(u,v) f(u,v+h) Figure 2. The area of f(s h ) is approximated by the area of P h. The area of P h is ( f(u+h, v) f(u, v)) ( f(u, v +h) f(u, v)). The rate that f changes area at (u, v) is given by area of P h ( f(u lim lim + h, v) f(u, v)) ( f(u, v + h) f(u, v)) h0 area of S h h0 h 2 f(u + h, v) f(u, lim v) h0 h f(u, v) u f(u, v) v. f(u, v + h) f(u, v) h Definition: The Jacobian for Functions from R 2 to R 2 or R 3 If f is a differentiable function from R 2 into R 2 or R 3, we define the Jacobian of f, denoted by f(u, v), to be

50 f(u, v) u f(u, v) v. It represents the rate that f changes area at (u, v). Of course, if f is a function from a subset of R 2 into R, then J f(u, v) 0. EXAMPLE 1: Let P be the polar transformation from rθ space defined by P (r, θ) (r cos θ, r sin θ). The partial derivatives of P are P(r, θ+h) (r, θ+h) (r, θ) (r+h, θ) P(r, θ) P P(r+h, θ) Figure 3. The square S in rθ space. The sector P (S) in xy space. P r (cos θ, sin θ) and P θ ( r sin θ, r cos θ). Thus the Jacobian of P is given by JP P (r, θ) r P θ r, and is the rate that P changes area at (r, θ). This fits the geometry of the function very well. Recall that the area of a sector of a circle of radius R is R2 ( θ) 2, where θ is the angle of the sector. If we consider the square S in rθ space with sides of length h and one vertex (r, θ) as in Figure 3, then P (S) is the portion of a sector as in Figure 3. The area of S is h 2, and the area of P (S) is (r+h)2 h (2rh+h 2 )h 2. Thus, [area of P (S)]/[ area of S] (2rh+h2 )h/2 h 2 2 r2 h 2 r + h 2. It follows that as h gets close to 0, then [area of P (S)]/[area of S] gets close to r. EXAMPLE 2: Find the Jacobian of h(u, v) (u, u 2 cos(v), u 2 sin(v)). Solution: h (u, v) (1, 2u cos(v), 2u sin(v)). u h v (u, v) ( 0, u 2 sin(v), u 2 cos(v) ). h u 2u h 3 cos 2 (v) + 2u 3 sin 2 (v) u 2 cos(v) v u 2 sin(v) 2u 3 u 2 cos(v) u 2 sin(v)

4.6 The Jacobian 51 J h(u, v) h u h v 4u 6 + u 4. Let D be a subset of R 3, and let f be a function from D into R 3. We are interested in how f changes volume at a point. We proceed exactly as we did with areas and functions from R 2 into R 2. Let B h be a square box with three adjacent edges the vectors [(u + h, v, w) (u, v, w)], [(u, v + h, w) (u, v, w)], and [(u, v, w + h) (u, v, w)]. Let P h be the parallelepiped with three adjacent edges [ f(u+h, v, w) f(u, v, w)], [ f(u, v+h, w) f(u, v, w)], and [ f(u, v, w+ h) f(u, v, w)]. We now approximate the volume of f(b h ) with the volume of P h (see Figure 4). (u, v, w+h) f(u, v+h, w) f(u+h, v, w) f(u, v, w) (u, v, w) f (u+h, v, w) (u, v+h, w) f(u, v, w+h) Figure 4. The volume of f(b h ) is approximated by the volume of a parallelepiped. The rate that f changes volume is given by volume of P h lim h0 volume of B h ( ) [( ) ( )] f(u + h, v, w) f(u, v, w) f(u, v + h, w) f(u, v, w) f(u, v, w + h) f(u, v, w) lim h0 h 3 ( ) ( ) ( ) f(u + h, v, w) f(u, v, w) f(u, v + h, w) f(u, v, w) f(u, v, w + h) f(u, v, w) lim h0 h h h f u ( f v f ) w.

52 Definition: The Jacobian for Functions from R 3 to R 3 If f: R 3 R 3 is a differentiable function, we define the Jacobian of f, denoted by Jf(u, v, w), to be ( f u f v f ) w. It represents the rate that f changes volume at (u, v, w). EXAMPLE 3: Let C z be the cylindrical transformation from rθz space into xyz space defined by C z (r, θ, z) (r cos θ, r sin θ, z). Thus and so C z (r, θ, z) r C z (r, θ, z) θ (cos θ, sin θ, 0), ( r sin θ, r cos θ, 0), C z (r, θ, z) (0, 0, 1), z JC z (r, θ, z) ( C z r C ) z θ c z r. EXAMPLE 4: Let S be the spherical transformation defined by ρ cos(θ) sin(φ) S(ρ, φ, θ) ρ sin(θ) sin(φ). ρ cos(φ) It is left as an exercise to show that J S(ρ, φ, θ) ρ 2 sin(φ). The following theorem follows from direct computation, and its proof is left as an exercise (Exercises 27 and 28). Theorem 1 If F is a differentiable function from a subset of R 2 into R 2 or from a subset of R 3 into R 3, and r is a point in R 2 or R 3, then JF ( r) det D( F ) r. If F is a differentiable function and G( r) F ( r) + r 0, then J F J G.

4.6 The Jacobian 53 The second part of Theorem 1 is not terribly surprising. A rigid motion such as a translation does not change area or volume. The next example leads to Theorem 2. EXAMPLE 5: Let T 1 (u, v) ( ) 2x 3y 1. Then x + y + 3 ( 5u + v v ), and let T 2 (x, y) T 2 T 1 (u, v) J( T 2 T 1 )(u, v) ( 2(5u + v) 3( v) 1 (5u + v) + ( v) + 3 ( ) 10u + 5v 1. 5u + 3 detdt 2 T 1 (u, v) ( ) 10 5 det 25. 5 0 ) Also, and Observe that ( ) JT 1 (u, v) 5 1 det 5 0 1 ( ) JT 2 (x, y) 2 3 det 5. 1 1 J( T 2 T 1 )(u, v) JT 2 (x, y)jt 2 (u, v). The observation in Example 5 is no coincidence. In general, the following theorem tells us that if f g is defined, then the rate that f g changes area (or volume) at r is the rate that g changes area (or volume) times the rate that f changes area (or volume) at g( r). Theorem 2 Suppose that F and H are functions from R 3 into R 3. Then J( F H)( r) JF ( H( r))j H(r). Suppose that F is a function from R 2 into R 2 and H is a function from R 2 into either R 2 or R 3. Then J( F H)( r) J F ( H( r))j H(r).

54 EXAMPLE 6: The polar function P changes area at a rate of r at the point (r, θ). The linear transformation T (x, y) (2x + 3y, x + y, x) changes area at the rate of 35. Then J( T P )(r, θ) 35 r. If we compute J( T P )(r, θ) directly, we have ( T P )(r, θ) 2r cos(θ) + 3r sin(θ) r cos(θ) + r sin(θ) r cos(θ), and J( T P )(r, θ) r ( T P )(r, θ) θ ( T P )(r, θ) 2 cos(θ) + 3 sin(θ) 2r sin(θ) + 3r cos(θ) cos(θ) + sin(θ) r sin(θ) + r cos(θ) cos(θ) r sin(θ) ( r, 3r, 5r) 35r 2 35 r. EXERCISES 4.6 In Exercises 1 4, find the Jacobian of the polar transformation P at the given point. (The angles are measured in radians.) 1. (2, 0). 2. (2, 2π). 3. (2, 2π). 4. ( 2, π). In Exercises 5 8, find the Jacobian of the cylindrical transformation C z at the given point. (The angles are measured in radians.) 5. (2, 0, 2). 6. (2, 2π, 2). 7. (2, 2π, 2). 8. ( 2, π, 1). In Exercises 9 12, find the Jacobian of the spherical transformation S(ρ, φ, θ) at the given point. (The angles are measured in radians.) 9. ( ) ( ) 2, 0, π 4. 10. 2, π, π 4. 11. ( 2, π 4, ) ( 7π 4. 12. 2, π 3, ) 7π 4.

4.6 The Jacobian 55 In Exercises 13 14, use Theorem 1 to find the Jacobian for the given function 13. S(r, φ, θ) + (2, 3, 4). 14. P (r, θ) + ( 1, 4). In Exercises 15 17, you are given J f and J g, find J( f g) at r 0. 15. J g(u, v) 3u 2v and J f(x, y) x 2 (y + 2) x. r 0 (1, 2) and g( r 0 ) ( 1, 4). 16. J g(u, v) u 2 + v 2 and J f(x, y) x 2 xy. r 0 ( 1, 2) and g( r 0 ) (2, 1). 17. J g(u, v, w) w u 2 + v 2 and Jf(x, y, z) x 2 xy + z. r 0 (1, 1, 2) and g( r 0 ) (2, 1, 1). Use Theorem 2 in Exercises 18 22 to find the Jacobian for f g. P is the polar function and S is the spherical transformation. 18. f P and g(u, v) (2u + 6v + 2, u). 19. f P and g(u, v) ( 2u + 6v + 2, u + v 1). 20. f S and g(u, v, w) ( 2u + 6v + 2w + 1, u + v + 2w 1, u + v 2w). 21. f(u, v) ( 2u + 6v + 2, u + v 1) and g P. 22. f(u, v, w) ( 2u + 6v + 2w + 1, u + v + 2w 1, u + v 2w) and f S. 23. Show that J( S(ρ, φ, θ)) ρ 2 sin(φ). 24. Calculate J( f(u, v)) for f(u, v) (2 cos(u) sin(v), 2 sin(u) sin(v), 2 cos(v)). 25. What is the image of f of Exercise 24? 26. Show that J F J( F + r 0 ). 27. Let F be a differentiable function from a subset of R 2 into R 2. Show that J( F (u, v)) det D(F ) (u,v). 28. Let F be a differentiable function from a subset of R 3 into R 3. Show that J( F (u, v, w)) det D(F ) (u,v,w). 29. Suppose that f is a function from R 3 into R 3 and g is a function from R 2 into R 3. What can be said about the relationship between J( f g) and the product J( f)j( g)?