Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics"

Transcription

1 Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy: economic growth, competition Finance: stochastic equations 1

2 Population models with interactions: dx dt dy dt = a(x) + b(x, y) = c(y) + d(x, y) a(x): self-growth rate of x b(x, y): effect of interaction on x c(y): self-growth rate of y d(x, y): effect of interaction on y 2

3 predator-prey: a: prey, b: predator a(x) > 0, b(x, y) < 0, c(y) < 0 and d(x, y) > 0 Example: Canadian lynx and hare, fox and rabbit,... competition: a(x) > 0, b(x, y) < 0, c(y) > 0 and d(x, y) < 0 Example: tiger and lion, rabbit and sheep, USA and Japan in WWII,... cooperative (mutualism): a(x) < 0, b(x, y) > 0, c(y) < 0 and d(x, y) > 0 Example: plant and pollinator 3

4 A system of nonlinear differential equations: dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy A competition model, with x and y both grow at a logistic rate, and the competing for the same food source slows down the growth of both species Questions to study: what are the qualitative behavior of the nonlinear system? Given an initial value, what is the fate of both species? 4

5 Studies of phase portrait (1): nullclines dx = f(x, y) dt dy = g(x, y) dt The sets f(x, y) = 0 and g(x, y) = 0 are curves on the phase portrait, and these curves are called nullclines. The set f(x, y) = 0 is the x-nullcline, where the vector field (f, g) is vertical. The set g(x, y) = 0 is the y-nullcline, where the vector field (f, g) is horizontal. 5

6 dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy x-nullcline: 2x(1 x 2 ) xy = 0 x(2 x y) = 0 so x = 0 or 2 x y = 0 y-nullcline: 3y(1 y 3 ) 2xy = 0 y(3 y 2x) = 0 so y = 0 or 3 y 2x = 0 6

7 The nullclines divide the phase portraits into regions, and in each region, the direction of vector field must be one of the following: north-east, south-east, north-west and south-west (So nullclines are where the vector field is exactly east, west, north and south) In each region, we use an arrow to indicate the direction. (In 1-d, we use only up-arrow and down-arrow in phase lines.) 7

8 dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy x-nullcline: x = 0 or 2 x y = 0 On x = 0, x = 0 and y = 3y y 2 On 2 x y = 0, x = 0 and y = 3y y 2 2y(2 y) = y + y 2 y-nullcline: y = 0 or 3 y 2x = 0 On y = 0, y = 0 and x = 2x x 2 On 3 y 2x = 0, y = 0 and x = 2x x 2 x(3 2x) = x + x 2 The nullclines divide the region x > 0 and y > 0 to four sub-regions. 8

9 Studies of phase portrait (2): equilibrium points dx = f(x, y) dt dy = g(x, y) dt Equilibrium points are points where f(x, y) = 0 and g(x, y) = 0. Equilibrium points are the intersection points of x-nullcline and y-nullcline. Equilibrium points are constant solutions of the system. Equilibrium points are also called steady state solutions, fixed points, etc. 9

10 dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy x-nullcline: x = 0 or 2 x y = 0 y-nullcline: so y = 0 or 3 y 2x = 0 Equilibrium points: (0, 0), (0, 3), (2, 0) and (1, 1) 10

11 Qualitative analysis from nullclines: Suppose that there is a solution from a point in one of the regions formed by nullclines, then there is only three possibilities for the orbit: A. tends to an equilibrium on the border of this region B. goes away to infinity C. enter other neighboring region following the arrow More information is needed for equilibrium points to further determination. 11

12 Studies of phase portrait (3): linearization at equilibrium dx = f(x, y) dt dy = g(x, y) dt Suppose that (x 0, y 0 ) is an equilibrium point. Near it, the behavior of the solutions is governed by the linearized equation dx dt = f(x 0, y 0 ) x dy dt = g(x 0, y 0 ) x Since f(x, y) f(x 0, y 0 ) + f(x 0, y 0 ) x and g(x, y) g(x 0, y 0 ) + g(x 0, y 0 ) x (x x 0 ) + f(x 0, y 0 ) (y y 0 ) y (x x 0 ) + g(x 0, y 0 ) (y y 0 ) y (x x 0 ) + f(x 0, y 0 ) (y y 0 ) y (x x 0 ) + g(x 0, y 0 ) (y y 0 ) y 12

