HÉNON HEILES HAMILTONIAN CHAOS IN 2 D
|
|
- Gabriel Allen
- 5 years ago
- Views:
Transcription
1 ABSTRACT HÉNON HEILES HAMILTONIAN CHAOS IN D MODELING CHAOS & COMPLEXITY 008 YOUVAL DAR PHYSICS DAR@PHYSICS.UCDAVIS.EDU Chaos in two degrees of freedom, demonstrated b using the Hénon Heiles Hamiltonian INTRODUCTION A star moving about a galactic center can be considered a point mass. The Hamiltonian governing this motion will have three degrees of freedom (si coordinates in phase space) and will look like H = p / m + V( ρ, z), p. In 964, while investigating the motion of stars around the galactic center, Michael Hénon and Carl Heiles wanted to simplif the problem, the reduced the degrees of freedom b arguing invariant energ l z, angular momentum (due to the clindrical smmetr) supported b the observed motion that was confined to two dimensions. The chose to use a Hamiltonian that is used for motion in a near elliptic equilibrium. In this paper I will eamine Hénon and Heiles Hamiltonian dnamic properties, to some etent. In particular, I will look at the chaotic nature of the motion described b the equation of motion that derived from this Hamiltonian. E PROJECT S OBJECTIVE To gain an understanding of the dnamics of this specific Hamiltonian, to locate the onset of chaos and to observe what general lesson we can learn from it. To do that using computational tools learned in the course and using Pthon. BASIC ANALYSIS The Hénon Heiles Hamiltonian: H p r p θ kr λ r sin( θ ) = (.) m mr In Cartesian coordinates p p ( ) = λ( ) (.) H k m m And converting, normalizing and making dimensionless constants
2 H ( = p + p + + ) + (.) This equation can be further generalized in the form H = ( p + p + A + B ) + D C (.4) Deriving the equation of motion using H q = p H p = q (.5) The equations of motion are = p = p p = p = (.6) For the fied points, we want = = p = p = 0 At the same time, so p p = 0 = 0 So = = 0, = 0, or 0, =, is a function of E (.7) We also have the condition of invariant energ which is reduced to ( ) E = + + (.8) So our fied points should compl with both (.7) and (.8)
3 The Jacobian: J = (.9) The Jacobian determinant: Det( J ) = Tr( J ) = 0 (.0) The Eigen values are: λ λ,,4 =± =± + (.) The sum of the Eigen values is zero, which is consistent with energ conservation. For chaos in torus we epect the Lapunov Characteristic Eponents (LCE) to be (+, 0, 0, ), to find them we can diagonalize the Jacobian and then find the Eigen values or use another method. v v,,4 ± ( ) ± ( ) = ± ( + ) ( + ) = (.) Using the conditions from (.7) Det J = + + with = Det J = ( ) 4 4, 0 ( ) Det J = with = E = ( ) 4, /6 (.) For = 0 Det( J ) 0, Det = for =, 0 E = / For = Det( J ) 0and Det = for =± E = /6
4 I did observe the epected Det = at the fied points (0,0),(0, ) but not at (, /).This solution is actuall a fied line along which onl point ields Det ( J ) =.I didn t figure out what is the meaning of that reasult. For the fied point = we get the energ level which is the boundar of the potential region that allows orbits and when = 0 we are out of that region. I want to note two interesting points in equation (.), E = /6is the limit of the energ above which we get no orbit. And the solution of E =, =0 indicates the range of possible values of and at those points we get elliptical orbit, elliptic in phase space oscillating through the center in space. GENERAL VISUALIZATION APPROACH M plan was to create graphical analsis tools and then probe the sstem. The graphical tool I wanted to build should have produced Poincaré plots for different intersecting planes, D phase space plots, capabilit to plot a single path and an evolution of arra of initial conditions, to plot D vs. plots, to calculate Lapunov eponents, Entrop rate h μ = λi λi > 0 and finall to evaluate what percentage of the solution space is chaotic. I started to work on that plan but had to significantl reduce m goals in order for this project to be completed on time. (The mid wa change of plan shows in the pthon program flow.) COMPUTATIONAL STRATEGY After choosing some initial conditions, I will use equations(.6), integrate them using Runga Kutta method and use equation (.) to make sure that we keep the condition of energ conservation, making sure not to accumulate error due to the numerical calculation. = p = p p = p = E = ( p + p + + ) + To get the Poincaré I will look at an intersection of the chosen D phase space solution with the plane (In the actual program I reduced it to the plane = 0 onl.) a + b cp d = 0 (.4) 4
5 p = E+ p p = E+ p (.5) The plan for error correction was to check the error of a certain number of iterations and randoml appl correction to one of the four phase space coordinates. I planned to do it this wa in order to avoid over correcting one coordinate. In practice, due to time limitations, I just monitored the numerical error and since it never went over e 6I never applied corrections (the correction functions, though not use are in the program). LOOKING AT OUR POTENTIAL Potential equipotentials lines, V(, ) = ( + + ) the straight lines are for E = and the 6 smaller closed curves are for a smaller E. Note the ± π /smmetr of the potential. Recall the results of the Jacobian determinant. Finding etremums b finding where the first derivatives are zero and then using the descriminant st derivatives : V = (+4), V = ( +- ) (.6) V = ( + 4 ), V =V =, V = ( - 4 ) (.7) The Descriminant: D = V V -V If D>0 and V (, )>0, the point is a relative minimum.. If D>0 and V (, )<0, the point is a relative maimum.. If D<0, the point is a saddle point. 4. If D=0, higher order tests must be used. (.8) The etremum points I found were: (0, 0), (0, ), ( /, -/), (- /, -/) 5
6 Using the descriminant point (0, 0) is a relative minimum and the rest are saddle points. So we can epect orbital motion to eit onl in the area within the triangle created b the points (0, ), ( /, -/), (- /, -/). When going out of that area the trajector just runs awa. (This can be verified b the program.) 6
7 CHAOTIC REGION The Poincaré plots together with the D phase space picture are ver useful in understanding what is going on. B plaing with the initial conditions it looks like the onset of chaos is around E = /9. When the energ was E < /9 we got non chaotic orbits, though some were ver comple. It looks like we have an elliptical orbit on the edge of possible solutions for the equations E = =0. (I did not investigate it enough to see, when we get period doubling). But, we could qualitativel observe how the chaotic region, in the solution space, is getting bigger with E when E > /9. FIGURE : POINCARÉ PLOT OF 5 INITIAL CONDITIONS DIMENTIONLESS ENERGY, E=/, NON CHAOTIC ORBIS FIGURE : E = /9. MORE COMPLEX PICTURE BUT STILL NO CHAOTIC REGIONS 7
8 FIGURE : CHAOTIC REGIONS STARTING TO APPEAR FIGURE 4: VERY CHAOTIC FIGURE 5: MOSTLY CHAOTIC 8
9 SCREEN SHOTS OF SOME D PHASE SPACE PLOTS FIGURE 6: E=/, = 0, = , p =,0 ELIPTICAL ORBIT IN PHASE SPACE, = 0.5 IS A BIFURACATION POINT THE COLOR INDICATES THE AGE OF THE CURVE FIGURE 7: E=/, = 0, = 0.4, p =,0 STILL NON CHAOTIC REGION THE COLOR INDICATES THE AGE OF THE CURVE 9
10 FIGURE 8: E=/6, = 0, = 0., p =,0 CHAOTIC ORBIT FIGURE 9: E=/6, = 0, = 0., p =,0 NON CHAOTIC ORBIT IN MOSTLY CHAOTIC LANDSCAPE PROJECT CONCLUSIONS The onset of chaos for this sstem was approimatel E = /9. I found it b looking at the Poincaré maps and tring different energies. Poincaré is a great tool to look at a sstem that has three independent coordinates; it is simple, and together with the D plots, ver helpful in getting some intuition and basic analsis for the sstem dnamics. When planning a project, if ou have more realistic goals ou can achieve more. I spent too much time thinking how to write the program to be more general and fleible, and at the end had to cut it in order to get results which left the program mess and left man interesting quantities not eplored. 0
11 PROGRAM ALGORITHMS I DID NOT INCLUDE IN THE PROGRAM Multi plane intersection The wa I thought to allow intersecting the D phase space with an plane was, to rotate the coordinate sstem in a wa that the intersecting plane will be at ' = 0, this, I thought, will make it much easier to determine when the curve crosses the intersection plane. The algorithm: The plane: a + b + cz = 0 The normal to the plane: nˆ = ( abc,, ) The normal to the plane in the final orientation: n = (0,0,) The two rotation angles are: θ = n nˆ i And φ =n ˆ i ( nˆ n ) Rotating to the new coordinate sstem Checking and collecting the points ( ', p', ') where ' changes sign Printing ( ', p' ) Lapunov Eponents Following the idea that the orthonormal vectors tangent to the curve are the quantities we need to follow Define all the derivative equations, d i, 4 D phase space 4 equations Initialize the vectors (,0,0,0...),(0,,0,0...)... Iterate the di equations and find the new vector length, di,(this is an indication of stretching or contracting in a certain direction) LCE ( ) i = LCEi + Log di Now use the Gramm Schmit process to make the new set of vectors orthonormal again. Repeat for as man iterations as desired and at the end do LCE = LCE /( iteration #) This process is taken from the eample pthon programs. BIBLIOGRAPHY. Jim Crutchfield lecture notes for the course Modeling Chaos & Compleit 008. Strogatz, Nonlinear Dnamics and Chaos. James D. Meiss, Differential Dnamical Sstems 4. Goldstein Poole & Safko, Classical Mechanics 5. HeilesEquation.html
4452 Mathematical Modeling Lecture 13: Chaos and Fractals
Math Modeling Lecture 13: Chaos and Fractals Page 1 442 Mathematical Modeling Lecture 13: Chaos and Fractals Introduction In our tetbook, the discussion on chaos and fractals covers less than 2 pages.
More information15. Eigenvalues, Eigenvectors
5 Eigenvalues, Eigenvectors Matri of a Linear Transformation Consider a linear ( transformation ) L : a b R 2 R 2 Suppose we know that L and L Then c d because of linearit, we can determine what L does
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationThe Force Table Introduction: Theory:
1 The Force Table Introduction: "The Force Table" is a simple tool for demonstrating Newton s First Law and the vector nature of forces. This tool is based on the principle of equilibrium. An object is
More informationEssential Question How can you solve a nonlinear system of equations?
.5 Solving Nonlinear Sstems Essential Question Essential Question How can ou solve a nonlinear sstem of equations? Solving Nonlinear Sstems of Equations Work with a partner. Match each sstem with its graph.
More informationWorksheet #1. A little review.
Worksheet #1. A little review. I. Set up BUT DO NOT EVALUATE definite integrals for each of the following. 1. The area between the curves = 1 and = 3. Solution. The first thing we should ask ourselves
More informationLocal Phase Portrait of Nonlinear Systems Near Equilibria
Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating
More informationN coupled oscillators
Waves 1 1 Waves 1 1. N coupled oscillators towards the continuous limit. Stretched string and the wave equation 3. The d Alembert solution 4. Sinusoidal waves, wave characteristics and notation T 1 T N
More informationPhysics Gravitational force. 2. Strong or color force. 3. Electroweak force
Phsics 360 Notes on Griffths - pluses and minuses No tetbook is perfect, and Griffithsisnoeception. Themajorplusisthat it is prett readable. For minuses, see below. Much of what G sas about the del operator
More informationENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 3 Solutions [Multiple Integration; Lines of Force]
ENGI 44 Advanced Calculus for Engineering Facult of Engineering and Applied Science Problem Set Solutions [Multiple Integration; Lines of Force]. Evaluate D da over the triangular region D that is bounded
More informationLab 5 Forces Part 1. Physics 211 Lab. You will be using Newton s 2 nd Law to help you examine the nature of these forces.
b Lab 5 Forces Part 1 Phsics 211 Lab Introduction This is the first week of a two part lab that deals with forces and related concepts. A force is a push or a pull on an object that can be caused b a variet
More informationMath 20 Spring 2005 Final Exam Practice Problems (Set 2)
Math 2 Spring 2 Final Eam Practice Problems (Set 2) 1. Find the etreme values of f(, ) = 2 2 + 3 2 4 on the region {(, ) 2 + 2 16}. 2. Allocation of Funds: A new editor has been allotted $6, to spend on
More informationLecture 1: A Preliminary to Nonlinear Dynamics and Chaos
Lecture 1: A Preliminary to Nonlinear Dynamics and Chaos Autonomous Systems A set of coupled autonomous 1st-order ODEs. Here "autonomous" means that the right hand side of the equations does not explicitly
More informationLaurie s Notes. Overview of Section 3.5
Overview of Section.5 Introduction Sstems of linear equations were solved in Algebra using substitution, elimination, and graphing. These same techniques are applied to nonlinear sstems in this lesson.
More informationMTHE 227 Problem Set 10 Solutions. (1 y2 +z 2., 0, 0), y 2 + z 2 < 4 0, Otherwise.
MTHE 7 Problem Set Solutions. (a) Sketch the cross-section of the (hollow) clinder + = in the -plane, as well as the vector field in this cross-section. ( +,, ), + < F(,, ) =, Otherwise. This is a simple
More informationCONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM
Chapter 4 CONSERVATION OF ANGULAR MOMENTUM FOR A CONTINUUM Figure 4.1: 4.1 Conservation of Angular Momentum Angular momentum is defined as the moment of the linear momentum about some spatial reference
More informationLESSON #42 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART 2 COMMON CORE ALGEBRA II
LESSON #4 - INVERSES OF FUNCTIONS AND FUNCTION NOTATION PART COMMON CORE ALGEBRA II You will recall from unit 1 that in order to find the inverse of a function, ou must switch and and solve for. Also,
More informationRELATIONS AND FUNCTIONS through
RELATIONS AND FUNCTIONS 11.1.2 through 11.1. Relations and Functions establish a correspondence between the input values (usuall ) and the output values (usuall ) according to the particular relation or
More informationMechanics Departmental Exam Last updated November 2013
Mechanics Departmental Eam Last updated November 213 1. Two satellites are moving about each other in circular orbits under the influence of their mutual gravitational attractions. The satellites have
More informationAlgebra II Notes Unit Six: Polynomials Syllabus Objectives: 6.2 The student will simplify polynomial expressions.
Algebra II Notes Unit Si: Polnomials Sllabus Objectives: 6. The student will simplif polnomial epressions. Review: Properties of Eponents (Allow students to come up with these on their own.) Let a and
More informationSome linear transformations on R 2 Math 130 Linear Algebra D Joyce, Fall 2013
Some linear transformations on R 2 Math 3 Linear Algebra D Joce, Fall 23 Let s look at some some linear transformations on the plane R 2. We ll look at several kinds of operators on R 2 including reflections,
More information4 The Cartesian Coordinate System- Pictures of Equations
The Cartesian Coordinate Sstem- Pictures of Equations Concepts: The Cartesian Coordinate Sstem Graphs of Equations in Two Variables -intercepts and -intercepts Distance in Two Dimensions and the Pthagorean
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More informationThe first change comes in how we associate operators with classical observables. In one dimension, we had. p p ˆ
VI. Angular momentum Up to this point, we have been dealing primaril with one dimensional sstems. In practice, of course, most of the sstems we deal with live in three dimensions and 1D quantum mechanics
More informationarxiv: v2 [nlin.cd] 24 Aug 2018
3D billiards: visualization of regular structures and trapping of chaotic trajectories arxiv:1805.06823v2 [nlin.cd] 24 Aug 2018 Markus Firmbach, 1, 2 Steffen Lange, 1 Roland Ketzmerick, 1, 2 and Arnd Bäcker
More information4.7. Newton s Method. Procedure for Newton s Method HISTORICAL BIOGRAPHY
4. Newton s Method 99 4. Newton s Method HISTORICAL BIOGRAPHY Niels Henrik Abel (18 189) One of the basic problems of mathematics is solving equations. Using the quadratic root formula, we know how to
More informationPhysics 207, Lecture 4, Sept. 15
Phsics 07, Lecture 4, Sept. 15 Goals for hapts.. 3 & 4 Perform vector algebra (addition & subtraction) graphicall or b, & z components Interconvert between artesian and Polar coordinates Distinguish position-time
More informationTriple Integrals. y x
Triple Integrals. (a) If is an solid (in space), what does the triple integral dv represent? Wh? (b) Suppose the shape of a solid object is described b the solid, and f(,, ) gives the densit of the object
More information9.1 VECTORS. A Geometric View of Vectors LEARNING OBJECTIVES. = a, b
vectors and POLAR COORDINATES LEARNING OBJECTIVES In this section, ou will: View vectors geometricall. Find magnitude and direction. Perform vector addition and scalar multiplication. Find the component
More informationMath Lesson 2-2 Properties of Exponents
Math-00 Lesson - Properties of Eponents Properties of Eponents What is a power? Power: An epression formed b repeated multiplication of the base. Coefficient Base Eponent The eponent applies to the number
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More informationDynamics of multiple pendula without gravity
Chaotic Modeling and Simulation (CMSIM) 1: 57 67, 014 Dnamics of multiple pendula without gravit Wojciech Szumiński Institute of Phsics, Universit of Zielona Góra, Poland (E-mail: uz88szuminski@gmail.com)
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More information3.1 Exponential Functions and Their Graphs
.1 Eponential Functions and Their Graphs Sllabus Objective: 9.1 The student will sketch the graph of a eponential, logistic, or logarithmic function. 9. The student will evaluate eponential or logarithmic
More information6.4 graphs OF logarithmic FUnCTIOnS
SECTION 6. graphs of logarithmic functions 9 9 learning ObjeCTIveS In this section, ou will: Identif the domain of a logarithmic function. Graph logarithmic functions. 6. graphs OF logarithmic FUnCTIOnS
More informationIncreasing and Decreasing Functions and the First Derivative Test. Increasing and Decreasing Functions. Video
SECTION and Decreasing Functions and the First Derivative Test 79 Section and Decreasing Functions and the First Derivative Test Determine intervals on which a unction is increasing or decreasing Appl
More informationAnalytic Trigonometry
CHAPTER 5 Analtic Trigonometr 5. Fundamental Identities 5. Proving Trigonometric Identities 5.3 Sum and Difference Identities 5.4 Multiple-Angle Identities 5.5 The Law of Sines 5.6 The Law of Cosines It
More informationAll parabolas through three non-collinear points
ALL PARABOLAS THROUGH THREE NON-COLLINEAR POINTS 03 All parabolas through three non-collinear points STANLEY R. HUDDY and MICHAEL A. JONES If no two of three non-collinear points share the same -coordinate,
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationSolutions to Two Interesting Problems
Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows
More informationSECTION 8-7 De Moivre s Theorem. De Moivre s Theorem, n a Natural Number nth-roots of z
8-7 De Moivre s Theorem 635 B eactl; compute the modulus and argument for part C to two decimal places. 9. (A) 3 i (B) 1 i (C) 5 6i 10. (A) 1 i 3 (B) 3i (C) 7 4i 11. (A) i 3 (B) 3 i (C) 8 5i 12. (A) 3
More information3D Schrödinger Eq. Today: Continue with hydrogen. Multi-electron atoms
3D Schrödinger Eq. Toda: Continue with hdrogen. Multi-electron atoms HWK 13 available online. Please fill out the online participation surve. Worth 10points on HWK 13. Final Eam is Monda Dec. 15 10:30A-1P
More informationD u f f x h f y k. Applying this theorem a second time, we have. f xx h f yx k h f xy h f yy k k. f xx h 2 2 f xy hk f yy k 2
93 CHAPTER 4 PARTIAL DERIVATIVES We close this section b giving a proof of the first part of the Second Derivatives Test. Part (b) has a similar proof. PROOF OF THEOREM 3, PART (A) We compute the second-order
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
More informationMath 123 Summary of Important Algebra & Trigonometry Concepts Chapter 1 & Appendix D, Stewart, Calculus Early Transcendentals
Math Summar of Important Algebra & Trigonometr Concepts Chapter & Appendi D, Stewart, Calculus Earl Transcendentals Function a rule that assigns to each element in a set D eactl one element, called f (
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationCh 3 Alg 2 Note Sheet.doc 3.1 Graphing Systems of Equations
Ch 3 Alg Note Sheet.doc 3.1 Graphing Sstems of Equations Sstems of Linear Equations A sstem of equations is a set of two or more equations that use the same variables. If the graph of each equation =.4
More informationf x, y x 2 y 2 2x 6y 14. Then
SECTION 11.7 MAXIMUM AND MINIMUM VALUES 645 absolute minimum FIGURE 1 local maimum local minimum absolute maimum Look at the hills and valles in the graph of f shown in Figure 1. There are two points a,
More informationSolving Linear Systems
1.4 Solving Linear Sstems Essential Question How can ou determine the number of solutions of a linear sstem? A linear sstem is consistent when it has at least one solution. A linear sstem is inconsistent
More informationLab 5 Forces Part 1. Physics 225 Lab. You will be using Newton s 2 nd Law to help you examine the nature of these forces.
b Lab 5 orces Part 1 Introduction his is the first week of a two part lab that deals with forces and related concepts. A force is a push or a pull on an object that can be caused b a variet of reasons.
More informationSystems of Linear Equations: Solving by Graphing
8.1 Sstems of Linear Equations: Solving b Graphing 8.1 OBJECTIVE 1. Find the solution(s) for a set of linear equations b graphing NOTE There is no other ordered pair that satisfies both equations. From
More information2 variables. is the same value as the solution of. 1 variable. You can use similar reasoning to solve quadratic equations. Work with a partner.
9. b Graphing Essential Question How can ou use a graph to solve a quadratic equation in one variable? Based on what ou learned about the -intercepts of a graph in Section., it follows that the -intercept
More informationAffine transformations. Brian Curless CSE 557 Fall 2014
Affine transformations Brian Curless CSE 557 Fall 2014 1 Reading Required: Shirle, Sec. 2.4, 2.7 Shirle, Ch. 5.1-5.3 Shirle, Ch. 6 Further reading: Fole, et al, Chapter 5.1-5.5. David F. Rogers and J.
More informationEigenvectors and Eigenvalues 1
Ma 2015 page 1 Eigenvectors and Eigenvalues 1 In this handout, we will eplore eigenvectors and eigenvalues. We will begin with an eploration, then provide some direct eplanation and worked eamples, and
More informationChapter 11. Correlation and Regression
Chapter 11 Correlation and Regression Correlation A relationship between two variables. The data can be represented b ordered pairs (, ) is the independent (or eplanator) variable is the dependent (or
More informationSecond-Order Linear Differential Equations C 2
C8 APPENDIX C Additional Topics in Differential Equations APPENDIX C. Second-Order Homogeneous Linear Equations Second-Order Linear Differential Equations Higher-Order Linear Differential Equations Application
More informationA Glance at the Standard Map
A Glance at the Standard Map by Ryan Tobin Abstract The Standard (Chirikov) Map is studied and various aspects of its intricate dynamics are discussed. Also, a brief discussion of the famous KAM theory
More informationnm nm
The Quantum Mechanical Model of the Atom You have seen how Bohr s model of the atom eplains the emission spectrum of hdrogen. The emission spectra of other atoms, however, posed a problem. A mercur atom,
More informationPolynomial and Rational Functions
Polnomial and Rational Functions Figure -mm film, once the standard for capturing photographic images, has been made largel obsolete b digital photograph. (credit film : modification of work b Horia Varlan;
More informationAccelerator School Transverse Beam Dynamics-2. V. S. Pandit
Accelerator School 8 Transverse Beam Dnamics- V. S. Pandit Equation of Motion Reference orbit is a single laner curve. Diole is used for bending and quadruole for focusing We use coordinates (r, θ, ) Diole
More informationResearch Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
Mathematical Problems in Engineering Volume, Article ID 88, pages doi:.//88 Research Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
More informationLESSON #48 - INTEGER EXPONENTS COMMON CORE ALGEBRA II
LESSON #8 - INTEGER EXPONENTS COMMON CORE ALGEBRA II We just finished our review of linear functions. Linear functions are those that grow b equal differences for equal intervals. In this unit we will
More informationKEY IDEAS. Chapter 1 Function Transformations. 1.1 Horizontal and Vertical Translations Pre-Calculus 12 Student Workbook MHR 1
Chapter Function Transformations. Horizontal and Vertical Translations A translation can move the graph of a function up or down (vertical translation) and right or left (horizontal translation). A translation
More informationNonlinear Systems Examples Sheet: Solutions
Nonlinear Sstems Eamples Sheet: Solutions Mark Cannon, Michaelmas Term 7 Equilibrium points. (a). Solving ẋ =sin 4 3 =for gives =as an equilibrium point. This is the onl equilibrium because there is onl
More information12.1 Systems of Linear equations: Substitution and Elimination
. Sstems of Linear equations: Substitution and Elimination Sstems of two linear equations in two variables A sstem of equations is a collection of two or more equations. A solution of a sstem in two variables
More information3 Polynomial and Rational Functions
3 Polnomial and Rational Functions 3.1 Quadratic Functions and Models 3.2 Polnomial Functions and Their Graphs 3.3 Dividing Polnomials 3.4 Real Zeros of Polnomials 3.5 Comple Zeros and the Fundamental
More informationChapter 8 Notes SN AA U2C8
Chapter 8 Notes SN AA U2C8 Name Period Section 8-: Eploring Eponential Models Section 8-2: Properties of Eponential Functions In Chapter 7, we used properties of eponents to determine roots and some of
More informationTangent Lines. Limits 1
Limits Tangent Lines The concept of the tangent line to a circle dates back at least to the earl das of Greek geometr, that is, at least 5 ears. The tangent line to a circle with centre O at a point A
More information14.1 Systems of Linear Equations in Two Variables
86 Chapter 1 Sstems of Equations and Matrices 1.1 Sstems of Linear Equations in Two Variables Use the method of substitution to solve sstems of equations in two variables. Use the method of elimination
More informationResearch Article Chaotic Attractor Generation via a Simple Linear Time-Varying System
Discrete Dnamics in Nature and Societ Volume, Article ID 836, 8 pages doi:.//836 Research Article Chaotic Attractor Generation via a Simple Linear Time-Varing Sstem Baiu Ou and Desheng Liu Department of
More informationFigure 1 Correlation diagram for BeH 2
Self-Stud Problems / Eam Preparation revise our computational chemistr workshop from last ear: http://www.huntresearchgroup.org.uk/teaching/ear1_lab_start.html o make sure ou have checked that the molecule
More information17.3. Parametric Curves. Introduction. Prerequisites. Learning Outcomes
Parametric Curves 7.3 Introduction In this Section we eamine et another wa of defining curves - the parametric description. We shall see that this is, in some was, far more useful than either the Cartesian
More informationFunctions. Introduction
Functions,00 P,000 00 0 970 97 980 98 990 99 000 00 00 Figure Standard and Poor s Inde with dividends reinvested (credit "bull": modification of work b Praitno Hadinata; credit "graph": modification of
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationUsing Intercept Form
8.5 Using Intercept Form Essential Question What are some of the characteristics of the graph of f () = a( p)( q)? Using Zeros to Write Functions Work with a partner. Each graph represents a function of
More informationSolving Linear-Quadratic Systems
36 LESSON Solving Linear-Quadratic Sstems UNDERSTAND A sstem of two or more equations can include linear and nonlinear equations. In a linear-quadratic sstem, there is one linear equation and one quadratic
More informationReview: critical point or equivalently f a,
Review: a b f f a b f a b critical point or equivalentl f a, b A point, is called a of if,, 0 A local ma or local min must be a critical point (but not conversel) 0 D iscriminant (or Hessian) f f D f f
More informationP 0 (x 0, y 0 ), P 1 (x 1, y 1 ), P 2 (x 2, y 2 ), P 3 (x 3, y 3 ) The goal is to determine a third degree polynomial of the form,
Bezier Curves While working for the Renault automobile compan in France, an engineer b the name of P. Bezier developed a sstem for designing car bodies based partl on some fairl straightforward mathematics.
More informationPeriodic Structures in FDTD
EE 5303 Electromagnetic Analsis Using Finite Difference Time Domain Lecture #19 Periodic Structures in FDTD Lecture 19 These notes ma contain coprighted material obtained under fair use rules. Distribution
More informationEssential Question What is the equation of a circle with center (h, k) and radius r in the coordinate plane?
10.7 Circles in the Coordinate Plane Essential Question What is the equation of a circle with center (h, k) and radius r in the coordinate plane? The Equation of a Circle with Center at the Origin Work
More informationQUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS
arxiv:1803.0591v1 [math.gm] QUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS Aleander Panfilov stable spiral det A 6 3 5 4 non stable spiral D=0 stable node center non stable node saddle 1 tr A QUALITATIVE
More informationMath-2. Lesson:1-2 Properties of Exponents
Math- Lesson:- Properties of Eponents Properties of Eponents What is a power? Power: An epression formed b repeated multiplication of the same factor. Coefficient Base Eponent The eponent applies to the
More informationFluid Mechanics II. Newton s second law applied to a control volume
Fluid Mechanics II Stead flow momentum equation Newton s second law applied to a control volume Fluids, either in a static or dnamic motion state, impose forces on immersed bodies and confining boundaries.
More informationA NOTE ON THE DYNAMICS AROUND THE L 1,2 LAGRANGE POINTS OF THE EARTH MOON SYSTEM IN A COMPLETE SOLAR SYSTEM MODEL
IAA-AAS-DCoSS1-8-8 A NOTE ON THE DYNAMICS AROUND THE L 1,2 LAGRANGE POINTS OF THE EARTH MOON SYSTEM IN A COMPLETE SOLAR SYSTEM MODEL Lian Yijun, Gerard Góme, Josep J. Masdemont, Tang Guojian INTRODUCTION
More informationProperties of Limits
33460_003qd //04 :3 PM Page 59 SECTION 3 Evaluating Limits Analticall 59 Section 3 Evaluating Limits Analticall Evaluate a it using properties of its Develop and use a strateg for finding its Evaluate
More informationUNCORRECTED SAMPLE PAGES. 3Quadratics. Chapter 3. Objectives
Chapter 3 3Quadratics Objectives To recognise and sketch the graphs of quadratic polnomials. To find the ke features of the graph of a quadratic polnomial: ais intercepts, turning point and ais of smmetr.
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More informationLESSON #24 - POWER FUNCTIONS COMMON CORE ALGEBRA II
1 LESSON #4 - POWER FUNCTIONS COMMON CORE ALGEBRA II Before we start to analze polnomials of degree higher than two (quadratics), we first will look at ver simple functions known as power functions. The
More information11.1 Inverses of Simple Quadratic and Cubic Functions
Locker LESSON 11.1 Inverses of Simple Quadratic and Cubic Functions Teas Math Standards The student is epected to: A..B Graph and write the inverse of a function using notation such as f (). Also A..A,
More informationCONSERVATION OF ENERGY FOR ACONTINUUM
Chapter 6 CONSERVATION OF ENERGY FOR ACONTINUUM Figure 6.1: 6.1 Conservation of Energ In order to define conservation of energ, we will follow a derivation similar to those in previous chapters, using
More informationMATRIX TRANSFORMATIONS
CHAPTER 5. MATRIX TRANSFORMATIONS INSTITIÚID TEICNEOLAÍOCHTA CHEATHARLACH INSTITUTE OF TECHNOLOGY CARLOW MATRIX TRANSFORMATIONS Matri Transformations Definition Let A and B be sets. A function f : A B
More informationUnit 10 - Graphing Quadratic Functions
Unit - Graphing Quadratic Functions PREREQUISITE SKILLS: students should be able to add, subtract and multipl polnomials students should be able to factor polnomials students should be able to identif
More informationSection 8.5 Parametric Equations
504 Chapter 8 Section 8.5 Parametric Equations Man shapes, even ones as simple as circles, cannot be represented as an equation where is a function of. Consider, for eample, the path a moon follows as
More informationPHY 5246: Theoretical Dynamics, Fall Assignment # 9, Solutions. y CM (θ = 0) = 2 ρ m
PHY 546: Theoretical Dnamics, Fall 5 Assignment # 9, Solutions Graded Problems Problem (.a) l l/ l/ CM θ x In order to find the equation of motion of the triangle, we need to write the Lagrangian, with
More information11.4 Polar Coordinates
11. Polar Coordinates 917 11. Polar Coordinates In Section 1.1, we introduced the Cartesian coordinates of a point in the plane as a means of assigning ordered pairs of numbers to points in the plane.
More informationC. Non-linear Difference and Differential Equations: Linearization and Phase Diagram Technique
C. Non-linear Difference and Differential Equations: Linearization and Phase Diaram Technique So far we have discussed methods of solvin linear difference and differential equations. Let us now discuss
More informationWe have examined power functions like f (x) = x 2. Interchanging x
CHAPTER 5 Eponential and Logarithmic Functions We have eamined power functions like f =. Interchanging and ields a different function f =. This new function is radicall different from a power function
More information1.5. Analyzing Graphs of Functions. The Graph of a Function. What you should learn. Why you should learn it. 54 Chapter 1 Functions and Their Graphs
0_005.qd /7/05 8: AM Page 5 5 Chapter Functions and Their Graphs.5 Analzing Graphs of Functions What ou should learn Use the Vertical Line Test for functions. Find the zeros of functions. Determine intervals
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More information