Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Size: px
Start display at page:

Download "Lecture 2 Lagrangian formulation of classical mechanics Mechanics"

Transcription

1 Lecture Lagrangian formulation of classical mechanics Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0, the initial positions an velocities r i t 0 ) an ṙ i t 0 ) for all point masses forming the system. Another formulation for the problem coul be to only give the positions of all point masses in the system, at the initial time t an final time t of the motion uner consieration. To see why this is mathematically equivalent, consier first t such that t = t t 0. Specifying r i t ) an r i t ) is then tantamount to specifying r i t ) an ṙ i t ). By continuity, for t sufficiently small, the motion remains unequivocally etermine if one specifies r i t ) an r i t ). As t grows arbitrarily, however, one can fin several possible motions with the same en points. The principle of stationary action, sometimes also calle less accurately the principle of least action, says that among all possible paths from r i t ) an r i t ), the physically realizable paths are the paths that extremize a functional calle the action S efine as follows Sq,..., q N ] = t Lq t),..., q N t), q t,..., q N, t) ) where L is calle the Lagrangian of the system. Note that the functions q i t) are generalize coorinates; they o not necessarily have to be Cartesian coorinates or other stanar coorinates. q i oes not have to have units of meter, an q i / oes not have to have the units of velocity. What the principle of stationary action is The principle of stationary action can be seen as a reformulation of the question of fining the proper evolution of a system in configuration space from a ifferential formulation to a variational/integral formulation. This has several avantages, which we briefly mentione in the last lecture, an which we will highlight as we start to use the variational formulation in practical examples. What the principle of stationary action is not The principle of stationary action oes not contain any new physics. All the physics is in fining appropriate Lagrangians o escribe a system of interest. We will soon learn a simple metho to construct Lagrangians for a large class of physical systems. However, there is no general metho to construct a Lagrangian for every system. In fact, in new fiels of physics, fining an appropriate Lagrangian to escribe the ynamics of a system can be a major accomplishment in its own right. Examples Before we learn how to construct Lagrangians for a certain class of physical systems, let us first see explicit variational formulations for the first two simple examples we looke at in the previous lecture. Newton s apple We saw that Newton s equations for Newton s apple are z = g ) The Newtonian ynamics for the apple can be written in a variational form as follows. A Lagrangian for the free-fall of an object suject to gravity is Lz, z, t) = ) z gz 3) where the mass of the object has been set to m = without any loss of generality. The action for the motion between times t an t thus is t Sz] = Lz, z t ) t, t) = z gz] t

2 As we have seen in our review of the calculus of variation, extremizing the action with fixe enpoints zt ) an zt ) leas to the following Euler-Lagrange equation for zt): Lz, z ), t) = Lz, z, t) 4) We have So Equation 4) can be rewritten as as expecte. The simple penulum A Lagrangian for the simple penulum is Minimizing the action Lz, z z, t) =, Lz, z, t) = g z = g Lθ, θ/, t) = ) θ l + gl cos θ 5) Sθ] = t t l ) θ + gl cos θ] leas to the E-L equations Lθ, θ/, t)) = Lθ, θ/, t) 6) We have Lθ, θ/, t) = l θ, Lθ, θ/, t) = gl sin θ so Equation 6) takes the form l g θ = gl sin θ θ + l sin θ = 0 which agrees with the secon-orer ODE we ha previously obtaine for that system. 3 Arbitrariness in the efinition of the Lagrangian Suppose a Lagrangian L escribes the ynamics of a system of point masses. Then, given the form of the E-L equation L) = L it is clear that the Lagrangian L = AL + B, where A an B are two constants, leas to the same equations of motion for the system. More interestingly, the equations of motion are also recovere if one replaces L with the Lagrangian L efine by L = L + ft, qt)) 7) Inee, if S is the action associate with the Lagrangian L an S is the action associate with the Lagrangian L, then Equation 7) implies t S q] = Sq] + t ft, qt)) = Sq] + ft, qt )) ft, qt )) 8) The principle of stationary action states that the physically relevant trajectory in configuration space is obtaine by extremals of the action holing the initial an final times an positions fixe. Since t, t, qt ) an qt ) are fixe, the extremals of S are the same as those of S, so both L an L can be use to escribe the motion of the system.