13 Quick Review of Multi-variable calculus: Linearization in 1-d: f(x) f(x 0 ) + f (x 0 )(x x 0 ) Linearization in 2-d: f(x, y) f(x 0, y 0 ) + f(x 0, y 0 ) (x x 0 ) + f(x 0, y 0 ) (y y 0 ) x y partial derivative: derivative of f w.r.t. x when y is fixed Notation: f(x 0, y 0 ) or f x (x 0, y 0 ) x Jacobian: ( all four partial ) derivatives of a vector field in a matrix fx (x 0, y 0 ) f y (x 0, y 0 ) g x (x 0, y 0 ) g y (x 0, y 0 ) 13

14 dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy Equilibrium points: (0, 0), (0, 3), (2, 0) and (1, 1) Jacobian J = ( 2 2x y x 2y 3 2y 2x ) J(0, 0) = ( ) (source) J(0, 3) = ( ) (sink) J(2, 0) = ( ) (sink) J(1, 1) = ( ) (saddle) 14

15 Classification of linear system: (Chap 3) Two real eigenvalues: 1. λ 1 > λ 2 > 0: source 2. λ 1 > λ 2 = 0: degenerate source 3. λ 1 > 0 > λ 2 : saddle 4. λ 1 = 0 > λ 2 : degenerate sink 5. 0 > λ 1 > λ 2 : sink Two complex eigenvalues: λ ± = a ± bi 1. a > 0: spiral source 2. a = 0: center 3. a < 0: spiral sink One real (repeated) eigenvalue: λ 1 = λ 2 = λ 1. λ > 0: star source or trying to spiral source 2. λ = 0: parallel lines or dumb system 3. λ < 0: star sink or trying to spiral sink 15

16 How to determine stability of equilibrium quickly (without solving eigenvalues)? ( a b c d ) ( x y ) = λ ( x y ), ( a λ b c d λ ) ( x y ) = ( 0 0 ), solve (a λ)(d λ) bc = 0 (characteristic equation) λ 2 (a + d)λ + (ad bc) = 0, or λ 2 T λ + D = 0 Trace=T = a + d = λ 1 + λ 2, Determinant=D = ad bc = λ 1 λ 2 (Routh-Hurwitz criteria for 2nd order systems): 1. Stable if T < 0 and D > 0 2. Unstable if D < 0, or T > 0 and D > 0 T 2 4D determines whether the eigenvalues are real or complex 16

17 Trace-determinant plane 1. unstable node, source; 2. unstable focus, spiral source; 3. stable focus, spiral sink; 4. stable node, sink; 5. saddle 17

18 Phase portrait of five types: 1. unstable node, source; 2. unstable focus, spiral source; 3. stable focus, spiral sink; 4. stable node, sink; 5. saddle 18

19 Suppose that P 0 = (x 0, y 0 ) is an equilibrium point, then a solution Y(t) is a stable orbit of P 0 if lim Y(t) = P 0 ; t and a solution Y(t) is a unstable orbit of P 0 if lim t Y(t) = P 0. If all solutions nearby are stable orbits, then the equilibrium point is stable. (sink and spiral sink) If an equilibrium point not stable, then it is unstable. (saddle, source and spiral source) If the equilibrium point is a saddle point, then there is only one pair of stable orbits, and one pair of unstable orbits. These orbits are called separatrices at the saddle point. (separatrix singular) When an equilibrium point is stable, the set of all its stable orbits is called the attraction basin of this equilibrium point. 19

20 dx dt = 2x(1 x 2 ) xy dy dt = 3y(1 y 3 ) 2xy At the saddle point (1, 1), there is a pair of stable orbits and a pair of unstable orbits, and they are the separatrices of this saddle point. Indeed, the stable orbits separates the whole region into two subregions. In the subregion above this separatrix, all solutions tend to the equilibrium point (0, 3) as t ; and in the subregion below the separatrix, all solutions tend to the equilibrium point (2, 0) as t. These two subregions are the attraction basins of (0, 3) and (2, 0) respectively. 20

