Maxima/minima with constraints

Similar documents
Laplace equation in polar coordinates

Chapter 4: Partial differentiation

Math 210, Final Exam, Fall 2010 Problem 1 Solution. v cosθ = u. v Since the magnitudes of the vectors are positive, the sign of the dot product will

DEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS. MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 2018

Solutions to Sample Questions for Final Exam

Marion and Thornton. Tyler Shendruk October 1, Hamilton s Principle - Lagrangian and Hamiltonian dynamics.

1MA6 Partial Differentiation and Multiple Integrals: I

Two dimensional oscillator and central forces

Taylor Series and stationary points

Math 210, Final Exam, Spring 2012 Problem 1 Solution. (a) Find an equation of the plane passing through the tips of u, v, and w.

Variational Principles

PHYS 705: Classical Mechanics. Hamiltonian Formulation & Canonical Transformation

Figure 25:Differentials of surface.

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

MATH H53 : Final exam

Notes 19 Gradient and Laplacian

Connection to Laplacian in spherical coordinates (Chapter 13)

Physics 5153 Classical Mechanics. Canonical Transformations-1

Lecture Notes for PHY 405 Classical Mechanics

Tutorial Exercises: Geometric Connections

A = (a + 1) 2 = a 2 + 2a + 1

Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011

Math221: HW# 7 solutions

Analytical Dynamics: Lagrange s Equation and its Application A Brief Introduction

Math 265H: Calculus III Practice Midterm II: Fall 2014

Expansion of 1/r potential in Legendre polynomials

Electromagnetism HW 1 math review

1. If the line l has symmetric equations. = y 3 = z+2 find a vector equation for the line l that contains the point (2, 1, 3) and is parallel to l.

Fermat s Principle. Fermat s Principle states that a ray of light in a medium will follow the path which takes the least amount of time.

Sketchy Notes on Lagrangian and Hamiltonian Mechanics

Constraints. Noninertial coordinate systems

CURRENT MATERIAL: Vector Calculus.

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

McGill University December Intermediate Calculus. Tuesday December 17, 2014 Time: 14:00-17:00

Math Review for Exam Compute the second degree Taylor polynomials about (0, 0) of the following functions: (a) f(x, y) = e 2x 3y.

Curves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,

MATH2321, Calculus III for Science and Engineering, Fall Name (Printed) Print your name, the date, and then sign the exam on the line

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

Introduction to the Calculus of Variations

Calculus III 2004 Summer Practice Final 8/3/2004

Brief introduction to groups and group theory

Classification of algebraic surfaces up to fourth degree. Vinogradova Anna

t 2 + 2t dt = (t + 1) dt + 1 = arctan t x + 6 x(x 3)(x + 2) = A x +

Lecture 5. Alexey Boyarsky. October 21, Legendre transformation and the Hamilton equations of motion

HW 3 Solution Key. Classical Mechanics I HW # 3 Solution Key

AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES

First Order Differential Equations

Review for the Final Test

Jim Lambers MAT 280 Summer Semester Practice Final Exam Solution. dy + xz dz = x(t)y(t) dt. t 3 (4t 3 ) + e t2 (2t) + t 7 (3t 2 ) dt

ESM 3124 Intermediate Dynamics 2012, HW6 Solutions. (1 + f (x) 2 ) We can first write the constraint y = f(x) in the form of a constraint

Engg. Math. II (Unit-IV) Numerical Analysis

the Cartesian coordinate system (which we normally use), in which we characterize points by two coordinates (x, y) and

Calculus of variations - Lecture 11

Introduction of Partial Differential Equations and Boundary Value Problems

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

Charge and current elements

1. For each function, find all of its critical points and then classify each point as a local extremum or saddle point.

ENGI Partial Differentiation Page y f x

Find the indicated derivative. 1) Find y(4) if y = 3 sin x. A) y(4) = 3 cos x B) y(4) = 3 sin x C) y(4) = - 3 cos x D) y(4) = - 3 sin x

McGill University April Calculus 3. Tuesday April 29, 2014 Solutions

Maxima and Minima. (a, b) of R if

Solution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3

Introduction and Vectors Lecture 1

CURRENT MATERIAL: Vector Calculus.

MATH 18.01, FALL PROBLEM SET # 6 SOLUTIONS

