Dynamics of a Delayed Competitive System Affected by Toxic Substances with Imprecise Biological Parameters
|
|
- Arnold Patrick
- 6 years ago
- Views:
Transcription
1 Filomat 31:16 17, Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: Dynamics of a Delayed Competitive System Affected by Toxic Substances with Imprecise Biological Parameters A. K. Pal a, P. Dolai b, G. P. Samanta b a Department of Mathematics, S. A. Jaipuria College, Kolkata-7, INDIA b Department of Mathematics, Indian Institute of Engineering Science and Technology, Shibpur, Howrah , INDIA Abstract. In this paper we have studied the dynamical behaviours of a delayed two-species competitive system affected by toxicant with imprecise biological parameters. We have proposed a method to handle these imprecise parameters by using parametric form of interval numbers. We have discussed the existence of various equilibrium points and stability of the system at these equilibrium points. In case of toxic stimulatory system, the delay model exhibits a stable limit cycle oscillation. Computer simulations are carried out to illustrate our analytical findings. 1. Introduction Environmental pollution is a burning issue for the industrialized society in the world today because it leads to damage the both terrestrial and aquatic environments. Various kind of industrial discharges and chemical wastes are polluting the air and contaminating the streams, lakes, rivers and oceans with varieties of toxicant and chemicals such as arsenic, lead, cadmium, zinc, copper, iron, mercury etc. [1, 8]. There are many species which becomes extinct and several others are at the verge of extinction because of uncontrolled effects of toxicant to the environment. Therefore, the study of the effects of the toxic substances on ecological communities are becoming quite important from both environmental and conservation point of view. The problem of estimating qualitatively the effect of toxic substances on the species by using mathematical models is a very effective way. The deterministic dynamic models with the effect of toxic substances on various ecosystems was analysed by Hallam et al.[16, 17], Hallam and De Luna [18], De Luna and Hallam [8], Freedman and Shukla [1], Ghosh et al. [11], He and Wang [], Das et al.[7] and many others. Another important observation made by researchers is that the increased population of one species might affect the growth of another species or of several other species through the production of allelopathic toxins or stimulators, thus influencing seasonal succession [38]. The toxin produced by the unicellular green alga, chlorella Vulgaris, is an autotoxin that limits the size of its own population [36, 37] as well as inhibits the growth of the planktonic algae Asterionella formosa and Nitzschia frustrum Bacillareae[39]. Several researchers have also observed toxic inhibition of phytoplankton by other phytoplankton as well as some of these algae produce auxin which stimulate the growth of the other algae [1]. Such allelopathic stimulators 1 Mathematics Subject Classification. Primary 9B; Secondary 9D, 9D3 Keywords. Competing species; Imprecise parameters; Interval number; Time delay; Hopf bifurcation. Received: 4 December ; Revised: 9 February 16; Accepted: May 16 Communicated by Miljana Jovanović addresses: akpal_@yahoo.co.in A. K. Pal, pralay_dolai@yahoo.in P. Dolai, g_p_samanta@yahoo.co.uk, gpsamanta@math.iiests.ac.in G. P. Samanta
2 A. K. Pal et al. / Filomat 31:16 17, and inhibitors certainly affect algal succession, blooms and pulses by causing stimulated inhibited species to have a selective advantage disadvantage in competition [, 38]. It has been suggested by Rice [38] that all meaningful functional ecological models will eventually have to include a category on allelopathic and allelochemic effects. In nature, time delays occur in almost every biological situation [6] and assume to be one of the causes of regular fluctuations on population biomass. In population dynamics, a time delay is introduced when the rate of change of population biomass is not only a function of the present population biomass but also depends on the past population biomass. Therefore, time delay can be incorporated in the mathematical population model due to various ways such as maturation time, capturing time and other reasons. Moreover, existence of time delays is frequently a source of instability in some way. Many researchers[4 6,, 33, 34, 4, 44] have introduced time delay in their respective biological models to make it more realistic. Most of the researches in theoretical ecology have considered models based on the assumption that the biological parameters are precisely known. But in real world, the values of all the biological parameters can not be known precisely for the lack of information, lack of data, mistakes done in the measurement process and determining the initial conditions. To overcome these difficulties, imprecise model is more realistic in the field of mathematical biology. The impreciseness of the bio-mathematical model is also occurring due to environmental fluctuations or due to imprecise biological phenomenon. There are several approaches to handle such models having imprecise parameters such as stochastic approach, fuzzy approach, fuzzy stochastic approach etc. In stochastic approach, the imprecise parameters are replaced by random variables with known probability distributions. In fuzzy approach, the imprecise parameters are replaced by fuzzy sets with known membership functions. In fuzzy stochastic approach, some parameters are as fuzzy in nature and rest of the parameters are taken as random variables. However it is very difficult to construct a suitable membership function or a suitable probability distribution for each of the imprecise biological parameters. Some researchers have introduced fuzzy models in predator-prey system such as Bassanezi et al. [3], Peixoto et al. [3], Guo et al. [14], Pal et al.[9 3], Sharma and Samanta [41] have presented an optimal harvesting predator-prey system with imprecise biological parameters and they have also exercised about bio-economic equilibrium and optimal harvesting policy. In this paper, we have considered a delayed two species competitive system which is affected by toxicant [7]. To make the delay model more realistic, we have considered imprecise biological parameters as interval numbers. We present the interval numbers in parametric function form and study the parametric model. The dynamical behaviour of the parametric model is investigated for different values of the parameter p [, 1]. In section, we discuss some basic definitions on interval numbers. In section 3, we present mathematical form of two species competitive system with time delay. Section 4 briefly describes the delayed two species competitive system with imprecise parameters. Then parametric form of the system is formulated for the study of different behaviours of the model. Section deals with the existence of equilibrium points of the system and stability. Our important analytical results are numerically verified in section 7. Finally, section 8 contains the conclusions of the paper.. Basic Definitions In this section we discuss some basic definitions of the interval number and interval-valued function which have been used to study the imprecise competition model. Definition 1 : Interval number An interval number A is represented by closed interval [a l, a u ] and defined by A = [a l, a u ] = {x : a l x a u, x R}, where R is the set of real numbers and a l, a u are the lower and upper limits of the interval number respectively. So, every real number can also be presented by the interval number [a, a], for all a R. Definition :Interval-valued function Let us consider the interval [a, b] where a >. The interval [a, b] can be expressed as a function φp = a 1 p b p for p [, 1]. This function is called interval-valued function.
3 A. K. Pal et al. / Filomat 31:16 17, Now we define some arithmetic operations on interval-valued functions. Let A = [a l, a u ] and B = [b l, b u ] be two interval numbers. Addition : A + B = [a l, a u ] + [b l, b u ] = [a l + b l, a u + b u ]. The interval-valued function for the interval number A + B is given by φp = a 1 p L a p U where a L = a l + b l, and a U = a u + b u. Subtraction : A B = [a l, a u ] [b l, b u ] = [a l b u, a u b l ], provided a l b u >. The interval-valued function for the interval A B is given by φp = b 1 p b p L U where b L = a l b u, and b U = a u b l. [αa l, αa u ], if α Scalar multiplication : αa = α[a l, a u ] = }, provided a l >. The interval-valued function for the interval αa is given by φp = 3. Basic Mathematical Model [αa u, αa l ], if α < c 1 p c p L U if α, where c L = αa l, c U = αa u, d L = α a l, d U = α a u. d p if α < U L d 1 p Maynerd-Smith [6] have considered a two species competing system : dn 1 dn = N 1 t[k 1 α 1 N 1 t β 1 N t γ 1 N 1 tn t] = N t[k α N t β 1 N 1 t γ N 1 tn t] 1 with initial data N 1 >, N >. Where N 1 t, N t denote the population biomass of two competing species at time t, having common food resources ; K 1, K, α 1, α, β 1, β 1, γ 1, γ are positive constants. Here K i is the intrinsic growth rate of species i, α i represents the intra-species competition coefficient of species i, β ij denote the inter-species competition rate of jth species upon the ith species and γ i denote the toxic inhibition rate of the ith species i j, i, j = 1,. Samanta [4] also analyzed dynamical behaviors of such model where a toxic substance is produced at a constant rate. When the toxic coefficients γ i > i = 1,, the model 1 represents toxic inhibited species system. If we assume γ i = γ i = 1,, where γ >, then the model 1 represents toxic stimulated system. i i We assume that each species produces a substance toxic to the other, but only when the other is present. Further, it is also reasonable to assume such model where a toxic substance is produced to the competing species will not be instantaneous, and takes discrete time lag which is regarded as maturity period of the species. Here we introduce the time lag τ for the maturity of the second species. Then the system 1 reduces to dn 1 = N 1 t[k 1 α 1 N 1 t β 1 N t τ γ 1 N 1 tn t τ] dn = N t[k α N t β 1 N 1 t γ N 1 tn t] The initial conditions for the model 3 take the form N 1 θ = ψ 1 θ, N θ = ψ θ, θ [ τ, ], ψ 1 >, ψ >, ψ 1 θ, ψ θ C[ τ, ], R +. 3
4 3.1. Positivity and Boundedness of the delayed system A. K. Pal et al. / Filomat 31:16 17, In theoretical ecology, positivity and boundedness of a system implies that the system is biologically well posed. The following proposition ensures the positivity and boundedness of the delayed system 3. Proposition 1 : Each component of the solution of the delayed system 3 is positive and bounded for all t. Proof : Since the right hand side of the system 3 is completely continuous and locally Lipschitzian on C, the solution N 1 t, N t of 3 exists and is unique on [, ξ, where < ξ []. From system 3, we have [ t ] N 1 t = N 1 exp K 1 α 1 N 1 s β 1 N s τ γ 1 N 1 sn s τds >, [ t ] N t = N exp K α N s β 1 N 1 s γ N 1 sn sds > Therefore, N 1 t >, N t >, t. Now we assume, Then, dw W = N 1 t + N t. α 1 N 1 K 1 α 1 α N K α + K 1 α 1 + K α K 1 N 1 K N µ νw, where µ = K 1 α 1 + K α, and ν = min {K 1, K }. Therefore, dw + νw µ. Applying a theorem on differential inequalities [4], we obtain and for t, < WN 1, N µ ν + WN 1, N e νt. < W µ ν. Thus, all solutions of the system 3 enter into the region B = {N 1, N : < W < µ } ν + ε for any ε >. This proves the theorem. 3.. Uniform Persistence of the system To prove the uniform persistence permanent of the system 1, we shall use the Average Lyapunov function[13, 43] Proposition : System 1 is uniform persistence if β 1 α < K 1 K < α 1 β 1 Proof: We consider the average Lyapunov function of the form VN 1, N = N β 1 1 Nβ 1, where β 1, β 1 are positive constants. In the interior of R+,we have V V = ψn 1, N = β 1 K 1 + β 1 K β 1 α 1 + β 1 N 1 β 1 β 1 + +α N β 1 γ 1 + β 1 γ N 1 N
5 A. K. Pal et al. / Filomat 31:16 17, To prove the uniform persistence of the system, we shall have to show that ψn 1, N > for all boundary equilibria of the system. The following conditions should be satisfied for equilibrium: E 1 k 1 α 1, ; K α 1 K 1 β 1 > C 1 E, k α : K 1 α K β 1 > C It can be easily shown that conditions C 1 and C are satisfied if the conditions stated in the proposition holds. 4. Imprecise Competition Model According to formation of the competition model 3, all the parameters are positive and precise. But due to lack of proper information of the data, the parameters are not always precise. Now if any of the parameters K i, α i, β ij, γ i i j, i, j = 1, are imprecise, i.e., if any parameter is interval number rather than a single value, then it become difficult to convert the equation to the standard form and analyse the dynamical behaviour of the system. For imprecise parameters, we present the system 3 with interval parameters as described below : 4.1. Competition model with interval parameters Case I : Toxic inhibition Let ˆK i, ˆα i, βˆ ij, ˆγ i i j, i, j = 1, be the interval counterparts of K i, α i, β ij, γ i respectively. Then the imprecise competition delay model 3 becomes : dn 1 dn = N 1 t[ Kˆ 1 αˆ 1 N 1 t βˆ 1 N t τ γˆ 1 N 1 tn t τ] = N t[ Kˆ αˆ N t βˆ 1 N 1 t γˆ N 1 tn t] 4 where ˆK i = [K il, K iu ], ˆα i = [α il, α iu ], j, i, j = 1,. Case II : Toxic stimulation ˆ β ij = [β ijl, β iju ], ˆγ i = [γ il, γ iu ], and K il >, α il >, β ijl >, γ il > i In this case, γˆ 1 = γ ˆ 1 and γˆ = γ ˆ where γ ˆ 1 = [γ 1l, γ 1u ] and γ ˆ = [γ l, γ u ] and γ ˆ 1l >, γ ˆ >, i.e., the last l term in each of the equations of system 4 gives a positive effect instead of negative effect as in case I. Then the imprecise competition delay model 3 becomes replacing γ ˆ 1 by γˆ 1 and γ ˆ by γˆ : dn 1 dn = N 1 t[ Kˆ 1 αˆ 1 N 1 t βˆ 1 N t τ + γˆ 1 N 1 tn t τ] = N t[ Kˆ αˆ N t βˆ 1 N 1 t + γˆ N 1 tn t]
6 A. K. Pal et al. / Filomat 31:16 17, Competition model with parametric interval parameters For fixed m, let us consider the interval-valued function φ m p = a 1 p m b p m, p [, 1] for an interval [a m, b m ]. Since φ m p is a strictly increasing and continuous function, the system 4 and can be written in the parametric form as follows : dn 1 t; p = N 1 [K 1l 1 p K 1u p α 1l 1 p α 1u p N 1 β 1l 1 p β 1u p N t τ γ 1l 1 p γ 1u p N 1 N t τ] dn t; p = N [K l 1 p K u p α l 1 p α u p N β 1l 1 p β 1u p N 1 γ l 1 p γ u p N 1 N ] 6 where p [, 1], and dn 1 t; p = N 1 [K 1l 1 p K 1u p α 1l 1 p α 1u p N 1 β 1l 1 p β 1u p N t τ + γ 1l 1 p γ 1u p N 1 N t τ] dn t; p = N [K l 1 p K u p α l 1 p α u p N β 1l 1 p β 1u p N 1 + γ l 1 p γ u p N 1 N ] 7 where p [, 1].. Equilibria and Local Stability of the Delayed System Case I : Toxic inhibition The system 6 has four positive steady states, namely i E,, the trivial equilibrium, ii E 1 N 1, and E, N, the axial equilibrium, where N 1 = K 1l 1 p K 1u p α 1l 1 p α 1u, N p = K l 1 p K u p, for all p [, 1] 8 α l 1 p α u p and iii E = N 1, N, the interior equilibrium, where N 1, N can be determined by a ij N i + b ij N i + c ij =, i j, i, j = 1, 9 where a ij = β ijl 1 p β iju p γ il 1 p γ iu p α il 1 p α iu p γ jl 1 p γ ju p, b ij = K il 1 p K iu p γ jl 1 p γ ju p K jl 1 p K ju p γ il 1 p γ iu p α il 1 p α iu p α jl 1 p α ju p + β ijl 1 p β iju p β jil 1 p β jiu p, 1 for all p [, 1]. Then c ij = K il 1 p K iu p α jl 1 p α ju p K jl 1 p K ju p β jl 1 p β ju p, N i = 1 b ij ± b a ij 4a ijc ij, i, j = 1, 11 ij
7 exists with the conditions A. K. Pal et al. / Filomat 31:16 17, a ij, b ij 4a ijc ij. 1 The variational matrix of the system 6 at E, is given by [ ] K1l 1 p K V E = 1u p K l 1 p K u p. Clearly, E, is a saddle node and hence unstable. The variational matrix of the system 6 at E 1 N 1, is given by [ α1l V E 1 1 p α = 1u p N 1 β 1l 1 p β 1u p N 1 γ 1l 1 p γ 1u p K l 1 p K u p β 1l 1 p β 1u p N 1 e λτ N 1 then the eigenvalues are λ 1 = α 1l 1 p α 1u p N 1 and λ = K l 1 p K u p β 1l 1 p β 1u p N 1. Therefore the equilibrium point E 1 N 1, is asymptotically stable if K 1l 1 p K 1u p K l 1 p K u > α 1l 1 p α 1u p p β 1l 1 p β 1u p and unstable if K 1l 1 p K 1u p K l 1 p K u < α 1l 1 p α 1u p p β 1l 1 p β 1u. 13 p Again, the variational matrix of the system 6 at E, N is given by [ ] K1l 1 p K V E = 1u p β 1l 1 p β 1u p N e λτ N {β 1l 1 p β 1u p + γ l 1 p γ u p N } α l 1 p α u p, then the characteristic equation of the system at E is of the form: {λ + α l 1 p α u p N }{λ K 1l 1 p K 1u p + β 1l 1 p β 1u p N e λτ } =. Here λ = α l 1 p α u p N is a eigenvalue. Now we consider the equation λ = K 1l 1 p K 1u p β 1l 1 p β 1u p N e λτ. 14 If τ =, then the equilibrium point E, N is locally asymptotically stable if K l 1 p K u p K 1l 1 p K 1u > α l 1 p α u p p β 1l 1 p β 1u and unstable if p K l 1 p K u p K 1l 1 p K 1u p < α l 1 p α u p β 1l 1 p β 1u p. Again by substituting λ = iη in equation 14 and equating real and imaginary parts, we obtain eliminating τ, we have K 1l 1 p K 1u p = β 1l 1 p β 1u p N cos ητ η = β 1l 1 p β 1u p N sin ητ, η = {β 1l 1 p β 1u p } N. 16 We know that the equation 16 has a positive root η + if K l 1 p K u p K 1l 1 p K 1u > α l 1 p α u p p β 1l 1 p β 1u. p Therefore, there is a positive constant τ + such that for τ > τ +, E, N is unstable. Hence combining 13 and, the condition required for the persistence of both the species is K il 1 p K iu p K jl 1 p K ju < α il 1 p α iu p i j; i, j = 1, for all p [, 1]. 17 p β jil 1 p β jiu p N ],
8 A. K. Pal et al. / Filomat 31:16 17, Condition 17 gives c ij > in equation 9. Then the system 6 has unique positive equilibrium if a ij < i.e., if α il 1 p α iu p β jil 1 p β jiu > γ il 1 p γ iu p i j; i, j = 1, for all p [, 1]. 18 p γ jl 1 p γ ju p So, combining equation 17 and 18, we have α il 1 p α iu p β jil 1 p β jiu > max{ γ il 1 p γ iu p p γ jl 1 p γ ju, K il 1 p K iu p } i j; i, j = 1, for all p [, 1], 19 p K jl 1 p K ju p as the condition of existence of unique positive interior equilibrium of the system 6. Now to investigate the local stability of the interior equilibrium E N 1, N, we linearize the system 6 by using the following transformations: N 1 = N 1 + n 1, N = N + n. Then the linear system is given by where and du = MUt + NUt τ, Ut = [n 1 n ] T, M = m ij, N = n ij, m 11 = N 1 [α 1l 1 p α 1u p + γ 1l 1 p γ 1u p N ], m 1 =, m 1 = N [β 1l 1 p β 1u p + γ l 1 p γ u p N ], m = N [α l 1 p α u p + γ l 1 p γ u p N 1 ], n 1 = N 1 [β 1l 1 p β 1u p + γ 1l 1 p γ 1u p N 1 ], and all other n ij =. 1 We now use the following theorem [1] to find the necessary and sufficient conditions for stability in the absence or presence of delay. Theorem 1 : A set of necessary and sufficient conditions for E to be asymptotically stable for τ is the following: a The real parts of all the roots of λ, = are negative. b For all real k and τ, ik, τ, where i = 1. Theorem : The unique interior equilibrium E of the system 6 with toxic inhibition is locally asymptotically stable for all τ. To prove this theorem, we assume a solution of the model 6 of the form Ut = ρe λt, ρ R. This leads to the following characteristic equation : λ, τ = λ + a 1 λ + a + a 3 e λτ =, where a 1 = m 11 + m ; a = m 11 m ; a 3 = m 1 n 1. When τ =, then equation becomes λ, = λ + a 1 λ + a + a 3 =. a 1 ± a 1 4a + a 3 λ =. Since, a 1 > and a + a 3 > [from equations 17 and 18], then both the roots of λ, = are negative. Hence the interior equilibrium E is asymptotically stable for τ =. The stability criteria of the system 6 for τ = will not necessarily ensure the stability of the same system with positive delay τ >. Now we assume τ.
9 A. K. Pal et al. / Filomat 31:16 17, It is well known that the signs of the real parts of the solution characterize the stability behavior of E. Therefore, by substituting λ = ξ + iη in, we obtain real and imaginary parts, respectively as ξ η + a 1 ξ + a + a 3 cos ητe ξτ = 3 ξη + a 1 η a 3 sin ητe ξτ = 4 A necessary condition for the change of stability of the equilibrium E is that the characteristic equation should have imaginary solutions. Hence to obtain the stability criterion, we take ξ = in 3 and 4. Then we have η + a + a 3 cos ητ = a 1 η a 3 sin ητ = 6 eliminating τ by squaring and adding and 6, we get the equation for determining η as η 4 + d 1 η + d = 7 where d 1 = a 1 a, and d = a a 3 As d 1 > and a a 3 >, then we must have a + a 3 >. Thus the equation 7 has no positive root η +. Hence the characteristic equation has no purely imaginary roots ±iη + for all values of the parameter p [, 1]. Therefore the interior equilibrium E of the toxic inhibited system is always asymptotically stable for all τ and hence the delay has no effect on the system. Case II : Toxic stimulation In this case, we study the system 7, which is same as the system 6 unless in the last term, where the toxic inhibition rate γ i i = 1, are negative i.e. either species produces auxin which stimulate the growth of the other species. Here the local stability analysis of the system 7 remains same as in case I for the trivial and axial equilibrium points and so as the criteria for persistence of both the species. Now we study the local stability of the system 7 for the interior equilibrium point E N which is determined by the equation a ij N i 1, N + b ij N i + c ij =, i j, i, j = 1, 8 where a ij = α il 1 p α iu p γ jl 1 p γ ju p β ijl 1 p β iju p γ il 1 p γ iu p, b ij = K jl 1 p K ju p γ il 1 p γ iu p K il 1 p K iu p γ jl 1 p γ ju p α il 1 p α iu p α jl 1 p α ju p + β ijl 1 p β iju p β jil 1 p β jiu p, 9 c ij = K il 1 p K iu p α jl 1 p α ju p K jl 1 p K ju p β jl 1 p β ju p, for all p [, 1]. Since, c ij > [from ], then the system 7 has unique positive interior equilibrium E N 1, N exists if a 1 >, a 1 <, and b 1 < 3 according to the persistence condition of both the species 17 together with the conditions 18 and 19. To study the local stability properties of E N 1, N, we linearize the system 7, then the relation 1 becomes m 11 = N 1 [α 1l 1 p α 1u p γ 1l 1 p γ 1u p N ], m 1 =, m 1 = N [β 1l 1 p β 1u p γ l 1 p γ u p N ], m = N [α l 1 p α u p γ l 1 p γ u p N 1 ],
10 A. K. Pal et al. / Filomat 31:16 17, n 1 = N 1 [β 1l 1 p β 1u p γ 1l 1 p γ 1u p N 1 ], and all other n ij =. 31 and the characteristic equation becomes λ, τ = λ + a 1 λ + a + a 3 e λτ = 3 where a 1 = m 11 + m ; a = m 11 m ; a 3 = m 1 n 1. From equation 3, for τ =, the interior equilibrium E is locally asymptotically stable if a + a 3 > i.e. m 11 m m 1 n 1 > 33 Now to study the stability behaviour for τ, we assume, λ = iη, η > is a root of the equation 3, then we have, η + ia 1 η + a + a 3 e iητ = 34 equating real and imaginary parts of equation 34, we have, and η + a a 3 cos ητ = 3 a 1 η a 3 sin ητ = 36 Eliminating τ by squaring and adding 3 and 36, we get the following fourth degree equation for determining η as η 4 + a 1 a η + a a 3 = 37 solving the above quadratic equation in η, we get, η = 1 a 1 a ± 1 a 1 a 4a a 3 Here we assume, a a 3 < i.e. a 3 < a < a 3, so that the equation 37 has only one positive solution η. Therefore the characteristic equation 3 has a purely imaginary roots ±iη for all values of the parameter p [, 1]. Putting the value of η in the equation 36 and solving for τ, we get, τ j = 1 η arcsin a 1 η a 3 From the above discussion we have the following lemma. Lemma 1 : for τ < τ, the equation 3 has a pair of imaginary roots ±iη for all p. + jπ, j =, 1,, η Theorem 3 : Let τ j be defined by 38 and the condition 33 is satisfied, then the equilibrium point E N 1, N of the system 7 is asymptotically stable for τ < τ and unstable for τ > τ. Further as τ increases through τ, E bifurcates into small amplitude periodic solutions, where τ = τ for j = and for all p [, 1]. j Proof : For τ =, the equilibrium point E N 1, N Hence by Butler s lemma [9], E N 1, N remains stable for τ < τ that dreλ dτ is asymptotically stable if the condition 33 is satisfied. for all p [, 1]. Now we have to show τ=τ,η=η >. This indicate that there exists at least one eigenvalue with positive real part for τ > τ for all p [, 1]. Also the condition of Hopf-bifurcation are then satisfied yielding the required periodic solution for all p [, 1]. Now differentiating 3 with respect to τ we get [ λ + a 1 τa 3 e λτ] dλ dτ = λa 3 e λτ 1 dλ = λ + a 1 dτ λa τ 3 e λτ λ 39
11 A. K. Pal et al. / Filomat 31:16 17, Since, e λτ = λ +a 1 λ+a a, the equation 39 becomes dλ dτ 1 λ+a = 1 3 λλ +a 1 λ+a τ λ. Thus [ ] [ dreλ dλ sign = sign dτ Re λλ + a 1 1 dτ = sign Re λ=iη λ λ + a λ=iη 1 λ + a τ ] λ λ=iη = 1 a 1 sign η η η a η a 1 η + a η = 1 η 4 + a 1 sign η a η a 1 η + a η [ ] As a 1 a is always positive, we have, dreλ dτ >. Therefore, the transversality condition holds η=η, τ=τ and hence Hopf-bifurcation occurs at η = η, τ = τ for all p [, 1]. Hence the theorem. 6. Direction and Stability of Hopf Bifurcation In this section our attention is focussed an investigation of the direction, stability and period of the periodic solution bifurcating from a stable equilibrium E = N 1, N. Following the ideas of Hassard et. al [19], we derive the explicit formula for determining the Hopf bifurcation at the critical value of τ j by using normal form and centre manifold theory. Without loss of generality, we denote any one of the critical values τ j where j =, 1,,... by τ at which equation 37 has a pair of purely imaginary roots ±iη and the system 7 undergoes a Hopf bifurcation at E. Let u 1 t = N 1 N 1, u t = N N, τ = µ + τ and t t τ where µ R, then at µ = the system 7 undergoes Hopf bifurcation. Linearizing system 7, we have du 1 = m 11 u 1 + n 1 u t τ du = m 1 u 1 + m u In space C[ 1, ], R, the system 4 is transformed into a functional differential equation as 4 ut = L µ u t + f µ 1 u t 41 where u t = u 1 t, u t T R and L µ : C R, f : R C R are represented by m L µ φ = µ + τ φ1 11 n m m + µ + τ φ 1 φ1 1 1, 4 φ 1 f µ, φ = µ + τ a11 φ 1 + a 1φ 1 φ 1 + a φ 1 φ 1 b 11 φ 1 φ + b 1 φ + b φ 1 φ, 43 φθ = φ 1 θ, φ θ T C ; the entries a ij and b ij are given as a 11 = αˆ 1 + γˆ 1 N, a1 = βˆ 1 + γˆ 1 N 1, a = γˆ 1, b 11 = βˆ 1 + γˆ N, b1 = αˆ + γˆ N 1, a = γˆ By Riesz representation theorem, there exist a function ηθ, µ of bounded variation for θ [ 1, ] such that In fact we can choose ηθ, µ = µ + τ L µ φ = m 11 m 1 m dηθ, φθ, f orφ C 44 1 δθ µ + τ n 1 δθ + 1, 4
12 A. K. Pal et al. / Filomat 31:16 17, where δ is a dirac delta function. For φ = C [ 1, ], R, define dφθ dθ, for θ [ 1, ; Aµφ = 46 dηµ, sφs, for θ =. 1 and The system 41 is equivalent to Rµφ = {, for θ [ 1, ; f µ, φ, for θ =. 47 u t = A µ u t + R µ u t, 48 where u θ = u t+θ for θ [ 1,. For ψ = C [, 1], R, where R is the -dimensional space of row vectors defined by dψs A ds, for s, 1]; ψs = 1 dηt t, ψ t, for s =. and a bilinear inner product ψs, φθ = ψφ = θ 1 ξ= 49 ψξ θdηθφξdξ, ηθ = ηθ, Then A and A are adjoint operators. Suppose that qθ and q θ are eigen vectors of A and A corresponding to iη τ and iη τ respectively. By direct computation, we have qθ = 1, m 11 + iη τ θ eiη where D = m 11 +iη n 1 e iη τ and q s, qθ = 1, q s, qθ =. Let u t be the solution of 48 when µ =. Define n 1 e iη τ m q s = D iη τ, 1 θ eiη m m 1 iη τ m m 1 + m 11 iη e iη τ 1 zt = q, u t, wt, θ = u t θ Reztqθ, on the centre manifold C, we have wt, θ = wzt, zt, θ, where wz, z, θ = w θ z + w 11θz z + w θ z + w 3θ z z and z are local coordinates for the centre manifold C in the direction of q and q. Note that w is real if u t is real and we consider only real solutions. For solution u t C of 48 since µ =, we have żt = iη τ z + q f, wz, z, + Rezqθ de f = iη τ z + q f z, z We rewrite this equation as żt = iη τ zt + z, z, 4
13 z where z, z = + z 11z z z z +... It follows form and 3 A. K. Pal et al. / Filomat 31:16 17, u t θ = wt, θ + Re{ztqt} = w θ z + w 11θz z + w θ z where q 1 = m 11 iη m 11 e iη τ θ. It follows together with 43 that +1, q 1 T e iη τ θ z + 1, q 1 T e iη τ θ z +..., z, z = q f z, z = q f, u t = τ D[{ q a 11 + a 1 q 1 e iη τ + a q 1 e iη τ + b 11 q 1 + b 1 q 1 + b q 1 }z +{ q a 11 + a 1 Req 1 + a Req 1 + b 11 Req 1 + b 1 q 1 +b q 1 }z z +{ q a 11 + a 1 q 1 e iη τ + a q 1 e iη τ + b 11 q 1 + b 1 q 1 + b q 1 } z 1 + { q a 11 w w 1 + a 1 1 q 1w 1 τ 1 + w 1eiη + q 1 w 1 τ 1 + w e iη + q a q 1 w τ 11 1e iη + w 1 q 1e iη τ + b 11 w 1 + w b 1 1 q 1w w 1 + q 1w 1 + w b w + w 11 q 1}z z] 6 where = m iη. Comparing the coefficients of z, m z, z z and z z, we have 1 = τ D q a 11 + a 1 q e iη τ + a q 1 e iη τ + b 11 q 1 + b 1 q 1 ; 11 = τ D q a 11 + a 1 q 1 e iη τ + a q 1 e iη τ + b 11 q 1 + b 1 q 1 + b q 1 ; = τ D q a 11 + a 1 Req 1 + a Req 1 + b 11 Req 1 + b 1 q 1 +b q 1 ; 1 = τ D q a 11 w 1 + w1 + a 1 q 1 w 1 τ 1 + w 1e iη 11 +q a q 1 w τ 11 1e iη + w 1 q 1e iη τ + b 11 w 1 + w1 +b 1 q 1 w + w q 1w 1 + w b w + w 11 q 1 7 Since there are w θ and w 11 θ in q 1, we need to find out their values at θ = and θ = 1. From the definition given by equation 48 and, we have 11 ẇ = u t żq ż q = { Aw Re{ q f qθ}, for θ [ 1, ; Aw Re{ q f q} + f, for θ =. = Aw + Hz, z, θ 8 where Hz, z, θ = H θ z + H 11θz z + H θ z Substituting the above series and computing the corresponding coefficients, we have A iη τ w = H θ, Aw 11 θ = Hw 11 θ,... 6 For θ [ 1,, we know that Hz, z, θ = q f qθ q f qθ = z, z, θqθ z, z, θ qθ 61
14 Comparing the coefficients, we have A. K. Pal et al. / Filomat 31:16 17, H θ = θ qθ θ qθandh 11 θ = 11 θq θ 11 θ q θ 6 From equation 6, 6 and the definition of A, it follows that We know that qθ = 1, q 1 T e iη τ θ. Hence ẇ = iη τ w θ + q θ qθ + q θ qθ. w θ = i τ θ η τ qeiη + i 3η q τ + E 1 e iη τ θ, 63 where E 1 = E 1 1, E 1 R is a constant vector. Similarly we can obtain where E = E 1, E R is a constant vector. w 11 θ = i 11 τ θ η τ qeiη + i 11 η qe τ iη τ θ + E, 64 Now we seek the appropriate value of E 1, E. From the definition of A and 6, we obtain where ηθ = η, θ. From equation 8, we have dηθw θ = η τ w H 6 1 and dηθw 11 θ = H 11, 66 1 H = q q + τ a11 + a 1 q 1 e iη τ + a q 1 e iη τ b 11 + b 1 q 1 + b q 1 H 11 = 11 q 11 q + τ Substituting 63 and 67 into 6 and noticing that iη τ I we obtain e iη τ θ dηθ 1 iη τ I q = and a 11 + a 1 Req 1 + a Req 1 b 11 Req 1 + b 1 q 1 + b q 1 e iη τ θ dηθ 1 iη τ I e iη τ θ dηθ q =, 1 E 1 = τ B1, 69 B where B 1 = a 11 + a 1 q 1 e iη τ + a q 1 e iη τ, B 1 = b 11 + b 1 q 1 + b q 1 Hence, equation 69 reduces to iη m 11 n 1 e iη τ By using Crammer s rule, we have, E 1 1 = A m 1 iη + m B 1 n τ 1 e iη B 1 iη + m, E 1 = A E 1 = B1 B iη m 11 B 1 m 1 B,
15 where A = A. K. Pal et al. / Filomat 31:16 17, iη m n τ 11 1 e iη m iη + m 1 Similarly, substituting 64 and 68 into 66, we have m n a a 1 Req 1 + a Req 1 m m E = 1 b 11 Req 1 + b 1 q 1 + b q 1 Hence we have E 1 1 = A E 1 = A a 11 + a 1 Req 1 + a Req 1 n 1 b 11 Req 1 + b 1 q 1 + b q 1 m m 11 a 11 + a 1 Req 1 + a Req 1 m 1 b 11 Req 1 + b 1 q 1 + b q 1,, where A = A θ=, η =. Then, we can find out w and w 11 from relation 63 and 64. Furthermore, we can determine 1 by the system parameters and delay in 7. Thus we can compute the following results: C 1 = i η τ , 7a µ = Re{C 1} Re{λ τ }, β = Re{C 1 }, 7b 7c τ = Im {C 1/} + µ Im{λ τ } τ η, 7d which determine the bifurcating periodic solution in the center manifold at the critical value τ, i.e., µ determines the direction of Hopf bifurcation: if µ > µ <, then the Hopf bifurcation is supercritical subcritical and the bifurcating periodic solution exists for τ > τ τ < τ. Here β determines the bifurcating periodic solutions: the bifurcating periodic solutions are stable unstable if β < β > and τ the period of periodic solution: the period increase decrease if τ > τ < Existence of switching stability The characteristic equation corresponding to the Jacobian matrix J at the interior equilibrium E N of the system with positive delay is λ + m 11 m λ + m 11 m = m 1 n 1 e λτ This equation can be written as Pλ + Qλe λτ =, where Pλ = λ + m 11 m λ and Qλ = m 1 n 1 Clearly, Pλ and Qλ are both analytic function in Reλ > how we have the following results: 1, N 1 i P + Q = m 11 m m 1 n 1
16 ii P in = A. K. Pal et al. / Filomat 31:16 17, = N 1 N [ αˆ 1αˆ + βˆ 1 γˆ 1 N 1 + βˆ 1 γˆ N βˆ 1 βˆ 1 + αˆ 1 γˆ N 1 + αˆ γˆ 1 N ] if αˆ 1 αˆ + βˆ 1 γˆ 1 N 1 + βˆ 1 γˆ N > βˆ 1 βˆ 1 + αˆ 1 γˆ N 1 + αˆ γˆ 1 N PiN, Q in = QiN Now FN = PiN QiN = N 4 + m 11 + m N + m 11 m m 1 n 1 which is a quadratic in N. Therefore FN = has at least one positive root if m 11 m < m 1 n. Hence by applying theorem 4.1 in 1 Kuang [], we see that the system 7 possesses at most finite number of stability switches. 7. Numerical Illustration Analytical studies can never be completed without numerical verification of the results. In this section we present computer simulation of some solutions of the system 6 and 7. Beside verification of our analytical findings, these numerical solutions are very important from practical point of view. We consider all fictitious set of data through analysis to establish our theoretical findings. Let us consider a set of imprecise biological values of parameters as follows in appropriate units : Kˆ 1 = [1.9,.], Kˆ = [.9, 1.1], αˆ 1 = [.,.8], αˆ = [.6,.9], βˆ 1 = [.4,.7], βˆ 1 = [.1,.], γˆ 1 = [.8,.1], γˆ = [.,.] and p [, 1]. Case I : Toxic Inhibition Using the parametric form of interval numbers and assuming the initial condition N 1, N = 1., 1., we find the conditions given in 19 are satisfied, which imply that the unique interior equilibrium point E N 1, N exist of the system 6 for all values of p [, 1]. From Table 1 we observe that for different values of the parameter p, the system 6 corresponds a unique positive equilibrium point which are locally asymptotically stable. We present the dynamics of the model for different values of p p =.,.3,.,.8, 1. in fig.1a-e. These figures show that the interior equilibrium E exist and asymptotically stable for all values of p [, 1]. But the values are different for different values of p. Fig.1f shows that both the species population decrease with increasing p, but N 1 species is decreasing rapidly where as N species decreasing slowly. However, in this case, analytically we have already shown that the system is always stable for all τ. Table 1 : Stable equilibrium points of toxic inhibited system for different p p Equilibrium Points , , , , ,.9463
17 A. K. Pal et al. / Filomat 31:16 17, p=. 3 p=.3 N 1 3 N 1 populationn 1,N populationn 1,N 1 1 N N timet Fig.1a timet Fig.1b 3 p=. 3 p=.8 N 1 N 1 populationn 1,N populationn 1,N 1 1 N N timet Fig.1c timet Fig.1d p=1. 3 N 1 3 N 1 populationn 1,N 1 population N 1,N 1 N N timet Fig.1e p Fig.1f Fig.1. Behavior of the toxic inhibited species with time by using the imprecise parameters values and initial conditions N 1, N = 1., 1. for p =.,.3,.,.8, 1. and Fig.1f shows the dynamical behavior of the two species population N 1, N with respect to p when the values of the other parameters are same
18 A. K. Pal et al. / Filomat 31:16 17, Case II : Toxic Stimulation Here also we consider the same parametric form of interval numbers and assuming the initial condition N 1, N = 1., 1., we find the conditions given in 3 are satisfied which imply that the unique interior equilibrium point E N 1, N exist of the system for all values of p [, 1]. Also when the delay value τ =, the condition 33 is satisfied and the interior equilibrium point E is locally asymptotically stable. From Table, we observe that for different values of the parameter p, the system 7 corresponds a unique positive equilibrium point which are locally asymptotically stable. we present the dynamics of the model for different values of p p =.,.3,.,.8, 1. in fig.a-e. These figures show that the interior equilibrium E exist and asymptotically stable for all values of p [, 1]. But values are different for different values of p. Fig.f shows that both the species population decrease with increasing p, but N species is decreasing rapidly where as N 1 species decreasing slowly. Table : Stable equilibrium points of toxic stimulated system for different p p Equilibrium Points , , , , , p=. 6 p= population, 4 3 population, timet Fig.a timet Fig.b 4 p=. 3 p= population, population, timet timet
19 A. K. Pal et al. / Filomat 31:16 17, Fig.c Fig.d 3 p= population, population, timet Fig.e p Fig.f Fig.. Behavior of the toxic stimulated species with time by using the imprecise parameters values and initial conditions N 1, N = 1., 1. for p =.,.3,.,.8, 1. and Fig.f shows the dynamical behavior of the two species population N 1, N with respect to p when the values of the other parameters are same Finally we investigate the dynamical behavior of the system 7 numerically in the presence of time delay using the same data. For the given values of the interval parameters the equation 37 has unique positive solution namely η which takes different value depending on the values of the parameter p given in Table 3. Now depending on these values of η for j =, from equation 38 we obtain the initial values of the bifurcation parameter time delay τ say τ given in Table 3. For a particular value of p we get a definite critical value of the delay τ as shown in Table 3. Now for a particular value of p, if value of τ is below the critical value τ then fig.3a-e shows that the interior equilibrium point E N 1, is asymptotically stable and both the species converge to their steady states in finite time. N Now if we gradually increase only the value of delay, the stability of the equilibrium point E N 1, N may switch over. By theorem as τ passes through the critical value τ as per Table 3, E losses its stability as shown in fig.3a-e, which is the case of Hopf-bifurcation. Also from fig.4a-e, the positive equilibrium point E N 1, N is unstable with a periodic orbit near E N 1, N as τ passes through τ. Table 3 : Values of η and critical values of time delay τ of toxic stimulated system for different p p η τ p=., τ=.7<τ * 3 p=., τ=.7> τ * p=., τ=.7>τ * 4 4 population, timet timet Fig.3a timet
20 A. K. Pal et al. / Filomat 31:16 17, p=.3, τ=.9< τ * 3 p=.3, τ=1.>τ * 14 p=.3, τ=1.>τ * population, timet timet Fig.3b timet 3 p=., τ=1.< τ * p=., τ=1.>τ * 11 p=., τ=1.>τ * population, timet timet Fig.3c timet 8 p=.8, τ=1.3< τ * p=.8, τ=1.4>τ * 9 p=.8, τ=1.4>τ * population, timet timet Fig.3d timet 7 p=1., τ=1.<τ * p=1., τ=1.6>τ * 8 p=1., τ=1.6>τ * population, timet timet Fig.3e timet Fig.3. Time course of two competitive species variation when τ < τ and τ > τ with same imprecise parameter values for p =.,.3,.,.8, 1.
21 A. K. Pal et al. / Filomat 31:16 17, p=., τ=.7>τ * 14 p=.3, τ=1.>τ * Fig.4a 1 3 Fig.4b 11 p=., τ=1.>τ * 9 p=.8, τ=1.4>τ * Fig.4c 1 1 Fig.4d 8 p=1., τ=1.6>τ * Fig.4e Fig.4. Phase portrait of species when τ > τ with same imprecise parameter values for p =.,.3,.,.8, 1.
22 A. K. Pal et al. / Filomat 31:16 17, Discussion In this paper, we have studied a delayed two species competition model with the effect of toxic substances. The equations of the model are obtained from the modified Lotka-Volterra type competition of two species which are affected by toxicant. These toxic substances may be toxic inhibitory or stimulatory to the other species, which are studied extensively in our paper. Most of the competition models are generally based on the assumption that the biological parameters are precisely known. But in real life situation, it is not possible to know the values of all biological parameters precisely. In this paper, we consider a two species competition model with delay in the second species using some imprecise parameters. Here we introduce the concept of interval numbers to the model by considering the biological parameters Kˆ 1, Kˆ, αˆ 1, αˆ, βˆ 1, βˆ 1, γˆ 1, γˆ which are imprecise in nature. We have used parametric functional form of the interval number to convert the imprecise competition model to the corresponding parametric competition delay model. We consider all fictitious set of data through analysis to establish our theoretical findings. We examine the dynamical behaviour of both the cases, toxic inhibitory and stimulatory model system in absence as well as presence of time delay for different values of the parameter p [, 1]. We have discussed the existence and stability of various equilibrium points of both the systems. We get different equilibrium level of the species and the critical value of the time delay depending upon values of the parameter p. Analytically it is obtained that the time delay does not affect the stability on the toxic inhibited system. But in case of toxic stimulatory system, it is shown that the system becomes unstable at different time delays for different values of the parameter p and leads to stable limit cycle periodic solutions through Hopf-bifurcation. All our important mathematical findings for the dynamical behaviour of the toxicant affected two species competition model are also numerically verified and graphical representation of a variety of solutions of system 6 and 7 are depicted by using MATLAB with some imprecise parameter values. These numerical results are very important to understand the system in both mathematical and ecological points of view. Finally, we conclude that impreciseness of biological parameters have great impact on the behaviour of the delay model. For the first time we use the concept of interval number to present imprecise delay competition model, which makes the situation more realistic as always it is not possible to know the parameter values precisely. Here we consider all the biological parameters are imprecise, except the delay parameter τ. The delay model can be made more realistic when incorporated with impreciseness in the delay term, it makes the model more interesting and is left for future work consideration. Acknowledgements: We are grateful to the anonymous referee and Prof. Miljana Jovanovic, Editor for their careful reading, valuable comments and helpful suggestions which have helped us to improve the presentation of this work significantly. The second author is thankful to the University Grants Commission, India for providing JRFRGNF. References [1] D. M. Anderson, Toxic algae blooms and red tides : a global perspective, Environmental Science and Toxicology, Elsevier, New York, USA, [] R. C. Bassanezi, L. C. Barros, A. Tonelli, Attractors and asymptotic stability for fuzzy dynamical systems, Fuzzy Sets Syst [3] H. Berglund, Stimulation of growth of two marine algae by organic substances excreted by Enteromorpha linza in unialgal and axeniccultures, Physicol. Plant [4] G. Birkhoff, G. C. Rota, Ordinary Differential Equations, Ginn, Boston, 198. [] C. Celik, The stability and Hopf bifurcation for a predator-prey system with time delay, Chaos, Solitons and Fractals [6] Y. Chen, J. Yu, C. Sun, Stability and Hopf bifurcation analysis in a three-level food chain system with delay, Chaos, Solitons and Fractals [7] T. Das, R. N. Mukherjee, K. S. Chaudhuri, Harvesting of a prey-predator fishery in the presence of toxicity, Appl. Math. Model,
23 A. K. Pal et al. / Filomat 31:16 17, [8] J. T. De Luna, T. G. Hallam, Effect of Toxicants on population : a qualitative approach IV. Resource-consumer-toxicants models, Ecol. Model [9] H. Freedman, V. S. H. Rao, The trade-off between mutual interference and tome lags in predator-prey systems, Bull. Math. Biol [1] H. I. Freedman, J. B. Shukla, Models for the effect of the toxicant in single species and predator-prey system, J. Math. Biol [11] M. Ghosh, P. Chandra, P. Sinha, A mathematical model to study the effect of toxic chemicals on a prey-predator type fishery, J. Biol. Syst [1] K. Gopalsamy, Stability and Oscillation in Delay-Differential Equations of Population Dynamics, Kluwer, Dordrecht, 199. [13] T. C. Gard, T. G. Hallam, persistence in food webs I: Lotka-Volterra food chains, Bull. Math. Bio [14] M. Guo, X. Xu, R. Li, Impulsive functional differential inclusions and fuzzy populations models, Fuzzy Sets Syst [] J. K. Hale, Theory of Functional Differential Equations, Springer-Verlag, New York, [16] T. G. Hallam, C. E. Clark, G. S. Jordan, Effects of toxicants on population : a qualitative approach II. First-order kinetics, J. Math. Biol [17] T. G. Hallam, C. E. Clark, R. R. Lassiter, Effects of toxicants on population :a qualitative approach I. Equilibrium environmental exposure, Ecol. Model [18] T. G. Hallam, J. T. De Luna, Effects of toxicants on population :a qualitative approach III. Environmental and food chain pathways, J. Theor. Biol [19] B. Hassard, D. Kazarinof, Y. Wan, Theory and application of Hopf bifurcation, Cambridge University Press, Cambridge [] J. He, K. Wang, The survival analysis for a single-species population model in a polluted environment, Appl. Math. Model [1] A. L. Jensen, J. S. Marshall, Application of a surplus production model to assess environmental impacts on exploited populations of Daphina pluex in the laboratory, Environ. Pollut. Series A [] Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Academic Press, New York, [3] W. M. Liu, Criteria of Hopf bifurcations without using eigen values, J. Math. Anal. Appl [4] A. Maiti, A. K. Pal, G. P. Samanta, Effect of time-delay on a food chain model, Appl. Math. Comp [] A. Maiti, A. K. Pal, G. P. Samanta, Usefulness of biocontrol of pests in tea : a mathematical model, Math. Model. Nat. Phenon [6] M. MacDonald, Biological delay system: linear stability theory, Cambridge University Press, Cambridge, [7] J. Maynard-Smith, Models in Ecology, Cambridge University Press, Cambridge, 197. [8] S. A. Nelson, The problem of oil pollution in the sea: F. S. Russell, M. Yonge Eds., Adv. in Marine Biol., Academic Press, London, [9] D. Pal, G. S. Mahapatra, G. P. Samanta, A proportional harvesting dynamical model with fuzzy intrinsic growth rate and harvesting quantity, Pac. Asian J. Math [3] D. Pal, G. S. Mahapatra, G. P. Samanta, Optimal harvesting of prey-predator system with interval biological parameters: A bioeconomic model, Math. Biosc [31] D. Pal, G. S. Mahapatra, G. P. Samanta, Bifurcation analysis of predator-prey model with time delay and harvesting efforts using interval parameter, Int. J. Dynam. Control 14. [3] D. Pal, G. S. Mahapatra, G. P. Samanta, Stability and bionomic analysis of fuzzy parameter based preypredator harvesting model using UFM, Nonlinear Dynamics [33] Sweta Pathak, A. Maiti, G. P. Samanta, Rich dynamics of a food chain model with HassellVarley type functional responses, App. Maths. and Comp [34] B. Patra, A. Maiti, G. P. Samanta, Effect of Time-Delay on a Ratio-Dependent Food Chain Model, Nonlinear Ana.: Modelling and Control [3] M. Peixoto, L. C. Barros, R. C. Bassanezi, Preor-prey fuzzy model, Ecol. Model [36] R. Pratt, Influence of the size of the inoculum on the growth of Chlorella vulgaris in freshly prepared culture medium, Am. J. Bot [37] R. Pratt, J. Fong, Studies on Chlorella vulgaris, II. Further evidence that chlorella cells form a growth-inhibiting substance, Am. J. Bot [38] E. L. Rice, Allelopathy, second ed., Academic Press, New York [39] T. R. Rice, Biotic influences affecting population growth of planktonic algae, US Fish Wildl. Serv., Fish Bull [4] G. P. Samanta, A two-species competitive system under the influence of toxic substances, Appl. Math. Comp [41] Swarnali Sharma, G. P. Samanta, Optimal harvesting of a two species competition model with imprecise biological parameters, Nonlinear Dynamics [4] Y. Song, M. Han, Y. Peng, Stability and Hopf bifurcations in a competitive Lotka-Volterra system with two delays, Chaos, Solitons and Fractals [43] W. Wang, Z. Ma, Permanance of population in a polluted enverironment, Math. Biosci [44] R. Xua, Q. Gan, Z. Ma, Stability and bifurcation analysis on a ratio-dependent predator-prey model with time delay, J. Comp. Appl. Maths
Direction and Stability of Hopf Bifurcation in a Delayed Model with Heterogeneous Fundamentalists
International Journal of Mathematical Analysis Vol 9, 2015, no 38, 1869-1875 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ijma201554135 Direction and Stability of Hopf Bifurcation in a Delayed Model
More informationHopf Bifurcation Analysis of a Dynamical Heart Model with Time Delay
Applied Mathematical Sciences, Vol 11, 2017, no 22, 1089-1095 HIKARI Ltd, wwwm-hikaricom https://doiorg/1012988/ams20177271 Hopf Bifurcation Analysis of a Dynamical Heart Model with Time Delay Luca Guerrini
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 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 informationHOPF BIFURCATION CONTROL WITH PD CONTROLLER
HOPF BIFURCATION CONTROL WITH PD CONTROLLER M. DARVISHI AND H.M. MOHAMMADINEJAD DEPARTMENT OF MATHEMATICS, FACULTY OF MATHEMATICS AND STATISTICS, UNIVERSITY OF BIRJAND, BIRJAND, IRAN E-MAILS: M.DARVISHI@BIRJAND.AC.IR,
More informationHOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY
Vol. 37 17 No. J. of Math. PRC HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY LI Zhen-wei, LI Bi-wen, LIU Wei, WANG Gan School of Mathematics and Statistics,
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 informationJournal of Mathematical Analysis and Applications
J Math Anal Appl 367 200) 249 259 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications wwwelseviercom/locate/jmaa Existence, uniqueness and stability analysis of
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 informationHARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL
ANZIAM J. 452004), 443 456 HARVESTING IN A TWO-PREY ONE-PREDATOR FISHERY: A BIOECONOMIC MODEL T. K. KAR 1 and K. S. CHAUDHURI 2 Received 22 June, 2001; revised 20 September, 2002) Abstract A multispecies
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 informationROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 6, Number 1, Winter 1997 ROLE OF TIME-DELAY IN AN ECOTOXICOLOGICAL PROBLEM J. CHATTOPADHYAY, E. BERETTA AND F. SOLIMANO ABSTRACT. The present paper deals with
More informationModeling and Simulation Study of Mutuality Interactions with Type II functional Response and Harvesting
American Journal of Applied Mathematics 201; 6(3): 109-116 http://www.sciencepublishinggroup.com/j/ajam doi: 10.116/j.ajam.2010603.12 ISSN: 2330-003 (Print); ISSN: 2330-006X (Online) Modeling and Simulation
More informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
More informationStability and Hopf Bifurcation Analysis of the Delay Logistic Equation
1 Stability and Hopf Bifurcation Analysis of the Delay Logistic Equation Milind M Rao # and Preetish K L * # Department of Electrical Engineering, IIT Madras, Chennai-600036, India * Department of Mechanical
More informationStability Analysis And Maximum Profit Of Logistic Population Model With Time Delay And Constant Effort Of Harvesting
Jurnal Vol 3, Matematika, No, 9-8, Juli Statistika 006 & Komputasi Vol No Juli 006 9-8, Juli 006 9 Stability Analysis And Maximum Profit Of Logistic Population Model With Time Delay And Constant Effort
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 information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More 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 informationDynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response
Liu Advances in Difference Equations 014, 014:314 R E S E A R C H Open Access Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response
More informationKey words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.
ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart
More informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationBehaviour of simple population models under ecological processes
J. Biosci., Vol. 19, Number 2, June 1994, pp 247 254. Printed in India. Behaviour of simple population models under ecological processes SOMDATTA SINHA* and S PARTHASARATHY Centre for Cellular and Molecular
More informationSTABILITY AND HOPF BIFURCATION OF A MODIFIED DELAY PREDATOR-PREY MODEL WITH STAGE STRUCTURE
Journal of Applied Analysis and Computation Volume 8, Number, April 018, 573 597 Website:http://jaac-online.com/ DOI:10.11948/018.573 STABILITY AND HOPF BIFURCATION OF A MODIFIED DELAY PREDATOR-PREY MODEL
More informationStability and nonlinear dynamics in a Solow model with pollution
Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 4, 565 577 565 http://dx.doi.org/10.15388/na.2014.4.3 Stability and nonlinear dynamics in a Solow model with pollution Massimiliano Ferrara
More informationLOTKA-VOLTERRA SYSTEMS WITH DELAY
870 1994 133-140 133 LOTKA-VOLTERRA SYSTEMS WITH DELAY Zhengyi LU and Yasuhiro TAKEUCHI Department of Applied Mathematics, Faculty of Engineering, Shizuoka University, Hamamatsu 432, JAPAN ABSTRACT Sufftcient
More informationChaos and adaptive control in two prey, one predator system with nonlinear feedback
Chaos and adaptive control in two prey, one predator system with nonlinear feedback Awad El-Gohary, a, and A.S. Al-Ruzaiza a a Department of Statistics and O.R., College of Science, King Saud University,
More informationGlobal Qualitative Analysis for a Ratio-Dependent Predator Prey Model with Delay 1
Journal of Mathematical Analysis and Applications 266, 401 419 (2002 doi:10.1006/jmaa.2001.7751, available online at http://www.idealibrary.com on Global Qualitative Analysis for a Ratio-Dependent Predator
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 informationPersistence of a generalized prey-predator model with prey reserved and herd behaviour of prey in
International Journal of Engineering Science Invention (IJESI) ISSN (Online): 239 6734, ISSN (Print): 239 6726 Volume 7 Issue 3 Ver. I March 208 PP. 40-48 Persistence of a Generalized Prey-Predator Model
More informationStability analysis of a prey-predator model with a reserved area
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 4, 5(3):93-3 ISSN: 976-86 CODEN (USA): AASRFC Stability analysis of a prey-predator model with a reserved area Neelima
More informationDYNAMICS OF A RATIO-DEPENDENT PREY-PREDATOR SYSTEM WITH SELECTIVE HARVESTING OF PREDATOR SPECIES. 1. Introduction
J. Appl. Math. & Computing Vol. 232007), No. 1-2, pp. 385-395 Website: http://jamc.net DYNAMICS OF A RATIO-DEPENDENT PREY-PREDATOR SYSTEM WITH SELECTIVE HARVESTING OF PREDATOR SPECIES TAPAN KUMAR KAR Abstract.
More informationBIFURCATION ANALYSIS ON A DELAYED SIS EPIDEMIC MODEL WITH STAGE STRUCTURE
Electronic Journal of Differential Equations, Vol. 27(27), No. 77, pp. 1 17. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) BIFURCATION
More informationOn a non-autonomous stochastic Lotka-Volterra competitive system
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 7), 399 38 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On a non-autonomous stochastic Lotka-Volterra
More informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo
More informationStability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks
International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in
More informationThe dynamics of disease transmission in a Prey Predator System with harvesting of prey
ISSN: 78 Volume, Issue, April The dynamics of disease transmission in a Prey Predator System with harvesting of prey, Kul Bhushan Agnihotri* Department of Applied Sciences and Humanties Shaheed Bhagat
More informationPrey Predator Model with Asymptotic. Non-Homogeneous Predation
Int. J. Contemp. Math. Sciences, Vol. 7, 0, no. 4, 975-987 Prey Predator Model with Asymptotic Non-Homogeneous Predation Harsha Mehta *, Neetu Trivedi **, Bijendra Singh * and B. K. Joshi *** *School of
More informationA Discrete Model of Three Species Prey- Predator System
ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj
More informationStability Analysis of Plankton Ecosystem Model. Affected by Oxygen Deficit
Applied Mathematical Sciences Vol 9 2015 no 81 4043-4052 HIKARI Ltd wwwm-hikaricom http://dxdoiorg/1012988/ams201553255 Stability Analysis of Plankton Ecosystem Model Affected by Oxygen Deficit Yuriska
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 information1 Introduction. ISSN , England, UK World Journal of Modelling and Simulation Vol. 11 (2015) No. 3, pp
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 11 (2015) No. 3, pp. 174-198 Stability, bifurcations and chaotic dynamics in a delayed hybrid tri-trophic food chain model with
More informationResearch Article The Stability of Gauss Model Having One-Prey and Two-Predators
Abstract and Applied Analysis Volume 2012, Article ID 219640, 9 pages doi:10.1155/2012/219640 Research Article The Stability of Gauss Model Having One-Prey and Two-Predators A. Farajzadeh, 1 M. H. Rahmani
More informationHopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationLOCAL AND GLOBAL STABILITY OF IMPULSIVE PEST MANAGEMENT MODEL WITH BIOLOGICAL HYBRID CONTROL
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 11 No. II (August, 2017), pp. 129-141 LOCAL AND GLOBAL STABILITY OF IMPULSIVE PEST MANAGEMENT MODEL WITH BIOLOGICAL HYBRID CONTROL
More informationThe effect of emigration and immigration on the dynamics of a discrete-generation population
J. Biosci., Vol. 20. Number 3, June 1995, pp 397 407. Printed in India. The effect of emigration and immigration on the dynamics of a discrete-generation population G D RUXTON Biomathematics and Statistics
More informationBifurcation Analysis of a Population Dynamics in a Critical State 1
Bifurcation Analysis of a Population Dynamics in a Critical State 1 Yong-Hui Xia 1, Valery G. Romanovski,3 1. Department of Mathematics, Zhejiang Normal University, Jinhua, China. Center for Applied Mathematics
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 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 informationANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS
Volume 1, No. 8, August 2013 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS Tarini
More informationA Comparison of Two Predator-Prey Models with Holling s Type I Functional Response
A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response ** Joint work with Mark Kot at the University of Washington ** Mathematical Biosciences 1 (8) 161-179 Presented by Gunog
More informationStability and Hopf bifurcation of a ratio-dependent predator prey model with time delay and stage structure
Electronic Journal of Qualitative Theory of Differential Equations 6, No 99, 3; doi: 43/ejqtde699 http://wwwmathu-szegedhu/ejqtde/ Stability and Hopf bifurcation of a ratio-dependent predator prey model
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 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 informationStability Analysis of a Population Dynamics Model with Allee Effect
Stability Analysis of a Population Dynamics Model with Allee Effect Canan Celik Abstract In this study, we focus on the stability analysis of equilibrium points of population dynamics with delay when the
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 informationPERSISTENCE AND EXTINCTION OF SINGLE POPULATION IN A POLLUTED ENVIRONMENT
Electronic Journal o Dierential Equations, Vol 2004(2004), No 108, pp 1 5 ISSN: 1072-6691 URL: http://ejdemathtxstateedu or http://ejdemathuntedu tp ejdemathtxstateedu (login: tp) PERSISTENCE AND EXTINCTION
More informationSPIRAL WAVE GENERATION IN A DIFFUSIVE PREDATOR-PREY MODEL WITH TWO TIME DELAYS
Bull. Korean Math. Soc. 52 (2015), No. 4, pp. 1113 1122 http://dx.doi.org/10.4134/bkms.2015.52.4.1113 SPIRAL WAVE GENERATION IN A DIFFUSIVE PREDATOR-PREY MODEL WITH TWO TIME DELAYS Wenzhen Gan and Peng
More informationEXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY. S. H. Saker
Nonlinear Funct. Anal. & Appl. Vol. 10 No. 005 pp. 311 34 EXISTENCE OF POSITIVE PERIODIC SOLUTIONS OF DISCRETE MODEL FOR THE INTERACTION OF DEMAND AND SUPPLY S. H. Saker Abstract. In this paper we derive
More informationPersistence and global stability in discrete models of Lotka Volterra type
J. Math. Anal. Appl. 330 2007 24 33 www.elsevier.com/locate/jmaa Persistence global stability in discrete models of Lotka Volterra type Yoshiaki Muroya 1 Department of Mathematical Sciences, Waseda University,
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 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 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 informationDelayed Dynamics in Heterogeneous Competition with Product Differentiation
Delayed Dynamics in Heterogeneous Competition with Product Differentiation Akio Matsumoto Department of Economics Chuo University Hachioji, Tokyo, 19-0393, Japan akiom@tamacc.chuo-u.ac.jp Ferenc Szidarovszky
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 Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses
Journal of Physics: Conference Series PAPER OPEN ACCESS Dynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses To cite this article: D Savitri 2018
More informationDelayed Dynamics in Heterogeneous Competition with Product Differentiation
Delayed Dynamics in Heterogeneous Competition with Product Differentiation Akio Matsumoto Department of Economics Chuo University Hachioji, Tokyo, 19-0393, Japan akiom@tamacc.chuo-u.ac.jp Ferenc Szidarovszky
More informationAnalysis of a Prey-Predator System with Modified Transmission Function
Research Paper American Journal of Engineering Research (AJER) e-issn : 2320-0847 p-issn : 2320-0936 Volume-3, Issue-9, pp-194-202 www.ajer.org Open Access Analysis of a Prey-Predator System with Modified
More informationOptimal harvesting policy of a stochastic delay predator-prey model with Lévy jumps
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 1 (217), 4222 423 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa Optimal harvesting policy of
More informationDynamical Analysis of a Harvested Predator-prey. Model with Ratio-dependent Response Function. and Prey Refuge
Applied Mathematical Sciences, Vol. 8, 214, no. 11, 527-537 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/12988/ams.214.4275 Dynamical Analysis of a Harvested Predator-prey Model with Ratio-dependent
More 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 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 informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationDelay SIR Model with Nonlinear Incident Rate and Varying Total Population
Delay SIR Model with Nonlinear Incident Rate Varying Total Population Rujira Ouncharoen, Salinthip Daengkongkho, Thongchai Dumrongpokaphan, Yongwimon Lenbury Abstract Recently, models describing the behavior
More informationMathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients
Nonlinear Analysis: Modelling and Control, 2007, Vol. 12, No. 4, 511 524 Mathematical Modeling and Analysis of Eutrophication of Water Bodies Caused by Nutrients A. K. Misra Department of Mathematics,
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 informationHopf Bifurcation Analysis of Pathogen-Immune Interaction Dynamics With Delay Kernel
Math. Model. Nat. Phenom. Vol. 2, No. 1, 27, pp. 44-61 Hopf Bifurcation Analysis of Pathogen-Immune Interaction Dynamics With Delay Kernel M. Neamţu a1, L. Buliga b, F.R. Horhat c, D. Opriş b a Department
More informationHopf bifurcation analysis for a model of plant virus propagation with two delays
Li et al. Advances in Difference Equations 08 08:59 https://doi.org/0.86/s366-08-74-8 R E S E A R C H Open Access Hopf bifurcation analysis for a model of plant virus propagation with two delays Qinglian
More informationPULSE-SEASONAL HARVESTING VIA NONLINEAR DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS IN FISHERY MANAGEMENT. Lev V. Idels
PULSE-SEASONAL HARVESTING VIA NONLINEAR DELAY DIFFERENTIAL EQUATIONS AND APPLICATIONS IN FISHERY MANAGEMENT Lev V. Idels University-College Professor Mathematics Department Malaspina University-College
More informationPERMANENCE IN LOGISTIC AND LOTKA-VOLTERRA SYSTEMS WITH DISPERSAL AND TIME DELAYS
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 60, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PERMANENCE
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 informationBIFURCATION ANALYSIS OF A MODEL OF THREE-LEVEL FOOD CHAIN IN A MANGROVE ECOSYSTEM
U.P.B. Sci. Bull., Series A, Vol. 78, Iss. 4, 2016 ISSN 1223-7027 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
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 informationA mathematical model to study the simultaneous effect of pollutants emitted from wood and non-wood based industries on forest resources
ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 12 2016 No. 1, pp. 34-47 A mathematical model to study the simultaneous effect of pollutants emitted from wood and non-wood based
More informationDynamical Systems. Pierre N.V. Tu. An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition.
Pierre N.V. Tu Dynamical Systems An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition With 105 Figures Springer-Verlag Berlin Heidelberg New York London Paris
More informationChaos Control and Bifurcation Behavior for a Sprott E System with Distributed Delay Feedback
International Journal of Automation and Computing 12(2, April 215, 182-191 DOI: 1.17/s11633-14-852-z Chaos Control and Bifurcation Behavior for a Sprott E System with Distributed Delay Feedback Chang-Jin
More informationIMPACT OF HARVESTING, NOISE AND DIFFUSION ON THE DYNAMICS OF A FOOD CHAIN MODEL WITH RATIO-DEPENDENT FUNCTIONAL RESPONSE III
ITALIAN JOURNAL OF PURE AND APPLIED MATHEMATICS N 39 8 69 89 69 IMPACT OF HARVESTING NOISE AND DIFFUSION ON THE DYNAMICS OF A FOOD CHAIN MODEL WITH RATIO-DEPENDENT FUNCTIONAL RESPONSE III G Basava Kumar
More informationBifurcation Analysis, Chaos and Control in the Burgers Mapping
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.4007 No.3,pp.171-185 Bifurcation Analysis, Chaos and Control in the Burgers Mapping E. M. ELabbasy, H. N. Agiza, H.
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 informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
More informationBifurcations in Delayed Lotka-Volterra Intraguild Predation Model
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 5-696 Vol. 37 Nos. - (4) pp. - Bifurcations in Delayed Lotka-Volterra Intraguild Predation Model Juancho A. Collera Department
More informationChaos suppression of uncertain gyros in a given finite time
Chin. Phys. B Vol. 1, No. 11 1 1155 Chaos suppression of uncertain gyros in a given finite time Mohammad Pourmahmood Aghababa a and Hasan Pourmahmood Aghababa bc a Electrical Engineering Department, Urmia
More informationDYNAMICS IN A COMPETITIVE LOTKA-VOLTERRA PREDATOR-PREY MODEL WITH MULTIPLE DELAYS. Changjin Xu 1. Yusen Wu. 1. Introduction
italian journal of pure and applied mathematics n. 36 2016 (231 244) 231 DYNAMICS IN A COMPETITIVE LOTKA-VOLTERRA PREDATOR-PREY MODEL WITH MULTIPLE DELAYS Changjin Xu 1 Guizhou Key Laboratory of Economics
More informationVariational iteration method for solving multispecies Lotka Volterra equations
Computers and Mathematics with Applications 54 27 93 99 www.elsevier.com/locate/camwa Variational iteration method for solving multispecies Lotka Volterra equations B. Batiha, M.S.M. Noorani, I. Hashim
More informationStability and optimal harvesting in a stage structure predator-prey switching model
Stability and optimal harvesting in a stage structure predator-prey switching model Q. J. A. Khan and Lakdere Benkherouf Department of Mathematics and Statistics, College of Science, Sultan Qaboos University,
More informationStability analysis of a continuous model of mutualism with delay dynamics
Stability analysis of a continuous model of mutualism with delay dynamics Item Type Article Authors Roberts, Jason A.; Joharjee, Najwa G. Citation Roberts, J. A., & Joharjee, N. G. (2016). Stability analysis
More informationTime-delay feedback control in a delayed dynamical chaos system and its applications
Time-delay feedback control in a delayed dynamical chaos system and its applications Ye Zhi-Yong( ), Yang Guang( ), and Deng Cun-Bing( ) School of Mathematics and Physics, Chongqing University of Technology,
More informationNumerical Solution of a Fractional-Order Predator-Prey Model with Prey Refuge and Additional Food for Predator
66 Numerical Solution of a Fractional-Order Predator-Prey Model with Prey Refuge Additional Food for Predator Rio Satriyantara, Agus Suryanto *, Noor Hidayat Department of Mathematics, Faculty of Mathematics
More information