21 Warning: If the linearized equation at an equilibrium point is a degenerate source, degenerate sink, center, parallel lines, then the nonlinear equation may have different phase plane near the equilibrium point. Linearization Theorem in 2-d: Suppose that (x 0, y 0 ) is an equilibrium point of x = f(x, y) and y = g(x, y), and the eigenvalues of Jacobian J(x 0, y 0 ) are λ 1 and λ 2. (1) λ 1 > λ 2 > 0, then the system is a source; (2) λ 1 > 0 > λ 2, then the system is a saddle; (3) 0 > λ 1 > λ 2, then the system is a sink; (4) λ 1,2 = a ± bi, a > 0, then the system is a spiral source; (5) λ 1,2 = a ± bi, a < 0, then the system is a spiral sink; (6) If the eigenvalues are other cases, then you need other information to determine the solution behavior near the equilibrium point. 21

22 Linearization Theorem in 1-d: Suppose that y = y 0 is an equilibrium point of y = f(y). if f (y 0 ) < 0, then y 0 is a sink; if f (y 0 ) > 0, then y 0 is a source; if f (y 0 ) = 0, then y 0 can be any type, but in addition if f (y 0 ) > 0 or f (y 0 ) < 0, then y 0 is a node. 22

23 When the linearization cannot provide sufficient information: dx dt = y (x2 + y 2 )x dy dt = x (x2 + y 2 )y In polar coordinates: dr dt = 2r2, dθ dt = 1

24 Job List for Phase plane analysis: 1. Find nullclines, and draw arrows on nullclines and bounded regions 2. Find equilibrium points, and linearization 3. Draw possible separatrices, and sketch the phase portrait 4. Explain the phase portrait Other possibilities: Periodic solutions Some points: 1. Nullclines may not be lines. (x 2 + y 2 = 2500) 2. Separatrices usually cannot be calculated explicitly, but near saddle point, they are tangent to the eigenvector direction. 23

25 Examples: (A) x = x(2 x y), y = y(y x 2 ). (x 0, y 0) (B) A bifurcation problem (more competition models) dx dy = x( x y + 70), = y( 2x y + a), (x 0, y 0) dt dt When 0 < a < 70, three equilibrium points, (0, 0), (70, 0) and (0, a). (70, 0) is a sink, and its attraction basin is the whole region x > 0, y > 0. When 70 < a < 140, four equilibrium points, (0, 0), (70, 0), (0, a) and (a 70, 140 a). Both (70, 0) and (0, a) are stable, their attraction basin splits the region. When a > 140, three equilibrium points, (0, 0), (70, 0) and (0, a). (0, a) is a sink, and its attraction basin is the whole region x > 0, y > 0. 24

26 Lotka-Volterra system of competition du 1 dτ = r 1U 1 (1 U 1 ) r 1α U 1 U 2, du 2 K 1 K 1 dτ = r 2U 2 (1 U 2 ) r 2β U 1 U 2 K 2 K 2 r 1, r 2 : growth rate K 1, K 2 : strength of the intraspecific competition α, β: strength of the interspecific competition Nondimensionalization: t = r 1 τ, u = U 1 K 1, v = U 2 K 2 New equation: du dt = u(1 u av), dv dt a = α K 2 K 1 b = β K 1 K 2 C = r 2 r 1 = cv(1 bu v) 25

27 Equilibrium points: (0, 0) (extinction), (1, 0) (u wins), (0, 1) (v wins), (u, v ) = ( 1 a 1 ab, 1 b 1 ab ) (coexistence) Case 1: a < 1 and b < 1, (u, v ) is stable Case 2: a > 1 and b > 1, (u, v ) is unstable, and (0, 1), (1, 0) are stable (bistable) Case 3: a < 1 < b, no (u, v ), (1, 0) stable, (0, 1) unstable (u is superior competitor) Case 4: a > 1 > b, no (u, v ), (0, 1) stable, (1, 0) unstable (v is superior competitor) 26

28 Principle of competition exclusion (Gause 1934) From textbook: if two species occupy the same ecological niche, then one of them will go extinct. An ecological niche is mode of existence that a species has within an ecosystem. Essentially it is the sum of all activities and relationships a species has while obtaining and using the resources needed to survive and reproduce. A species niche includes: (a). Habitat - where it lives in the ecosystem; (b). Relationships - all interactions with other species in the ecosystem; (c). Nutrition - its method of obtaining food. 27