3 4 Hamilton s principle 4. Statement of the principle An important class of problems in classical mechanics is the class of problems that involves conservative forces only. Conservative forces are forces that can be expresse as the graient of a scalar function one often calls a potential, that epens only on time an the positions of the particles: F = V t, rt)). For these situations, the construction of a Lagrangian is straightforwar, an given by Hamilton s principle, which we state below: A system of point masses for which the forces are erive from a potential energy that is inepenent of velocity evolves along a path q for which the action Sq] = t t Lt, qt), q ) is stationary with respect to variations of the path q that leave the enpoints fixe, where L = T V is the ifference between the kinetic energy T an the potential energy V. 4. Illustration The examples provie in section are illustrations of Hamilton s principle. Let us stuy one more important example. Newton s law of gravitation tells us that two objects with masses m an M separate by a istance R exert a force with magnitue F G = G mm R on each other, an the force is an attractive force, irecte along the line joining the centers of gravity of the two objects. G N.m.kg is calle the gravitational constant. We will consier the particular case in which M m, so the object of mass M is essentially at rest, an the lighter object, with mass m, orbits aroun the larger object. This is a very goo approximation for the motion of the Earth aroun the sun, or satellites, natural the Moon) or not, orbiting aroun Earth. It easy to show that the motion of the object of mass m is containe in a plane. If we take the center of gravity of the massive object as the origin of our coorinate system, we can write that the force acting on the light object is F G = GmM r r 3 Let us compute the quantity r r ṙ) = r r = r GM r 3 = 0 where the secon equality was obtaine by applying Newton s law. We conclue that Lt) = r ṙ = Cst = L0) = r0) ṙ0) 9) This means that the motion of the object of mass m is at all times in the plane perpenicular to L0). The analysis of the motion of the object can thus be reuce to stuying a system with two egrees of freeom. Let us consier the polar coorinate system r, θ) in the plane perpenicular to L0), with the origin fixe at the center of gravity of the massive object. In this coorinate system, the force F G can be written as F G = V where V = GmM 0) r is the gravitational potential. We can therefore use Hamilton s principle to write a Lagrangian for the object of mass m. The kinetic energy of the object is T = m r + r θ ) ] 3

4 so a Lagrangian for the object is Lr, θ, r, θ, t) = T V = m r + r θ ) ] + GmM r ) The equations for the motion of the object are given by the E-L equations. For the r coorinate, the E-L equation is 3Lr, θ, r, θ, t) = Lr, θ, r, θ ), t) r θ = r GM r ) The E-L equation for the θ coorinate is 4Lr, θ, r, θ, t) = Lr, θ, r, θ, t) r θ ) = 0 3) Equation 3) tells us that the quantity r θ is a conserve quantity of the motion. This quantity is calle the angular momentum. It is the magnitue of the vector L introuce in Equation 9). We will soon see that this quantity is conserve for all potentials V such that V = V r ), sometimes calle central potentials the associate forces are calle central forces). Let us write r θ = Γ = Cst Plugging Γ into Equation ), we have r = Γ r 3 GM r Multiplying this equation by r/, one fins ) r + Γ r GM r ] = 0 We conclue that ) r + Γ r GM r = Cst = E 4) The total energy of the system is conserve: the first two terms correspon to the kinetic energy, the last term is the potential energy. Equation 4) can be seen as the energy equation for a system with one egree of freeom, with an effective potential energy V eff r) = Γ r GM r This potential is ominate by the attractive gravitational part at large istances, but ominate by the repulsive angular momentum part at short istances. V eff r) is shown in Figure for the particular case Γ = GM =. Since 4) can be rewritten as ) r = E V eff r) point masses with E < 0 have boune orbits: there exists r in the attractive part of the potential such that E = V eff r ), an this r is a turning point of the trajectory. On the other han, point masses with E > 0 have unboune orbits. We will soon see that there are very efficient ways of fining out from the outset, from the functional epenence of the Lagrangian, that energy an angular momentum are conserve in the system stuie here. We coul have save a bit of algebra. Still, it oes not hurt too much to get practice with the umb ways at first. 4

