2D-Volterra-Lotka Modeling For 2 Species

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "2D-Volterra-Lotka Modeling For 2 Species"

Transcription

1 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 of this paper is to model multi-species interactions using Volterra Lotka equations in two dimensions which is a family of first order autonomous ordinary differential equations. The model for a pair of interacting populations of the predator and prey is studied and the stability points have been viewed. Lyapunov function is used to view the stable points. The population dynamics of the resulting systems are analyzed in terms of stability around equilibrium points and within invariant surfaces. Introduction The Volterra lotka (VL) problem, originally introduced in 192 by A. Lotka and later applied by V. Volterra to model the prey predator interactions. It is determine the predator-prey system with continuous of non-linear predator-prey interaction, the predicts oscillatory behaviour of populations with constant amplitudes dependent on the initial conditions [1], [3]. The discussion starts by presenting the basic exponential growth 2D VL. of differential equations and analyzing it in terms of stability of stationary points [5]. In this paper we examine some linear and non-linear two dimensional systems that have been used as mathematical model of the growth of the species showing a common environment [4]. In the first unit, which treats only a single species variation mathematical a summations on the growth rate are discussed. In the second unit, the simplest types of equations that model predator-prey ecology are investigated [2], [3]. The object is to find out the long-run qualitative behavior of these dynamical systems. 1 Alhashmi Darah, PO Box 675 Zliten Libya, << >>. 63