29 Another version of Principle of competition exclusion: If the coexistence (u, v ) satisfies u 1 + v B is not possible. 1, then coexistence Evolution strategy: how to gain competition advantage r 1, r 2 : growth rate K 1, K 2 : strength of the intraspecific competition α, β: strength of the interspecific competition a = α K 2 b = β K 1 C = r 2 K 1 K 2 r 1 K-strategy: increase your carrying-capacity r-strategy: increase your growth rate 28

30 Ecological conclusions: 1. Similar species in the same habitat will not coexist. 2. Coexistence will only occur when the competition is less intense (for competition in same niche, ab = αβ = 1) 3. This principle may be not very useful in practice, as we can always find a way in which two species differ. 4. If two competitors coexist, there is some reason for them to do so. Or to make two competitors coexist, create difference in niche. 5. Species may distributed in patchy way, so competition can be minimized; or time-sharing arrangement with succession of species or seasonal variability can effect similar result. Gradual evolution of differing genetic traits (character displacement) to minimize competition, and more complex multi-species interactions in which predation mediates competition. 29

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

Math 232, Final Test, 20 March 2007

Math 232, Final Test, 20 March 2007 Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.

More information

Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 Solutions Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

More information

6.3. Nonlinear Systems of Equations

6.3. Nonlinear Systems of Equations G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: Two-Dimensional Nonlinear Systems. ritical Points and Linearization. The Hartman-Grobman Theorem. Part Two: ompeting

More information

3.5 Competition Models: Principle of Competitive Exclusion

3.5 Competition Models: Principle of Competitive Exclusion 94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless

More information

Lecture 20/Lab 21: Systems of Nonlinear ODEs

Lecture 20/Lab 21: Systems of Nonlinear ODEs Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species

More information

Models Involving Interactions between Predator and Prey Populations

Models Involving Interactions between Predator and Prey Populations Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate

More information

1 The pendulum equation

1 The pendulum equation Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

More information

MATH 215/255 Solutions to Additional Practice Problems April dy dt

MATH 215/255 Solutions to Additional Practice Problems April dy dt . For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the

More information

Lab 5: Nonlinear Systems

Lab 5: Nonlinear Systems Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting

More information

Is chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2?

Is chaos possible in 1d? - yes - no - I don t know. What is the long term behavior for the following system if x(0) = π/2? Is chaos possible in 1d? - yes - no - I don t know What is the long term behavior for the following system if x(0) = π/2? In the insect outbreak problem, what kind of bifurcation occurs at fixed value

More information

Ordinary Differential Equations

Ordinary Differential Equations Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations

More information

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12 Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider

More information

Lotka Volterra Predator-Prey Model with a Predating Scavenger

Lotka Volterra Predator-Prey Model with a Predating Scavenger Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

REUNotes08-ODEs May 30, Chapter Two. Differential Equations

REUNotes08-ODEs May 30, Chapter Two. Differential Equations Chapter Two Differential Equations 4 CHAPTER 2 2.1 LINEAR EQUATIONS We consider ordinary differential equations of the following form dx = f(t, x). dt where x is either a real variable or a vector of several

More information

Math Lecture 46

Math Lecture 46 Math 2280 - Lecture 46 Dylan Zwick Fall 2013 Today we re going to use the tools we ve developed in the last two lectures to analyze some systems of nonlinear differential equations that arise in simple

More information

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland

FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material

More information

Systems of Ordinary Differential Equations

Systems of Ordinary Differential Equations Systems of Ordinary Differential Equations Scott A. McKinley October 22, 2013 In these notes, which replace the material in your textbook, we will learn a modern view of analyzing systems of differential

More information

Math 266: Phase Plane Portrait

Math 266: Phase Plane Portrait Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions

More information

Math 312 Lecture Notes Linearization

Math 312 Lecture Notes Linearization Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior

More information

Nonlinear Autonomous Systems of Differential

Nonlinear Autonomous Systems of Differential Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

More information

Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm.

Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm. 1 competing species Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm. This section and the next deal with the subject of population biology. You will already have seen examples of this. Most calculus

More information

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy

1.2. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy .. Direction Fields: Graphical Representation of the ODE and its Solution Let us consider a first order differential equation of the form dy = f(x, y). In this section we aim to understand the solution