Complex Differentials and the Stokes, Goursat and Cauchy Theorems

Figure 21:The polar and Cartesian coordinate systems.

Review for the First Midterm Exam

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

= π + sin π = π + 0 = π, so the object is moving at a speed of π feet per second after π seconds. (c) How far does it go in π seconds?

Module Two: Differential Calculus(continued) synopsis of results and problems (student copy)

Practice Midterm 2 Math 2153

Physics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics

LESSON 25: LAGRANGE MULTIPLIERS OCTOBER 30, 2017

Math 223 Final. July 24, 2014

Math 201 Solutions to Assignment 1. 2ydy = x 2 dx. y = C 1 3 x3

Gradient, Divergence and Curl in Curvilinear Coordinates

Hertz potentials in curvilinear coordinates

Attempt QUESTIONS 1 and 2, and THREE other questions. penalised if you attempt additional questions.

Math 234 Final Exam (with answers) Spring 2017

Electrodynamics PHY712. Lecture 4 Electrostatic potentials and fields. Reference: Chap. 1 & 2 in J. D. Jackson s textbook.

2 Use the factor theorem to show that a + b c is a factor of. (a + b + c) 3 6(a + b + c)(a 2 + b 2 + c 2 )+8(a 3 + b 3 + c 3 ).

SOME PROBLEMS YOU SHOULD BE ABLE TO DO

Physics 3323, Fall 2016 Problem Set 2 due Sep 9, 2016

Errata Instructor s Solutions Manual Introduction to Electrodynamics, 3rd ed Author: David Griffiths Date: September 1, 2004

Lecture 2 : Curvilinear Coordinates

S12.1 SOLUTIONS TO PROBLEMS 12 (ODD NUMBERS)

Mathematical Tripos Part IA Lent Term Example Sheet 1. Calculate its tangent vector dr/du at each point and hence find its total length.

Math 265 (Butler) Practice Midterm III B (Solutions)

Solutions to old Exam 3 problems

Notes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.

MTH4101 CALCULUS II REVISION NOTES. 1. COMPLEX NUMBERS (Thomas Appendix 7 + lecture notes) ax 2 + bx + c = 0. x = b ± b 2 4ac 2a. i = 1.

Problem 1 (Equations with the dependent variable missing) By means of the substitutions. v = dy dt, dv

From quantum to classical statistical mechanics. Polyatomic ideal gas.

1 + x 2 d dx (sec 1 x) =

Review problems for the final exam Calculus III Fall 2003

Integration by Substitution

ENGI Partial Differentiation Page y f x

Curvilinear Coordinates

Math 268: Homework. Assignment 2: Assigned Wed 09/09. Due Wed 09/16

Transcription:

Maxima/minima with constraints Very often we want to find maxima/minima but subject to some constraint Example: A wire is bent to a shape y = 1 x 2. If a string is stretched from the origin to the wire, at what point along the wire is the length of the string minimized? We want to minimize d 2 = x 2 + y 2. We can eliminate the y 2 using y 2 = (1 x 2 ) 2, so we minimize the function f (x) = x 2 + (1 x 2 ) 2 = x 4 x 2 + 1 We find df dx = 4x 3 2x. The minimum occurs at 2x 2 1 = 0 so x = ± 1/2. We also find a local maxima at x = 0. Check with second derivatives!

Another approach... Let s do the same problem, but starting from f (x, y) = x 2 + y 2, and then the differential Or we can write as, df = 2xdx + 2ydy df dx = 2x + 2y dy dx Then we can obtain dy dx y = 1 x 2 from the equation of constraint dy dx = 2x

Another approach continued... Then we substitute into df dx = 0, df = 2x 4xy = 0 dx We can also solve df = 0 since dx is arbitrary, so 2x 4xy = 0 We still get 2x 4x(1 x 2 1 ) = 0 so that x = ± 2 or x = 0 as before

Method of Lagrange multipliers The approaches above work, either by substituting and eliminating a variable, or by finding dy dx However, these approaches can often lead to inconvenient algrebra We note that we can write the constraint φ(x, y) = constant (sometimes we write it so that φ(x, y) = 0) Then we have, for minimization of f (x, y) with constraint φ(x, y) = constant, df = dx + x y dy = 0 Note that this is what we solve since f is really a function of one variable due to the constraint dφ = φ φ dx + x y dy = 0

