A first order system of ordinary di erential equations is a set of equations of the form .. (24) x 1 (t 0 )=x 10 x 2 (t 0 )=x 20.
|
|
- Paula Patterson
- 5 years ago
- Views:
Transcription
1 Lecture notes for Numerical Analysis 4 Systems of ODEs Topics: 1 Problem statement 2 Two motivating examples 3 Connection between systems and higher order equations 4 Numerical solutions 5 Phase plane analysis 6 Notes and further reading 41 Problem statement A first order system of ordinary di erential equations is a set of equations of the form x 0 1(t) =f 1 (t, x 1 (t),,x n (t)) x 0 2(t) =f 2 (t, x 1 (t),,x n (t)) (24) x 0 n(t) =f n (t, x 1 (t),,x n (t)), for some set of functions f i, i =1 n Note that each f i depends on all the variables x 1,,x n, and thus that solving the system thus corresponds to finding trajectories x i (t) such that all n equations are simultaneously true In general, the solution to the system (24) is not unique, as arbitrary constants are lost in the di erentiation process However, if initial values are given for each of the x i, x 1 (t 0 )=x 10 x 2 (t 0 )=x 20 x n (0) = x n0, (25) then the subsequent trajectories of each of the x i are generally uniquely determined Solving the system (24) with a particular set of initial values (25) is referred to as solving an initial value problem 42 Two motivating examples Applied science is rife with systems of ordinary di erential equations Here were present two simple examples designed to illustrate how such systmes might materialize in practice Example 41 Consider an ecosystem in which at time t the number of sheep is given by s(t) and the number of wolves by w(t) Left to its own devices, the sheep population will grow, but it is subject to predation by wolves Assuming sheep growth is exponential, and that predation is proportional to both the number of sheep and the number of wolves, we are led to the equation s 0 (t) = s(t) s(t)w(t) (26) On the other hand, it is reasonable to assume that the rate at which the wolf population grows will be proportional to level of predation, while the rate at which wolves die o is proportional to the total number of wolves Together, these assumptions lead to the following wolf growth model: w 0 (t) = s(t)w(t) w(t) (27) 21
2 43 Connection between systems and higher order equations Lecture notes for Numerical Analysis Taken together, equations (26) and (27) represent a first order system of the form s 0 (t) =f(t, s(t),w(t)) w 0 (t) =g(t, s(t),w(t)) where f(t,s(t),w(t)) = s(t) s(t)w(t) and g(t,s(t),w(t)) = s(t)w(t) w(t) If the number of sheep and wolves at time t 0 specified as s 0 and w 0, respectively, this system becomes an initial value problem with initial conditions s 0 and w 0 Example 42 Let x(t) denote the vertical position of a mass thrown straight up Let v(t) denote the vertical velocity of this mass Since velocity is defined as the rate of change of position, and the rate of change is given by the derivative, we have x 0 (t) =v(t) (28) On the other hand, Newton s Law states that the force on an object is equal to its mass times its acceleration, F (t) =m a(t) Since the gravitational force on an object is equal to the object s mass times the gravitational constant of 98m/s 2, and the acceleration of an object is the derivative of velocity, we have Together, (28) and (29) represent a first order system of the form v 0 (t) = 98 (29) x 0 (t) =f(t, x(t),v(t)) v 0 (t) =g(t, x(t),v(t)), where f(t,x(t),v(t)) = v(t) and g(t,s(t),w(t)) = 98 If the position and velocity at some particular time t 0 are specified as x 0 and v 0, respectively, this system becomes an initial value problem with initial values x 0 and v 0 43 Connection between systems and higher order equations A higher order di erential equation in variable is an equation of the form f(t, x(t),x 0 (t),,x (n) (t)) = 0 (30) for some function f : R n+1! R Note that the left hand side is a rather compact way of writing any combination of x(t) and its higher order derivatives up to x (n) Since t shows up as an argument, functions of t can be coe cients or free standing terms in this expression Higher order di erential equations show up often in applied science Example 43 Above, we saw an example that used Newton s Law Note that since acceleration is the second derivative of position, Newton s Law can be restated as mx 00 = F, ie x 00 (t) = F (t) m, where x(t) is position and F (t) is the force at time t This is an equation of the form where f(t, x, x 0,x 00 )=x 00 (t) F (t) m f(t, x, x 0,x 00 )=0, In Example 42, we saw how a one-dimension example involving Newton s Law could be written as a first order system with two equations In Example (43), we saw how the same equation could be written as a single second order equation As the following theorem shows, this feature of being able to swap out higher derivatives for more first order equations is not an accident, but rather a very general property Indeed, this property is one of the reasons that the study of first order systems is so important 22
3 44 Numerical solutions Lecture notes for Numerical Analysis Theorem 41 Suppose that for some nth order di erential equation, the function f in (30) has the form f(t, x, x 0,,x (n) )=x (n) (t)+g(t, x(t),x 0 (t),,x (n 1) (t)) (31) for some function g Then this nth order equation can be expressed in the form (24), ie as a first order system of n equations in n unknowns Proof Given any equation of the form define the following variables: f(t, x, x 0,,x (n) )=0, y 1 (t) =x(t) y 2 (t) =x 0 (t) y 3 (t) =x 00 (t) y n (t) =x (n 1) (t) Then note that with this notation we have the following di erential equations: y1(t) 0 =y 2 (t) y2(t) 0 =y 3 (t) (32) yn 0 1(t) =y n (t) f(t, y 1 (t),y 2 (t), y n (t),yn(t)) 0 = 0 By (31), the last of these equations can be written as yn(t) 0 = g(t, y(t),y 0 (t),,y (n 1) (t)), which, together with the first n unknowns 1 equations in (32), represents a first order system of n equations in n It is worth noting that even when f does not have the form (31), the equation f(t, y 1 (t),y 2 (t), y n (t),y 0 n(t)) generally defines y 0 n(t) implicitly In other words, it is generally true that any nth order ODE can be written as a system of n first order ODEs in n variables 44 Numerical solutions Numerical solutions for systems are not substantively di erent than numerical solutions for single variable problems As with the single variable case, we distinguish between explicit and implicit schemes, and as with the single variable case, these schemes come in various flavors To develop concrete numerical algorithms for systems, the easiest approach is to write the system in vector notation, and proceed to develop finite di erence approximations exactly as in the single variable case To this end, define 2 3 and write the system (24) as X(t) = 6 4 x 1 (t) x n (t) 7 5 d X(t) =F (t, X(t)), (33) dt 23
4 45 Phase plane analysis Lecture notes for Numerical Analysis where F : R n+1! R n is the vector valued function whose components are just the functions f i, f 1 (t, X(t)) f 1 (t, x 1 (t),,x n (t)) F (t, X(t)) = 4 5 = 4 5 f n (t, X(t)) f n (t, x 1 (t),,x n (t)) To develop numerical schemes, we fix a time interval [t 0,t 1 ], let t 0,t 1,,t n denote uniformly spaced samples within that interval, and let X i denote X(t i ), i =0,,n To develop an explicit scheme, we let h denote the time step t i t i 1, and approximate (33) with a finite di erence of the form X i+1 h X i = F (t i,x i ) (explicit scheme) (34) To develop an implicit scheme, we replace the right hand side of (34) with F (t i+1,x i+1 ), ie X i+1 h X i = F (t i+1,x i+1 ) (implicit scheme) (35) Note that both (34) and (35) employ = signs, not signs: both are approximations of (33), but the equality is needed because each defines a numerical scheme, ie the succession of points X i Note that just as in the single variable case, the explicit scheme is easy to implement Indeed, once X i is known, X i+1 is given by X i+1 = X i + h F (t i,x i ) (explicit update step) Implicit schemes are harder to implement, since to find X i+1 we need to solve an equation of the form X i+1 h F (t i+1,x i+1 ) X i = 0 (implicit update step) To solve this, we need to know how to find roots of non-linear algebraic systems This is the topic of the next unit 45 Phase plane analysis Consider the simple sheep-wolf system described in Example 41 If we set = = = = 1 and suppose that the initial sheep and wolf populations are s 0 =19 and w 0 = 1, respectively, then we have a well defined initial value problem that we can solve using whatever numerical method we wish Below are some plots that illustrate the solutions: Figure 1: This figure illustrates how the populations of sheep and wolves evolve over time The first figure plots the respective sheep and wolve populations over time Note that these plots reveal a cyclic pattern: populations rise and fall, and there is no natural steady state 24
5 46 Notes and further reading Lecture notes for Numerical Analysis Figure 2: This figure is a phase-space diagram It shows which sheep-wolf pairs are possible The second figure is called a phase plane portrait Each point on that graph represents the population of both sheep and wolves at a particular time The phase plane portrait is a (deformed) circle, again indicating cyclicity Note that when the wolf population is maximal, the sheep population is small, and when the sheep population is maximal, the wolf population is small These findings jump out from a phase plane portrait in ways that they might not with a standard time-series plot For this reason, phase plane portraits form a basic tool in the numerical analysis of di erential systems 46 Notes and further reading Most of this material closely followed our class text The text provides a detailed analysis of a simplified sheep-wolf system In these notes I have suppressed that analysis, as it is rather specific to one rather unrealistic system, and I wanted to focus on big ideas However, the interested student should consult the text to get ideas for how such analysis might proceed In general, the mathematical analysis of di erential systems is rather ad hoc: we try to say something meaningful about the system, and the tricks we need to employ to do so will vary according to the system The fact that an nth order ODE can be expressed a first order system of n equations in n unknowns is discussed in Boyce and Deprima This reduction is one of the main reasons that the study of di erential systems is so important A variety of software exists to perform phase plan analysis One open source tool that you might consult is called pplane it has a couple of canned algorithms for calculating phase planes, and allows you to experiment easily with changing parameter values Sources: 1 Our textbook 2 PPLANE Copyright , John C Polking, Rice University ~dfield/dfpphtml 3 Boyce and DePrima, Elementary Di erential Equations, Wiley and Sons 25
Mathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationMATH 320, WEEK 2: Slope Fields, Uniqueness of Solutions, Initial Value Problems, Separable Equations
MATH 320, WEEK 2: Slope Fields, Uniqueness of Solutions, Initial Value Problems, Separable Equations 1 Slope Fields We have seen what a differential equation is (relationship with a function and its derivatives),
More informationD1.3 Separable Differential Equations
Section 5.3 Separable Differential Equations D.3 Separable Differential Equations Sketching solutions of a differential equation using its direction field is a powerful technique, and it provides a wealth
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
More informationCHAPTER 1. First-Order Differential Equations and Their Applications. 1.1 Introduction to Ordinary Differential Equations
CHAPTER 1 First-Order Differential Equations and Their Applications 1.1 Introduction to Ordinary Differential Equations Differential equations are found in many areas of mathematics, science, and engineering.
More informationMath 215/255 Final Exam, December 2013
Math 215/255 Final Exam, December 2013 Last Name: Student Number: First Name: Signature: Instructions. The exam lasts 2.5 hours. No calculators or electronic devices of any kind are permitted. A formula
More informationLecture 9. Systems of Two First Order Linear ODEs
Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predator - Prey Model Trajectories and the nonlinear conservation law James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline
More information(1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min.
CHAPTER 1 Introduction 1. Bacground Models of physical situations from Calculus (1) Rate of change: A swimming pool is emptying at a constant rate of 90 gal/min. With V = volume in gallons and t = time
More informationUsing Abel's Theorem to Explain Repeated Roots of the Characteristic Equation
CODEE Journal Volume 8 Article 3 7-26-20 Using Abel's Theorem to Explain Repeated Roots of the Characteristic Equation William Green Follow this and additional works at: http://scholarship.claremont.edu/codee
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations MATH 365 Ordinary Differential Equations J Robert Buchanan Department of Mathematics Fall 2018 Objectives Many physical problems involve a number of separate
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
More information(1 + 2y)y = x. ( x. The right-hand side is a standard integral, so in the end we have the implicit solution. y(x) + y 2 (x) = x2 2 +C.
Midterm 1 33B-1 015 October 1 Find the exact solution of the initial value problem. Indicate the interval of existence. y = x, y( 1) = 0. 1 + y Solution. We observe that the equation is separable, and
More informationLecture Notes in Mathematics. Arkansas Tech University Department of Mathematics
Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics Introductory Notes in Ordinary Differential Equations for Physical Sciences and Engineering Marcel B. Finan c All Rights
More informationMODULE 12. Topics: Ordinary differential equations
Topics: Ordinary differential equations MODULE 12 We shall be concerned with mappings = functions = operators which map a function space into a function space. We shall write, sort of generically and in
More informationV 1 V 2. r 3. r 6 r 4. Math 2250 Lab 12 Due Date : 4/25/2017 at 6:00pm
Math 50 Lab 1 Name: Due Date : 4/5/017 at 6:00pm 1. In the previous lab you considered the input-output model below with pure water flowing into the system, C 1 = C 5 =0. r 1, C 1 r 5, C 5 r r V 1 V r
More information4.9 Anti-derivatives. Definition. An anti-derivative of a function f is a function F such that F (x) = f (x) for all x.
4.9 Anti-derivatives Anti-differentiation is exactly what it sounds like: the opposite of differentiation. That is, given a function f, can we find a function F whose derivative is f. Definition. An anti-derivative
More information1 The relation between a second order linear ode and a system of two rst order linear odes
Math 1280 Spring, 2010 1 The relation between a second order linear ode and a system of two rst order linear odes In Chapter 3 of the text you learn to solve some second order linear ode's, such as x 00
More informationp(t)dt a p(τ)dτ , c R.
11 3. Solutions of first order linear ODEs 3.1. Homogeneous and inhomogeneous; superposition. A first order linear equation is homogeneous if the right hand side is zero: (1) ẋ + p(t)x = 0. Homogeneous
More informationOrdinary Di erential Equations Lecture notes for Math 133A. Slobodan N. Simić
Ordinary Di erential Equations Lecture notes for Math 133A Slobodan N. Simić c Slobodan N. Simić 2016 Contents Chapter 1. First order di erential equations 1 1.1. What is a di erential equation? 1 1.2.
More informationWithout fully opening the exam, check that you have pages 1 through 11.
Name: Section: Recitation Instructor: INSTRUCTIONS Fill in your name, etc. on this first page. Without fully opening the exam, check that you have pages through. Show all your work on the standard response
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationMath 353, Fall 2017 Study guide for the final
Math 353, Fall 2017 Study guide for the final December 12, 2017 Last updated 12/4/17. If something is unclear or incorrect, please let me know so I can update the documents. Updates 12/4: De-emphasized
More informationSolutions Definition 2: a solution
Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential
More informationENGI9496 Lecture Notes State-Space Equation Generation
ENGI9496 Lecture Notes State-Space Equation Generation. State Equations and Variables - Definitions The end goal of model formulation is to simulate a system s behaviour on a computer. A set of coherent
More informationAntiderivatives and Indefinite Integrals
Antiderivatives and Indefinite Integrals MATH 151 Calculus for Management J. Robert Buchanan Department of Mathematics Fall 2018 Objectives After completing this lesson we will be able to use the definition
More information1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.
G. NAGY ODE January 9, 2018 1 1.2. Introduction to Modeling Section Objective(s): Review of Exponential Growth. The Logistic Population Model. Competing Species Model. Overview of Mathematical Models.
More informationOrdinary Differential Equation Introduction and Preliminaries
Ordinary Differential Equation Introduction and Preliminaries There are many branches of science and engineering where differential equations arise naturally. Now days, it finds applications in many areas
More information1.1 Differential Equation Models. Jiwen He
1.1 Math 3331 Differential Equations 1.1 Differential Equation Models Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math3331 Jiwen He, University of
More informationLaplace Transforms Chapter 3
Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important
More informationFirst Order Systems of Linear Equations. or ODEs of Arbitrary Order
First Order Systems of Linear Equations or ODEs of Arbitrary Order Systems of Equations Relate Quantities Examples Predator-Prey Relationships r 0 = r (100 f) f 0 = f (r 50) (Lokta-Volterra Model) Systems
More informationLecture Notes of EE 714
Lecture Notes of EE 714 Lecture 1 Motivation Systems theory that we have studied so far deals with the notion of specified input and output spaces. But there are systems which do not have a clear demarcation
More informationFirst Order Differential Equations Lecture 3
First Order Differential Equations Lecture 3 Dibyajyoti Deb 3.1. Outline of Lecture Differences Between Linear and Nonlinear Equations Exact Equations and Integrating Factors 3.. Differences between Linear
More informationMotion Along a Straight Line
Chapter 2 Motion Along a Straight Line PowerPoint Lectures for University Physics, Twelfth Edition Hugh D. Young and Roger A. Freedman Lectures by James Pazun Copyright 2008 Pearson Education Inc., publishing
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationMATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES
MATH 353 LECTURE NOTES: WEEK 1 FIRST ORDER ODES J. WONG (FALL 2017) What did we cover this week? Basic definitions: DEs, linear operators, homogeneous (linear) ODEs. Solution techniques for some classes
More informationProjects and Investigations to accompany Mathematics for Electrical Engineering and Computing
Projects and Investigations to accompany Mathematics for Electrical Engineering and Computing Investigation 1 An investigation of the number 'e' Introduction An important equation which describes many
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationLecture 2. Introduction to Differential Equations. Roman Kitsela. October 1, Roman Kitsela Lecture 2 October 1, / 25
Lecture 2 Introduction to Differential Equations Roman Kitsela October 1, 2018 Roman Kitsela Lecture 2 October 1, 2018 1 / 25 Quick announcements URL for the class website: http://www.math.ucsd.edu/~rkitsela/20d/
More information37. f(t) sin 2t cos 2t 38. f(t) cos 2 t. 39. f(t) sin(4t 5) 40.
28 CHAPTER 7 THE LAPLACE TRANSFORM EXERCISES 7 In Problems 8 use Definition 7 to find {f(t)} 2 3 4 5 6 7 8 9 f (t),, f (t) 4,, f (t) t,, f (t) 2t,, f (t) sin t,, f (t), cos t, t t t 2 t 2 t t t t t t t
More informationDIFFERENTIAL EQUATIONS COURSE NOTES, LECTURE 2: TYPES OF DIFFERENTIAL EQUATIONS, SOLVING SEPARABLE ODES.
DIFFERENTIAL EQUATIONS COURSE NOTES, LECTURE 2: TYPES OF DIFFERENTIAL EQUATIONS, SOLVING SEPARABLE ODES. ANDREW SALCH. PDEs and ODEs, order, and linearity. Differential equations come in so many different
More informationPhysics 141. Lecture 3. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 03, Page 1
Physics 141. Lecture 3. Frank L. H. Wolfs Department of Physics and Astronomy, University of Rochester, Lecture 03, Page 1 Physics 141. Lecture 3. Today's Topics: Course Information: Laboratories - software.
More informationThese notes are based mostly on [3]. They also rely on [2] and [1], though to a lesser extent.
Chapter 1 Introduction These notes are based mostly on [3]. They also rely on [2] and [1], though to a lesser extent. 1.1 Definitions and Terminology 1.1.1 Background and Definitions The words "differential
More informationMATH 319, WEEK 2: Initial Value Problems, Existence/Uniqueness, First-Order Linear DEs
MATH 319, WEEK 2: Initial Value Problems, Existence/Uniqueness, First-Order Linear DEs 1 Initial-Value Problems We have seen that differential equations can, in general, given rise to multiple solutions.
More informationMath 310 Introduction to Ordinary Differential Equations Final Examination August 9, Instructor: John Stockie
Make sure this exam has 15 pages. Math 310 Introduction to Ordinary Differential Equations inal Examination August 9, 2006 Instructor: John Stockie Name: (Please Print) Student Number: Special Instructions
More informationIERCU. Dynamic Monopoly with Demand Delay
IERCU Institute of Economic Research, Chuo University 50th Anniversary Special Issues Discussion Paper No.15 Dynamic Monopoly with Demand Delay Akio Matsumoto Chuo University Ferenc Szidarovszky University
More informationFirst-Order ODEs. Chapter Separable Equations. We consider in this chapter differential equations of the form dy (1.1)
Chapter 1 First-Order ODEs We consider in this chapter differential equations of the form dy (1.1) = F (t, y), where F (t, y) is a known smooth function. We wish to solve for y(t). Equation (1.1) is called
More informationLS.5 Theory of Linear Systems
LS.5 Theory of Linear Systems 1. General linear ODE systems and independent solutions. We have studied the homogeneous system of ODE s with constant coefficients, (1) x = Ax, where A is an n n matrix of
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More information! 1.1 Definitions and Terminology
! 1.1 Definitions and Terminology 1. Introduction: At times, mathematics aims to describe a physical phenomenon (take the population of bacteria in a petri dish for example). We want to find a function
More informationUsing web-based Java pplane applet to graph solutions of systems of differential equations
Using web-based Java pplane applet to graph solutions of systems of differential equations Our class project for MA 341 involves using computer tools to analyse solutions of differential equations. This
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationAPPPHYS 217 Tuesday 6 April 2010
APPPHYS 7 Tuesday 6 April Stability and input-output performance: second-order systems Here we present a detailed example to draw connections between today s topics and our prior review of linear algebra
More informationPHYSICS 107. Lecture 5 Newton s Laws of Motion
PHYSICS 107 Lecture 5 Newton s Laws of Motion First Law We saw that the type of motion which was most difficult for Aristotle to explain was horizontal motion of nonliving objects, particularly after they've
More informationPhysics 1010: The Physics of Everyday Life. TODAY More Acceleration Newton s Second Law Gravity
Physics 11: The Physics of Everyday Life TODAY More Acceleration Newton s Second Law Gravity 1 Help, Office and Tutorial Hours TUTORIALS (G1B75/77/79) Isidoros: 12:32:3 Thursday (same as office hours)
More informationDonald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles
More informationComputer lab for MAN460
Computer lab for MAN460 (version 20th April 2006, corrected 20 May) Prerequisites Matlab is rather user friendly, and to do the first exercises, it is enough to write an m-file consisting of two lines,
More informationVector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave
Vector Fields and Solutions to Ordinary Differential Equations using MATLAB/Octave Andreas Stahel 5th December 27 Contents Vector field for the logistic equation 2 Solutions of ordinary differential equations
More informationDifferential Equations & Separation of Variables
Differential Equations & Separation of Variables SUGGESTED REFERENCE MATERIAL: As you work through the problems listed below, you should reference Chapter 8. of the recommended textbook (or the equivalent
More informationLecture 01: Mathematical Modeling and Physics
Massachusetts Institute of Technology MITES 2017 Physics III Lecture 01: Mathematical Modeling and Physics In these notes, we define physics and discuss how the properties of physical theories suggest
More informationf(s)ds, i.e. to write down the general solution in
1. First order ODEs. Simplest examples. Existence and uniqueness theorem Example 1.1. Consider the equation (1.1) x (t) = 1 This is an equation because there is an unknown: the function x(t). This is a
More informationMATH3203 Lecture 1 Mathematical Modelling and ODEs
MATH3203 Lecture 1 Mathematical Modelling and ODEs Dion Weatherley Earth Systems Science Computational Centre, University of Queensland February 27, 2006 Abstract Contents 1 Mathematical Modelling 2 1.1
More informationGoals: To develop skills needed to find the appropriate differential equations to use as mathematical models.
Unit #16 : Differential Equations Goals: To develop skills needed to find the appropriate differential equations to use as mathematical models. Differential Equation Modelling - 1 Differential Equation
More informationPHYS 211 Lecture 9 - Examples of 3D motion 9-1
PHYS 211 Lecture 9 - Examples of 3D motion 9-1 Lecture 9 - Examples of 3D motion Text: Fowles and Cassiday, Chap. 4 In one dimension, the equations of motion to be solved are functions of only one position
More information1 Differential. Equations. A differential equation is any equation that involves a derivative. For example, Newton s second law F ma
1 Differential Equations b The strange attractor for a Sprott system consisting of three quadratic differential equations. 1 A differential equation is any equation that involves a derivative. For example,
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationMIDTERM 1 PRACTICE PROBLEM SOLUTIONS
MIDTERM 1 PRACTICE PROBLEM SOLUTIONS Problem 1. Give an example of: (a) an ODE of the form y (t) = f(y) such that all solutions with y(0) > 0 satisfy y(t) = +. lim t + (b) an ODE of the form y (t) = f(y)
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 1: Vectors, Representations Algebra and Linear Algebra Algebra: numbers and operations on numbers 2 + 3 = 5 3 7 = 21 Linear Algebra: tuples, triples... of numbers
More informationMA 125 CALCULUS I SPRING 2007 April 27, 2007 FINAL EXAM. Name (Print last name first):... Student ID Number (last four digits):...
CALCULUS I, FINAL EXAM 1 MA 125 CALCULUS I SPRING 2007 April 27, 2007 FINAL EXAM Name (Print last name first):............................................. Student ID Number (last four digits):........................
More informationMath 127: Course Summary
Math 27: Course Summary Rich Schwartz October 27, 2009 General Information: M27 is a course in functional analysis. Functional analysis deals with normed, infinite dimensional vector spaces. Usually, these
More informationPhysics 101 Discussion Week 3 Explanation (2011)
Physics 101 Discussion Week 3 Explanation (2011) D3-1. Velocity and Acceleration A. 1: average velocity. Q1. What is the definition of the average velocity v? Let r(t) be the total displacement vector
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationParticle Motion. Typically, if a particle is moving along the x-axis at any time, t, x()
Typically, if a particle is moving along the x-axis at any time, t, x() t represents the position of the particle; along the y-axis, yt () is often used; along another straight line, st () is often used.
More informationChapter1. Ordinary Differential Equations
Chapter1. Ordinary Differential Equations In the sciences and engineering, mathematical models are developed to aid in the understanding of physical phenomena. These models often yield an equation that
More informationLecture 8: Ordinary Differential Equations
MIT-WHOI Joint Program Summer Math Review Summer 2015 Lecture 8: Ordinary Differential Equations Lecturer: Isabela Le Bras Date: 31 July 2015 Disclaimer: These notes are for the purposes of this review
More informationMATH 1231 MATHEMATICS 1B Calculus Section 3A: - First order ODEs.
MATH 1231 MATHEMATICS 1B 2010. For use in Dr Chris Tisdell s lectures. Calculus Section 3A: - First order ODEs. Created and compiled by Chris Tisdell S1: What is an ODE? S2: Motivation S3: Types and orders
More informationSection 9.8 Higher Order Linear Equations
Section 9.8 Higher Order Linear Equations Key Terms: Higher order linear equations Equivalent linear systems for higher order equations Companion matrix Characteristic polynomial and equation A linear
More information5 Finding roots of equations
Lecture notes for Numerical Analysis 5 Finding roots of equations Topics:. Problem statement. Bisection Method 3. Newton s Method 4. Fixed Point Iterations 5. Systems of equations 6. Notes and further
More informationUse separation of variables to solve the following differential equations with given initial conditions. y 1 1 y ). y(y 1) = 1
Chapter 11 Differential Equations 11.1 Use separation of variables to solve the following differential equations with given initial conditions. (a) = 2ty, y(0) = 10 (b) = y(1 y), y(0) = 0.5, (Hint: 1 y(y
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationLinear Systems of ODE: Nullclines, Eigenvector lines and trajectories
Linear Systems of ODE: Nullclines, Eigenvector lines and trajectories James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 6, 2013 Outline
More informationPhysics 6010, Fall Relevant Sections in Text: Introduction
Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational
More informationGetting Started With The Predator - Prey Model: Nullclines
Getting Started With The Predator - Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator
More informationUniversity of Alabama Department of Physics and Astronomy. PH 105 LeClair Summer Problem Set 3 Solutions
University of Alabama Department of Physics and Astronomy PH 105 LeClair Summer 2012 Instructions: Problem Set 3 Solutions 1. Answer all questions below. All questions have equal weight. 2. Show your work
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationCourse Name : Physics I Course # PHY 107. Note - 3 : Motion in One Dimension
Course Name : Physics I Course # PHY 107 Note - 3 : Motion in One Dimension Abu Mohammad Khan Department of Mathematics and Physics North South University https://abukhan.weebly.com Copyright: It is unlawful
More informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More information3. Solutions of first order linear ODEs
10 3. Solutions of first order linear ODEs 3.1. The homogeneous equation and its solutions. A first order linear equation is homogeneous if the right hand side is zero: (1) ẋ + p(t)x = 0. Homogeneous linear
More informationFirst Year Physics: Prelims CP1. Classical Mechanics: Prof. Neville Harnew. Problem Set III : Projectiles, rocket motion and motion in E & B fields
HT017 First Year Physics: Prelims CP1 Classical Mechanics: Prof Neville Harnew Problem Set III : Projectiles, rocket motion and motion in E & B fields Questions 1-10 are standard examples Questions 11-1
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationChapter 9: Impulse and Momentum
Midterm: covers everything in chapters 1-8 - three problems, each worth 10 points. - first problem is actually five short-answer (1 line) questions (definitions, F = this, a = that, what is m?) - second
More informationInverse and Implicit Mapping Theorems (Sections III.3-III.4)
MATH W8 Daily Notes Preamble As an executive decision, I am skipping Section III.2. It is something like an 8-page lemma, with a number of concepts and results introduced at this stage mainly for the purpose
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationLecture 1, August 21, 2017
Engineering Mathematics 1 Fall 2017 Lecture 1, August 21, 2017 What is a differential equation? A differential equation is an equation relating a function (known sometimes as the unknown) to some of its
More information