More information

A Stability Analysis on Models of Cooperative and Competitive Species

A Stability Analysis on Models of Cooperative and Competitive Species Research Journal of Mathematical and Statistical Sciences ISSN 2320 6047 A Stability Analysis on Models of Cooperative and Competitive Species Abstract Gideon Kwadzo Gogovi 1, Justice Kwame Appati 1 and

More information

Field experiments on competition. Field experiments on competition. Field experiments on competition

Field experiments on competition. Field experiments on competition. Field experiments on competition INTERACTIONS BETWEEN SPECIES Type of interaction species 1 species 2 competition consumer-resource (pred, herb, para) mutualism detritivore-detritus (food is dead) Field experiments on competition Example

More information

8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds

8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds c Dr Igor Zelenko, Spring 2017 1 8. Qualitative analysis of autonomous equations on the line/population dynamics models, phase line, and stability of equilibrium points (corresponds to section 2.5) 1.

More information

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point

More information

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic

More information

Stable Coexistence of a Predator-Prey model with a Scavenger

Stable Coexistence of a Predator-Prey model with a Scavenger ANTON DE KOM UNIVERSITEIT VAN SURINAME INSTITUTE FOR GRADUATE STUDIES AND RESEARCH Stable Coexistence of a Predator-Prey model with a Scavenger Thesis in partial fulfillment of the requirements for the

More information

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt

Lokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period

More information

Continuous time population models

Continuous time population models Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations

More information

The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor

The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor , pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo

More information

D. Correct! Allelopathy is a form of interference competition in plants. Therefore this answer is correct.

D. Correct! Allelopathy is a form of interference competition in plants. Therefore this answer is correct. Ecology Problem Drill 18: Competition in Ecology Question No. 1 of 10 Question 1. The concept of allelopathy focuses on which of the following: (A) Carrying capacity (B) Limiting resource (C) Law of the

More information

BIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences

BIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences BIOS 5970: Plant-Herbivore Interactions Dr. Stephen Malcolm, Department of Biological Sciences D. POPULATION & COMMUNITY DYNAMICS Week 10. Population models 1: Lecture summary: Distribution and abundance

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point

More information

Mathematical Modeling I

Mathematical Modeling I Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex

More information

. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point,

. For each initial condition y(0) = y 0, there exists a. unique solution. In fact, given any point (x, y), there is a unique curve through this point, 1.2. Direction Fields: Graphical Representation of the ODE and its Solution Section Objective(s): Constructing Direction Fields. Interpreting Direction Fields. Definition 1.2.1. A first order ODE of the

More information

Chapter 7. Nonlinear Systems. 7.1 Introduction

Chapter 7. Nonlinear Systems. 7.1 Introduction Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)

More information

Gerardo Zavala. Math 388. Predator-Prey Models

Gerardo Zavala. Math 388. Predator-Prey Models Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked

More information

Introduction to Dynamical Systems

Introduction to Dynamical Systems Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology Rabbit-Fox System Trout System Trout System

More information

Key words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.

Key words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability. ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart

More information

Global Stability Analysis on a Predator-Prey Model with Omnivores

Global Stability Analysis on a Predator-Prey Model with Omnivores Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani

More information

Dynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory

More information

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium

More information

4 Second-Order Systems

4 Second-Order Systems 4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization

More information

Lecture 3. Dynamical Systems in Continuous Time

Lecture 3. Dynamical Systems in Continuous Time Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a

More information

Dynamical Systems: Ecological Modeling

Dynamical Systems: Ecological Modeling Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied

More information

ANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1

ANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1 ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate

More information

arxiv: v1 [math.ds] 14 Mar 2019

arxiv: v1 [math.ds] 14 Mar 2019 Resilience Analysis for Competing Populations Artur César Fassoni Instituto de Matemática e Computação, Universidade Federal de Itajubá, Avenida BPS 1303, Pinheirinho, CEP 37.500 903, Itajubá, MG, Brazil.

More information

Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations

Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

Section 5. Graphing Systems

Section 5. Graphing Systems Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each

More information

Matrices, Linearization, and the Jacobi matrix. y f x g y g J = dy/dt = g(x, y) Theoretical Biology, Utrecht University