5 .5 V eff r) r Figure : V eff r) as a function of r for Γ = GM = A little igression regaring the classical -boy problem It turns out that any system consisting of two boies with masses m an m, at positions r t) an r t) an each exerting the force F = ±G m m r r 3 r r ) on the other boy can be stuie within the framework we presente above, even if the masses m an m are comparable an none of the boies can be consiere static. Here is the reason why. Newton s equations for the boies are Let us efine the position vectors m r = G m m r r 3 r r ) 5) m r = G m m r r 3 r r ) 6) R = m r + m r m + m r = r r R is the position of the center of mass of the system. Aing Equations 5) an 6), one fins the following equation for R: R = 0 7) This is expecte. The system of the two boies is isolate, so in the absence of external forces, the center of mass of the whole system has a constant velocity. Subtracting 6)/m from 5)/m, one fins a secon orer ODE for r: µ r = G m m r 3 r 8) with µ = m m /m + m ), which is often calle the reuce mass. Equation 8) escribes the motion of a particle of mass µ an position vector r with respect to a fixe center of force. This is precisely the situation we have stuie with the Lagrangian formulation above. 5

6 The separate motion of the two boies is obtaine by solving for r using either Newton s approach or a variational approach), an reconstructing r t) an r t) with the relations m r = R + r m + m m r = R r m + m 6

Lagrangian and Hamiltonian Mechanics

Lagrangian and Hamiltonian Mechanics Lagrangian an Hamiltonian Mechanics.G. Simpson, Ph.. epartment of Physical Sciences an Engineering Prince George s Community College ecember 5, 007 Introuction In this course we have been stuying classical

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

and from it produce the action integral whose variation we set to zero:

and from it produce the action integral whose variation we set to zero: Lagrange Multipliers Monay, 6 September 01 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const.

G4003 Advanced Mechanics 1. We already saw that if q is a cyclic variable, the associated conjugate momentum is conserved, L = const. G4003 Avance Mechanics 1 The Noether theorem We alreay saw that if q is a cyclic variable, the associate conjugate momentum is conserve, q = 0 p q = const. (1) This is the simplest incarnation of Noether

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

221A Lecture Notes Notes on Classica Mechanics I

221A Lecture Notes Notes on Classica Mechanics I 1A Lecture Notes Notes on Classica Mechanics I 1 Precursor: Fermat s Principle in Geometric Optics In geometric optics, you talk about how light rays go. In homogeneous meiums, the light rays go straight.

More information

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion

Lecture 27: Generalized Coordinates and Lagrange s Equations of Motion Lecture 27: Generalize Coorinates an Lagrange s Equations of Motion Calculating T an V in terms of generalize coorinates. Example: Penulum attache to a movable support 6 Cartesian Coorinates: (X, Y, Z)

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012

Free rotation of a rigid body 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 Free rotation of a rigi boy 1 D. E. Soper 2 University of Oregon Physics 611, Theoretical Mechanics 5 November 2012 1 Introuction In this section, we escribe the motion of a rigi boy that is free to rotate

More information

G j dq i + G j. q i. = a jt. and

G j dq i + G j. q i. = a jt. and Lagrange Multipliers Wenesay, 8 September 011 Sometimes it is convenient to use reunant coorinates, an to effect the variation of the action consistent with the constraints via the metho of Lagrange unetermine

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

4. Important theorems in quantum mechanics

4. Important theorems in quantum mechanics TFY4215 Kjemisk fysikk og kvantemekanikk - Tillegg 4 1 TILLEGG 4 4. Important theorems in quantum mechanics Before attacking three-imensional potentials in the next chapter, we shall in chapter 4 of this

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes

12.5. Differentiation of vectors. Introduction. Prerequisites. Learning Outcomes Differentiation of vectors 12.5 Introuction The area known as vector calculus is use to moel mathematically a vast range of engineering phenomena incluing electrostatics, electromagnetic fiels, air flow

More information

MATH2231-Differentiation (2)

MATH2231-Differentiation (2) -Differentiation () The Beginnings of Calculus The prime occasion from which arose my iscovery of the metho of the Characteristic Triangle, an other things of the same sort, happene at a time when I ha

More information

u t v t v t c a u t b a v t u t v t b a

