ARTICLE IN PRESS. Nonlinear Analysis: Real World Applications. Contents lists available at ScienceDirect
|
|
- Adrian Davidson
- 5 years ago
- Views:
Transcription
1 Nonlinear Analsis: Real World Applications ( Contents lists available at ScienceDirect Nonlinear Analsis: Real World Applications journal homepage: Qualitative analsis of a predator pre model with constant-rate pre harvesting incorporating a constant pre refuge Lili Ji, Chengqiang Wu College of Mathematics and Computer Science, Fuzhou Universit, Fuzhou, Fujian, 3000, PR China a r t i c l e i n f o a b s t r a c t Article histor: Received Ma 009 Accepted Jul 009 Kewords: Predator pre model Pre refuge Limit ccle Harvesting Global stabilit We consider a predator pre model with Holling tpe II functional response incorporating a constant pre refuge and a constant-rate pre harvesting. Depending on the constant pre refuge m, which provides a condition for protecting m of pre from predation, and the constant-rate pre harvesting, some sufficient conditions for the instabilit and global stabilit of the equilibria, and the eistence and uniqueness of limit ccles of the model are obtained. We also show the influences of pre refuge and harvesting efforts on equilibrium densit values. Numerical simulations are carried out to illustrate the feasibilit of the obtained results and the dependence of the namic behavior on the harvesting efforts or pre refuge. 009 Elsevier Ltd. All rights reserved.. Introduction The stu of the consequences of hiding behavior of pre on the namics of predator pre interactions can be recognized as a major issue in applied mathematics and theoretical ecolog [ 3]. Some of the empirical and theoretical work have investigated the effect of pre refuges and drawn a conclusion that the refuges used b pre have a stabilizing effect on the considered interactions and pre etinction can be prevented b the addition of refuges [ ]. Gonzalez-Olivares and Ramos-Jiliberto [9] studied the namic consequences of the following predator pre sstems with constant number of pre using refuges, which protects m of pre from predation: = α ( k β( m + a( m, = + cβ( m + a( m. The showed the effects of the two models of refuge on the namic propert, such as the local stabilit of equilibria and the eistence of limit ccle. Chen [0] gave the complete qualitative analsis for model (.. Motivated b the paper of Gonzalez-Olivares and Ramos-Jiliberto [9], we consider the following predator pre model with Holling tpe II functional response incorporating a positive constant pre refuge and a nonzero constant-rate pre harvesting = α ( k = + cβ( m + a( m, β( m + a( m h, (. (. Corresponding author. addresses: lilms@3.com (L. Ji, cqw00@fzu.edu.cn (C. Wu. -/$ see front matter 009 Elsevier Ltd. All rights reserved. doi:0.0/j.nonrwa
2 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( where and denote the pre and predator population sizes, respectivel, at an time t; α > 0 represents the intrinsic growth rate of the pre; k is the carring capacit of the pre in the absence of predator and harvesting; m is a constant number of pre using refuges, which protects m of pre from predation; the term β/(+a denotes the functional response of the predator, which is known as Holling tpe II response function []; c > 0 is the conversion factor denoting the number of newl born predators for each captured pre; d > 0 is the death rate of the predator; h > 0 is the rate of pre harvesting. From the point of view of human needs, the eploitation of biological resources and the harvest of population are commonl practiced in fisher, forestr and wildlife management. Concerning the conservation for the long-term benefits of humanit, there is a wide range of interest in the use of bioeconomic modeling to gain insight into the scientific management of renewable resources like fisheries and forestries. The problem of predator pre interactions under constant rate of harvesting of either species or both species simultaneousl have been studied b some authors. For eample, Brauer and Soudak [ ] and Xiao [,] studied a class of predator pre models under constant rate of harvesting and under constant quota of harvesting of both species simultaneousl. The showed how to classif the possibilities of the quantitative behavior of the solutions to locate the set of initial values in which the trajectories of the solutions approach to either an asmptotic stable equilibrium or an asmptoticall stable limit ccle. Recentl, Dai and Tang [9] studied a predator pre model in which two ecological interacting species are harvested independentl with constant rates. The showed that the harvested predator pre sstem ma ehibit ver complicated namics such as a spontaneous appearance of a homoclinic orbit and multiple limit ccles. The rest of the paper is organized as follows. Basic properties such as the eistence, stabilit and instabilit of the equilibria of the model are given in Section. In Section 3, sufficient conditions for the global stabilit of the unique positive equilibrium are obtained. In Section, we derive the eistence and uniqueness of limit ccle. In Section, we analze the influence of pre refuge and harvesting efforts and give numerical stimulations. This paper ends b a brief conclusion.. Basic properties of the model Let Ω 0 = {(, > m, > 0}. For practical biological meaning, we simpl stu sstem (. in Ω 0 or in Ω 0. From the first equation of sstem (., it is eas to derive lim sup t + (t k. Lemma.. The solution ((t, (t of sstem (. with the initial values (0 > m, (0 > 0 is positive and bounded for all t 0. Proof. Obviousl, the solution ((t, (t of sstem (. with the initial values (0 > m, (0 > 0 is positive for all t 0. Define the function ω(, = c +. Given an ε > 0, (t k + ε for t sufficientl large. From sstem (., it follows that ( ( dω = cα k ch cα k then ω cα(k+ε, for t sufficientl large, which completes the proof. min{d,α} Lemma. ([0]. Assume that stable (unstable limit ccle. = X(,, min{d, α}ω + cα(k + ε, = Y(,. Along a limit ccle L, if div(x, Y < 0(> 0 holds, then L is a L Lemma.3 ([0]. Assume that = X(,, ε, = Y(,, ε, X(,, ε, Y(,, ε structures a rotation vector field in the wider sense. And furthermore suppose that the set of equations admit a semi-stable limit ccle L ε0 for ε = ε 0. Then, when parameter ε changes according to appropriate direction, L ε0 at least decomposes to a stable and an unstable limit ccle, and the la in inboard and outboard of L ε0, respectivel; while ε changes in opposite direction, L ε0 disappears. For simplicit, we take the following scaling: = m, = β, = ( + a( mdτ and then sstem (. takes the following form (still denote,, τ as,, t = b( + a[ + (k m + km m u] P(,, = (cβ a Q (,, where b = α k, u = h b. Clearl, Ω 0 transforms into Ω = (, > 0, > 0 and sstem (. is bounded. Considering the function Φ( = + (k m + km m u, the discriminant of Φ( = 0 takes the following form = (k m + (km m u = k u. If = k u > 0 and Φ(0 = km m u > 0 hold, Φ( has onl one positive root 0 = k m+ k u. (.
3 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( 3 So, if 0 < u < km m holds, sstem (. has at most two equilibria P 0 ( 0, 0, P (,, where d = cβ ad, = b( + a Φ(. P (, is a unique positive equilibrium if and onl if cβ ad > 0, 0 < m < k+ k u d. cβ ad In what follows, we onl consider the case that 0 < u < km m. Lemma.. If m > k+ k u d cβ ad holds, then P 0 is a stable node point; if 0 < m < k+ k u d cβ ad holds, then P 0 is a saddle point. Furthermore assume that m = [(k + a ( + a + k u]/, then (I P is locall asmptoticall stable for all 0 < m < k+ k u d, if ak ad3 u(cβ ad 3 holds; cβ ad (cβ add (II P is locall asmptoticall stable for m < m < k+ k u d and P cβ ad is locall unstable for 0 < m < m, if ak > ad3 u(cβ ad 3 holds. (cβ add Proof. The Jacobian matri of sstem (. for the equilibrium point P 0 is given b ( ρ J(P 0 = 0, 0 ρ where ρ = b( + a 0 ( 0 + k m and ρ = (cβ ad 0 d. If m > k+ k u d holds, the eigenvalues cβ ad of matri J(P 0 are ρ, ρ which are negative, hence, P 0 is locall asmptoticall stable and furthermore it is a stable node point. If 0 < m < k+ k u d holds, one of the eigenvalues of matri J(P cβ ad 0 is ρ which is positive, hence, P 0 is locall unstable and furthermore P 0 it is a saddle point. The Jacobian matri of sstem (. for P is given b ( P J(P =, ω 0 where P = abφ( + b( + a φ ( and ω = cβ ad. Note that tr(j(p = b [m (k + a m ( ak a 3 + u]. Considering the function R( = (k + a ( ak a 3 + u, the discriminant of R( = 0 takes the following form = (k + a + [( ak a 3 + u] = + k u > 0, therefore, R( has two roots m = [(k + a ( + a + k u]/, m = [(k + a + ( + a + k u]/. It is eas to see that m > k. (I If ak ad3 u(cβ ad 3 (cβ add holds, then m 0 and for all m < < m, R( < 0, hence, for all 0 < m < k+ k u d cβ ad < 0, P is locall asmptoticall stable. (II If ak > ad3 u(cβ ad 3 holds, then m (cβ add > 0 and R( < 0 for all m < < m, R( > 0 for all 0 < < m. Notice that > ( + a, then m < (k+a (+a < k+ k u d, hence, if m cβ ad < m < k+ k u d holds, cβ ad then tr(j(p < 0 and P is locall asmptoticall stable; if 0 < m < m holds, then tr(j(p > 0 and P is locall unstable. This completes the proof. Theorem.. If m > k+ k u d cβ ad hold, then P 0 is globall asmptoticall stable. Proof. Notice that m > k+ k u d.p cβ ad 0( 0, 0 is the unique equilibrium of sstem (.. If sstem (. eists a closed orbit in Ω, then there must eist an equilibrium in the interior of the closed orbit, which is impossible. Hence, sstem (. does not eist limit ccle for all m > k+ k u d. From the boundedness of sstem (., the stable node point P cβ ad 0 is globall asmptoticall stable.
4 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( 3. Global stabilit for the positive equilibrium Theorem 3.. Assume that 0 < u < km m, and one of the following conditions holds: (H ak min{ ad3 u(cβ ad 3, am}, 0 < m < k+ k u d ; (cβ add cβ ad (H ad3 u(cβ ad 3 < ak am, m (cβ add < m < k+ k u d. cβ ad Then the positive equilibrium P (, of sstem (. is globall asmptoticall stable. Proof. Define a Dulac function B(, =, from sstem (., we have D = (BP + (BQ = b [a 3 + ( ak + am + km m u] b ϕ(. Then ϕ ( = [3a + ( ak + am] = 0 has two roots = 0, ak + am =. 3a If ak am, then 0. For all 0, ϕ ( 0 and ϕ ( = 0 holds if and onl if = 0. Hence, ϕ( is increasing in [0, +. From ϕ(0 = km m u > 0, we derive ϕ( > 0, for all > 0, therefore, for all 0, D < 0, sstem (. does not eist limit ccle. If one of (H, (H holds, P is locall asmptoticall stable and sstem (. does not eist limit ccle in Ω. From the boundedness of sstem (., P is globall asmptoticall stable.. Eistence and uniqueness of limit ccle Theorem.. Suppose that 0 < u < km m, ak > ad3 u(cβ ad 3, and 0 < m < m. Sstem (. admits at least one limit ccle in Ω. (cβ add Proof. For sstem (., constructing a Bendison ring OA BCD including P. Define OA as a length of line L = = 0, AB as a length of line L = k m+ k u = 0. Define = (a 0, = [ d + (cβ ad], where a 0 = ma 0 b( + a[ + (k m + km m u]/. The orbit of sstem (. with the initial value point B( k m+ k u, a 0 intersects with line = and the intersection point is C(,, then we derive the curve BC. Define CD as a length of line L 3 = = 0, DO as a length of line L = = 0. Since dl (. = k m+ k u < 0( > 0, dl 3 (. = (cβ a < 0(0 < <, dl (. = b(km m u > 0, the orbits of sstem (. go through into the interior of the Bendison ring from the outer of AB, CD and DO; on BC, compared sstem (. with sstem (., we have (. < (. < 0 and (. = (. > 0, then the orbits of sstem (. go through into the interior of the Bendison ring from the outer of BC. On the other hand, under the assumptions of Theorem., P (, is an unstable equilibrium. B Poincare Bendison Theorem, Sstem (. admits at least one limit ccle in region OA BCD Ω. This completes the proof of Theorem.. Theorem.. Suppose that 0 < u < km m, ak > ad3 u(cβ ad 3, and 0 < m < m. Sstem (. admits at most one limit ccle which is globall asmptoticall stable in Ω. (cβ add Proof. In order to prove Theorem., we take the following change of variables v =, u = ln, dτ =, still denote v, u, τ as,, t, then sstem (. takes the following form = b( + a[ + (k m + km m u] e +, d + (cβ ad =, and the positive equilibrium P (, P (, ln. (. (.
5 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( Taking translation transformation = +, = + ln, then P transform to the origin O(0, 0, still denote, as,, then sstem (. takes the following form where = e + b( + a( + [ ( + + (k m( + + km m u] + φ( F(( >, = d + (cβ ad( + g(, + F( = b( + a( + [ ( + + (k m( + + km m u] +, f ( = F ( = b[a( ( ak + am( + + km m u] ( +. Notice that under the assumptions of Theorem., P is an unstable equilibrium and tr(j(p > 0, thus f (0 = tr(j(p < 0 and g(0 = 0, when 0, g( > 0. Under the assumptions of Theorem., O(0, 0 is an unstable equilibrium of sstem (.3. Thus, if there eist closed orbits of sstem (.3, then define L as the nearest orbit to the equilibrium O and L is inboard stable. Define the parameter equation of the closed orbit L : = (t, = (t. B Lemma., [ ( φ( F( div(x, Y = L L = F ( L = f ( (t < 0, L + (g( then L f ( (t 0. Furthermore suppose that there eists a closed orbit L L of sstem (.3, where L is net to L and the parameter equation of the closed orbit L is = (t, = (t. Thereinafter, we depart the proof of Theorem. to two steps: (I f ( L (t > f ( L (t; (II L is impossible to be a semi-stable limit ccle. Therefore, if (II holds, then L must be stable; b (I, there eists no other closed orbits ecept L, which completes the proof of uniqueness. ] (.3 Proof (I. Constructing a function f ( = f ( + δg(. Suppose that Q, P are the least and most -coordinate points in L respectivel, and the must la on the curve φ( F( = 0. Let δ = f ( Q g( Q, then f (0 = f (0 < 0, f ( Q = 0. Doing a perpendicular line passing through Q, the perpendicular line intersects with -ais and the foot of a perpendicular M, then f ( M = f ( Q = 0. We shall show: (T f ( must intersect with positive -semi-ais and the -coordinate of the intersection point N satisfing N P ; (T For all ( M, N, f ( < 0; (T 3 For all (, M ( N, +, function f ( is non-decreasing. g( Proof (T. If (T does not hold, then on region L, f ( 0, and therefore f ( (t = f ( (t δ g( (t L L L = f ( (t δ L L = f ( (t < 0, L which contradicts f ( L (t 0. This completes the proof of (T. From the proof of (T, we derive f L ( i (t i > f ( L i i (t, i =,. Thus, (I can be rewrite as the following form (I f ( (t > f ( (t. L L
6 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( Proof (T. f ( = f ( + δg ( { } b[a( + 3 (km m u] = + δd ( + + ψ( ( +, ψ ( = b[a( (km m u] > 0, >. + Clearl, ψ( is increasing in (, +. If ψ( 0, then f ( 0, which contradicts f ( M = 0, f (0 < 0; if ψ( 0 then f ( 0, which contradicts f (0 < 0, f ( N = 0. Hence, there must eist a 0 ( M, N so that ψ( 0 = 0 and < 0, ψ( < 0, > 0, ψ( > 0, then f ( 0 = 0 and < 0, f ( < 0, > 0, f ( > 0. These combined with f ( M = f ( N = 0, impl that for all ( M, N, f ( < 0. Proof (T 3. where f ( g( = f ( g( bψ ( [(cβ ad( + d]( +, Ψ ( = (cβ ad( + [a( + 3 (km m u] d[a( ( ak + am( + (km m u]. From tr(j(p > 0 and d = (cβ ad, we derive Ψ ( = (cβ ad[a( + 3 a ( + ( ak + am( + (km m u], Ψ ( = (cβ ad[a( + ( ak + am]. Since tr(j(p > 0, ak + am < 0 and for all >, Ψ ( > 0, so, Ψ ( is increasing for all >. Since Ψ ( = (cβ ad(km m u < 0, Ψ (0 = (cβ ad[ a 3 ( ak + am (km m u] > 0, there eists a 0 (, 0 such that Ψ ( 0 = 0, where Ψ ( 0 = (cβ ad[a( a ( 0 + ( ak + am( 0 + (km m u] = 0, and hence, when < < 0, Ψ ( < 0; when > 0, Ψ ( > 0, so Ψ ( takes the least value at = 0. Notice that 0 < 0 and then a( + 3 a ( + (km m u = ( ak + am( +, Ψ ( 0 = (cβ ad( + [a( + 3 (km m u] (cβ ad [a( ( ak + am( + (km m u] = (cβ ad(m km + u + (cβ ad( + [a( (m km + u a ( + ( ak + am( + ] = (cβ ad(m km + u 0 a 0 (cβ ad( + 3 > 0. We derive Ψ ( > 0 for all >, which means for all (, M ( N, +, the function f ( is non-decreasing. B (T, (T, (T 3 and Green formula (the regions G, G, G 3, G refer to the figure in Zhang et al. [0], page, we have f ( (t f ( (t = ( ( f ( f ( + L L G G 3 φ( F( G G g( = f ( ( e (φ( + F( + f ( > 0, g( which completes the proof of (I. G G 3 Proof. (II. Assume that L is a semi-stable limit ccle and line = 3 < 0 intersects with L. Constructing a new function F( = F( εr(, where 0 < ε, { 0, 3, r( = (g( g( 3, < 3. G G g(
7 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( It is obvious that (I f ( > L L f (, holds, where f ( = F (. On the other hand, for the parameters 0 < ε < ε, we have φ( F( + ε r( g( φ( F( + ε r( g( = r(g((ε ε 0. Therefore, ( φ( F(, g( structures a rotation vector field in the wider sense on parameter ε. Considering = φ( F(, = g(, (. when ε = 0, sstem (. is transformed into sstem (.3. B Lemma.3, when 0 < ε, the semi-stable limit ccle L of sstem (.3 decomposes into at least two rings L L of sstem (. and L is inboard stable, L is outboard unstable, then f ( 0, f ( 0, L L which contradicts (I. This completes the proof of (II and then Theorem. follows.. The influence of pre refuge and harvesting efforts and simulations I. The influence of pre refuge on model (. Notice that the scaling from sstem (. to sstem (.: = m, = β, = (+a( mdτ, then the coordinates of the positive equilibrium P (, of sstem (. take = + m, and = cα (k ch kd d = cα ( ( d d m + k m ch kd cβ ad cβ ad d. d cβ ad Notice that is a continuous differential function of parameter m and dm = > 0. The above inequalit shows that is the strictl increasing function of parameter m and that increasing the pre refuge leads to the increasing of the densit of pre species. One could see that is also a continuous differential function of parameter m. Simple computation shows that dm = cα kd ( k m d cβ ad. We will discuss (. in two cases: Case : Assume that the inequalit k d 0 holds, then < 0 for all m > 0, thus cβ ad dm is the strictl decreasing function of parameter m. That is, increasing the amount of pre refuge can decrease the densit of predator species. In this cα case, reaches the maimum value k(cβ ad (k d ch at m = 0. cβ ad d Case : Assume that the inequalit k d > 0 holds, denoting cβ ad m = k d, then it follows that > 0 for all cβ ad dm m (0, m and < 0 for all m dm (m, k. Thus, there eists a threshold m = m such that is the strictl increasing function of parameter m for all m (0, m, while is the strictl decreasing function of parameter m for all m (m, k. attains its maimum value c(kα h at m = m. d The above analsis shows that in the interval (0, m, increasing the amount of pre refuge can increase the densit of predator species and this is happened due to predator species still have enough food for predation with m being small, while as m large than m and the pre refuge becomes large enough, then increasing the amount of pre refuge can decrease the densit of predator species and this is happened due to the loss of food for predator species. Letα =, k = 0, d =, a = c = 0., β =, h =. B simple computation, we have m =., u =, m = 0, k+ k u d = 9.993, ak =, ad3 u(cβ ad 3 =.. Then cβ ad (cβ add (. (.
8 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( Fig.. There is a stable limit ccle surrounding P (,. with m = Fig.. P (0,.33 is globall asmptoticall stable with m = 0. ( If 0 < m <., the conditions of Theorems. and. hold, so model (. admits eactl one limit ccle; ( If. < m < 9.999, (H holds, then P is globall stable; (3 If < m < 0, the conditions of Theorem. hold, so P 0 is globall stable. Figs. show that the densit of pre species is increasing while m is increasing and Figs. and show that the densit of predator species increase as m < 0 and Figs. 3 and show that the densit of predator species decrease as m > 0. Figs. also show that the dependence of the namic behavior of sstem (. on the pre refuge m. When m is small, there is a stable limit ccle surrounding the unique positive equilibrium, and when m is large enough, the limit ccle is broken and both the pre and predator population converge to their equilibrium values respectivel, which means that if we change the value of m, it is possible to prevent the cclic behavior of the predator pre sstem and to drive it to a required stable state. II. The influence of harvesting efforts on model (. We will discuss this for an fied parameter m. It is eas to calculate that: dh = c d < 0. The above analsis shows that increasing the harvesting effort h can decrease predator densit, which is obvious, because the increasing of h must lead to the loss of food for predator species. Furthermore, for fied m, we notice that the densit of pre species still have no change as t + and alwas keep the value +m, though increase harvesting pre species. Obviousl, the decreasing of densit of predator species can lead to the final immovabilit of densit of pre species. Let α =, k = 0, d =, a = c = 0., β =, m =. Figs. and shows that, the densities of pre species keep unchangeable when there eists the positive equilibrium, while the densities of predator species decrease with an increase in h. d cβ ad (.3
9 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( Fig. 3. P (3,. is globall asmptoticall stable with m = Fig.. P (3,.9 is globall asmptoticall stable with m = Fig.. P (,. is globall asmptoticall stable with h =.. Conclusion In this paper we consider a predator pre model with Holling tpe II functional response incorporating a constant pre refuge and constant-rate pre harvesting. When the assumption 0 < h < α(km m /k holds, we show that the constant-
10 0 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( Fig.. There is a stable limit ccle surrounding P (, 3. with h = 0. rate pre harvesting has no influence on the limit ccle. The sstem also admits a unique stable limit ccle when the positive equilibrium is unstable. Differentl from the conclusion of Chen [0] increasing the amount of refuge can increase pre densities and lead to population outbreaks, this paper verifies that increasing the amount of refuge can increase pre densities, and when the assumption k d 0 holds, increasing the amount of pre refuge can decrease the predator densities; when the assumption k cβ ad d cβ ad > 0 holds, there eists a threshold m, such that for the pre refuge smaller than this threshold, increasing the amount of pre refuge can increase the predator densities and if the pre refuge is larger than the threshold, increasing the amount of pre refuge can decrease the predator densities. The analsis of the influence of harvesting effort on equilibrium densit values indicates that for fied pre refuge, harvesting has no influence on the final densit of the pre species, while the densit of predator species is decreasing with the increasing of harvesting effort on pre species. Our results and numerical simulation also indicate that namic behavior of the model ver much depends on the pre refuge parameter m or harvesting effort h. Hence, it is possible to control the sstem in such a wa that the sstem approaches a required state, using the effort h or pre refuge m as controls. Acknowledgements This work is supported b the Natural Science Foundation of Fujian Province (00J009, and the Foundation of Science and Technolog of Fuzhou Universit (00-XQ-0. References [] C.S. Holling, Some characteristics of simple 0 tpes of predation and parasitism, The Canada Entomological 9 ( [] M.P. Hassell, The Dnamics of Arthropod Predator Pre Sstems, Princeton Univ. Press, Princeton, NJ, 9. [3] M.A. Ho, Almonds (California, in: W. Helle, M.W. Sabelis (Eds., Spider Mites: Their Biolog, Natural Enemies and Control, World Crop Pests, vol. B, Elsevier, Amsterdam, 9. [] A. Sih, Pre refuges and predator pre stabilit, Theoretical Population Biolog 3 (9. [] V. Krivan, Effects of optimal antipredator behavior of pre on predator pre namics:the role of refuges, Theoretical Population Biolog 3 (99 3. [] T.K. Kar, Stabilit analsis of a pre predator model incorporating a pre refuge, Communications in Nonlinear Science and Numerical Simulation 0 (00 9. [] W. Ko, K. Ru, Qualitative analsis of a predator pre model with Holling tpe II functional response incorporating a pre refuge, Journal of Differential Equations 3 ( [] J.B. Collings, Bifurcation and stabilit analsis of a temperature-dependent mite predator pre interaction model incorporating a pre refuge, Bulletin of Mathematical Biolog (99 3. [9] E. Gonzalez-Olivares, R. Ramos-Jiliberto, Dnamic consequences of pre refuges in a simple model sstem: More pre, fewer predators and enhanced stabilit, Ecological Modelling ( [0] L. Chen, F. Chen, Qualitative analsis of a predator pre model with Holling tpe II functional response incorporating a constant pre refuge, Nonlinear Analsis: Real World Applications (00 doi:0.0/j.nonrwa [] C.S. Holling, The functional response of predators to pre densit and its role in mimicr and population regulations, Memoirs of the Entomological Societ of Canada ( [] F. Brauer, A.C. Soudack, Stabilit regions and transition phenomena for harvested predator pre sstems, Journal of Mathematical Biolog ( [3] F. Brauer, A.C. Soudack, Stabilit regions in predator pre sstems with constant rate pre harvesting, Journal of Mathematical Biolog (99. [] F. Brauer, A.C. Soudack, Constant-rate stocking of predator pre sstems, Journal of Mathematical Biolog (9. [] F. Brauer, A.C. Soudack, Coeistence properties of some predator pre sstems under constant rate harvesting and stocking, Journal of Mathematical Biolog (9 0. [] F. Brauer, A.C. Soudack, On constant effort harvesting and stocking in a class of predator pre sstems, Journal of Theoretical Biolog 9 ( (9.
11 L. Ji, C. Wu / Nonlinear Analsis: Real World Applications ( [] D. Xiao, S. Ruan, Bogdanov Takens bifurcations in predator pre sstems with constant rate harvesting, Fields Institute Communications ( [] D. Xiao, L. Jennings, Bifurcations of a ratio-dependent predator pre sstem with constant rate harvesting, SIAM Applied Mathematics ( [9] G. Dai, M. Tang, Coeistence region and global namics of a harvested predator pre sstem, SIAM Journal on Applied Mathematics ( ( [0] Z.F. Zhang, T.R. Ding, W.Z. Huang, Z.X. Dong, Qualitative theor of differential equations, first ed., Science Publishing Compan, Beijing, 9.
A PREY-PREDATOR MODEL WITH HOLLING II FUNCTIONAL RESPONSE AND THE CARRYING CAPACITY OF PREDATOR DEPENDING ON ITS PREY
Journal of Applied Analsis and Computation Volume 8, Number 5, October 218, 1464 1474 Website:http://jaac-online.com/ DOI:1.11948/218.1464 A PREY-PREDATOR MODEL WITH HOLLING II FUNCTIONAL RESPONSE AND
More information* τσ σκ. Supporting Text. A. Stability Analysis of System 2
Supporting Tet A. Stabilit Analsis of Sstem In this Appendi, we stud the stabilit of the equilibria of sstem. If we redefine the sstem as, T when -, then there are at most three equilibria: E,, E κ -,,
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
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 informationNew approach to study the van der Pol equation for large damping
Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for
More informationHeteroclinic Bifurcation and Multistability in a Ratio-dependent Predator-Prey System with Michaelis-Menten Type Harvesting Rate
Proceedings of the World Congress on ngineering Vol I WC Jul 4-6 London U.K. Heteroclinic Bifurcation and Multistabilit in a Ratio-dependent Predator-Pre Sstem with Michaelis-Menten Tpe Harvesting Rate
More information7.7 LOTKA-VOLTERRA M ODELS
77 LOTKA-VOLTERRA M ODELS sstems, such as the undamped pendulum, enabled us to draw the phase plane for this sstem and view it globall In that case we were able to understand the motions for all initial
More informationLocal Phase Portrait of Nonlinear Systems Near Equilibria
Local Phase Portrait of Nonlinear Sstems Near Equilibria [1] Consider 1 = 6 1 1 3 1, = 3 1. ( ) (a) Find all equilibrium solutions of the sstem ( ). (b) For each equilibrium point, give the linear approimating
More informationInterspecific Segregation and Phase Transition in a Lattice Ecosystem with Intraspecific Competition
Interspecific Segregation and Phase Transition in a Lattice Ecosstem with Intraspecific Competition K. Tainaka a, M. Kushida a, Y. Ito a and J. Yoshimura a,b,c a Department of Sstems Engineering, Shizuoka
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 informationChapter 8 Equilibria in Nonlinear Systems
Chapter 8 Equilibria in Nonlinear Sstems Recall linearization for Nonlinear dnamical sstems in R n : X 0 = F (X) : if X 0 is an equilibrium, i.e., F (X 0 ) = 0; then its linearization is U 0 = AU; A =
More informationQualitative analysis of a harvested predator-prey system with Holling-type III functional response incorporating a prey refuge
Wang and Pan Advances in Difference Equations 2012, 2012:96 R E S E A R C H Open Access Qualitative analysis of a harvested pretor-prey system with Holling-type III functional response incorporating a
More informationEE222: Solutions to Homework 2
Spring 7 EE: Solutions to Homework 1. Draw the phase portrait of a reaction-diffusion sstem ẋ 1 = ( 1 )+ 1 (1 1 ) ẋ = ( 1 )+ (1 ). List the equilibria and their tpes. Does the sstem have limit ccles? (Hint:
More informationSTABILITY AND NEIMARK-SACKER BIFURCATION OF A SEMI-DISCRETE POPULATION MODEL
Journal of Applied Analsis and Computation Website:http://jaac-online.com/ Volume 4, Number 4, November 14 pp. 419 45 STABILITY AND NEIMARK-SACKER BIFURCATION OF A SEMI-DISCRETE POPULATION MODEL Cheng
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 informationThe approximation of piecewise linear membership functions and Łukasiewicz operators
Fuzz Sets and Sstems 54 25 275 286 www.elsevier.com/locate/fss The approimation of piecewise linear membership functions and Łukasiewicz operators József Dombi, Zsolt Gera Universit of Szeged, Institute
More informationLotka Volterra Model with Time Delay
International Journal of Mathematics Research. ISSN 976-584 Volume 6, Number (4), pp. 5- International Research Publication House http://www.irphouse.com Lotka Volterra Model with Time Dela Tapas Kumar
More informationMath 2930 Worksheet Equilibria and Stability
Math 2930 Worksheet Equilibria and Stabilit Week 3 September 7, 2017 Question 1. (a) Let C be the temperature (in Fahrenheit) of a cup of coffee that is cooling off to room temperature. Which of the following
More informationThe Hopf Bifurcation Theorem: Abstract. Introduction. Transversality condition; the eigenvalues cross the imginary axis with non-zero speed
Supercritical and Subcritical Hopf Bifurcations in Non Linear Maps Tarini Kumar Dutta, Department of Mathematics, Gauhati Universit Pramila Kumari Prajapati, Department of Mathematics, Gauhati Universit
More informationQUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS
arxiv:1803.0591v1 [math.gm] QUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS Aleander Panfilov stable spiral det A 6 3 5 4 non stable spiral D=0 stable node center non stable node saddle 1 tr A QUALITATIVE
More informationNonlinear Systems Examples Sheet: Solutions
Nonlinear Sstems Eamples Sheet: Solutions Mark Cannon, Michaelmas Term 7 Equilibrium points. (a). Solving ẋ =sin 4 3 =for gives =as an equilibrium point. This is the onl equilibrium because there is onl
More informationFunctional Response to Predators Holling type II, as a Function Refuge for Preys in Lotka-Volterra Model
Applied Mathematical Sciences, Vol. 9, 2015, no. 136, 6773-6781 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.53266 Functional Response to Predators Holling type II, as a Function Refuge
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 informationStability Analysis for Linear Systems under State Constraints
Stabilit Analsis for Linear Sstems under State Constraints Haijun Fang Abstract This paper revisits the problem of stabilit analsis for linear sstems under state constraints New and less conservative sufficient
More informationMULTIPLE MIXED-TYPE ATTRACTORS IN A COMPETITION MODEL. 1. Introduction
MULTIPLE MIXED-TYPE ATTRACTORS IN A COMPETITION MODEL J. M. CUSHING, SHANDELLE M. HENSON, AND CHANTEL C. BLACKBURN Abstract. We show that a discrete time, two species competition model with Ricker (eponential)
More informationCh 2.5: Autonomous Equations and Population Dynamics
Ch 2.5: Autonomous Equations and Population Dnamics In this section we examine equations of the form d/dt = f (), called autonomous equations, where the independent variable t does not appear explicitl.
More informationControl Schemes to Reduce Risk of Extinction in the Lotka-Volterra Predator-Prey Model
Journal of Applied Mathematics and Phsics, 2014, 2, 644-652 Published Online June 2014 in SciRes. http://www.scirp.org/journal/jamp http://d.doi.org/10.4236/jamp.2014.27071 Control Schemes to Reduce Risk
More informationProblem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form
Problem set 6 Math 207A, Fall 2011 s 1 A two-dimensional gradient sstem has the form x t = W (x,, x t = W (x, where W (x, is a given function (a If W is a quadratic function W (x, = 1 2 ax2 + bx + 1 2
More informationSolutions to Two Interesting Problems
Solutions to Two Interesting Problems Save the Lemming On each square of an n n chessboard is an arrow pointing to one of its eight neighbors (or off the board, if it s an edge square). However, arrows
More informationUsing MatContM in the study of a nonlinear map in economics
Journal of Phsics: Conference Series PAPER OPEN ACCESS Using MatContM in the stud of a nonlinear map in economics To cite this article: Niels Neirnck et al 016 J. Phs.: Conf. Ser. 69 0101 Related content
More informationResearch Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
Mathematical Problems in Engineering Volume, Article ID 88, pages doi:.//88 Research Article A New Four-Scroll Chaotic Attractor Consisted of Two-Scroll Transient Chaotic and Two-Scroll Ultimate Chaotic
More informationFebruary 11, JEL Classification: C72, D43, D44 Keywords: Discontinuous games, Bertrand game, Toehold, First- and Second- Price Auctions
Eistence of Mied Strateg Equilibria in a Class of Discontinuous Games with Unbounded Strateg Sets. Gian Luigi Albano and Aleander Matros Department of Economics and ELSE Universit College London Februar
More informationTwo Limit Cycles in a Two-Species Reaction
Two Limit Ccles in a Two-Species Reaction Brigita Ferčec 1 Ilona Nag Valer Romanovski 3 Gábor Szederkéni 4 and János Tóth 5 The Facult of Energ Technolog Krško 1 Center for Applied Mathematics and Theoretical
More informationResearch Article Chaotic Attractor Generation via a Simple Linear Time-Varying System
Discrete Dnamics in Nature and Societ Volume, Article ID 836, 8 pages doi:.//836 Research Article Chaotic Attractor Generation via a Simple Linear Time-Varing Sstem Baiu Ou and Desheng Liu Department of
More informationChaos Control and Synchronization of a Fractional-order Autonomous System
Chaos Control and Snchronization of a Fractional-order Autonomous Sstem WANG HONGWU Tianjin Universit, School of Management Weijin Street 9, 37 Tianjin Tianjin Universit of Science and Technolog College
More informationSubband Coding and Wavelets. National Chiao Tung University Chun-Jen Tsai 12/04/2014
Subband Coding and Wavelets National Chiao Tung Universit Chun-Jen Tsai /4/4 Concept of Subband Coding In transform coding, we use N (or N N) samples as the data transform unit Transform coefficients are
More informationx y plane is the plane in which the stresses act, yy xy xy Figure 3.5.1: non-zero stress components acting in the x y plane
3.5 Plane Stress This section is concerned with a special two-dimensional state of stress called plane stress. It is important for two reasons: () it arises in real components (particularl in thin components
More information(4p) See Theorem on pages in the course book. 3. Consider the following system of ODE: z 2. n j 1
MATEMATIK Datum -8- Tid eftermiddag GU, Chalmers Hjälpmedel inga A.Heintz Telefonvakt Aleei Heintz Tel. 76-786. Tenta i ODE och matematisk modellering, MMG, MVE6 Answer rst those questions that look simpler,
More informationGlobal dynamics and bifurcation of a tri-trophic food chain model
ISSN 746-7, England, UK World Journal of Modelling and Simulation Vol. 8 ) No., pp. 66-8 Global dnamics and bifurcation of a tri-trophic food chain model Tapan Kumar Kar, Abhijit Ghorai, Ashim Batabal
More informationChapter 5: Systems of Equations
Chapter : Sstems of Equations Section.: Sstems in Two Variables... 0 Section. Eercises... 9 Section.: Sstems in Three Variables... Section. Eercises... Section.: Linear Inequalities... Section.: Eercises.
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 informationSection 1.2: A Catalog of Functions
Section 1.: A Catalog of Functions As we discussed in the last section, in the sciences, we often tr to find an equation which models some given phenomenon in the real world - for eample, temperature as
More informationFirst-Order Differential Equations
First-Order Differential Equations. Solution Curves Without a Solution.. Direction Fields.. Autonomous First-Order DEs. Separable Equations.3 Linear Equations. Eact Equations.5 Solutions b Substitutions.6
More informationThe Sigmoidal Approximation of Łukasiewicz Operators
The Sigmoidal Approimation of Łukasiewicz Operators József Dombi, Zsolt Gera Universit of Szeged, Institute of Informatics e-mail: {dombi gera}@inf.u-szeged.hu Abstract: In this paper we propose an approimation
More informationStability Analysis of Mathematical Model for New Hydraulic Bilateral Rolling Shear
SJ nternational, Vol. 56 (06), SJ nternational, No. Vol. 56 (06), No., pp. 88 9 Stabilit Analsis of Mathematical Model for New Hdraulic Bilateral Rolling Shear QingXue HUANG, ) Jia L, ) HongZhou L, ) HeYong
More informationResearch Article Dynamics of a Predator-Prey System with a Mate-Finding Allee Effect
Abstract and Applied Analsis, Article ID 673424, 14 pages http://d.doi.org/1.1155/214/673424 Research Article Dnamics of a Predator-Pre Sstem with a Mate-Finding Allee Effect Ruiwen Wu and Xiuiang Liu
More informationJones & Bartlett Learning, LLC NOT FOR SALE OR DISTRIBUTION. Jones & Bartlett Learning, LLC NOT FOR SALE OR DISTRIBUTION
FIRST-ORDER DIFFERENTIAL EQUATIONS 2 Chapter Contents 2. Solution Curves without a Solution 2.. Direction Fields 2..2 Jones Autonomous & Bartlett First-Order Learning, DEs LLC 2.2 Separable Equations 2.3
More informationFrom bell-shaped solitary wave to W/M-shaped solitary wave solutions in an integrable nonlinear wave equation
PRAMANA c Indian Academ of Sciences Vol. 74, No. journal of Januar 00 phsics pp. 9 6 From bell-shaped solitar wave to W/M-shaped solitar wave solutions in an integrable nonlinear wave equation AIYONG CHEN,,,
More informationGlobal Analysis of an Epidemic Model with Nonmonotone Incidence Rate
Global Analysis of an Epidemic Model with Nonmonotone Incidence Rate Dongmei Xiao Department of Mathematics, Shanghai Jiaotong University, Shanghai 00030, China E-mail: xiaodm@sjtu.edu.cn and Shigui Ruan
More informationSpace frames. b) R z φ z. R x. Figure 1 Sign convention: a) Displacements; b) Reactions
Lecture notes: Structural Analsis II Space frames I. asic concepts. The design of a building is generall accomplished b considering the structure as an assemblage of planar frames, each of which is designed
More informationMAT 1275: Introduction to Mathematical Analysis. Graphs and Simplest Equations for Basic Trigonometric Functions. y=sin( x) Function
MAT 275: Introduction to Mathematical Analsis Dr. A. Rozenblum Graphs and Simplest Equations for Basic Trigonometric Functions We consider here three basic functions: sine, cosine and tangent. For them,
More informationJones & Bartlett Learning, LLC, an Ascend Learning Company. NOT FOR SALE OR DISTRIBUTION
stefanel/shutterstock, Inc. CHAPTER 2 First-Order Differential Equations We begin our stu of differential equations with first-order equations. In this chapter we illustrate the three different was differential
More informationTHE most studied topics in theoretical ecology is the
Engineering Letters :3 EL 3_5 Stabilit Nonlinear Oscillations Bifurcation in a Dela-Induced Predator-Pre Sstem with Harvesting Debaldev Jana Swapan Charabort Nadulal Bairagi Abstract A harvested predator-pre
More informationLecture 7: Advanced Coupling & Mixing Time via Eigenvalues
Counting and Sampling Fall 07 Lecture 7: Advanced Coupling & Miing Time via Eigenvalues Lecturer: Shaan Oveis Gharan October 8 Disclaimer: These notes have not been subjected to the usual scrutin reserved
More informationResearch Article Chaos and Control of Game Model Based on Heterogeneous Expectations in Electric Power Triopoly
Discrete Dnamics in Nature and Societ Volume 29, Article ID 469564, 8 pages doi:.55/29/469564 Research Article Chaos and Control of Game Model Based on Heterogeneous Epectations in Electric Power Triopol
More informationEstimators in simple random sampling: Searls approach
Songklanakarin J Sci Technol 35 (6 749-760 Nov - Dec 013 http://wwwsjstpsuacth Original Article Estimators in simple random sampling: Searls approach Jirawan Jitthavech* and Vichit Lorchirachoonkul School
More informationMATHEMATICAL MODELING AND THE STABILITY STUDY OF SOME CHEMICAL PHENOMENA
HENRI COND IR FORCE CDEM ROMNI INTERNTIONL CONFERENCE of SCIENTIFIC PPER FSES rasov -6 Ma GENERL M.R. STEFNIK RMED FORCES CDEM SLOVK REPULIC MTHEMTICL MODELING ND THE STILIT STUD OF SOME CHEMICL PHENOMEN
More informationOPTIMAL CONTROL FOR A PARABOLIC SYSTEM MODELLING CHEMOTAXIS
Trends in Mathematics Information Center for Mathematical Sciences Volume 6, Number 1, June, 23, Pages 45 49 OPTIMAL CONTROL FOR A PARABOLIC SYSTEM MODELLING CHEMOTAXIS SANG UK RYU Abstract. We stud the
More informationQUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS
Journal of the Operations Research Societ of Japan 008, Vol. 51, No., 191-01 QUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS Tomonari Kitahara Shinji Mizuno Kazuhide Nakata Toko Institute of Technolog
More informationDynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems
. Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak
More informationHeteroclinic Bifurcation of a Predator-Prey System with Hassell-Varley Functional Response and Allee Effect
International Journal of ngineering Research And Management (IJRM) ISSN: 49-058 Volume-05 Issue-0 October 08 Heteroclinic Bifurcation of a Predator-Prey System with Hassell-Varley Functional Resonse and
More informationFIRST- AND SECOND-ORDER IVPS The problem given in (1) is also called an nth-order initial-value problem. For example, Solve: Solve:
.2 INITIAL-VALUE PROBLEMS 3.2 INITIAL-VALUE PROBLEMS REVIEW MATERIAL Normal form of a DE Solution of a DE Famil of solutions INTRODUCTION We are often interested in problems in which we seek a solution
More informationEffects of birefringence on Fizeau interferometry that uses polarization phase shifting technique
Effects of birefringence on Fizeau interferometr that uses polarization phase shifting technique Chunu Zhao, Dongel Kang and James H. Burge College of Optical Sciences, the Universit of Arizona 1630 E.
More informationReview of Essential Skills and Knowledge
Review of Essential Skills and Knowledge R Eponent Laws...50 R Epanding and Simplifing Polnomial Epressions...5 R 3 Factoring Polnomial Epressions...5 R Working with Rational Epressions...55 R 5 Slope
More informationLinear and Nonlinear Systems of Equations. The Method of Substitution. Equation 1 Equation 2. Check (2, 1) in Equation 1 and Equation 2: 2x y 5?
3330_070.qd 96 /5/05 Chapter 7 7. 9:39 AM Page 96 Sstems of Equations and Inequalities Linear and Nonlinear Sstems of Equations What ou should learn Use the method of substitution to solve sstems of linear
More information= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :
1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationHeterogeneous Strategies in Nonlinear Duopoly with Product Differentiation
Heterogeneous Strategies in Nonlinear Duopol with Product Differentiation Akio Matsumoto Department of Economics Chuo Universit, Toko, Japan akiom@tamacc.chuo-u.ac.jp Tamotsu Onozaki Facult of Management
More informationFractal dimension of the controlled Julia sets of the output duopoly competing evolution model
Available online at www.isr-publications.com/jmcs J. Math. Computer Sci. 1 (1) 1 71 Research Article Fractal dimension of the controlled Julia sets of the output duopol competing evolution model Zhaoqing
More informationPhase Portraits of 1-D Autonomous Equations
Phase Portraits of 1-D Autonomous Equations In each of the following problems [1]-[5]: (a) find all equilibrium solutions; (b) determine whether each of the equilibrium solutions is stable, asmptoticall
More informationMathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes
Mathematics 309 Conic sections and their applicationsn Chapter 2. Quadric figures In this chapter want to outline quickl how to decide what figure associated in 2D and 3D to quadratic equations look like.
More informationGlobal dynamics of a predator-prey system with Holling type II functional response
4 Nonlinear Analysis: Modelling and Control, 011, Vol. 16, No., 4 53 Global dynamics of a predator-prey system with Holling type II functional response Xiaohong Tian, Rui Xu Institute of Applied Mathematics
More informationA route to computational chaos revisited: noninvertibility and the breakup of an invariant circle
A route to computational chaos revisited: noninvertibilit and the breakup of an invariant circle Christos E. Frouzakis Combustion Research Laborator, Paul Scherrer Institute CH-5232, Villigen, Switzerland
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics http://jipam.vu.edu.au/ Volume 2, Issue 2, Article 15, 21 ON SOME FUNDAMENTAL INTEGRAL INEQUALITIES AND THEIR DISCRETE ANALOGUES B.G. PACHPATTE DEPARTMENT
More information14.1 Systems of Linear Equations in Two Variables
86 Chapter 1 Sstems of Equations and Matrices 1.1 Sstems of Linear Equations in Two Variables Use the method of substitution to solve sstems of equations in two variables. Use the method of elimination
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
More informationStochastic Interpretation for the Arimoto Algorithm
Stochastic Interpretation for the Arimoto Algorithm Serge Tridenski EE - Sstems Department Tel Aviv Universit Tel Aviv, Israel Email: sergetr@post.tau.ac.il Ram Zamir EE - Sstems Department Tel Aviv Universit
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationYou don't have to be a mathematician to have a feel for numbers. John Forbes Nash, Jr.
Course Title: Real Analsis Course Code: MTH3 Course instructor: Dr. Atiq ur Rehman Class: MSc-II Course URL: www.mathcit.org/atiq/fa5-mth3 You don't have to be a mathematician to have a feel for numbers.
More informationSection 1.5 Formal definitions of limits
Section.5 Formal definitions of limits (3/908) Overview: The definitions of the various tpes of limits in previous sections involve phrases such as arbitraril close, sufficientl close, arbitraril large,
More information8.4. If we let x denote the number of gallons pumped, then the price y in dollars can $ $1.70 $ $1.70 $ $1.70 $ $1.
8.4 An Introduction to Functions: Linear Functions, Applications, and Models We often describe one quantit in terms of another; for eample, the growth of a plant is related to the amount of light it receives,
More informationChapter 1 Coordinates, points and lines
Cambridge Universit Press 978--36-6000-7 Cambridge International AS and A Level Mathematics: Pure Mathematics Coursebook Hugh Neill, Douglas Quadling, Julian Gilbe Ecerpt Chapter Coordinates, points and
More informationTHE collective control of multi-agent systems is a rapidly
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 52, NO. 5, MAY 27 1 Stabilization of Planar Collective Motion: All-to-all Communication Rodolphe Sepulchre, Member, IEEE, Derek A. Pale, Student Member, IEEE,
More informationFlows and Connectivity
Chapter 4 Flows and Connectivit 4. Network Flows Capacit Networks A network N is a connected weighted loopless graph (G,w) with two specified vertices and, called the source and the sink, respectivel,
More informationAN EXTENSION OF THE ANTOCI-DEI-GALEOTTI EVOLUTIONARY MODEL FOR ENVIRONMENT PROTECTION THROUGH FINANCIAL INSTRUMENTS
ISSN 1974-4110 WP-EMS Working Papers Series in Economics, Mathematics and Statistics AN EXTENSION OF THE ANTOCI-DEI-GALEOTTI EVOLUTIONARY MODEL FOR ENVIRONMENT PROTECTION THROUGH FINANCIAL INSTRUMENTS
More informationTwo Dimensional Linear Systems of ODEs
34 CHAPTER 3 Two Dimensional Linear Sstems of ODEs A first-der, autonomous, homogeneous linear sstem of two ODEs has the fm x t ax + b, t cx + d where a, b, c, d are real constants The matrix fm is 31
More informationInternational Examinations. Advanced Level Mathematics Pure Mathematics 1 Hugh Neill and Douglas Quadling
International Eaminations Advanced Level Mathematics Pure Mathematics Hugh Neill and Douglas Quadling PUBLISHED BY THE PRESS SYNDICATE OF THE UNIVERSITY OF CAMBRIDGE The Pitt Building, Trumpington Street,
More informationAdditional Material On Recursive Sequences
Penn State Altoona MATH 141 Additional Material On Recursive Sequences 1. Graphical Analsis Cobweb Diagrams Consider a generic recursive sequence { an+1 = f(a n ), n = 1,, 3,..., = Given initial value.
More informationMath Notes on sections 7.8,9.1, and 9.3. Derivation of a solution in the repeated roots case: 3 4 A = 1 1. x =e t : + e t w 2.
Math 7 Notes on sections 7.8,9., and 9.3. Derivation of a solution in the repeated roots case We consider the eample = A where 3 4 A = The onl eigenvalue is = ; and there is onl one linearl independent
More informationAutonomous Equations / Stability of Equilibrium Solutions. y = f (y).
Autonomous Equations / Stabilit of Equilibrium Solutions First order autonomous equations, Equilibrium solutions, Stabilit, Longterm behavior of solutions, direction fields, Population dnamics and logistic
More informationCHAPTER 2. First-Order Differential Equations CHAPTER CONTENTS. Jones & Bartlett Learning, LLC NOT FOR SALE OR DISTRIBUTION
CHAPTER First-Order Differential Equations CHAPTER CONTENTS.1 Solution Curves Without NOT a FOR Solution SALE OR DISTRIBUTION.1.1 Direction Fields.1. Autonomous First-Order DEs. Separable Equations.3 Linear
More informationF and G have continuous second-order derivatives. Assume Equation (1.1) possesses an equilibrium point (x*,y*) so that
3.1.6. Characterizing solutions (continued) C4. Stability analysis for nonlinear dynamical systems In many economic problems, the equations describing the evolution of a system are non-linear. The behaviour
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More information8.7 Systems of Non-Linear Equations and Inequalities
8.7 Sstems of Non-Linear Equations and Inequalities 67 8.7 Sstems of Non-Linear Equations and Inequalities In this section, we stud sstems of non-linear equations and inequalities. Unlike the sstems of
More informationTHE collective control of multi-agent systems is a rapidly
SUBMITTED TO IEEE TRASACTIOS O AUTOMATIC COTROL 1 Stabilization of Planar Collective Motion: All-to-all Communication Rodolphe Sepulchre, Member, IEEE, Derek A. Pale, Student Member, IEEE, and aomi Ehrich
More informationPermanence and global attractivity of a nonautonomous modified Leslie-Gower predator-prey model with Holling-type II schemes and a prey refuge
Xie et al. Advances in Difference Equations 016 016:184 DOI 10.1186/s1366-016-089-5 R E S E A R C H Open Access Permanence and global attractivity of a nonautonomous modified Leslie-Gower predator-prey
More informationDynamical aspects of multi-round horseshoe-shaped homoclinic orbits in the RTBP
Dnamical aspects of multi-round horseshoe-shaped homoclinic orbits in the RTBP E. Barrabés, J.M. Mondelo, M. Ollé December, 28 Abstract We consider the planar Restricted Three-Bod problem and the collinear
More informationC. Non-linear Difference and Differential Equations: Linearization and Phase Diagram Technique
C. Non-linear Difference and Differential Equations: Linearization and Phase Diaram Technique So far we have discussed methods of solvin linear difference and differential equations. Let us now discuss
More informationIntermediate Math Circles Wednesday November Inequalities and Linear Optimization
WWW.CEMC.UWATERLOO.CA The CENTRE for EDUCATION in MATHEMATICS and COMPUTING Intermediate Math Circles Wednesda November 21 2012 Inequalities and Linear Optimization Review: Our goal is to solve sstems
More information4 Inverse function theorem
Tel Aviv Universit, 2013/14 Analsis-III,IV 53 4 Inverse function theorem 4a What is the problem................ 53 4b Simple observations before the theorem..... 54 4c The theorem.....................
More information