Dynamics of a Population Model Controlling the Spread of Plague in Prairie Dogs
|
|
- Cody White
- 5 years ago
- Views:
Transcription
1 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 Catalin.Georgescu@usd.edu Dan Van eursem The University of South Dakota Department of Mathematical Sciences 414 East Clark Street, Vermillion, SD USA dpeursem@usd.edu Abstract: The infestation of prairie dog population with the deadly Yersinia estis has been recently the object of investigation of numerous biologists due to the potentially high risk of contamination of the human population. In the present paper we construct and analyze a model describing the interaction between the prairie dog population as the host population) and the flea population as the vector population) and derive the existence under certain conditions impose to the main parameter) of an equilibrium that plays the role of a global attractor. Key Words: Stability of Equilibria, Global attractor, Bifurcation. 1 Background. In recent years, prairie dog populations have been hit hard by the plague caused by Yersinia pestis. The literature suggests that total colonies can be wiped out upon infestation and the disease is rapidly spread from one colony to the next see [1],[2],[6]). There are reasons to expect that the spread of plague from face to face contact or through respiratory methods is not the predominate mode of transferring the disease as Hoogland et al. have reported see [5]) that upon killing fleas with the insecticide yraperm the spread of the disease was stopped for a time. The agreed upon mode of transmission of the vector is through infected fleas hosting on the prairie dogs or other mammals. Seery et al. and Hoogland et. al have shown that by using insecticides on an infected colony one can stop the spread of an epizootic and can be used to save the colony. We present a mathematical model that will allow one to show analytically what kind of control one needs on the fleas in order to prevent an epizootic once the plague has been manifested in a colony. The model shows that when one starts with healthy and infected prairie dogs and fleas, there are three equilibrium points that can develop. We will consider α, the ratio between the normal death rates of flea population and prairie dog population, as the working parameter of the problem. We will show that if this parameter stays smaller than a certain value α 0, the dynamics of the system is governed by three equilibria: there is the case where healthy and infected prairie dogs and fleas can coexist and it is an asymptotically stable equilibrium, there is the case where all infected fleas and infected prairie dogs are wiped out, and there is the case where only the healthy prairie dogs exist, these last two equilibria being unstable. If the parameter exceeds the bifurcation point α 0, one witness a collapse of equilibria to only healthy prairie dogs existing and the change in stability of the last two mentioned above equilibria. 2 The Origins of the Model. The model will consist of a typical host-vector model with the prairie dogs acting as the host population and the flea population being the vector population. Let be the total prairie dog population, S being the number of infected or sick prairie dogs, F be the total flea population, and D be the number of diseased or infected fleas. The governing equations are: ) d dt = r p 1 kp d r, 1) df dt = r f 1 F ) F d f F, 2) m ) ds S dt = t p D c 1 d r S, 3) dd dt = t S f F D) c 2d f D. 4) ISSN:
2 The overall populations of the prairie dogs and fleas are governed by a logistic growth model where r p and r f are the growth rates of the prairie dog and flea populations respectively, k p, is the carrying capacity of the prairie dog population, m is the maximum number of fleas that would be found on a prairie dog, thus giving m as the carrying capacity for the fleas, and d r and d f being the normal death rate of the overall prairie dog and flea populations respectively. The number of sick prairies dogs grows at a rate that is proportional to the product of the proportion of healthy prairie dogs and the number of diseased fleas with t p being the proportionality constant. This is similar to the interactions among predators and prey in the standard Lotka-Volterra equationssee [7]). We make a simplifying assumption here that with the infected prairie dogs, the deathrate c 1 d r is a constant multiple of the over all death rate. The infected flea population is very similar in that it grows at a rate proportional to the product of the proportion of diseased prairie dogs and the number of healthy fleas. The proportionality constant is t f and again we make the simplifying assumption that the death rate, c 2 d f, is a constant multiple of the overall death rate for the total flea population.this system is very similar to the governing equations of [3] where they were looking at ticks as being the vector population on deers. There are however some simplifying assumptions we made here that will allow for full analytical analysis of the system. One simplification involves the elimination of the interaction between the infected and healthy hosts due to the evidence that face to face and airborne transmissions seemed to be minor as mentioned earlier. We also eliminate the interaction between infected and healthy fleas since a lack of competitiveness has been noticed at this type of population. In scaling the equations we will let = /k p, t = t d r, F = F/mk p ), S = S/k p, and D = D/mk p ). These are the natural scalings where we scale the populations by their carrying capacity and time is scaled by the reciprocal of the prairie dog death rate. Introducing the positive dimensionless parameters of k 1 = r p /d r, k 2 = r f /d r, k 3 = mt p /d r, k 4 = t f /d r, and α = d f /d r, upon dropping the bars the scaled equations are as follows: d dt = k 11 ), 5) df dt = k 21 F )F αf, 6) ds dt = k 31 S )D c 1S, 7) dd dt = k S 4 F D) c 2αD. 8) We will use α as a control parameter for the system as it is essentially the death rate of the vector population. We will investigate controlling the population by varying this parameter. It is noted that the only physical solutions are for, F, S, and D to all be positive. The parameter α will play a particulary important role, since bifurcation occurs with respect to this parameter and it is essentially our control on how effective we are in killing the fleas. 3 Dynamical Analysis of the Model. Since the first two equations 5) and 6) of the original system are decoupled from the other two, it is natural to start with the analysis of the phase portrait in the first quadrant of the plane, F ). An immediate computation shows that in order to have equilibria in the first quadrant, the following two conditions have to be true: k 1 > 1 9) 0 < α < k 2, 10) and in this case the two nonzero equilibria will be, 0) and, F 0 ), where: ) = 1 1k1, F 0 = 1 α ). 11) k 2 Linearizing the system at an arbitrary point, F ) leads to the Jacobian: J, F ) = k 1 1) 2k 1 k 2 F 2 2 k 2 α) 2k 2F Under the assumptions 9) and 10) the equilibrium, F 0 ) will be a sink since the eigenvalues of the linearized system will be λ 1 = 1 k 1 < 0 and λ 2 = α k 2 < 0 and the equilibrium, 0) is a saddle with λ 1 = 1 k 1 < 0 and λ 2 = k 2 α > 0. Notice also that F = 0 and = are two invariant lines and on the line F = 0 we have d dt > 0 if < ISSN:
3 F Figure 1: hase ortrait in F-plane,, F 0. The values k 1 = 2, k 2 = 1 and α = 0.5 were used. and d dt < 0 if >, so the phase portrait will look as in Figure 1. Remark 1. When α = k 2, the two equilibria will coincide with, 0) which becomes a sink in fact an attractor for all orbits in the first quadrant since df < 0 and the invariant line F = 0 stops the orbits dt from exiting the first quadrant). For α > k 2 the equilibrium, F 0 ) moves out of the first quadrant with, 0) continuing to be a sink. In order to analyze the general case, we will show that if one restricts themselves to the plane =, F = F 0 in the space of coordinates, F, S, D), all orbits in the first quadrant) starting in this plane will converge to an equilibrium S 0, D 0 ) in this plane. Since projections on the, F ) plane of all orbits except those starting with F = 0) converge to, F 0 ) as we have seen from above, we will have that almost all orbits will converge to a global attractor due to the following result: Lemma 2. Consider the following system of differential equations: dx dt = fx) dy = gx, y), dt where x R n and y R m, and the functions f and g are considered of class C 1. Assume that there is 1.0 an x 0 such that fx 0 ) = 0 and any solution xt) of the system dx/dt = fx) has the property that lim t xt) = x 0. Assume also that if x 0, yt)) is a solution of the above system, then lim t yt) = y 0, for some y 0 R m. Then for any solution xt), yt)) of the system we have that lim t xt), yt)) = x 0, y 0 ). We will apply this lemma for x =, F ) and y = S, D). We now proceed to analyze the general system. First we determine the equilibria. Consider = in ds/dt = dd/dt = 0 equations and eliminate the term S D from the last two equations with =. This leads to: D = βs, where β = c 1 k 4 + k 3 k 4 F k 3 k 4 + k 3 c 2 α 12) Notice that 12) implies that S = 0 is equivalent to D = 0 due to positivity of coefficients. The implication D = 0 S = 0 is particulary important from observations in the field since this shows no transfer of disease from prairie dog to prairie dog, as prior research suggests. If S 0, then we obtain the equilibrium: S 0 = 1 c ) 1, D 0 = β 1 c ) 1 βk 3 βk 3 13) Hence we distinguish two cases: Case 1. If F 0 = 0, then 12) shows that 1 c 1 = βk 3 αc 2 /k 4 and hence the equilibrium S 0, D 0 ) has both components negative so it is not of interest for our model. Thus in this case in the hyperplane =, F = 0 only the equilibrium S = 0, D = 0) is of interest. The linearization at this point in the above hyperplane gives the Jacobian: J0, 0) = c1 k 3 0 c 2 α and hence the two eigenvalues are negative which shows this is a sink. Since the hyperplane F = 0 is invariant to the flow, in the 4 dimensional picture this equilibrium will have only one repelling direction, the line connecting, 0) with, F 0 ) in the S = 0, D = 0 plane. This reflects the fact that if there are no fleas, the prairie dog population tends to a constant level. Case 2. t) =, F t) = F 0. This is the interesting ) ISSN:
4 case since the plane t) =, F t) = F 0 is invariant and an attractor. We will denote by f 3 and f 4 the third and fourth components of the vector field in this plane, hence: f 3 S, D) = k 3 D k 3 SD c 1 S f 4 S, D) = k 4F 0 S k 4 SD c 2 αd Note that since the divergence f 3 S + f 4 is negative D in the first quadrant, the system will have no periodic orbits in the plane t) =, F t) = F 0. Before we analyze the nature of equilibria an important observation is necessary. Notice that from 13) we must have 1 c 1 /βk 3 > 0. Using 12) we obtain: 1 c 1 = 1 c 1 k 3 k 4 + c 2 α) βk 3 k 3 k 4 c 1 + k 3 1 α/k 2 )) = k 3k 4 αc 1 c 2 + k 3 k 4 /k 2 ) k 4 c 1 + k 3 1 α/k 2 )) so the parameter α should additionally satisfy: k 3 k 4 0 < α < α 0 =: = 14) k 3 k 4 /k 2 + c 1 c 2 = k c 1c 2 k 2 k 3 k 4 < k 2. We can see that if α satisfies 14) it will automatically satisfy 10). Unless otherwise stated, we will assume from now on that α < α 0. There are two equilibria in this case: 0, 0) and S 0, D 0 ). For the first one the Jacobian is: J0, 0) = c 1 k 3 k 4 F 0 15) c 2 α The trace of this matrix is obviously negative and the determinant, using 11), is D = αc 1 c 2 + k 3 k 4 /k 2 ) k 3 k 4 and hence from 14) we have D < 0 so this equilibrium is a saddle for the linear case and from Hartman-Grobman theorem it is a saddle for the original nonlinear) system. At the other equilibrium S 0, D 0 ), the Jacobian will be after some rearrangements): JS 0, D 0 ) = ) k3 D 0 +c 1 )/ k 3 S 0 )/ k 4 F 0 D 0 )/ k 4 S 0 +c 2 α )/ Again the trace of this matrix is obviously negative, but the sign of the determinant is less obvious. We clearly have: = {k 3 D 0 + c 1 )k 4 S 0 + c 2 α ) k 3 k 4 S 0 )F 0 D 0 )} 16) Using 11) and 13)we will express in 16) D 0, S 0 and in terms of. An immediate computation leads to: = βk 3 k 4 +c 2 k 3 α) 1 c 1 k 4 c 1k 4 α + c2 1 k ) 4 β k 2 k 3 This implies that: > 0 β 2 > c 1 + c 1 α 1k 4 k 3 k 2 = k 3 k 4 + αc 2 = c 1k 4 k 2 k3 2 k2k 3 αk 3 + c 1 k 2 17) k 4 + αc 2 The expression of β when F = F 0 is derived from 12): β = k 4 k 2 c 1 + k 2 k 3 αk 3 k 2 k 3 k 4 + c 2 α After some rearrangements, the inequality in 17) becomes: k 3 k 4 α < = α 0. k 3 k 4 /k 2 + c 1 c 2 This shows that in the plane = and F = F 0 the equilibrium S = S 0 and D = D 0 is a sink if α < α 0 which is also the condition for which this nontrivial) equilibrium exists. We will now show that in the plane =, F = F 0 there is a connecting orbit between the two equilibria 0, 0) and S 0, D 0 ). Consider λ the positive eigenvalue of the linearized system at the equilibrium 0, 0). If 1, y) gives the direction of the unstable space, we will have that J0, 0) 1, y) t = λ1, y) t. The first equation of this system yields: c 1 + k 3 y = λ which implies y > 0. Hence the unstable manifold which will be tangent at the origin to this direction will enter in the first quadrant. The flow enters into a box situated in the first quadrant having two sides along the coordinate axes as it can also be seen from figure 2), so the ω-limit set of this unstable trajectories is inside this box. From the oincaré-bendixon criterion, since the system has no periodic orbits as we proved above, this orbit will have the equilibrium S 0, D 0 ) as a limit point. The arguments above can be summarized by the following: ISSN:
5 d S Figure 2: hase ortrait in SD-plane, S, D 0. The values k 1 = 2, k 2 = 1, α = 0.5 and k 3 = k 4 = c 1 = c 2 = 1 were used. roposition 3. Let us consider k 1 > 1. If 0 < α < α 0, then there are only three equilibria of interest: E 1 =, 0, 0, 0), E 2 =, F 0, 0, 0) and E 3 =, F 0, S 0, D 0 ). Any solution of interest having initial condition, F 0, S 0, D 0 ) with F 0 strictly positive either tends to the equilibrium E 2 if S 0 = D 0 = 0) and hence if initially there are no diseased fleas or prairie dogs, the two populations coexist or to the equilibrium E 3 if S 0 > 0 and D 0 > 0). Therefore, except on a set of measure zero, solutions of interest tend to E 3, the equilibrium which plays the role of a global attractor. If F 0 = 0 then all solutions will tend toward E 1. For α = α 0 the equilibrium E 3 coincides with E 2 and all solutions of interest with F 0 > 0 will tend to this equilibrium. For α 0 < α < k 2 the equilibrium E 3 will have the last two coordinates negative so it is no longer of interest. In fact, the components f 3 and f 4 of the vector field will become negative for S > 0 and D > 0 so solutions of interest with F 0 > 0 will continue to approach E 2. They can not go below S = 0 and D = 0 since they would intersect then the solutions in the, F ) plane contradicting the uniqueness theorem which obviously holds for this system. When α gets closer to k 2 the equilibrium E 2 gets closer to E 1. When α = k 2, as explained in Remark 1, there will be only one equilibrium with all solutions of interest going toward E 1 = E 2 = E 3 for F 0 > 0. Solutions in the, F ) plane also go toward this equilibrium. We conclude with the following proposition. roposition 4. Consider the system 1) with k 1 > 1. If α [α 0, k 1 ], all solutions with initial condition having F 0 > 0 tend toward E 2 and those with F 0 = 0 tend toward E 1. In other words, the populations St) and Dt) become eventually depleted. 4 Conclusion The model presented in this paper generates a four dimensional system of differential equations depending on several parameters. Rescaling and fixing some proportionality constants leads to a system with one parameter and the bifurcation problem is described. The idea behind managing such a system was to take advantage of the fact that the first two equations are decoupled. We proved the existence of an attractor for all solutions starting in the first quadrant for the first two equations and using this we then prove the existence of an equilibrium that is a global attractor for solutions having positive initial conditions. This shows that for k 1 > 1 and 0 < α < α 0, we will have coexistence with healthy prairie dogs and fleas if there are no initial diseased fleas or prairie dogs. With initial diseased prairie dogs and fleas both prairie dogs and fleas will coexist with disease present in both populations. For α 0 α < k 2 only healthy prairie dogs and fleas will exist. Finally for α k 2 only healthy prairie dogs will exist. Acknowledgements: The research was partially supported by a University of South Dakota Travel Grant. References [1] Barnes,A., A Review of plague and its relevance to prairie dog population and the black-footed ferret, Washington, DC: United States Fish and Wildlife Service, [2] Cully,F.,Williams,S., Interspecific comparisons of sylvatic plague in prairie dogs, Journal of Mammology, 82, , [3] Gaff,H.D., Gross,L.J., Modeling Tick-Borne Disease: A Metapopulation Model, Bulletin of Mathematical Biology, 69, , [4] Guckenheimer,J., Holmes,., Nonlinear Oscillations, Dynamical Systems, and Bifurcation of Vector Fields, Springer-Verlag, ISSN:
6 [5] Hoogland,J.,Davis,S.,Benson- Amram,S.,Labruna,D.,Goossens,B.,Hoogland,M., yraperm Kills Fleas and Halts lague Among Utah rairie Dogs, The Southwestern Naturalist, 493), ; [6] Seery,D.,Biggins,D.,Montenieri,J., Enscore,R.,Tanda,D.,Gage,K., Treatment of Blacked-Tailed rairie Dog burrows with Deltamethrin to Control Fleas Insecta: Siphonaptera) and lague, Journal of Medical Entomology, Vol.40, no.5, [7] Takeuchi,Y.,Global Dynamical roperties of Lotka-Volterra Systems Singapore, River Edge NJ: World Scientific, ISSN:
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 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 information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More 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 informationPhysics: 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 informationAn 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 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 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:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationNonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.
Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical
More informationThe 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 informationApplications in Biology
11 Applications in Biology In this chapter we make use of the techniques developed in the previous few chapters to examine some nonlinear systems that have been used as mathematical models for a variety
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 informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationWe 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 information6.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 informationMA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total
MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all
More information3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
More informationOutline. Learning Objectives. References. Lecture 2: Second-order Systems
Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!
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 informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
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 Rabbit-Fox System Trout System Trout System
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 informationBifurcation 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 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 informationDynamics 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 informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
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 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 informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationFundamentals 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 informationMath 312 Lecture Notes Linearization
Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior
More informationResearch Article Global Dynamics of a Competitive System of Rational Difference Equations in the Plane
Hindawi Publishing Corporation Advances in Difference Equations Volume 009 Article ID 1380 30 pages doi:101155/009/1380 Research Article Global Dynamics of a Competitive System of Rational Difference Equations
More informationEpidemics in Two Competing Species
Epidemics in Two Competing Species Litao Han 1 School of Information, Renmin University of China, Beijing, 100872 P. R. China Andrea Pugliese 2 Department of Mathematics, University of Trento, Trento,
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 informationA 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 informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More 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 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 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 informationLotka 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 informationMath Lecture 46
Math 2280 - Lecture 46 Dylan Zwick Fall 2013 Today we re going to use the tools we ve developed in the last two lectures to analyze some systems of nonlinear differential equations that arise in simple
More informationStable 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 informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More 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 (Springer-Verlag
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 information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More 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 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 Massachusetts-Dartmouth Coupled ODEs: Species
More informationPOPULATION 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 informationLogistic 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 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 LIH-ING W. ROEGER* and ROGER
More information10 Back to planar nonlinear systems
10 Back to planar nonlinear sstems 10.1 Near the equilibria Recall that I started talking about the Lotka Volterra model as a motivation to stud sstems of two first order autonomous equations of the form
More 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 first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More 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 informationPermanency and Asymptotic Behavior of The Generalized Lotka-Volterra Food Chain System
CJMS. 5(1)(2016), 1-5 Caspian Journal of Mathematical Sciences (CJMS) University of Mazandaran, Iran http://cjms.journals.umz.ac.ir ISSN: 1735-0611 Permanency and Asymptotic Behavior of The Generalized
More 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 informationPredator-Prey Model with Ratio-dependent Food
University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013 Table of Contents
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 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 informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationTopic 30 Notes Jeremy Orloff
Toic 30 Notes Jeremy Orloff 30 Alications to oulation biology 30.1 Modeling examles 30.1.1 Volterra redator-rey model The Volterra redator-rey system models the oulations of two secies with a redatorrey
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 Predator-prey models are used to show the intricate
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 informationLAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC
LAW OF LARGE NUMBERS FOR THE SIRS EPIDEMIC R. G. DOLGOARSHINNYKH Abstract. We establish law of large numbers for SIRS stochastic epidemic processes: as the population size increases the paths of SIRS epidemic
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 informationKey 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 informationANALYSIS OF A MODEL OF TWO PARALLEL FOOD CHAINS. Sze-Bi Hsu. Christopher A. Klausmeier. Chiu-Ju Lin
DISCRETE AND CONTINUOUS doi:10.3934/dcdsb.2009.12.337 DYNAMICAL SYSTEMS SERIES B Volume 12, Number 2, September 2009 pp. 337 359 ANALYSIS OF A MODEL OF TWO PARALLEL FOOD CHAINS Sze-Bi Hsu Department of
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationMathematical Modeling I
Mathematical Modeling I Dr. Zachariah Sinkala Department of Mathematical Sciences Middle Tennessee State University Murfreesboro Tennessee 37132, USA November 5, 2011 1d systems To understand more complex
More informationQualitative Analysis of Tumor-Immune ODE System
of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System
More informationGerardo 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 informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More informationLecture 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 information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
More informationExam 2 Study Guide: MATH 2080: Summer I 2016
Exam Study Guide: MATH 080: Summer I 016 Dr. Peterson June 7 016 First Order Problems Solve the following IVP s by inspection (i.e. guessing). Sketch a careful graph of each solution. (a) u u; u(0) 0.
More informationFeedback control for a chemostat with two organisms
Feedback control for a chemostat with two organisms Patrick De Leenheer and Hal Smith Arizona State University Department of Mathematics and Statistics Tempe, AZ 85287 email: leenheer@math.la.asu.edu,
More informationModeling the Immune System W9. Ordinary Differential Equations as Macroscopic Modeling Tool
Modeling the Immune System W9 Ordinary Differential Equations as Macroscopic Modeling Tool 1 Lecture Notes for ODE Models We use the lecture notes Theoretical Fysiology 2006 by Rob de Boer, U. Utrecht
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
More informationDynamics of a predator-prey system with prey subject to Allee effects and disease
Dynamics of a predator-prey system with prey subject to Allee effects and disease 3 Yun Kang, Sourav Kumar Sasmal,, Amiya anjan Bhowmick, 3, Joydev Chattopadhyay, 4 4 5 6 Abstract n this article, we propose
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 informationStability Implications of Bendixson s Criterion
Wilfrid Laurier University Scholars Commons @ Laurier Mathematics Faculty Publications Mathematics 1998 Stability Implications of Bendixson s Criterion C. Connell McCluskey Wilfrid Laurier University,
More informationMA 138: Calculus II for the Life Sciences
MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationMath 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 informationDynamical Analysis of a Harvested Predator-prey. Model with Ratio-dependent Response Function. and Prey Refuge
Applied Mathematical Sciences, Vol. 8, 214, no. 11, 527-537 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/12988/ams.214.4275 Dynamical Analysis of a Harvested Predator-prey Model with Ratio-dependent
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: Vol.8, No.7, pp , 2015
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304 Vol.8, No.7, pp 99-, 05 Lotka-Volterra Two-Species Mutualistic Biology Models and Their Ecological Monitoring Sundarapandian
More informationSpotlight on Modeling: The Possum Plague
70 Spotlight on Modeling: The Possum Plague Reference: Sections 2.6, 7.2 and 7.3. The ecological balance in New Zealand has been disturbed by the introduction of the Australian possum, a marsupial the
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 informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationNon-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 informationOn the stabilizing effect of specialist predators on founder-controlled communities
On the stabilizing effect of specialist predators on founder-controlled communities Sebastian J. Schreiber Department of Mathematics Western Washington University Bellingham, WA 98225 May 2, 2003 Appeared
More informationStability analysis of a prey-predator model with a reserved area
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 4, 5(3):93-3 ISSN: 976-86 CODEN (USA): AASRFC Stability analysis of a prey-predator model with a reserved area Neelima
More information