2DVolterraLotka Modeling For 2 Species


 Naomi Pitts
 1 years ago
 Views:
Transcription
1 Majalat AlUlum AlInsaniya wat  Tatbiqiya 2DVolterraLotka 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 multispecies 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 predatorprey system with continuous of nonlinear predatorprey 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 nonlinear 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 predatorprey ecology are investigated [2], [3]. The object is to find out the longrun qualitative behavior of these dynamical systems. 1 Alhashmi Darah, PO Box 675 Zliten Libya, << >>. 63
2 2DVolterraLotka 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 birthrate 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 AlUlum AlInsaniya 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 2DVolterraLotka modeling for 2 species 2 VolterraLotka System We consider the PredatorPrey 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 predatorprey 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 AlUlum AlInsaniya wat  Tatbiqiya n Theorem (Principle of Linearzed Stabilit: Let F C 1 ( U, ) for n U with F ( P ). then the nonlinear 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 2DVolterraLotka 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 AlUlum AlInsaniya 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 VolterraLotka is closed orbit except the equilibrium points, and then there is no limit cycle. The phase portrait is shown in Figure 2. 69
8 2DVolterraLotka 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 AlUlum AlInsaniya wat  Tatbiqiya Solution: We used the RungKutta 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 2DVolterraLotka 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 Predatorprey 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 AlUlum AlInsaniya 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 2DVolterraLotka 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 AlUlum AlInsaniya 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 2DVolterraLotka modeling for 2 species where the growth rates M and N are C 1 functions of nonnegative 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 nonnegative C map f:[,a] R, 1 such that f () a. Then M > below the curve μ and M < above it. Fig. 6a. 1 In the same way, the set v N () is smooth curve of the form {(: 1 x=g(}, where g: [,b] R is nonnegative C map with 1 g () b. The function N is positive to the left of v and negative to the right. Fig. 9a, b. 76
15 Majalat AlUlum AlInsaniya wat  Tatbiqiya μ ν M < N > M > N > Fig. 9a. Fig. 9b. 4.1Case: 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 2DVolterraLotka 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 AlUlum AlInsaniya 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 twodimensional VolterraLotka system exhibits stable periodic behavior for all nonzero 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 2DVolterraLotka modeling for 2 species References [1] Farkas H. and Noszyiczius Z., Generalized LotkaVolterra 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 LotkaVolterraBrusselator (LVB) Model. Act Physica., Hungary 66, 2322, [3] Guckenheimer J. and Holmes P., Oscillations, Dynamical Systems and Bifurcations of Vector Fields. SpringerVerlag, 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., SpringerVerlag. New York, Inc.,
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 informationGetting Started With The Predator  Prey Model: Nullclines
Getting Started With The Predator  Prey Model: Nullclines James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline The Predator
More informationPhysics: springmass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: springmass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationLab 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 informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS114.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 informationProblem 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 informationLotka Volterra PredatorPrey Model with a Predating Scavenger
Lotka Volterra PredatorPrey 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 informationMATH 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 informationIntroduction to Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology RabbitFox System Trout System Trout System
More informationFundamentals of Dynamical Systems / DiscreteTime Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / DiscreteTime Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationEG4321/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 informationA 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 informationEAD 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 informationAPPM 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 informationLesson 9: PredatorPrey and ode45
Lesson 9: PredatorPrey 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 informationGerardo Zavala. Math 388. PredatorPrey Models
Gerardo Zavala Math 388 PredatorPrey 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 informationNonlinear 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 informationMath 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 informationLecture 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 MassachusettsDartmouth Coupled ODEs: Species
More information154 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 informationAPPPHYS217 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 (SpringerVerlag
More informationComplex 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 informationCoordinate 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 informationChapter 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é (18541912)
More informationStable Coexistence of a PredatorPrey model with a Scavenger
ANTON DE KOM UNIVERSITEIT VAN SURINAME INSTITUTE FOR GRADUATE STUDIES AND RESEARCH Stable Coexistence of a PredatorPrey model with a Scavenger Thesis in partial fulfillment of the requirements for the
More information6.3. Nonlinear Systems of Equations
G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: TwoDimensional Nonlinear Systems. ritical Points and Linearization. The HartmanGrobman Theorem. Part Two: ompeting
More informationE209A: 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 informationDynamics 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 informationSection 5. Graphing Systems
Section 5. Graphing Systems 5A. The Phase Plane 5A1. Find the critical points of each of the following nonlinear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A2. Write each
More informationMath 312 Lecture Notes Linear Twodimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Twodimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
20081107 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 informationThe 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 informationNUMERICAL SIMULATION DYNAMICAL MODEL OF THREESPECIES FOOD CHAIN WITH LOTKAVOLTERRA LINEAR FUNCTIONAL RESPONSE
Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 4450 ISSN: 18238556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREESPECIES FOOD
More information1 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 informationStability 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 informationMotivation 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 informationContinuous 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 informationPOPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL. If they coexist 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 informationB5.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 informationSymmetry 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 informationStability and bifurcation in a two species predatorprey model with quintic interactions
Chaotic Modeling and Simulation (CMSIM) 4: 631 635, 2013 Stability and bifurcation in a two species predatorprey model with quintic interactions I. Kusbeyzi Aybar 1 and I. acinliyan 2 1 Department of
More informationLogistic Map, Euler & RungeKutta Method and LotkaVolterra Equations
Logistic Map, Euler & RungeKutta Method and LotkaVolterra 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 informationM469, 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:005:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationBIFURCATION 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 informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 3546 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 informationWe have two possible solutions (intersections of nullclines. 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. LotkaVolterra predatorprey model: phase plane analysis. Earlier we introduced the system of equations for prey
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver YueXian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationNonlinear 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 informationOrdinary 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 informationMath 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 informationBifurcation and Stability Analysis of a Preypredator System with a Reserved Area
ISSN 746733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 859 Bifurcation and Stability Analysis of a Preypredator System with a Reserved Area Debasis Mukherjee Department
More informationChapter 4: Firstorder differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: Firstorder differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: Firstorder differential equations Phase portrait Singular point Separatrix
More information3.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 informationDynamical 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 informationMATH 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 informationDynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBASIJENS Vol:13 No:05 55 Dynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationCopyright (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 informationDynamical 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 informationFIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRSTORDER 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 information1 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 informationLECTURE 8: DYNAMICAL SYSTEMS 7
15382 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 informationANSWERS 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 information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More information8 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 informationA Discrete Model of Three Species Prey Predator System
ISSN(Online): 398753 ISSN (Print): 347670 (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 informationMA 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 (LotkaVolterra PredatorPrey
More informationModels 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 Predatorprey models are used to show the intricate
More informationMath 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 informationLotkaVolterra Models Nizar Ezroura M53
LotkaVolterra Models Nizar Ezroura M53 The LotkaVolterra equations are a pair of coupled firstorder ODEs that are used to describe the evolution of two elements under some mutual interaction pattern.
More informationAn Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems
An Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The HartmanGrobman and PoincaréBendixon Theorems
More informationDynamical 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 informationDynamical 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 informationIntroduction 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 informationLyapunov functions and stability problems
Lyapunov functions and stability problems Gunnar Söderbacka, Workshop Ghana, 29.510.5, 2013 1 Introduction In these notes we explain the power of Lyapunov functions in determining stability of equilibria
More informationSolutions 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 informationIs 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 informationPredator  Prey Model Trajectories and the nonlinear conservation law
Predator  Prey Model Trajectories and the nonlinear conservation law James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University October 28, 2013 Outline
More informationHalf 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 informationLinearization 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 informationPreservation 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 LIHING W. ROEGER* and ROGER
More informationChapter 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 firstorder system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationStability 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 PredatorPrey
More informationApplied 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 information2015 Holl ISU MSM Ames, Iowa. A Few Good ODEs: An Introduction to Modeling and Computation
2015 Holl MiniConference @ 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 informationNonLinear Models. NonLinear Models Cont d
Focus on more sophistiated interaction models between systems. These lead to nonlinear, rather than linear, DEs; often not soluble exactly in analytical form so use PhasePlane Analysis. This is a method
More informationPhenomenon: Canadian lynx and snowshoe hares
Outline Outline of Topics Shan He School for Computational Science University of Birmingham Module 0623836: Computational Modelling with MATLAB Phenomenon: Canadian lynx and snowshoe hares All began with
More information1.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 informationTHE 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 informationSTUDY 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 informationGLOBAL 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 LotkaVolterra 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 informationENGI Linear Approximation (2) Page Linear Approximation to a System of NonLinear ODEs (2)
ENGI 940 4.06  Linear Approximation () Page 4. 4.06 Linear Approximation to a System of NonLinear ODEs () From sections 4.0 and 4.0, the nonlinear system dx dy = x = P( x, y), = y = Q( x, y) () with
More informationWorkshop on Theoretical Ecology and Global Change March 2009
20223 Workshop on Theoretical Ecology and Global Change 218 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with onedimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
More informationA ProducerConsumer Model With Stoichiometry
A ProducerConsumer 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 informationFinal 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 informationInteractions. 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 informationB. 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 informationAutonomous systems. Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous.
Autonomous equations Autonomous systems Ordinary differential equations which do not contain the independent variable explicitly are said to be autonomous. i f i(x 1, x 2,..., x n ) for i 1,..., n As you
More information2.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 informationCalculus 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