u t v t v t c a u t b a v t u t v t b a Nonlinear Dynamical Systems In orer to iscuss nonlinear ynamical systems, we must first consier linear ynamical systems. Linear ynamical systems are just systems of linear equations like we have been stuying

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v

Math 11 Fall 2016 Section 1 Monday, September 19, Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v Math Fall 06 Section Monay, September 9, 06 First, some important points from the last class: Definition: A vector parametric equation for the line parallel to vector v = x v, y v, z v passing through

More information

The Ehrenfest Theorems

The Ehrenfest Theorems The Ehrenfest Theorems Robert Gilmore Classical Preliminaries A classical system with n egrees of freeom is escribe by n secon orer orinary ifferential equations on the configuration space (n inepenent

More information

Vectors in two dimensions

Vectors in two dimensions Vectors in two imensions Until now, we have been working in one imension only The main reason for this is to become familiar with the main physical ieas like Newton s secon law, without the aitional complication

More information

Calculus of variations - Lecture 11

Calculus of variations - Lecture 11 Calculus of variations - Lecture 11 1 Introuction It is easiest to formulate the problem with a specific example. The classical problem of the brachistochrone (1696 Johann Bernoulli) is the search to fin

More information

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations

Lecture XII. where Φ is called the potential function. Let us introduce spherical coordinates defined through the relations Lecture XII Abstract We introuce the Laplace equation in spherical coorinates an apply the metho of separation of variables to solve it. This will generate three linear orinary secon orer ifferential equations:

More information

Quantum mechanical approaches to the virial

Quantum mechanical approaches to the virial Quantum mechanical approaches to the virial S.LeBohec Department of Physics an Astronomy, University of Utah, Salt Lae City, UT 84112, USA Date: June 30 th 2015 In this note, we approach the virial from

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1

Physics 5153 Classical Mechanics. The Virial Theorem and The Poisson Bracket-1 Physics 5153 Classical Mechanics The Virial Theorem an The Poisson Bracket 1 Introuction In this lecture we will consier two applications of the Hamiltonian. The first, the Virial Theorem, applies to systems

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form

LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1. where a(x) and b(x) are functions. Observe that this class of equations includes equations of the form LINEAR DIFFERENTIAL EQUATIONS OF ORDER 1 We consier ifferential equations of the form y + a()y = b(), (1) y( 0 ) = y 0, where a() an b() are functions. Observe that this class of equations inclues equations

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation

More information

Short Intro to Coordinate Transformation

Short Intro to Coordinate Transformation Short Intro to Coorinate Transformation 1 A Vector A vector can basically be seen as an arrow in space pointing in a specific irection with a specific length. The following problem arises: How o we represent

More information

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations

Designing Information Devices and Systems II Fall 2017 Note Theorem: Existence and Uniqueness of Solutions to Differential Equations EECS 6B Designing Information Devices an Systems II Fall 07 Note 3 Secon Orer Differential Equations Secon orer ifferential equations appear everywhere in the real worl. In this note, we will walk through

More information

Switching Time Optimization in Discretized Hybrid Dynamical Systems

Switching Time Optimization in Discretized Hybrid Dynamical Systems Switching Time Optimization in Discretize Hybri Dynamical Systems Kathrin Flaßkamp, To Murphey, an Sina Ober-Blöbaum Abstract Switching time optimization (STO) arises in systems that have a finite set

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

arxiv: v1 [physics.class-ph] 20 Dec 2017

arxiv: v1 [physics.class-ph] 20 Dec 2017 arxiv:1712.07328v1 [physics.class-ph] 20 Dec 2017 Demystifying the constancy of the Ermakov-Lewis invariant for a time epenent oscillator T. Pamanabhan IUCAA, Post Bag 4, Ganeshkhin, Pune - 411 007, Inia.

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere

Objective: To introduce the equations of motion and describe the forces that act upon the Atmosphere Objective: To introuce the equations of motion an escribe the forces that act upon the Atmosphere Reaing: Rea pp 18 6 in Chapter 1 of Houghton & Hakim Problems: Work 1.1, 1.8, an 1.9 on pp. 6 & 7 at the

More information

Topic 2.3: The Geometry of Derivatives of Vector Functions

Topic 2.3: The Geometry of Derivatives of Vector Functions BSU Math 275 Notes Topic 2.3: The Geometry of Derivatives of Vector Functions Textbook Sections: 13.2 From the Toolbox (what you nee from previous classes): Be able to compute erivatives scalar-value functions

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

The Principle of Least Action and Designing Fiber Optics

The Principle of Least Action and Designing Fiber Optics University of Southampton Department of Physics & Astronomy Year 2 Theory Labs The Principle of Least Action an Designing Fiber Optics 1 Purpose of this Moule We will be intereste in esigning fiber optic

More information

Multivariable Calculus: Chapter 13: Topic Guide and Formulas (pgs ) * line segment notation above a variable indicates vector

Multivariable Calculus: Chapter 13: Topic Guide and Formulas (pgs ) * line segment notation above a variable indicates vector Multivariable Calculus: Chapter 13: Topic Guie an Formulas (pgs 800 851) * line segment notation above a variable inicates vector The 3D Coorinate System: Distance Formula: (x 2 x ) 2 1 + ( y ) ) 2 y 2