Matrices, Linearization, and the Jacobi matrix. y f x g y g J = dy/dt = g(x, y) Theoretical Biology, Utrecht University Matrices, Linearization, and the Jacobi matrix { dx/dt f(x, y dy/dt g(x, y ( x f J y f x g y g λ 1, tr ± tr 4 det Theoretical Biology, Utrecht University i c Utrecht University, 018 Ebook publically available

More information

2.10 Saddles, Nodes, Foci and Centers

2.10 Saddles, Nodes, Foci and Centers 2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one

More information

Non-Dimensionalization of Differential Equations. September 9, 2014

Non-Dimensionalization of Differential Equations. September 9, 2014 Non-Dimensionalization of Differential Equations September 9, 2014 1 Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side

More information

MathQuest: Differential Equations

MathQuest: Differential Equations MathQuest: Differential Equations Geometry of Systems 1. The differential equation d Y dt = A Y has two straight line solutions corresponding to [ ] [ ] 1 1 eigenvectors v 1 = and v 2 2 = that are shown

More information

We have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),

We have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v), Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey

More information

BIO S380T Page 1 Summer 2005: Exam 2

BIO S380T Page 1 Summer 2005: Exam 2 BIO S380T Page 1 Part I: Definitions. [5 points for each term] For each term, provide a brief definition that also indicates why the term is important in ecology or evolutionary biology. Where I ve provided

More information

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total

MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all

More information

Non-Linear Models. Non-Linear Models Cont d

Non-Linear Models. Non-Linear Models Cont d Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method

More information

E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions

E209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions E09A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions Michael Vitus Stanford University Winter 007 Problem : Planar phase portraits. Part a Figure : Problem a This phase portrait is correct.

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

x 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1

x 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1 ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability

More information

Global dynamics of a mutualism competition model with one resource and multiple consumers

Global dynamics of a mutualism competition model with one resource and multiple consumers University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln USGS Staff -- Published Research US Geological Survey 8-5-08 Global dynamics of a mutualism competition model with one resource

More information

Section 9.3 Phase Plane Portraits (for Planar Systems)

Section 9.3 Phase Plane Portraits (for Planar Systems) Section 9.3 Phase Plane Portraits (for Planar Systems) Key Terms: Equilibrium point of planer system yꞌ = Ay o Equilibrium solution Exponential solutions o Half-line solutions Unstable solution Stable

More information

M469, Fall 2010, Practice Problems for the Final

M469, Fall 2010, Practice Problems for the Final M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems

More information

MATH 650 : Mathematical Modeling Fall, Written Assignment #3

MATH 650 : Mathematical Modeling Fall, Written Assignment #3 1 Instructions: MATH 650 : Mathematical Modeling Fall, 2017 - Written Assignment #3 Due by 11:59 p.m. EST on Sunday, December 3th, 2017 The problems on this assignment involve concepts, solution methods,

More information

THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS

THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS Maria P. Skhosana and Stephan V. Joubert, Tshwane University

More information

Matrices, Linearization, and the Jacobi matrix

Matrices, Linearization, and the Jacobi matrix Matrices, Linearization, and the Jacobi matrix { dx/dt f(x, y dy/dt g(x, y ( x f J y f x g y g λ 1, tr ± tr 4 det Alexander V. Panfilov, Kirsten H.W.J. ten Tusscher & Rob J. de Boer Theoretical Biology,

More information

Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area

Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

More information

The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z).

The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). Review of Multi-variable calculus: The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). The spatial variable represents the environment where the species

More information

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington

Interactions. Yuan Gao. Spring Applied Mathematics University of Washington Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,

More information

A Producer-Consumer Model With Stoichiometry

A Producer-Consumer Model With Stoichiometry A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann

More information

1.2. Introduction to Modeling. P (t) = r P (t) (b) When r > 0 this is the exponential growth equation.

1.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 information

Problem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J)

Problem Set 2. 1 Competitive and uncompetitive inhibition (12 points) Systems Biology (7.32/7.81J/8.591J) Problem Set 2 1 Competitive and uncompetitive inhibition (12 points) a. Reversible enzyme inhibitors can bind enzymes reversibly, and slowing down or halting enzymatic reactions. If an inhibitor occupies

More information

Math 2930 Worksheet Introduction to Differential Equations. What is a Differential Equation and what are Solutions?

Math 2930 Worksheet Introduction to Differential Equations. What is a Differential Equation and what are Solutions? Math 2930 Worksheet Introduction to Differential Equations Week 1 January 25, 2019 What is a Differential Equation and what are Solutions? A differential equation is an equation that relates an unknown

More information

Effect of Species 2 on Species 1 Competition - - Predator-Prey + - Parasite-Host + -

Effect of Species 2 on Species 1 Competition - - Predator-Prey + - Parasite-Host + - Community Ecology Community - a group of organisms, of different species, living in the same area Community ecology is the study of the interactions between species The presence of one species may affect

More information

Community Ecology. Classification of types of interspecific interactions: Effect of Species 1 on Species 2

Community Ecology. Classification of types of interspecific interactions: Effect of Species 1 on Species 2 Community Ecology Community - a group of organisms, of different species, living in the same area Community ecology is the study of the interactions between species The presence of one species may affect

More information

EAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science

EAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Transient Response of a Chemical Reactor Concentration of a substance in a chemical reactor

More information

Interspecific Competition

Interspecific Competition Use Interspecific Competition 0.8 0.6 0.4 0. 0 0 0.5.5.5 3 Resource The niche and interspecific competition Species A Use Use 0.8 0.6 0.4 0. Species B 0 0 0.5.5.5 3 0.8 0.6 0.4 0. 0 0 0.5.5.5 3 Resource

More information

Modelling and Mathematical Methods in Process and Chemical Engineering

Modelling and Mathematical Methods in Process and Chemical Engineering Modelling and Mathematical Methods in Process and Chemical Engineering Solution Series 3 1. Population dynamics: Gendercide The system admits two steady states The Jacobi matrix is ẋ = (1 p)xy k 1 x ẏ

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

Practice Problems for Final Exam

Practice Problems for Final Exam Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,

More information

Differential Equations Handout A

Differential Equations Handout A Differential Equations Handout A Math 1b November 11, 2005 These are problems for the differential equations unit of the course, some topics of which are not covered in the main textbook. Some of these

More information

7 Planar systems of linear ODE

7 Planar systems of linear ODE 7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution

More information

Qualitative Analysis of Tumor-Immune ODE System

Qualitative Analysis of Tumor-Immune ODE System of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

A Discrete Model of Three Species Prey- Predator System

A Discrete Model of Three Species Prey- Predator System ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj

More information

Permanency and Asymptotic Behavior of The Generalized Lotka-Volterra Food Chain System

Permanency and Asymptotic Behavior of The Generalized Lotka-Volterra Food Chain System CJMS. 5(1)(2016), 1-5 Caspian Journal of Mathematical Sciences (CJMS) University of Mazandaran, Iran http://cjms.journals.umz.ac.ir ISSN: 1735-0611 Permanency and Asymptotic Behavior of The Generalized

More information

7.7 LOTKA-VOLTERRA M ODELS

7.7 LOTKA-VOLTERRA M ODELS 77 LOTKA-VOLTERRA M ODELS sstems, such as the undamped pendulum, enabled us to draw the phase plane for this sstem and view it globall In that case we were able to understand the motions for all initial

More information

Multiple choice 2 pts each): x 2 = 18) Essay (pre-prepared) / 15 points. 19) Short Answer: / 2 points. 20) Short Answer / 5 points

Multiple choice 2 pts each): x 2 = 18) Essay (pre-prepared) / 15 points. 19) Short Answer: / 2 points. 20) Short Answer / 5 points P 1 Biology 217: Ecology Second Exam Fall 2004 There should be 7 ps in this exam - take a moment and count them now. Put your name on the first p of the exam, and on each of the ps with short answer questions.

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

EECS 700: Exam # 1 Tuesday, October 21, 2014

EECS 700: Exam # 1 Tuesday, October 21, 2014 EECS 700: Exam # 1 Tuesday, October 21, 2014 Print Name and Signature The rules for this exam are as follows: Write your name on the front page of the exam booklet. Initial each of the remaining pages

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

Dynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs

Dynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs Dynamics of a opulation Model Controlling the Spread of lague in rairie Dogs Catalin Georgescu The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA

More information