BIFURCATION ANALYSIS OF A MODEL OF THREE-LEVEL FOOD CHAIN IN A MANGROVE ECOSYSTEM
|
|
- Darren Hardy
- 5 years ago
- Views:
Transcription
1 U.P.B. Sci. Bull., Series A, Vol. 78, Iss. 4, 2016 ISSN BIFURCATION ANALYSIS OF A MODEL OF THREE-LEVEL FOOD CHAIN IN A MANGROVE ECOSYSTEM Cristina Bercia 1 and Romeo Bercia 2 In this paper, we consider a three-dimensional non-linear dynamical system of predator-prey type with 12 parameters which models a food-chain in a mangrove ecosystem. Our goal is to study the dynamical properties of the model with constant rate harvesting on the top-predator. It is shown that transcritical, saddle-node and supercritical Hopf ifurcations occur when one parameter is varied. By numerical integration of the system, we plot the phase portraits for the important types of dynamics and we show the presence of a stale limit cycle. We deduce the controlling role of the harvesting rate upon the stale equlirium or periodic state of coexistence of the species. Keywords: Predator-prey system, Staility, Limit cycle, Bifurcation. 1. Introduction We consider a three dimensional ordinary differential autonomous system which models a food chain in a mangrove ecosystem with detritus recycling and with constant rate harvesting of top predators. The three levels of the food chain are the detritus (x(t) denotes its density at time t) which consists of algal species and leaves of the mangrove plants, then detritivores (y(t))- unicellular animals, crustacean, amphiopod and others (see [4]) which depend on rich detritus-ase and the predators of detritivores (z(t))-fish and prawn. Some of the predators have commercial value and undergo harvesting. The system has the following form: x (t) = x 0 ax βxye x + γz y (t) = β 1 xye x d 1 y cyz k + y z (t) = c 1yz k + y d 2z hz with the initial conditions x(0), y(0), z(0) 0. (1) 1 Lecturer, Faculty of Applied Sciences, University Politehnica of Bucharest, Romania, cristina ercia@yahoo.co.uk 2 Lecturer, Faculty of Applied Sciences, University Politehnica of Bucharest, Romania 105
2 106 Cristina Bercia, Romeo Bercia The main assumptions of the model are the following: 1. The detritus has a constant input rate (x 0 ) maintained y the decomposed mangrove leaves. [4] Also the detritus is washout from the system with a certain rate (a) and is supplied with dead organic matter converted into detritus y the action of the micro-organisms. 2. The interaction etween detritus and detritivores is different from the classical models of predator-prey in that the detritivores response function is not a monotone incresing function of prey density, ut rather is only monotone increasing until some critical density and then ecomes monotone decreasing. (as it is descried in [3] etween phytoplankton and zooplankton). This phenomenon is called group defense in which predation (in our case y detritivores) decreases when the density of the prey population is sufficiently large, which is also related to the nutrient uptake inhiition phenomenon in chemical kinetics (see [1]). One of the nonmonotone functional response introduced (see [6]) is p(x) = βxe x where 1 > 0 is the density of the prey at which predation reaches its maximum. 3. The amount of detritivores consumed y the top predator is assumed to follow a Holling type-ii functional response. It exhiits saturation effect when y-population is aundent and k is the half-saturation constant. The top predators (z(t)) are harvested with a constant rate, h. 4. The detritus-detritivores conversion rate is less then the detritus uptake rate, so β 1 < β. 5. For the same reason, the detritivore-predator conversion rate is c 1 < c. 6. The detritus recycle rate due to the death of predators, γ is less than d 2 + h. First we study the oundedness of the solutions of the system (1). Proposition 1.1. The first octant R 3 + is positively invariant under the flow generated y the system (1). Proof. Let e (v 1, v 2, v 3 ) the vector field which defines the differential system (1). We study the vector-field on the oundaries of the first octant. In Oxzplane, v 2 (x, 0, z) = 0, therefore all trajectories which initiate in this plane, remain in it, for any t 0, so the plane y = 0 is an invariant set for the system. In the same situation is the plane z = 0 ecause v 3 (x, y, 0) = 0. But trajectories which start in Oyz-plane, y, z > 0 are directed towards the interior of the first octant, since v 1 (0, y, z) = x 0 + γz > 0. In consequence all solutions with x(0), y(0), z(0) 0, remain in the first octant. Proposition 1.2. With the hypothesis upon the parameters 4,5 and 6, every solution of the system (1) with positive initial values, is ounded. Proof. Let (x(t), y(t), z(t)) e a solution of the system with positive initial conditions. We deduce that (x+y+z) (t) < x 0 ax(t)+(γ d 2 h)z(t) d 1 y(t), t 0. If M = min(a, d 2 + h γ, d 1 ), then (x + y + z) < x 0 M(x + y + z). Next we denote u(t) = (x + y + z)(t). By multiplying the last inequality with
3 Bifurcation Analysis of a Model of Three-level Food Chain in a Mangrove Ecosystem 107 e Mt, we get that u verifies (ue Mt ) x 0 e Mt < 0, t 0, or, equivalently from integration, u(t) (u(0) x 0 M )e Mt + x 0, that is u(t) max(u(0), x 0 ), t 0. M M Therefore, any solution of the system (1) starting in R 3 +, is ounded. Let us nondimensionalize the system (1) with the following scaling: t at; y β cβ y; z z and the system ecomes a a 2 x (t) = x xye x + µz (2) y (t) = xye x y yz K + y z (t) = ryz K + y Hz where = x 0 a ; = β 1 ; D a 1 = d 1 a ; K = kβ ; r = c 1 a a ; H = d 2+h ; µ = γa are positive a cβ constants. For simplicity, we keep the ecological implications of parameters, µ,,, K, r and H the same as x 0, γ, β 1, d 1, k, c 1 and h, respectively. In consequence, in the following study, we assume that µ < H, < 1 and r < Staility and ifurcation analysis Now we determine the location and the existence criteria for the equiliria of the system (2) in the first octant. Proposition 2.1. a) There is an axial equlirium E 0 (, 0, 0), for any values of the parameters; ) If e < < e and > 1, the system has two oundary equiliria E i (x i, yi, 0), i = 1, 2, where x 1 < x 2 are the two solutions of the equation xe x = (3) and yi = ( x i ),i = 1, 2; c) If > e, there exists only one oundary equilirium E 1 (x 1, y1, 0) with x 1 < 1 and y 1 given in ); d) If < e, there are no oundary equiliria in R 3 +. Proof. ) Un equilirium with z = 0 has its first component, solution of the equation xe x =. A simple classical analysis shows that the equation has solutions iff e and in the case with strict inequality there are two solutions x 1 < x 2 which are smaller than when e <, > 1. c) The condition > e is equivalent with x 1 < < x 2 and in consequence, only E 1 is in the first octant. We take the case when ( )M such that xe x x=m >, which, from the second equation in system (2), means that otherwise y could not survive on the prey at any density in the asence of z. In the following study, e and this is the case when the equation (3) posses at least one solution.
4 108 Cristina Bercia, Romeo Bercia Next we give a sufficient existence criterion for the interior equilirium E 3 (x 3, y 3, z 3 ). A simple calculation shows that x 3 is a solution of the equation: g(x) := K H µr r H xe x + x + µrk r H = 0, (4) y 3 = HK r H and z 3 = (x 3 e x 3 ) rk r H. Remark 2.1. We assume that r > H from now on, that is the detritivorepredator conversion rate has to e greater than the harvesting rate, so that there exists un interior equilirium which is the case with iological relevance. Proposition 2.2. i) If < e and > 1 (i.e. the system has two oundary equiliria) and H 2 < H < H 1, then it is at least one interior equilirium, E 3, where H i := r( x i ) i = 1, 2; (5) K + ( x i ), ii) If > e (i.e. there exists only one oundary equilirium) and H < H 1, then it is also at least one interior equilirium, E 3. iii) When the system has two oundary equiliria, H 2 < H 1. Proof. First of all note that z 3 > 0 is equivalent to x 3 e x 3 > or x 3 (x 1, x 2), where x 1 < x 2 are the two solutions of the equation xe x =. Then notice that in equation (4) the coefficient H µr > 0, due to the assumptions (4,5,6) of the model. A sufficient condition for equation (4) to have a solution x 3 (x 1, x 2) is g(x 1)g(x 2) < 0 or in the form x 1 < KH (r H) < x 2 which in case i) ecomes H 2 < H < H 1 and in case ii), H < H 1. We shall see that for H = H 1, E 3 = E 1 and in case i), for H = H 2, E 3 = E 2. Corollary 2.1. If, in addition to the conditions in Proposition 2.2 we require Kµ + e2 H < r =: H K + e 2 0, (6) it follows that the equation g(x) = 0 has only one solution x 3 (x 1, x 2), that is the equilirium E 3 exists and it is unique. Otherwise, g(x) = 0 may have maximum three positive solutions, so there are maximum three interior equiliria. Proof. g (x) = 1 + K H µr (1 r H x)e x attains its minimum m = 1 K(H µr) e 2 (r H) at x = 2. Because g (0) > 0 and lim g (x) > 0, it turns out that if m > 0 x (the condition (6)), then g is increasing. If m < 0, then g (x) = 0 has two solutions and in consequence, the equation g(x) = 0 may have three positive solutions.
5 Bifurcation Analysis of a Model of Three-level Food Chain in a Mangrove Ecosystem 109 Next we discuss the dynamics of the system (2) in the neighorhood of each equilirium. The Jacoian matrix of the linearization of the system is Dv D(x, y, z) = Evaluating 1 ye x (1 x) xe x µ ye x (1 x) xe x Kz (y+k) 2 Dv D(x,y,z) Proposition 2.3. i) If 0 rkz (y+k) 2 y y+k ry y+k H at each equilirium, we get the following results: < e (7), the equilirium E 0 is a hyperolic stale node. If > e, E 0 is a hyperolic saddle. In any case, it is attractive in two directions in the plane y = 0; ii) The oundary equilirium E 1 is locally asymptotically stale for H > H 1 and it is a hyperolic saddle for H < H 1 ; iii) E 2 is a hyperolic saddle for the values of the parameters which ensure its existence, namely < e, > 1 and for any H. Proof. i) The eigenvalues of the Jacoian matrix evaluated at E 0 are λ 1 = 1, λ 2 = H, λ 3 = e and the eigenvectors corresponding to λ 1,2 < 0 are u 1 = (1, 0, 0), u 2 = (µ, 0, 1 H) which imply i). ii),iii) For E i, i = 1, 2, the eigenvalues are λ i 1 = r( x i ) H and K +( x i ) λi 2,3 are the roots of the equation λ 2 + (1 + p i )λ + p i = 0 with p i = ( x i )(1 x i ). x i Hence, p 1 > 0 and in consequence, Re(λ 1 2,3) < 0, ut p 2 < 0 so, λ 2 2λ 2 3 < 0 and E 2 ecomes a saddle. The eigenvalues of the Jacoian matrix at the interior equilirium E 3 (x 3, y 3, z 3 ) verify the characteristic equation where A 1 = H r λ 3 + A 1 λ 2 + A 2 λ + A 3 = 0 (8) H r(r H) e x 3 [Kr x 3 ((r H) + Kr)]; A 2 = H r (x 3e x 3 )[r H 1 KH r H e x 3 (1 x 3 )] + KH r H e 2x 3 x 3 (1 x 3 ); (9) A 3 = (x 3 e x 3 )(r H)[ KH r H e x 3 (1 x 3 )( H r µ) + H r ]. The necessary and sufficient condition for E 3 to e asymptotically stale is given y the Routh-Hurwitz criterion, i.e. A 1 A 2 > A 3 and A 1, A 3 > 0;
6 110 Cristina Bercia, Romeo Bercia The third component of E 3 has to e strictly positive, or equivalently x 3 e x 3 > 0. Then we notice that A 1 > 0 if x 3 < rk (r H) + Kr. (10) This implies x 3 < 1 and since H > µr, we have A 3 > 0. With simple algera, we deduce the following Proposition 2.4. If E 3 exists, the condition (10) holds and ( H r (x 3e x 3 ) 1 KH ) r H e x 3 (1 x 3 ) ( x3 e x 3 ( 1 + KH ) r r H e x 3 (1 x 3 ) K ) r H e 2x 3 x 3 (1 x 3 ) (11) ( ) +(x 3 e x 3 ) µke x 3 H(r H) (1 x 3 ) (x r 2 3 e x 3 ) > 0 then the equilirium E 3 is locally asymptotically stale. The previous propositions are sintethized in the following Theorem 2.1. Let > 1, H 1 < H 0 and e < < e. Then: i) If H 2 > µ, for H (µ, H 2 ), the system (2) has at least three equiliria: E 0 as attractive node, E 1 and E 2 hyperolic saddles; ii) For max(µ, H 2 ) < H < H 1, the system has the equilirium points E 0, E 1, E 2 as in the previous case and E 3 which is locally asymptotically stale (l.a.s.) if conditions (10) and (11) hold; iii) For H = H 2, E 2 and E 3 concide; iv) For H = H 1, E 1 and E 3 concide; v) For H 1 < H < r, E 3 may e unphysical, E 1 ecomes l.a.s., E 0 remains attractive node and E 2 hyperolic saddle. Theorem 2.2. Let H 1 < H 0 and > e. Then: i) For µ < H < H 1, the system has the equilirium points E 0, E 1 hyperolic saddles and E 3 which is l.a.s. if conditions (10) and (11) hold; ii) For H = H 1, E 1 and E 3 concide; iii) For H 1 < H < r, E 3 may e unphysical, E 1 ecomes l.a.s., E 0 remains hyperolic saddle. 3. Bifurcation analysis First we take as the control parameter for the system (2). Theorem 3.1. If H er( 1), X e( 1)+K 0 > 1, for = ed 1, the oundary equiliria E 1 and E 2 appear through a saddle-node ifurcation. For < e, these two equiliria do not exist.
7 Bifurcation Analysis of a Model of Three-level Food Chain in a Mangrove Ecosystem 111 Proof. We fix = e, then x 1 = x 2 = 1 < X 0 and y1 = y2, so E 1 = E 2. The eigenvalues corresponding to E 1 are λ 1 = er( 1) H, λ e( 1)+K 2 = 0, λ 3 = 1, so we need λ 1 0 in order to have the codimension 1 ifurcation of saddle-node type. Let e DF Dv (E 1, e ) the Jacoian matrix of the system written in the form v = F (v, ), with the ifurcation parameter and v = (x, y, z), evaluated in E 1 for = e. We use a theorem (Sotomayor [2]) which gives necessary and sufficient conditions for a saddle-node ifurcation. These are: (SN1) DF Dv (E 1, e ) has a single zero eigenvalue which is fulfilled, due to the hypothesis of the theorem. Let e u = (1, e, 0) T its right eigenvector and w = (0, r, 1) T its left eigenvector. (SN2) w F (E 1, e ) = r( 1) 0. (SN3) w D vv F (E 1, e )(u, u) = 2 e ( 1)r 0. The nonzero conditions (SN2) and (SN3) imply that SN = { = e, H er( 1), X e( 1)+K 0 > 1 } is a saddle-node ifurcation surface in the parameters space. When the parameters pass from one side of the surface to the other side, the numer of the oundary equiliria changes from zero when < e to two hyperolic equiliria E 1,2 when > e, in the neighorhood of = e. These equiliria are connected y an orit that is asymptotic to E 1 for t and to E 2 for t. Now we take H as a control parameter for the system. Theorem 3.2. If > 1, > e and µ > K ( x 1) e 2, (12) ( x 1) + K then a) the equiliria E 1 and E 3 coincide for H = r( x 1 ) =: H K +( x 1 ) 1 at a point of transcritical ifurcation; ) E 2 and E 3 coincide for H = r( x 2 ) =: H K +( x 2 ) 2 also at a point of transcritical ifurcation, if, in addition < e. Proof. a) i) We fix H = H 1. The condition (12) is equivalent to H 1 < r Kµ+e2 K+e 2 which (see Proposition 2.2 and corrolary) implies that the solution x 3 of the equation (4), g(x) = 0, exists and it is unique. Then H = H 1 x 1 = KH. But (r H) x 1 < 1 and verifies xe x =, in consequence x 1e x 1 + x 1 = µrk r H 1 and x 1 = x 3. Also y 3 = KH 1 r H 1 = y1, z 3 = 0 K H 1 µr r H 1 which imply E 1 = E 3. From the proof of Proposition 2.3, we have: ii) For H = H 1, the eigenvalues for E 1 are λ 1 = 0, Re(λ 2,3 ) < 0, so only one eigenvalue is zero. iii) On the other hand, the equilirium E 1, from unstale (when H < H 1 ) ecomes stale for H > H 1. iv) We evaluate the matrix ( DF Dv F H ) at the ifurcation point (E1, H 1 ) and we
8 112 Cristina Bercia, Romeo Bercia find that F = (0, 0, H 0)T. It implies rank ( DF ) F Dv H = 2. From i)-iv) (see [5]), we deduce that in R 3 + {H > 0} the ranches of equiliria E 1 and E 3 intersect through a transcritical ifurcation at H = H 1 and change their stale manifold when we pass the critical parameter value. ) There are similar arguments. The condition < e ensures the existence of the equilirium E 2. Then H = H 2 x 3 = x 2 > 1. Also only one eigenvalue (see proof of Proposition 2.3) corresponding to E 2 is zero for H = H 1. These values of the parameters H = H 1, H = H 2, = e delimitate strata on parameters space induced y topological equivalence of the phase portraits. We are now investigating dynamic ifurcations. The equilirium E 3 is the only one which may experience Hopf ifurcation ecause only the interior equilirium can have a pair of purely imaginary eigenvalues. This necessary condition for E 3 to undergo a Hopf ifurcation is equivalent to A 1 A 2 = A 3, A 2 > 0 (13) (with A i given y (9)) together with the sufficient conditions for the existence of the interior equilirium (see Proposition 2.2). Also we need that the third eigenvalue of E 3 to e non-zero, i.e. A 3 0, which is satisfied if x 3 < 1 with g(x 3 ) = 0. In this case, we have A 3 > 0 and so the characteristic equation (8) admits the roots λ 1,2 = ±iω, ω > 0 and λ 3 = A 1 < Numerical results We take H as a control parameter. We fix the parameters in order to e e in the case e < < D 1 (see Theorem 2.1): K = 1; = 1.2; r = 0.3; = 0.66; = 0.2; µ = 0.1; = 1. This case corresponds with the existence of two oundary equiliria E i (x i, yi, 0), i = 1, 2. With a program in MAPLE, LamertW (0, D1 ) = we find the solutions of the equation xe x =, x 1 = ,x 2 = LamertW ( 1, ) = , then y 1 = , y2 = The static ifurcation parameters are H 2 = < H 1 = We solve numerically the system A 2 > 0, x 3 < 1, x 3e x 3 > 0 A 1 A 2 = A 3 g(x 3 ) = 0 (14) which are the necessary conditions for Hopf ifurcation of E 3 and we find the ifurcation parameter value H = H cr = and x 3 = H cr (H 2, H 1 ), H 1 < H 0 and we are in the hypothesis of Proposition 2.2 and its corollary. There is only one interior equilirium for any H (H 2, H 1 ).
9 Bifurcation Analysis of a Model of Three-level Food Chain in a Mangrove Ecosystem 113 Then, we investigate the appearance of a limit cycle when H is in the neighorhood of the critical parameter value. We find that for H < H cr an asymptotically stale limit cycle appears. For example, when H = 0.11 (H 2, H cr ), the solution of equation g(x 3 ) = 0, x 3 e x 3 > is x 3 = Then y 3 = ; z 3 = With initial conditions close to E 3, we integrate numerically the system. With a program in MATLAB we otain the phase portrait. The trajectories come together close to the saddle connection etween E 1 and E 2, towards E 1 and then tend to a stale limit cycle, while E 3 is repelling. (see Figure 1) Figure 1. Left: A stale limit cycle when K = 1; = 1.2; r = 0.3; = 0.66; = 0.2; µ = 0.1; = 1; H = 0.11 < H cr. Right: Time oscillations of populations corresponding to one of the trajectories, with initial values (0.6228; ; ). Figure 2. Left: Trajectories tending to the attractive focus E 3. The modified parameter is H = 0.13 > H cr. Right: Time evolution of the three populations corresponding to one of the trajectories, with initial values (0.7666; ; ).
10 114 Cristina Bercia, Romeo Bercia For H = 0.13 (H cr, H 1 ), trajectories initiated in a small neighorhood of E 3, come together very close to the saddle connection E 2 E 1 and tend for t to E 3 as a focus. (see Figure 2) In consequence, (E 3, H cr ) is a point of supercritical Hopf ifurcation ecause, when H varies and passes the critical value, from a stale equilirium E 3 for H > H cr, it appears a stale limit cycle for H < H cr, while E 3 loses its staility. 5. Conclusions The harvesting rate (H) of top predator can control the stale equilirium point or periodic state of coexistence of the three species. First of all, the model illustrates the phenomenon of iological overharvesting when H > H 1 and in this case the top-predator goes to extinction. If the detritus-detritivores conversion rate is not high, e<< e, the scenario of extinction of oth detritivores and top-predator is possile, for any value of the harvesting rate, depending on initial population levels. When we decrease H, for H (H cr, H 1 ), all the populations coexist in a form of a stale equilirium, under certain initial conditions. When H passes through a critical value, H cr, the system undergoes a Hopf ifurcation, namely for H (max(µ, H 2 ), H cr ), different populations of the system will start to oscillate with a finite period around the equilirium point of coexistence. Since H = d 2+h, it turns out that also for increasing large values of a washout rate (a) of detritus, the system changes from a stale state to an unstale state and the populations will survive through periodic fluctuations. Mathematically, we would like to point out here that our analysis of the model is a first look at the local ifurcations, ut it is not complete. For example a study of codimension 2 ifurcations will reveal richer dynamics. R E F E R E N C E S [1] H. I. Freedman, G. S. K. Wolkowicz, Predator-prey systems with group defence: The paradox of enrichment revisited, Bull. Math. Biol., 48(1986), [2] J. Guckenheimer, P. Holmes, Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields, Springer Verlag New York, 1983 [3] B. Mukhopadhyay, R. Bhattacharyya, Modelling Phytoplankton Allelopathy in a Nutrient-Plankton Model with Spatial Heterogeneity, Ecological Modelling, 198(2006), 198, p [4] B. Mukhopadhyay, R. Bhattacharyya, Diffusion Induced Shift of Bifurcation Point in a Mangrove Ecosystem Food-chain Model with Harvesting, Natural Resource Modeling, 22(2009), No 3, p [5] A. H. Nayfeh, B. Balachandran, Applied Nonlinear Dynamics, Analytical, Computational and Experimental Methods, Wiley Series, [6] S. Ruan, H. I. Freedman, Persistence in the three species food chain models with group defence, Math. Biosci., 107(1991), p
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 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 informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationEQUILIBRIA AND STABILITY ANALYSIS OF A BRANCHED METABOLIC NETWORK WITH FEEDBACK INHIBITION. F. Grognard Y. Chitour G. Bastin
EQUILIBRIA AND STABILITY ANALYSIS OF A BRANCHED METABOLIC NETWORK WITH FEEDBACK INHIBITION F. Grognard Y. Chitour G. Bastin Projet COMORE. INRIA Sophia-Antipolis. BP 93 06902 Sophia-Antipolis Cedex, France
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 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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationMath 216 Second Midterm 28 March, 2013
Math 26 Second Midterm 28 March, 23 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
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 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 informationHopf bifurcation and Turing instability in the reaction diffusion Holling Tanner predator prey model
IMA Journal of Applied Mathematics 203 78, 287 306 doi:0.093/imamat/hxr050 Advance Access pulication on Novemer 7, 20 Hopf ifurcation and Turing instaility in the reaction diffusion Holling Tanner predator
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 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 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 informationAndronov Hopf and Bautin bifurcation in a tritrophic food chain model with Holling functional response types IV and II
Electronic Journal of Qualitative Theory of Differential Equations 018 No 78 1 7; https://doiorg/10143/ejqtde018178 wwwmathu-szegedhu/ejqtde/ Andronov Hopf and Bautin bifurcation in a tritrophic food chain
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 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 informationSTABILITY. 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 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 informationWorkshop 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 informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationarxiv: v1 [nlin.ao] 16 Jul 2017
arxiv:1707.05676v1 [nlin.ao] 16 Jul 017 Andronov Hopf ifurcation with and without parameter in a cuic memristor oscillator with a line of equiliria Ivan A. Korneev 1 1, a) and Vladimir V. Semenov Department
More informationDifference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay
Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,
More informationNonlinear 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 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 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 informationStability Domain of a Linear Differential Equation with Two Delays
ELSEVIER An International Journal Availale online at www.sciencedirect.com computers &.c,..c. ~--~c,..c.. mathematics with applications Computers and Mathematics with Applications 51 (2006) 153-159 www.elsevier.com/locate/camwa
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 informationStability of Ecosystem Induced by Mutual Interference between Predators
Available online at www.sciencedirect.com Procedia Environmental Sciences (00) 4 48 International Society for Environmental Information Sciences 00 Annual Conference (ISEIS) Stability of Ecosystem Induced
More informationNUMERICAL 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 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 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 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 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 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 informationDynamics of an HIV pathogenesis model with CTL immune response and two saturated rates
ISSN 746-7233 England UK World Journal of Modelling and Simulation Vol. (24) No. 3 pp. 25-223 Dynamics of an HIV pathogenesis model with CTL immune response and two saturated rates Youssef Tait Khalid
More informationHopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models
Vol.37 (SUComS 6), pp.57-66 http://dx.doi.org/.457/astl.6.37. Hopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models Shuang Guo, Xiuli Li, Jian Yu, Xiu Ren Department of Mathematics, Teacher
More informationHopf bifurcations in an extended Lorenz system
Zhou et al. Advances in Difference Equations 2017 2017:28 DOI 10.1186/s13662-017-1083-8 R E S E A R C H Open Access Hopf ifurcations in an extended Lorenz system Zhiming Zhou 1,2*, Gheorghe Tigan 3 and
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 information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
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 informationMATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation.
MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation. Bifurcation theory The object of bifurcation theory is to study changes that maps undergo
More informationStability and bifurcation of a simple neural network with multiple time delays.
Fields Institute Communications Volume, 999 Stability and bifurcation of a simple neural network with multiple time delays. Sue Ann Campbell Department of Applied Mathematics University of Waterloo Waterloo
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 informationLinearized geometric dynamics of Tobin-Benhabib-Miyao economic flow
Linearized geometric dynamics of Tobin-Benhabib-Miyao economic flow Constantin Udrişte and Armando Ciancio Abstract The aim of this paper is to study the influence of the Euclidean-Lagrangian structure
More informationMultiple focus and Hopf bifurcations in a predator prey system with nonmonotonic functional response
Multiple focus and Hopf bifurcations in a predator prey system with nonmonotonic functional response Dongmei Xiao Department of Mathematics Shanghai Jiao Tong University, Shanghai 200030, China Email:
More informationA Cournot duopoly model with complementary goods: multistability and complex structures in the basins of attraction.
Fernando Bignami, Anna Agliari A Cournot duopoly model with complementary goods: multistaility and complex structures in the asins of attraction. In this paper we study a Cournot duopoly, proposed y Matsumoto
More informationHarvesting Model for Fishery Resource with Reserve Area and Modified Effort Function
Malaya J. Mat. 4(2)(2016) 255 262 Harvesting Model for Fishery Resource with Reserve Area and Modified Effort Function Bhanu Gupta and Amit Sharma P.G. Department of Mathematics, JC DAV College, Dasuya
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 informationDroop models of nutrient-plankton interaction with intratrophic predation
Droop models of nutrient-plankton interaction with intratrophic predation S. R.-J. Jang 1, J. Baglama 2 1. Department of Mathematics, University of Louisiana at Lafayette, Lafayette, LA 754-11 2. Department
More informationJOURNAL OF MATH SCIENCES -JMS- Url: Jl. Pemuda, No. 339 Kolaka Southeast Sulawesi, Indonesia
JOURNAL OF MATH SCIENCES -JMS- Url: http://usnsj.com/index.php/jms Jl. Pemuda, No. 339 Kolaka Southeast Sulawesi, Indonesia STABILITY ANALYSIS OF ONE PREY TWO PREDATOR MODEL WITH HOLLING TYPE III FUNCTIONAL
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 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 informationDynamical Systems Solutions to Exercises
Dynamical Systems Part 5-6 Dr G Bowtell Dynamical Systems Solutions to Exercises. Figure : Phase diagrams for i, ii and iii respectively. Only fixed point is at the origin since the equations are linear
More informationDYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS. Wei Feng
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 7 pp. 36 37 DYNAMICS IN 3-SPECIES PREDATOR-PREY MODELS WITH TIME DELAYS Wei Feng Mathematics and Statistics Department
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationDepartment of Mathematics IIT Guwahati
Stability of Linear Systems in R 2 Department of Mathematics IIT Guwahati A system of first order differential equations is called autonomous if the system can be written in the form dx 1 dt = g 1(x 1,
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 informationSharp estimates of bounded solutions to some semilinear second order dissipative equations
Sharp estimates of ounded solutions to some semilinear second order dissipative equations Cyrine Fitouri & Alain Haraux Astract. Let H, V e two real Hilert spaces such that V H with continuous and dense
More informationEN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015
EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions
More informationStudy of a detritus-based ecosystem model for generalized functional responses
Study of a detritus-based ecosystem model for generalized functional responses M.R. Mandala), N.H. Gazib) a) Department of Mathematics, Siliguri College, Siliguri, West Bengal 7341, India email: mdreduanur.mandal@gmail.com
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 informationINTRATROPHIC PREDATION IN A SIMPLE FOOD CHAIN WITH FLUCTUATING NUTRIENT. S. R.-J. Jang. J. Baglama. P. Seshaiyer. (Communicated by Shigui Ruan)
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 5, Number 2, May 25 pp. 335 352 INTRATROPHIC PREDATION IN A SIMPLE FOOD CHAIN WITH FLUCTUATING NUTRIENT S. R.-J.
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 informationBifurcation Analysis of Prey-Predator Model with Harvested Predator
International Journal of Engineering Research and Development e-issn: 78-67X, p-issn: 78-8X, www.ijerd.com Volume, Issue 6 (June 4), PP.4-5 Bifurcation Analysis of Prey-Predator Model with Harvested Predator
More informationWeak Allee effect in a predator-prey system involving distributed delays
Volume 30, N. 3, pp. 675 699, 2011 Copyright 2011 SBMAC ISSN 0101-8205 www.scielo.br/cam Weak Allee effect in a predator-prey system involving distributed delays PAULO C.C. TABARES 1, JOCIREI D. FERREIRA
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 informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
More informationCHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM
CHARACTERIZATION OF SADDLE-NODE EQUILIBRIUM POINTS ON THE STABILITY BOUNDARY OF NONLINEAR AUTONOMOUS DYNAMICAL SYSTEM Fabíolo Moraes Amaral, Josaphat Ricardo Ribeiro Gouveia Júnior, Luís Fernando Costa
More informationAdditive resonances of a controlled van der Pol-Duffing oscillator
Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University
More informationCusp-Bodanov-Takens Bifurcation in a Predator-Prey Type of Dynamical System with Time-Periodic Perturbation
Cusp-Bodanov-Takens Bifurcation in a Predator-Prey Type of Dynamical System with Time- Kus Prihantoso Johan Matheus October 24 th, 2011 ITB-Bandung Indonesia Predator-Prey System Cusp-Bogdanov-Takens Bifurcation
More informationGrowth models for cells in the chemostat
Growth models for cells in the chemostat V. Lemesle, J-L. Gouzé COMORE Project, INRIA Sophia Antipolis BP93, 06902 Sophia Antipolis, FRANCE Valerie.Lemesle, Jean-Luc.Gouze@sophia.inria.fr Abstract A chemostat
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
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 informationLocal Stability Analysis of a Mathematical Model of the Interaction of Two Populations of Differential Equations (Host-Parasitoid)
Biology Medicine & Natural Product Chemistry ISSN: 089-6514 Volume 5 Number 1 016 Pages: 9-14 DOI: 10.1441/biomedich.016.51.9-14 Local Stability Analysis of a Mathematical Model of the Interaction of Two
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 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 informationUNIVERSITY of LIMERICK
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics & Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS08 SEMESTER: Spring 0 MODULE TITLE: Dynamical
More informationQUARTERLY OF APPLIED MATHEMATICS
QUARTERLY OF APPLIED MATHEMATICS Volume LIV December 1996 Number 4 DECEMBER 1996, PAGES 601-607 ON EXISTENCE OF PERIODIC ORBITS FOR THE FITZHUGH NERVE SYSTEM By S. A. TRESKOV and E. P. VOLOKITIN Institute
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 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 informationApplication demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems
Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics
More informationComplex dynamics of a three species food-chain model with Holling type IV functional response
Nonlinear Analysis: Modelling and Control, 211, Vol. 16, No. 3, 353 374 353 Complex dynamics of a three species food-chain model with Holling type IV functional response Ranjit Kumar Upadhyay, Sharada
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationLong food chains are in general chaotic
OIKOS 109: 135/144, 2005 Long food chains are in general chaotic Thilo Gross, Wolfgang Ebenhöh and Ulrike Feudel Gross, T., Ebenhöh, W. and Feudel, U. 2005. Long food chains are in general chaotic. / Oikos
More informationImpulsive Stabilization and Application to a Population Growth Model*
Nonlinear Dynamics and Systems Theory, 2(2) (2002) 173 184 Impulsive Stabilization and Application to a Population Growth Model* Xinzhi Liu 1 and Xuemin Shen 2 1 Department of Applied Mathematics, University
More informationA New Modified Hyperchaotic Finance System and its Control
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.8(2009) No.1,pp.59-66 A New Modified Hyperchaotic Finance System and its Control Juan Ding, Weiguo Yang, Hongxing
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 informationSTABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 20, Number 4, Winter 2012 STABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS XIHUI LIN AND HAO WANG ABSTRACT. We use an algebraic
More informationConverse Lyapunov theorem and Input-to-State Stability
Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts
More informationExact Free Vibration of Webs Moving Axially at High Speed
Eact Free Viration of Wes Moving Aially at High Speed S. HATAMI *, M. AZHARI, MM. SAADATPOUR, P. MEMARZADEH *Department of Engineering, Yasouj University, Yasouj Department of Civil Engineering, Isfahan
More informationNONSTANDARD NUMERICAL METHODS FOR A CLASS OF PREDATOR-PREY MODELS WITH PREDATOR INTERFERENCE
Sixth Mississippi State Conference on Differential Equations and Computational Simulations, Electronic Journal of Differential Equations, Conference 15 (2007), pp. 67 75. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
More informationMULTIPLE FOCUS AND HOPF BIFURCATIONS IN A PREDATOR-PREY SYSTEM WITH NONMONOTONIC FUNCTIONAL RESPONSE
SIAM J. APPL. MATH. Vol. 66, No. 3, pp. 802 89 c 2006 Society for Industrial and Applied Mathematics MULTIPLE FOCUS AND HOPF BIFURCATIONS IN A PREDATOR-PREY SYSTEM WITH NONMONOTONIC FUNCTIONAL RESPONSE
More informationOBSERVABILITY AND OBSERVERS IN A FOOD WEB
OBSERVABILITY AND OBSERVERS IN A FOOD WEB I. LÓPEZ, M. GÁMEZ, AND S. MOLNÁR Abstract. The problem of the possibility to recover the time-dependent state of a whole population system out of the observation
More informationThe Invariant Curve in a Planar System of Difference Equations
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 13, Number 1, pp. 59 71 2018 http://campus.mst.edu/adsa The Invariant Curve in a Planar System of Difference Equations Senada Kalabušić,
More informationOn predator-prey systems and small-gain theorems
On predator-prey systems and small-gain theorems Patrick De Leenheer, David Angeli and Eduardo D. Sontag December 4, Abstract This paper deals with an almost global attractivity result for Lotka-Volterra
More information