More information

Quantum Mechanics in Three Dimensions

Quantum Mechanics in Three Dimensions Physics 342 Lecture 20 Quantum Mechanics in Three Dimensions Lecture 20 Physics 342 Quantum Mechanics I Monay, March 24th, 2008 We begin our spherical solutions with the simplest possible case zero potential.

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

The continuity equation

The continuity equation Chapter 6 The continuity equation 61 The equation of continuity It is evient that in a certain region of space the matter entering it must be equal to the matter leaving it Let us consier an infinitesimal

More information

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1

Angles-Only Orbit Determination Copyright 2006 Michel Santos Page 1 Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close

More information

Goldstein Chapter 1 Exercises

Goldstein Chapter 1 Exercises Golstein Chapter 1 Exercises Michael Goo July 17, 2004 1 Exercises 11. Consier a uniform thin isk that rolls without slipping on a horizontal plane. A horizontal force is applie to the center of the isk

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

Numerical Integrator. Graphics

Numerical Integrator. Graphics 1 Introuction CS229 Dynamics Hanout The question of the week is how owe write a ynamic simulator for particles, rigi boies, or an articulate character such as a human figure?" In their SIGGRPH course notes,

More information

Statics. There are four fundamental quantities which occur in mechanics:

Statics. There are four fundamental quantities which occur in mechanics: Statics Mechanics isabranchofphysicsinwhichwestuythestate of rest or motion of boies subject to the action of forces. It can be ivie into two logical parts: statics, where we investigate the equilibrium

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing

Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Examining Geometric Integration for Propagating Orbit Trajectories with Non-Conservative Forcing Course Project for CDS 05 - Geometric Mechanics John M. Carson III California Institute of Technology June

More information

THERE ARE BASICALLY TWO approaches to conservative mechanics, due principally

THERE ARE BASICALLY TWO approaches to conservative mechanics, due principally CHAPTER 2 LAGRANGE S EQUATIONS THERE ARE BASICALLY TWO approaches to conservative mechanics, ue principally to Lagrange an Hamilton respectively Each has avantages over the other, although the Lagrangian

More information

Lecture XVI: Symmetrical spacetimes

Lecture XVI: Symmetrical spacetimes Lecture XVI: Symmetrical spacetimes Christopher M. Hirata Caltech M/C 350-17, Pasaena CA 91125, USA (Date: January 4, 2012) I. OVERVIEW Our principal concern this term will be symmetrical solutions of

More information

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes Implicit Differentiation 11.7 Introuction This Section introuces implicit ifferentiation which is use to ifferentiate functions expresse in implicit form (where the variables are foun together). Examples

More information

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9

Prep 1. Oregon State University PH 213 Spring Term Suggested finish date: Monday, April 9 Oregon State University PH 213 Spring Term 2018 Prep 1 Suggeste finish ate: Monay, April 9 The formats (type, length, scope) of these Prep problems have been purposely create to closely parallel those

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Math 1271 Solutions for Fall 2005 Final Exam

Math 1271 Solutions for Fall 2005 Final Exam Math 7 Solutions for Fall 5 Final Eam ) Since the equation + y = e y cannot be rearrange algebraically in orer to write y as an eplicit function of, we must instea ifferentiate this relation implicitly

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

1 Heisenberg Representation