2 2D-Volterra-Lotka modeling for 2 species n Definition: LetM, the differentiable function φ : M M is called dynamical system if: 1- φ(, x) x. 2- φ( t, φ( s, x)) φ( t s, x) for every t, s, x M. Definition: An equilibrium solution of the system: dx x ( t) f ( x( t)), (1) dt n x : a vector valued function, is a point x such that f ( x). where Definition: x(t) is said to be stable if for any given, there exists n ( ), such that, for any other solution y( t) of (1) satisfying x t ) y( ), then x t ) y( ) for all t t, t. ( t ( t Definition: x(t) is said to be asymptotically stable if it is stable and if there exists a constant b such that, if x( t ) y( ) b, then lim x ( t ) y( t ). t 1. One Specie t The birth-rate of human population is usually given in terms of the number of births per thousand in one year [1]. Similarly, one can define death rate and the growth rate = (birth rate death rate)/ total population. The growth rate is thus the net change in population per unit of time divided by the total population at the beginning of the time period. The growth rate is: y y( t) lim. yt y( t) t 64

3 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya If the growth rate is a constant =, then: Thus, y y d (ln. dt t y( t) e y, y y(, is the predator population at the time t where ) exponential growth rate. If a ( ), a,, we have an then the growth rate is positive if, and negative if, while in case, the growth rate equal zero. The solution of the differential equation: dy dt a( ) y( t) (1) is: y ( t ) exp[( a ( ) t ]. If we assume that is the limit of the population, then we assume that the growth rate is proportional to ( y ): c(, c, Then differential equation (1) takes the form: dy c( y, c, (2) dt The equilibrium of (2) accumulates at y = and at y =, the equilibrium at is asymptotically stable since the derivative of cy( - at is -c which is negative. A more realistic mold of single species is y ym (, where M denotes an arbitrary function. (3) 65

4 2D-Volterra-Lotka modeling for 2 species 2- Volterra-Lotka System We consider the Predator-Prey Interaction of two species; the predator y and its prey x. The prey population is the total food supply for the predator at any given time, then: y a( x ) y, a,, which can be written in the form: y ( Cx D) y, C, D. We investigate the growth rate of the prey. At each small time, number of prey is eaten [6]. The prey species is assumed to have a constant per a capita food supply available sufficient to increase its population in the absence of predators. Therefore the prey is subject to differential equation of the form: x ( A B A, B. Now we arrive to the predator-prey equation of Volterra and Lotka: x ( A B ( 4) y ( Cx D) y, where x(t), y(t) are the prey and predator populations at time t respectively, A, D are the natural growth and decay coefficients, and A, B, C and D are positive numbers. The line terms A -Dy model the natural growth and decay of the prey and predator respectively, the quadratic terms Bxy and Cxy model the effects of interaction on the rates of change of the two species. This system has stationary points at P1 ( (, ) and at P ( ( D C, A B). 2 To determine the stability of the stationary points we linearize the system by taking the partial derivative of the system. The stability of points in linear system of equation 4 can be determined by finding the Eigenvalues of the matrix for the linear system at the points and applying the following theorem: 66

5 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya n Theorem (Principle of Linearzed Stabilit: Let F C 1 ( U, ) for n U with F ( P ). then the non-linear system a F(a) the following is true: 1- Re( ( F( P ))) P is asymptotically stable. 2- P is stable Re( ( F( P ))). According to the linear part of the vector field (4): A By Bx F ( x,. Cy Cx D The Jacobian at the point (, ): A F (, ), D The point (,) has Eigenvalues A and C, and we have A > and C > so the Eigenvalue is A, then the point (,) is saddle point, figure.1, so the system is unstable at the point (,). This is not surprising which stases the prey population increases in the absence of predators. Fig. 1 The Jacobian at the point ( D C, A B ) is: 67

6 2D-Volterra-Lotka modeling for 2 species BA B BD C A B CD C A F( D C, A B), C D BD C. AC B The trace of this point P is zero, the matrix has the Eigenvalues i AD 2, the real part of these values is zero, so the point could be either stable or instable. Another method of analysis is needed to find out more. Theorem 2: (Liapunov's theorem): Let ~ x W be an equilibrium for the system x f (x). Let V : U be a continuous function defined on a neighborhood U W of differentiable on U ~ x, such that: 1- V (x ~ ) and V (x) if x ~ x ; 2- V in U ~ x. Then ~ x is stable. Furthermore, if also 3- V inu ~ x, Then ~ x is asymptotically stable. Suppose that the Liapunov function is in the following form: v( F( G(. Then,. d v( v[ x( t), y( t)], dt df dg x y, dx dy df dg ( A B ( Cx D). dx dy Let v, then: df x y dg = constant. dx Cx D By A dy 68

7 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya If we put this constant = 1, we get: By integrating: df dx dg dy C D x, B A y. F( x) Gx D log G( By Alog y. Then: v ( Cx Dlog x By Alog y, x, y, By take the derivative with respect to x ( to y ) we get: d v y C D x, dx d vd C, A B, dx d v y B A y, dy d vd C, A B. dy Thus P is absolute minimum of v, then v is Liapunov function and P is 2 2 stable equilibrium point. Since v, therefore, every trajectory of the Volterra-Lotka is closed orbit except the equilibrium points, and then there is no limit cycle. The phase portrait is shown in Figure 2. 69

8 2D-Volterra-Lotka modeling for 2 species y Therefore, for any given initial population (x(), y()) with x() and y() other than P 2, the population of predator and prey will oscillate cyclically, no matter what the number of prey and predator are, neither species will die out, nor will grow indefinitely. On the other hand except for the state P 2, which is improbable, the population will not remain constant. Example: Assume that the initial number of foxes and rabbits are = 1 and, = 5 respectively, and that the coefficients, A 2, B. 1, D. 8, C. 2 are used to form the system of D. E.'s dx dt Fig. 2. Ax Bxy 2 x. 1xy, and x dy dt Cxy Dy. 2xy. 8y Solve the system of D. E.'s for x(t)and y(t) over the interval a t b. 7

9 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya Solution: We used the Rung-Kutta numerical method with initial data: x 5, y 1, a., b 12., n 6. The graph of the solution, over a larger time interval, the graph traces over itself, because the solution is periodic. Fig 3. Fig. 3. Fig. 4. The number of rabbits. 71

10 2D-Volterra-Lotka modeling for 2 species Fig. 5. The number of foxes. The periodic behavior is clearly present in the latter two graphs. Notice the phase lag for the number of prey. 3- Predator-prey equation of species with limited (Overcrowding) This model is more realistic from the biological point of view. Two species may compete for a resource in short supply. One can model for competitive interaction with overcrowding is: x ( A By x) y ( Cx D y, ( 4) where λ, μ, A, B, C and D are positive numbers, and x >, y >. We divide the positive quadratic Q (x >, y > ) into sectors by two lines L and M: L : A By x M : Cx D y. We get two cases according to the position of L and M. Case1: (they are not intersection) We can see that in the figure (6), 72

11 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya A B L M x x y y x y A D C We obtain two of equilibrium points P,, P A,. of (4) is A By 2x F( Cy F( P ), 1 Fig By, Cx D 2y The linear part A, D then the point P is a saddle as in the previous section. 1 The linearization at the point P : 2 A AB F ( P 2 ), AC ( D) then the point P is stable. Note, from the figure 4, one can be seen that, 2 the point P ( A, ) is globally asymptotically stable, because every 2 trajectory coming from the position quadrant ends to this point. 73

12 2D-Volterra-Lotka modeling for 2 species y x y x y x y A / λ Fig. 7. x Case 2: (L intersects M) In this case we add the point z ( which is the intersection of L and M to the equilibrium points P 1 and P 2, figure 8. Then the Jacobian at the point z is: A B y 2x B x F (. C y C x D 2y We have from L and M that: then A By x, Cx D y, 74

13 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya x B x F ( x,, C y y Therefore, z ( is asymptotically stable (stable node). It is not easy to know whether there is any limited cycle, but if we assume that there is a rectangle whose corners are (,), (P,), (P,Q) and (,Q) with: P A, Q A, and ( P, Q) M. Every trajectory at boundary points of either enters is positively invariant or is a part of the boundary. Therefore, is positively invariant. Figure 8. (,Q) (P,Q) Г L M z D/C Fig. 8. A/λ (P,) 4- Competing Species We consider two species x and y, which compete for some thing (common food supply for example).the equations of the growth of the two species are formed as: x M ( y N( y, ( 5) 75

14 2D-Volterra-Lotka modeling for 2 species where the growth rates M and N are C 1 functions of non-negative variables x and y with the following three assumptions: (1) If either species is increasing, the growth rate of the other goes down. Hence, M y N x, and. (2) If either population is very large, neither species can multiply. Hence there exists K, such that: M (, and N(, if K x or K y (3) In the absence either species, the other has positive growth rate up to a certain population, and negative growth rate beyond it. 1 By 1 and 3 each vertical line {x}r meats the set M () exactly once if x a and not at all if x > a. By (1) and the implicit function 1 theorem, μ is the growth rate of non-negative C map f:[,a] R, 1 such that f () a. Then M > below the curve μ and M < above it. Fig. 6-a. 1 In the same way, the set v N () is smooth curve of the form {(: 1 x=g(}, where g: [,b] R is non-negative C map with 1 g () b. The function N is positive to the left of v and negative to the right. Fig. 9-a, b. 76

15 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya μ ν M < N > M > N > Fig. 9-a. Fig. 9-b. 4.1-Case: 1 ( and do not intersect and below ) With this assumption we have three equilibrium points: (, ), ( a, ) and (, b ) The equilibrium point (, b ) is asymptotically stable, figure 1. 77

16 2D-Volterra-Lotka modeling for 2 species (,b) x y x y Fig Case: 2 ( and are intersections) x y (a,) is a finite set, then the Suppose that and are intersection and curves and, and the coordinate axes bound a finite number of connected open sets in the positive quadrant, these sets where x, and y, figure 11. (1) Points in the set (vertices). (2)Points on or on but not on both and not on the coordinate axes. (Ordinary points). (3)Points on the axes. A vertex is an equilibrium, (let us denote by p, q), and the other equilibrium points are (,), (a,), (,b). At an ordinary point B, the vector ( y ) is either vertical ( if ), or horizontal ( if v ), it points either into or out of B since μ has no 78

17 Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya vertical tangents and ν has no horizontal tangents. We call ω an inward or out ward point of B, accordingly. (,b) x y q x y x y p x y Figure 11. x y (a,) Conclusion The two-dimensional Volterra-Lotka system exhibits stable periodic behavior for all non-zero initial conditions. These trajectories run along closed paths around the stationary point ( D C, A B), which is non asymptotically stable. The other stationary point is at (; ), for which both populations are extinct. This point is instable. If only the predator population y is extinct, then the prey population x grows without bound. If only the prey population x is extinct, then the predator population y approaches extinction. Constants A;B;C;D are positive by dentition, so no alteration of the constants changes this behavior. Periodic stability is present for all possible combinations of variables. 79

18 2D-Volterra-Lotka modeling for 2 species References [1] Farkas H. and Noszyiczius Z., Generalized Lotka-Volterra of Two Dimensional Exploratory Cores and their Liapunov Bifurcations., J. Chem. Soc., Faraday Trans. 2, 1985, 81, [2] Farkas H. and Noszticzius Z., Two Dimensional Exploratory. 2. Global Analysis of Lotka-Volterra-Brusselator (LVB) Model. Act Physica., Hungary 66, 23-22, [3] Guckenheimer J. and Holmes P., Oscillations, Dynamical Systems and Bifurcations of Vector Fields. Springer-Verlag, New York, [4] Hirsch M. and Smale S., differential Equations. Dynamical System, and Linear Algebra, Academic Press, INC., London, [5] Krasnov M. L., Kiselyov A. I and Makarenko G. I., A Book of Problems in ordinary Differential Equations., Mir Publishers, Mosco, [6] Wiggins S., Introduction to Applied Nonlinear Dynamical System and Chaos., Springer-Verlag. New York, Inc.,

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

Getting Started With The Predator - Prey Model: Nullclines

Getting 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 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

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

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:

More information

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

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

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

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

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

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt

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

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

APPM 2360 Lab #3: The Predator Prey Model

APPM 2360 Lab #3: The Predator Prey Model APPM 2360 Lab #3: The Predator Prey Model 1 Instructions Labs may be done in groups of 3 or less. One report must be turned in for each group and must be in PDF format. Labs must include each student s:

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

Lesson 9: Predator-Prey and ode45

Lesson 9: Predator-Prey and ode45 Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and

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 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

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

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

APPPHYS217 Tuesday 25 May 2010

APPPHYS217 Tuesday 25 May 2010 APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

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

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

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

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

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

Coordinate Curves for Trajectories

Coordinate Curves for Trajectories 43 The material on linearizations and Jacobian matrices developed in the last chapter certainly expanded our ability to deal with nonlinear systems of differential equations Unfortunately, those tools

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

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of

More information

Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo

Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo 2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical

More information

NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE

NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 44-50 ISSN: 1823-8556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD

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

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

Stability of Dynamical systems

Stability of Dynamical systems Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u

More information

Motivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt

Motivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt Motivation and Goals Modelling with ODEs 24.10.01 Motivation: Ordinary Differential Equations (ODEs) are very important in all branches of Science and Engineering ODEs form the basis for the simulation

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

Symmetry Properties of Confined Convective States

Symmetry Properties of Confined Convective States Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

Stability and bifurcation in a two species predator-prey model with quintic interactions

Stability and bifurcation in a two species predator-prey model with quintic interactions Chaotic Modeling and Simulation (CMSIM) 4: 631 635, 2013 Stability and bifurcation in a two species predator-prey model with quintic interactions I. Kusbeyzi Aybar 1 and I. acinliyan 2 1 Department of

More information

The Liapunov Method for Determining Stability (DRAFT)

The Liapunov Method for Determining Stability (DRAFT) 44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level

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

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

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

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

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

Nonlinear Autonomous Dynamical systems of two dimensions. Part A

Nonlinear Autonomous Dynamical systems of two dimensions. Part A Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous

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

POPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL. If they co-exist in the same environment:

POPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL. If they co-exist in the same environment: POPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL Next logical step: consider dynamics of more than one species. We start with models of 2 interacting species. We consider,

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

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix

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

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

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd

More information

MATH 250 Homework 4: Due May 4, 2017

MATH 250 Homework 4: Due May 4, 2017 Due May 4, 17 Answer the following questions to the best of your ability. Solutions should be typed. Any plots or graphs should be included with the question (please include the questions in your solutions).

More information

Copyright (c) 2006 Warren Weckesser

Copyright (c) 2006 Warren Weckesser 2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and

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

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

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

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

1 Functions and Graphs

1 Functions and Graphs 1 Functions and Graphs 1.1 Functions Cartesian Coordinate System A Cartesian or rectangular coordinate system is formed by the intersection of a horizontal real number line, usually called the x axis,

More information

LECTURE 8: DYNAMICAL SYSTEMS 7

LECTURE 8: DYNAMICAL SYSTEMS 7 15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin

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

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

28. 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) 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 information

8 Ecosystem stability

8 Ecosystem stability 8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single

More information

Lotka-Volterra Models Nizar Ezroura M53

Lotka-Volterra Models Nizar Ezroura M53 Lotka-Volterra Models Nizar Ezroura M53 The Lotka-Volterra equations are a pair of coupled first-order ODEs that are used to describe the evolution of two elements under some mutual interaction pattern.

More information

MA 138 Calculus 2 with Life Science Applications Handout

MA 138 Calculus 2 with Life Science Applications Handout .. MA 138 Calculus 2 with Life Science Applications Handout Alberto Corso alberto.corso@uky.edu Department of Mathematics University of Kentucky February 17, 2017 . Example 4 (Lotka-Volterra Predator-Prey

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

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems

More information

Dynamical Systems in Biology

Dynamical Systems in Biology Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems

More information

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0

More information

Solutions to Math 53 Math 53 Practice Final

Solutions to Math 53 Math 53 Practice Final Solutions to Math 5 Math 5 Practice Final 20 points Consider the initial value problem y t 4yt = te t with y 0 = and y0 = 0 a 8 points Find the Laplace transform of the solution of this IVP b 8 points

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

Predator - Prey Model Trajectories and the nonlinear conservation law

Predator - 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

Preservation of local dynamics when applying central difference methods: application to SIR model

Preservation of local dynamics when applying central difference methods: application to SIR model Journal of Difference Equations and Applications, Vol., No. 4, April 2007, 40 Preservation of local dynamics when applying central difference methods application to SIR model LIH-ING W. ROEGER* and ROGER

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

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

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

Applied Dynamical Systems

Applied Dynamical Systems Applied Dynamical Systems Recommended Reading: (1) Morris W. Hirsch, Stephen Smale, and Robert L. Devaney. Differential equations, dynamical systems, and an introduction to chaos. Elsevier/Academic Press,

More information

2015 Holl ISU MSM Ames, Iowa. A Few Good ODEs: An Introduction to Modeling and Computation

2015 Holl ISU MSM Ames, Iowa. A Few Good ODEs: An Introduction to Modeling and Computation 2015 Holl Mini-Conference @ ISU MSM Ames, Iowa A Few Good ODEs: An Introduction to Modeling and Computation James A. Rossmanith Department of Mathematics Iowa State University June 20 th, 2015 J.A. Rossmanith

More information

Lyapunov functions and stability problems

Lyapunov functions and stability problems Lyapunov functions and stability problems Gunnar Söderbacka, Workshop Ghana, 29.5-10.5, 2013 1 Introduction In these notes we explain the power of Lyapunov functions in determining stability of equilibria

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

Phenomenon: Canadian lynx and snowshoe hares

Phenomenon: Canadian lynx and snowshoe hares Outline Outline of Topics Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Phenomenon: Canadian lynx and snowshoe hares All began with

More information

1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.

1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos. Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based

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

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

STUDY OF THE DYNAMICAL MODEL OF HIV

STUDY OF THE DYNAMICAL MODEL OF HIV STUDY OF THE DYNAMICAL MODEL OF HIV M.A. Lapshova, E.A. Shchepakina Samara National Research University, Samara, Russia Abstract. The paper is devoted to the study of the dynamical model of HIV. An application

More information

Introduction to the Phase Plane

Introduction to the Phase Plane Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called

More information

One Dimensional Dynamical Systems

One Dimensional Dynamical Systems 16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

More information

Final Exam December 20, 2011

Final Exam December 20, 2011 Final Exam December 20, 2011 Math 420 - Ordinary Differential Equations No credit will be given for answers without mathematical or logical justification. Simplify answers as much as possible. Leave solutions

More information

GLOBAL DYNAMICS OF A LOTKA VOLTERRA MODEL WITH TWO PREDATORS COMPETING FOR ONE PREY JAUME LLIBRE AND DONGMEI XIAO

GLOBAL DYNAMICS OF A LOTKA VOLTERRA MODEL WITH TWO PREDATORS COMPETING FOR ONE PREY JAUME LLIBRE AND DONGMEI XIAO This is a preprint of: Global dynamics of a Lotka-Volterra model with two predators competing for one prey, Jaume Llibre, Dongmei Xiao, SIAM J Appl Math, vol 742), 434 453, 214 DOI: [11137/1392397] GLOBAL

More information

Workshop on Theoretical Ecology and Global Change March 2009

Workshop on Theoretical Ecology and Global Change March 2009 2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological

More information

Linearization of Differential Equation Models

Linearization of Differential Equation Models Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking

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

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

B. Differential Equations A differential equation is an equation of the form

B. Differential Equations A differential equation is an equation of the form B Differential Equations A differential equation is an equation of the form ( n) F[ t; x(, xʹ (, x ʹ ʹ (, x ( ; α] = 0 dx d x ( n) d x where x ʹ ( =, x ʹ ʹ ( =,, x ( = n A differential equation describes

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

Calculus C (ordinary differential equations)

Calculus C (ordinary differential equations) Calculus C (ordinary differential equations) Lesson 1: Ordinary first order differential equations Separable differential equations Autonomous differential equations Logistic growth Frank Hansen Institute

More information

Stability Analysis of Predator- Prey Models via the Liapunov Method

Stability Analysis of Predator- Prey Models via the Liapunov Method Stability Analysis of Predator- Prey Models via the Liapunov Method Gatto, M. and Rinaldi, S. IIASA Research Memorandum October 1975 Gatto, M. and Rinaldi, S. (1975) Stability Analysis of Predator-Prey

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