Lagrange multipliers continued... We can obviously add df + λdφ = 0, where λ is a Lagrange muliplier We will find a λ such that we can solve separate equations, x + λ φ x = 0 y + λ φ y = 0 We can do this because λ is purely arbitrary! We can choose however we want. Back to example above, we have f (x, y) = x 2 + y 2 and φ(x, y) = y + x 2 = 1, we create F (x, y) to minimize F (x, y) = f (x, y) + λφ(x, y) = x 2 + y 2 + λ(y + x 2 )

Example of Lagrange mulitpliers, continued Next we solve the following two equation, x = 2x(1 + λ) = 0 y = 2y + λ = 0 We could have λ = 1, then y = 1 2 and from y = 1 x 2 we get x = ± 1/2 Alternately we could have x = 0, and then y = 1, and λ = 2 Same points as before! Still need to check with second derivatives whether this is a local minima or a maxima

Langrange multipliers in three dimensions and beyond We might want to minimize f (x, y, z) subject to constraint φ(x, y, z) = constant Due to constraint, only two of the variables can be considered independent, for example x,y We construct the function to minimize F (x, y, z) Then df = df + λdφ is df = F (x, y, z) = f (x, y, z) + λφ(x, y, z) ( ) ( ) ( ) x + λ φ dx + x y + λ φ dy + y z + λ φ dz z Then we choose λ such that z + λ φ z = 0

Lagrange multipliers in more dimensions Then, since dx and dy are independent, we can always take either dx = 0 or dy = 0, and then we have x + λ φ x = 0 y + λ φ y = 0

Lagrange multipliers with multiple constraint equations Say we want to minimize f (x, y, z, w) with constraints φ 1 (x, y, z, w) = constant and φ 2 (x, y, z, w) = constant Now we introduce two multipliers λ 1 and λ 2 There are two independent variables, so we can choose, for example, dx = 0 or dy = 0 We minimize then F (x, y, z, w) = f (x, y, z, w) + λ 1 φ 1 (x, y, z, w) + λ 2 φ 2 (x, y, z, w) We wind up with equations to solve, w + λ φ 1 1 w + λ φ 2 2 w = 0 And we can write identical equations for partials with respect to z, and then x and y.

Change of variables Sometimes problems are more easily solved in a different set of variables Best example is Cartesian, polar, cylindrical, spherical, etc. coordinate systems Example: Write the wave equation in terms of r = x + vt and s = x vt 2 F x 2 1 v 2 2 F t 2 = 0 df t = r dr t + s ds t Then dr t = r x dx = dx, and ds t = s x dx = dx x = r + s We can then see x = r + s, and t = v ( r ) s

Change of variables, wave equation 2 F r s = 0 Solutions F = f (s) + g(r) = f (x vt) + g(x + vt) One solution is f (x vt) = e ik(x vt) and g(x + vt) = e ik(x+vt) for right/left traveling waves

Laplace equation in polar coordinates The Laplace equation is given by 2 F x 2 + 2 F y 2 = 0 We have x = r cos θ, y = r sin θ, and also r 2 = x 2 + y 2, tan θ = y/x We have for the partials with respect to x and y, x = r r x + θ θ x Then 2rdr = 2xdx + 2ydy, and d tan θ = (1 + sin 2 θ/ cos 2 θ)dθ = y dx + 1 x 2 x dy The first relation shows r x = x/r = cos θ The second relation shows θ x = sin θ r

Laplace equation in polar coordinates, continued So finally we get for x, and also y x y = cos θ r sin θ r θ = sin θ r + cos θ r θ We can repeat this process, taking y of the above results Finally we obtain Laplace equation in polar coordinates, 1 r r x and ( r ) + 1 2 F r r 2 2 θ = 0

Another example of change of variables, Hamiltonian and Lagrangian Legendre transformations are also a way to change the independent variables In classical mechanics, we can work with L(q, q) or H(q, p) Since L(q, q) depends on independent variables q and q, we can find dl It turns out that dl = L q L dq + q d q is give by, dl = ṗdq + pd q Define H = qp L (A Legendre transformation!) dh = pd q + qdp dl = ṗdq + qdp Therefore H(q, p) (function of independent variables q and p Constructing H(q, p) is the usual starting point for quantum mechanics Chapter 4, section 11, problem 11