1 Heisenberg Representation 1 Heisenberg Representation What we have been ealing with so far is calle the Schröinger representation. In this representation, operators are constants an all the time epenence is carrie by the states.

More information

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b.

Lecture 1b. Differential operators and orthogonal coordinates. Partial derivatives. Divergence and divergence theorem. Gradient. A y. + A y y dy. 1b. b. Partial erivatives Lecture b Differential operators an orthogonal coorinates Recall from our calculus courses that the erivative of a function can be efine as f ()=lim 0 or using the central ifference

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Math 115 Section 018 Course Note

Math 115 Section 018 Course Note Course Note 1 General Functions Definition 1.1. A function is a rule that takes certain numbers as inputs an assigns to each a efinite output number. The set of all input numbers is calle the omain of

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Gyroscopic matrices of the right beams and the discs

Gyroscopic matrices of the right beams and the discs Titre : Matrice gyroscopique es poutres roites et es i[...] Date : 15/07/2014 Page : 1/16 Gyroscopic matrices of the right beams an the iscs Summary: This ocument presents the formulation of the matrices

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

GEODESIC BOUNDARY VALUE PROBLEMS WITH SYMMETRY. Colin J. Cotter. Darryl D. Holm. (Communicated by the associate editor name)

GEODESIC BOUNDARY VALUE PROBLEMS WITH SYMMETRY. Colin J. Cotter. Darryl D. Holm. (Communicated by the associate editor name) Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX GEODESIC BOUNDARY VALUE PROBLEMS WITH SYMMETRY Colin J. Cotter Department of Aeronautics Imperial

More information

Mathematical Review Problems

Mathematical Review Problems Fall 6 Louis Scuiero Mathematical Review Problems I. Polynomial Equations an Graphs (Barrante--Chap. ). First egree equation an graph y f() x mx b where m is the slope of the line an b is the line's intercept

More information

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann

Gravitation as the result of the reintegration of migrated electrons and positrons to their atomic nuclei. Osvaldo Domann Gravitation as the result of the reintegration of migrate electrons an positrons to their atomic nuclei. Osvalo Domann oomann@yahoo.com (This paper is an extract of [6] liste in section Bibliography.)

More information

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems

1 M3-4-5A16 Assessed Problems # 1: Do 4 out of 5 problems D. D. Holm M3-4-5A16 Assesse Problems # 1 Due 1 Nov 2012 1 1 M3-4-5A16 Assesse Problems # 1: Do 4 out of 5 problems Exercise 1.1 (Poisson brackets for the Hopf map) Figure 1: The Hopf map. In coorinates

More information

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes

3.7 Implicit Differentiation -- A Brief Introduction -- Student Notes Fin these erivatives of these functions: y.7 Implicit Differentiation -- A Brief Introuction -- Stuent Notes tan y sin tan = sin y e = e = Write the inverses of these functions: y tan y sin How woul we

More information

2.20 Marine Hydrodynamics Lecture 3

2.20 Marine Hydrodynamics Lecture 3 2.20 Marine Hyroynamics, Fall 2018 Lecture 3 Copyright c 2018 MIT - Department of Mechanical Engineering, All rights reserve. 1.7 Stress Tensor 2.20 Marine Hyroynamics Lecture 3 1.7.1 Stress Tensor τ ij

More information

Linear and quadratic approximation

Linear and quadratic approximation Linear an quaratic approximation November 11, 2013 Definition: Suppose f is a function that is ifferentiable on an interval I containing the point a. The linear approximation to f at a is the linear function

More information

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask

5.4 Fundamental Theorem of Calculus Calculus. Do you remember the Fundamental Theorem of Algebra? Just thought I'd ask 5.4 FUNDAMENTAL THEOREM OF CALCULUS Do you remember the Funamental Theorem of Algebra? Just thought I' ask The Funamental Theorem of Calculus has two parts. These two parts tie together the concept of

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Physics 115C Homework 4

Physics 115C Homework 4 Physics 115C Homework 4 Problem 1 a In the Heisenberg picture, the ynamical equation is the Heisenberg equation of motion: for any operator Q H, we have Q H = 1 t i [Q H,H]+ Q H t where the partial erivative

More information

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003

arxiv:physics/ v2 [physics.ed-ph] 23 Sep 2003 Mass reistribution in variable mass systems Célia A. e Sousa an Vítor H. Rorigues Departamento e Física a Universiae e Coimbra, P-3004-516 Coimbra, Portugal arxiv:physics/0211075v2 [physics.e-ph] 23 Sep

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y Ph195a lecture notes, 1/3/01 Density operators for spin- 1 ensembles So far in our iscussion of spin- 1 systems, we have restricte our attention to the case of pure states an Hamiltonian evolution. Toay

More information

Partial Differential Equations

Partial Differential Equations Chapter Partial Differential Equations. Introuction Have solve orinary ifferential equations, i.e. ones where there is one inepenent an one epenent variable. Only orinary ifferentiation is therefore involve.

More information

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs

Problem Sheet 2: Eigenvalues and eigenvectors and their use in solving linear ODEs Problem Sheet 2: Eigenvalues an eigenvectors an their use in solving linear ODEs If you fin any typos/errors in this problem sheet please email jk28@icacuk The material in this problem sheet is not examinable

More information

Discrete Hamilton Jacobi Theory and Discrete Optimal Control

Discrete Hamilton Jacobi Theory and Discrete Optimal Control 49th IEEE Conference on Decision an Control December 15-17, 2010 Hilton Atlanta Hotel, Atlanta, GA, USA Discrete Hamilton Jacobi Theory an Discrete Optimal Control Tomoi Ohsawa, Anthony M. Bloch, an Melvin

More information

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7.

1 dx. where is a large constant, i.e., 1, (7.6) and Px is of the order of unity. Indeed, if px is given by (7.5), the inequality (7. Lectures Nine an Ten The WKB Approximation The WKB metho is a powerful tool to obtain solutions for many physical problems It is generally applicable to problems of wave propagation in which the frequency

More information

Basic Differentiation Rules and Rates of Change. The Constant Rule

Basic Differentiation Rules and Rates of Change. The Constant Rule 460_00.q //04 4:04 PM Page 07 SECTION. Basic Differentiation Rules an Rates of Change 07 Section. The slope of a horizontal line is 0. Basic Differentiation Rules an Rates of Change Fin the erivative of

More information

Physics 170 Week 7, Lecture 2

Physics 170 Week 7, Lecture 2 Physics 170 Week 7, Lecture 2 http://www.phas.ubc.ca/ goronws/170 Physics 170 203 Week 7, Lecture 2 1 Textbook Chapter 12:Section 12.2-3 Physics 170 203 Week 7, Lecture 2 2 Learning Goals: Learn about

More information

1 Lecture 20: Implicit differentiation

1 Lecture 20: Implicit differentiation Lecture 20: Implicit ifferentiation. Outline The technique of implicit ifferentiation Tangent lines to a circle Derivatives of inverse functions by implicit ifferentiation Examples.2 Implicit ifferentiation

More information

Pure Further Mathematics 1. Revision Notes

Pure Further Mathematics 1. Revision Notes Pure Further Mathematics Revision Notes June 20 2 FP JUNE 20 SDB Further Pure Complex Numbers... 3 Definitions an arithmetical operations... 3 Complex conjugate... 3 Properties... 3 Complex number plane,

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

arxiv: v2 [math.ds] 26 Apr 2017

arxiv: v2 [math.ds] 26 Apr 2017 GAUSSIAN CURVATURE AND GYROSCOPES GRAHAM COX AND MARK LEVI arxiv:1607.03217v2 [math.ds] 26 Apr 2017 Abstract. We relate Gaussian curvature to the gyroscopic force, thus giving a mechanical interpretation

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A

AN INTRODUCTION TO AIRCRAFT WING FLUTTER Revision A AN INTRODUCTION TO AIRCRAFT WIN FLUTTER Revision A By Tom Irvine Email: tomirvine@aol.com January 8, 000 Introuction Certain aircraft wings have experience violent oscillations uring high spee flight.

More information

6 Wave equation in spherical polar coordinates

6 Wave equation in spherical polar coordinates 6 Wave equation in spherical polar coorinates We now look at solving problems involving the Laplacian in spherical polar coorinates. The angular epenence of the solutions will be escribe by spherical harmonics.

More information