Research article

Global stability of a tridiagonal competition model with seasonal succession

  • Received: 03 December 2022 Revised: 05 January 2023 Accepted: 16 January 2023 Published: 18 January 2023
  • In this paper, we investigate a tridiagonal three-species competition model with seasonal succession. The Floquet multipliers of all nonnegative periodic solutions of such a time-periodic system are estimated via the stability analysis of equilibria. Together with the Brouwer degree theory, sufficient conditions for existence and uniqueness of the positive periodic solution are given. We further obtain the global dynamics of coexistence and extinction for three competing species in this periodically forced environment. Finally, some numerical examples are presented to illustrate the effectiveness of our theoretical results.

    Citation: Xizhuang Xie, Meixiang Chen. Global stability of a tridiagonal competition model with seasonal succession[J]. Mathematical Biosciences and Engineering, 2023, 20(4): 6062-6083. doi: 10.3934/mbe.2023262

    Related Papers:

    [1] Ranjit Kumar Upadhyay, Swati Mishra, Yueping Dong, Yasuhiro Takeuchi . Exploring the dynamics of a tritrophic food chain model with multiple gestation periods. Mathematical Biosciences and Engineering, 2019, 16(5): 4660-4691. doi: 10.3934/mbe.2019234
    [2] Luca Galbusera, Sara Pasquali, Gianni Gilioli . Stability and optimal control for some classes of tritrophic systems. Mathematical Biosciences and Engineering, 2014, 11(2): 257-283. doi: 10.3934/mbe.2014.11.257
    [3] Bruno Buonomo , Francesco Giannino , Stéphanie Saussure , Ezio Venturino . Effects of limited volatiles release by plants in tritrophic interactions. Mathematical Biosciences and Engineering, 2019, 16(5): 3331-3344. doi: 10.3934/mbe.2019166
    [4] Manoj K. Singh, Brajesh K. Singh, Poonam, Carlo Cattani . Under nonlinear prey-harvesting, effect of strong Allee effect on the dynamics of a modified Leslie-Gower predator-prey model. Mathematical Biosciences and Engineering, 2023, 20(6): 9625-9644. doi: 10.3934/mbe.2023422
    [5] Hongqiuxue Wu, Zhong Li, Mengxin He . Dynamic analysis of a Leslie-Gower predator-prey model with the fear effect and nonlinear harvesting. Mathematical Biosciences and Engineering, 2023, 20(10): 18592-18629. doi: 10.3934/mbe.2023825
    [6] Saheb Pal, Nikhil Pal, Sudip Samanta, Joydev Chattopadhyay . Fear effect in prey and hunting cooperation among predators in a Leslie-Gower model. Mathematical Biosciences and Engineering, 2019, 16(5): 5146-5179. doi: 10.3934/mbe.2019258
    [7] Shuangte Wang, Hengguo Yu . Stability and bifurcation analysis of the Bazykin's predator-prey ecosystem with Holling type Ⅱ functional response. Mathematical Biosciences and Engineering, 2021, 18(6): 7877-7918. doi: 10.3934/mbe.2021391
    [8] Christian Cortés García, Jasmidt Vera Cuenca . Impact of alternative food on predator diet in a Leslie-Gower model with prey refuge and Holling Ⅱ functional response. Mathematical Biosciences and Engineering, 2023, 20(8): 13681-13703. doi: 10.3934/mbe.2023610
    [9] Eduardo González-Olivares, Claudio Arancibia-Ibarra, Alejandro Rojas-Palma, Betsabé González-Yañez . Bifurcations and multistability on the May-Holling-Tanner predation model considering alternative food for the predators. Mathematical Biosciences and Engineering, 2019, 16(5): 4274-4298. doi: 10.3934/mbe.2019213
    [10] Mengyun Xing, Mengxin He, Zhong Li . Dynamics of a modified Leslie-Gower predator-prey model with double Allee effects. Mathematical Biosciences and Engineering, 2024, 21(1): 792-831. doi: 10.3934/mbe.2024034
  • In this paper, we investigate a tridiagonal three-species competition model with seasonal succession. The Floquet multipliers of all nonnegative periodic solutions of such a time-periodic system are estimated via the stability analysis of equilibria. Together with the Brouwer degree theory, sufficient conditions for existence and uniqueness of the positive periodic solution are given. We further obtain the global dynamics of coexistence and extinction for three competing species in this periodically forced environment. Finally, some numerical examples are presented to illustrate the effectiveness of our theoretical results.



    Heavy metals are a major constituent in wastewater being released from glass, mining, metal plating, battery, and paints industries. Zn, Ni, Fe, Cu, Cr, and Cd are some of the most common heavy metals which contaminate effluent streams and are not biodegradable. These metals are contaminate the receiving water bodies which then lead to them being bio accumulated in the living beings including marine life, human being and flora and fauna around which then leads to health issues like cancer, kidney problems, ulcers, renal failure as well as stomach damage in rodents [1]. Problems arising from heavy metal pollution are extremely serious and require immediate action [2]. Efforts to remove metal ions have increased over time and newer and cheaper methods are being looked at. Recent efforts have focused on using cost-effective and easily available polymers to remove metal ions. Poly(methyl methacrylates) have some proper- ties such as chemical stability and mechanical strength, and they have been used for the development of application in separation techniques. Polymethyl acrylate with hydrazine chains contains carbonyl groups to pro- vide interaction between metals and polymers. In this study the complex of hydrazine and PMA were prepared, and used in batch equilibration technique for removal of Cu(Ⅱ) ions from aqueous solutions.

    This study focuses on removal of Cu2+ ions using poly methyl acrylate (PMA) as an adsorbent. Parameters such as temperature, initial Cu2+ ion concentration, stirring speed, and adsorbent dose were investigated. Langmuir and Freundlich models were used to analyze the isotherm data. The pseudo-first-order and pseudo-second-order models were employed to examine the kinetic data. Scanning Electron Microscopy (SEM) was used to examine the pore structure in the synthesized and activated PMAs, and the I.R test was employed to examine the functional groups.

    Emulsion polymerization is the most convenient method for preparation as bulk polymerization is impractical as there is high rate and heat evolution during acrylate polymerization. Water acts as both a moderator as well as a solvent for the catalyst system. Bromate ions react with hydrogen sulfite to produce hydrogen bisulfite radicals (HSO3) and hydroxyl radicals (OH) as catalysts. Poly(methyl acrylate) is an addition polymer. The reaction involved is:

    Methyl acrylate is polymerized in an aqueous emulsion with a free-radical catalyst being present. The polymer then coagulates as a white product in a concentrated solution of sodium chloride. Acetone solution is used to produce a film.

    A solution of poly-methyl acrylate (2 gm) in di-methyl sulfoxide (60 mL) was treated with 99% hydrazine hydrate (20 mL) and boiled under reflux with constant stirring on a hot plate for 20 h. The polymer was added with constant stirring to a large excess volume of methanol (400 mL). The precipitated poly methyl acrylic acid hydrazide polymer was filtered and dried in a vacuum pump. The obtained polymer was purified by dissolution in lest amount of water, then reprecipitated by absolute methanol. The pure polymer was then filtered dried in a vacuum desiccator for ten days over concentrated sulfuric acid.

    JSM-5600 model of SEM was employed to analyze the morphological features of the synthesized and activated PMA.

    10 mL (0.05 M, 0.1 M, 0.25 M, 0.5 M) of CuSO4 was measured and added to 100 mL dry and clean bottles which was followed by addition of the modified PMA (0.05, 0.1, 0.2 g). The conical flask was then placed on a rotary shaker and shaken at various speeds (50, 100, 150, 200 rpm) at different temperatures (25, 35, 45, 50 ℃) for 60 min to attain equilibrium. Whatman filter paper was employed to filter the suspension, followed by the determination of the metal ion concentration in the filtrate. Eq 1 was used to calculate the quantity of ions adsorbed:

    qt=coctw×v (1)

    The adsorption capacity was determined using the following mass balance equation (Eq 2) [3].

    qe=cocew×v (2)

    where Co and Ct refer to the initial and final concentrations (in mg/L) of the metal ions in wastewater before and after adsorption, for time t, respectively. Ce is the concentration of the cupric ion in CuSO4 (mg/L) when equilibrium was attained. V is the volume of CuSO4 (mL). And W is the mass (g) of the adsorbent used. The below Eq 3 was used to determine the percentage of metal ions removed:

    %R=CoCtCo×100 (3)

    where (%R) denotes the ratio between the difference in metal concentrations before and after adsorption [4].

    Various concentrations of CuSO4 (0.05, 0.25, 0.5, 0.1 M) solution was mixed with 0.1 g of the modified PMA in conical flasks which were then sealed and agitated using a shaker for 1 h at 150 rpm, 25 ℃ and a pH of 4.6. The contents were then filtered out and analyzed.

    Different adsorbent doses (0.05–0.1–0.2 g) were added to 10 mL CuSO4 solution in a conical flask which was then sealed and agitated in a shaker for 1 h at an rpm of 150 at 25 ℃. The flask contents were filtered and analyzed [5].

    Different shaking speeds were considered to analyze their effects when 1 g of the adsorbent was added to 10 mL CuSO4 and agitated at different speeds (50, 100, 150, 200 rpm) for 1 h at 25 ℃. The contents were filtered and analyzed [6].

    This was studied by agitating 0.1 g of adsorbent in 10 mL CuSO4 solution for 1 h at 150 rpm and different temperatures (25, 35, 45, 50 ℃). Contents were filtered and analyzed.

    The usual method of carrying out isotherm studies for solid–liquid interfaces involves varying the CuSO4 solution concentration while keeping the other parameters like volume, size, pH of solution, rpm and temperature constant. Performance of this adsorbent was determined by the uptake. Langmuir, and Freundlich isotherm models were used to carry out these adsorption isotherm studies. The correlation coefficient, R2 value, was used to investigate and compare the applicability of the above mentioned isotherm models.

    Langmuir isotherm: Eq 4 has been used to compute the Langmuir isotherm [7,8]:

    ceQe=1Q0b+1Q0Ce (4)

    Here, Q0 is the maximum quantity of Cu2+ ions that can be taken up as “per unit mass polymer” (mg/g) The Langmuir constant, b, is associated with energy of sorption and expressed as b.

    Freundlich isotherm: For non-ideal adsorption with heterogeneous desorption Freundlich isotherm is often used. Eq 5 expresses the formula [7,9]:

    qe=KFC1/ne (5)

    The above equation can be modified to a linear form by using logarithm on both sides as:

    logqe=logKF+1nlogCe (6)

    KF, Freundlich constant, and 1/n both describe adsorption capacity and adsorption intensity. Plotting log of qe vs logCe from Eq 6 will give the values of KF and n from the slope and intercept of the graph. The feasibility of the adsorption can be judged from the magnitude of 1/n exponent. When the magnitude is less than 1, adsorption is inadequate; 1-2 means that there is moderate adsorption and 2-10 stands for god adsorption [10].

    Kinetic Analysis: The adsorption kinetics was studied by analyzing the Cu(Ⅱ) uptake from CuSO4 at different intervals of time during adsorption. For Cu(Ⅱ) adsorption by activated PMA kinetic study, pseudo-first and pseudo-second order reaction models were applied. The appropriate model for the adsorption kinetics can be found out when linearity of both models is checked via graphs. Eqs 7 and 8 are used to describe pseudo-first-order systems:

    dqtdt=k(qeqt) (7)

    qe and qt respectively (mg/g) describe the sorption capacities at equilibrium and at time t. k describes the pseudo-first order rate constant (L min-1). The equation in integrated with the application of the following boundary conditions, from t = 0 to t = t and qt = 0 to qt = qt, to yield [11]:

    log(qeqt)=logqek2.303t (8)

    The following equation is used for pseudo-second order model Eq 9 [12]:

    tqt=1k21q2e+1qet (9)

    qe and qt denote adsorption capacities during equilibrium and at a fixed time t (mg/g) and k (g/mg min) describes the rate constant in a pseudo-second order reaction.

    Scanning Electron Microscopy (SEM) images for a selected sample of poly-methyl acrylate and its hydrazide derivative were taken in the faculty of science, Alexandria University. The purpose of this study is to understand the microstructure relationship of the sample matrix. SEM serves as a qualitative or semi-quantitative technique which is used for studying morphology and observed the pore structure of PMA. The pores present in PMA act as active sites and the adsorption takes place in those sites. Figure 1a presents images of PMA in which active sites are not visible. The modified PMA was produced by dissolving it in dimethyl sulfoxide (DMSO) and became chemically active. Figure 1b depicts the image of chemically active hydrazide derivative of PMA where active sites are clearly visible (see Figure 1).

    Figure 1.  (a) SEM images for poly-methyl acrylate, and (b) hydrazide derivative of poly-methyl acrylate.

    The IR spectra of the poly-methyl acrylate and its hydrazide derivative, which were collected and analyzed by Perkin-Elmer BX-II FT-IR spectrometer are shown in Figure 2 (see Figure 2).

    Figure 2.  Infrared spectra for poly-methyl acrylate and hydrazine derivatives of poly-methyl acrylate.

    The spectra of poly-methyl acrylate showed the peak of carbonyl stretching at 1716 cm-1 while the absorptions at 2854 cm-1 and 2927 cm-1 are characteristics of CH2 groups in alkyl chains. The spectra of acid hydrazide showed a series of bands characteristic to the hydrazide (–CO–NH–NH2) residue. A number of bands observed due to (NH2) symmetric stretching vibration and (NH2) asymmetric stretching vibration in the range of 3400–3100 cm-1. The peak at 3200 cm-1 attributes to (NH) of the secondary NH group. Similarly, the hydrazide group of amides showed the spectra of acid hydrazide, and acid hydrazone also showed the IR bands for both in-plane deformation and out-of-plane deformation of C = O at 670 cm-1 And 510 cm-1, respectively.

    Figure 3 show how the initial concentration of Cu(Ⅱ) ions affect the residual concentration and percentage removal, respectively. Our findings showed that the increase in the initial concentration of Cu(Ⅱ) led to an increase in the residual concentration and a decrease in percentage removal. Also, the number of sites available for adsorption decreased with the increase in concentration. This can be explained by the competitive diffusion process of Cu(Ⅱ) ions through the micro-channel and pores in PMA When metal ions are adsorbed, the inlet of channel on the surface gets locked. This does not allow the metal ions to pass deeply inside the PMA since the adsorption process occurs only on the surfaces. Previous studies have also showed similar results for natural zeolite and olive cake [13,14]. They found that the maximum adsorption of Cu2+ ions achieved at 30 ℃ with a particle size of 45 mm using natural zeolite [15]. In another study, the maximum removal of Cu2+ occurred at 28 ℃ with an initial pH of the solution of 4.5 by using olive cake as absorbent [15]. Besides, Karthikeyan et al. did the experiment of adsorbing Fe3+ ions using chitin and Rao et al. used carbon absorbents to remove the color of wastewater in the dye manufacturing industry [16,17] (see Figure 3).

    Figure 3.  Effect of concentration on the percentage removal at 150 rpm, 0.1 g of polymer, Temperature = 25 ℃, and pH = 4.6.

    The effluences of PMA masses on percentage removal of Cu(Ⅱ) are shown in Figure 4. It can be observed that the PMA amount increases with the decrease in residual concentration and simultaneously the percentage removal of Cu (Ⅱ) increases. The higher concentration of PMA signifies that a larger surface area is available, forming a large number of adsorption sites which in turn increases the capacity of adsorption [18] (see Figure 4).

    Figure 4.  Effect of adsorbent dosage the percentage removal at 150 rpm, 0.25 M CuSO4, Temperature = 25 ℃, and pH = 4.6.

    Shaking speed is a crucial parameter on the adsorption. Figure 5 show how shaking speed affects the residual concentration and percentage removal of Cu(Ⅱ). An increase in the shaking speed decreased the residual concentration and increased the percentage removal. Increased shaking speed increased the stirring rate, which further increased the percentage removal of metal ions due to the diffusion of the ions to the adsorbent’s surface. This also reduced the film boundary layer around the adsorbent [19]. In fact, the efficiency increased because of increased turbulence and decreased external mass transfer resistance thickness around the adsorbent particles [20] (see Figure 5).

    Figure 5.  Effect of rpm on the percentage removal at constant 0.1 g of polymer, 0.25 M CuSO4, Temperature = 25 ℃, and pH = 4.6.

    Figure 6 denote the temperature effects on the residual concentration and percentage removal of Cu(Ⅱ), respectively. Our findings showed that the increase in temperature decreased the percentage removal decreases due to decreased adsorption or the weak attractive forces between the metal ions and adsorbent surface. The metal ions have an increased tendency to escape from the surface of the biomass in the solution phase at high temperatures. This results in decreased thickness of the boundary layer and decreased adsorption (Aksu and Kutsal, 1991). In addition, the increased temperature increases the loading capacity and decreases the initial adsorbate concentration. This indicates the significance of the rate of desorption compared to the rate of adsorption; thus, this implies that the nature of the reaction process is exothermic [21] (see Figure 6).

    Figure 6.  Effect of temperature on the percentage removal at 0.1 g of polymer, 0.25 M CuSO4, 150 rpm, and pH = 4.6.

    Langmuir isotherm: The Langmuir adsorption is commonly described as the equilibrium between an adsorbate and adsorbent system in which the adsorbate behaves like an ideal gas at the isothermal condition. A pressure relationship called Langmuir Isotherm describes the release of adsorbed gas. The isothermal Langmuir adsorption assumes the gas attaches to the surface and covers as a single layer. The results obtained from the Langmuir model for the removal of Cu2+ by using poly methyl acrylate as absorbent are shown in Figure 7a. The experiments showed strong positive evidence on the adsorption of Cu2+ onto poly methyl acrylate by following the Langmuir isotherm. The high correlation coefficient, R2 > 0.99 proved the effectiveness and applicability of the linear forms of Langmuir models, indicating its feasibility in the sorption system. Meanwhile, the results obtained by using the Freundlich isotherm in which R2 = 0.9169 are shown in Figure 7b. By comparing the above two models, we conclude that Langmuir isotherm is the best model to fit the data as R2 value is higher here (see Figure 7).

    Figure 7.  (a) Langmuir isotherm model, and (b) Freundlich isotherm model.

    Adsorption kinetics must be considered for the accurate design of the adsorption method [22]. It is used in calculate adsorption mechanism, adsorption features and the interaction between adsorbent and adsorbate. Both pseudo-first-order and pseudo-second-order equations have been applied here to determine the adsorption mechanism based on the present experimental data [11,23]. Figure 8a presents the plot of pseudo-first-order kinetics for the adsorption of Cu2+ on PMA in which qe and k values were calculated from the slope. In the plot, the intercept of the line was used to calculate the correlation coefficient, R2. Figure 8b depicts the plot of pseudo-second-order kinetics for the adsorption of Cu2+ on PMA. Similarly, the values of qt and k2 were calculated from the slope and obtained R2 from the intercept of the line. Comparing the two plots presented in Figures 8, it is seen that R2 values is much higher in pseudo-first-order model than in pseudo-second-order model, which confirms that pseudo-first-order model is the best model to fit the data (see Figure 8).

    Figure 8.  (a) Pseudo-first-order model, and (b) Pseudo-second-order model.

    In conclusion, the PMA can be prepared and treated easily. PMA has been proven to be effective in cupric ion (Cu2+) removal from effluent system. The correlation coefficient, R2, value suggests that the adsorption model by Langmuir correctly describes the PMA adsorption system at equilibrium. Pseudo-first-order model is the best model which represents the behavior of adsorption process with time and to identify the adsorption rate at room temperature. Further studies may consider the following recommendations: (a) effect of using the lowest available concentration; (b) decreasing the temperature during the process of adsorption; (c) increasing higher speed for agitation (rpm); (d) increasing the quantity of the polymer being used, i.e., PMA.

    The authors declare that there is no conflict of interest regarding the publication of this paper.



    [1] H. L. Smith, Monotone Dynamical Systems. An Introduction to the Theory of Competitive and Cooperative Systems, American Mathematical Society, 1995. https://doi.org/10.1090/surv/041
    [2] J. D. Murray, Mathematical Biology I: An Introduction. 3rd edition, Springer-Verlag, Berlin, 2002.
    [3] M. W. Hirsch, H. L. Smith, Monotone dynamical systems, in Handbook of Differential Equations: Ordinary Differential Equations, Elsevier, Amsterdam, 2005. https://doi.org/10.1016/S1874-5725(05)80006-9
    [4] X. Q. Zhao, Dynamical systems in population biology. 2nd edition, Springer, New York, 2017. https://doi.org/10.1007/978-3-319-56433-3
    [5] J. Smillie, Competitive and cooperative tridiagonal systems of differential equations, SIAM J. Math. Anal., 15 (1984), 530–534. https://doi.org/10.1137/0515040 doi: 10.1137/0515040
    [6] H. L. Smith, Periodic tridiagonal competitive and cooperative systems of differential equations, SIAM J. Math. Anal., 22(1991), 1102–1109. https://doi.org/10.1137/0522071 doi: 10.1137/0522071
    [7] M. Gyllenberg, Y. Wang, Periodic tridiagonal systems modeling competitive-cooperative ecological interactions, Discrete Contin. Dyn. Syst. B, 5 (2005), 289–298. https://doi.org/10.3934/dcdsb.2005.5.289 doi: 10.3934/dcdsb.2005.5.289
    [8] C. Fang, M. Gyllenberg, Y. Wang, Floquet bundles for tridiagonal competitive-cooperative systems and the dynamics of time-recurrent systems, SIAM J. Math. Anal., 45 (2013), 2477–2498. https://doi.org/10.1137/120878021 doi: 10.1137/120878021
    [9] C. A. Klausmeier, Successional state dynamics: A novel approach to modeling nonequilibrium foodweb dynamics, J. Theor. Biol., 262 (2010), 584–595. https://doi.org/10.1016/j.jtbi.2009.10.018 doi: 10.1016/j.jtbi.2009.10.018
    [10] C. F. Steiner, A. S. Schwaderer, V. Huber, C. A. Klausmeier, E. Litchman, Periodically forced food chain dynamics: model predictions and experimental validation, Ecology, 90 (2009), 3099–3107. https://doi.org/10.1890/08-2377.1 doi: 10.1890/08-2377.1
    [11] S. B. Hsu, X. Q. Zhao, A Lotka-Volterra competition model with seasonal succession, J. Math. Biol., 64 (2012), 109–130. https://doi.org/10.1007/s00285-011-0408-6 doi: 10.1007/s00285-011-0408-6
    [12] L. Niu, Y. Wang, X. Xie, Carrying simplex in the Lotka-Volterra competition model with seasonal succession with applications, Discrete Contin. Dyn. Syst. B, 26 (2021), 2161–2172. https://doi.org/10.3934/dcdsb.2021014 doi: 10.3934/dcdsb.2021014
    [13] X. Xie, L. Niu, Global stability in a three-species Lotka-Volterra cooperation model with seasonal succession, Math. Meth. Appl. Sci., 44 (2021), 14807–14822. https://doi.org/10.1002/mma.7744 doi: 10.1002/mma.7744
    [14] Y. Zhang, X. Q. Zhao, Bistable travelling waves for reaction and diffusion model with seasonal succession, Nonlinearity, 26 (2013), 691–709. https://doi.org/10.1088/0951-7715/26/3/691 doi: 10.1088/0951-7715/26/3/691
    [15] D. Xiao, Dynamics and bifurcations on a class of population model with seasonal constant-yield harvesting, Discrete Contin. Dyn. Syst. B, 21 (2016), 699–719. https://doi.org/10.3934/dcdsb.2016.21.699 doi: 10.3934/dcdsb.2016.21.699
    [16] P. G. Barrientos, J. A. Rodriguez, A. Ruiz-Herrera, Chaotic dynamics in the seasonally forced SIR epidemic model, J. Math. Biol., 75 (2017), 1655–1668. https://doi.org/10.1007/s00285-017-1130-9 doi: 10.1007/s00285-017-1130-9
    [17] M. Han, X. Hou, L. Sheng, C. Wang, Theory of rotated equations and applications to a population model, Discrete Contin. Dyn. Syst., 38 (2018), 2171–2185. https://doi.org/10.3934/dcds.2018089 doi: 10.3934/dcds.2018089
    [18] Y. Tang, D. Xiao, W. Zhang, D. Zhu, Dynamics of epidemic models with asymptomatic infection and seasonal succession, Math. Biosci. Eng., 14 (2017), 1407–1424. https://doi.org/10.3934/mbe.2017073 doi: 10.3934/mbe.2017073
    [19] H. Amann, Fixed point equations and nonlinear eigenvalue problems in Ordered Banach Spaces, SIAM Rev., 18 (1976), 620–709. https://doi.org/10.1137/1018114 doi: 10.1137/1018114
    [20] X. Liang, J. Jiang, On the finite-dimensional dynamical systems with limited competition, Trans. Am. Math. Soc., 354 (2002), 3535–3554. https://doi.org/10.1090/S0002-9947-02-03032-5 doi: 10.1090/S0002-9947-02-03032-5
    [21] E. C. Balreira, S. Elaydi, R. Luís, Global stability of higher dimensional monotone maps, J. Differ. Equation Appl., 23 (2017), 2037–2071. https://doi.org/10.1080/10236198.2017.1388375 doi: 10.1080/10236198.2017.1388375
    [22] A. Ruiz-Herrera, Exclusion and dominance in discrete population models via the carrying simplex, J. Differ. Equation Appl., 19 (2013), 96–113. https://doi.org/10.1080/10236198.2011.628663 doi: 10.1080/10236198.2011.628663
    [23] A. Ruiz-Herrera, Topological criteria of global attraction with applications in population dynamics, Nonlinearity, 25 (2012), 2823–2841. https://doi.org/10.1088/0951-7715/25/10/2823 doi: 10.1088/0951-7715/25/10/2823
    [24] H. L. Smith, Competing subcommunities of mutualists and a generalized Kamke theorem, SIAM J. Math. Anal., 46 (1986), 856–874. https://doi.org/10.1137/0146052 doi: 10.1137/0146052
    [25] H. L. Smith, P. Waltman, The Theory of the Chemostat, Cambridge University Press, 1995. https://doi.org/10.1017/CBO9780511530043
    [26] Y. Wang, J. Jiang, Uniqueness and attractivity of the carrying simplex for discrete-time competitive dynamical systems, J. Differ. Equations, 186 (2002), 611–632. https://doi.org/10.1016/S0022-0396(02)00025-6 doi: 10.1016/S0022-0396(02)00025-6
  • This article has been cited by:

    1. Victor Moreno, Baltazar Espinoza, Kamal Barley, Marlio Paredes, Derdei Bichara, Anuj Mubayi, Carlos Castillo-Chavez, The role of mobility and health disparities on the transmission dynamics of Tuberculosis, 2017, 14, 1742-4682, 10.1186/s12976-017-0049-6
    2. Quentin Richard, Samuel Alizon, Marc Choisy, Mircea T. Sofonea, Ramsès Djidjou-Demasse, Gabor Balazsi, Age-structured non-pharmaceutical interventions for optimal control of COVID-19 epidemic, 2021, 17, 1553-7358, e1008776, 10.1371/journal.pcbi.1008776
    3. Hao Kang, Xi Huo, Shigui Ruan, On first-order hyperbolic partial differential equations with two internal variables modeling population dynamics of two physiological structures, 2021, 200, 0373-3114, 403, 10.1007/s10231-020-01001-5
    4. Biao Tang, Yanni Xiao, Jianhong Wu, Implication of vaccination against dengue for Zika outbreak, 2016, 6, 2045-2322, 10.1038/srep35623
    5. P. J. Dodd, C. Pretorius, B. G. Williams, 2019, Chapter 3, 978-3-030-29107-5, 25, 10.1007/978-3-030-29108-2_3
    6. Sergey Kabanikhin, Olga Krivorotko, Victoriya Kashtanova, Varvara Latyshenko, 2017, Identification the mathematical model of the transmission TB/HIV co-infection in endemic areas, 978-1-5386-1596-6, 77, 10.1109/SIBIRCON.2017.8109841
    7. Yu Mu, Tsz-Lik Chan, Hsiang-Yu Yuan, Wing-Cheong Lo, Transmission Dynamics of Tuberculosis with Age-specific Disease Progression, 2022, 84, 0092-8240, 10.1007/s11538-022-01032-4
    8. Ai Lei, Gopal Chaudhary, Knowledge Transfer Analysis and Management of Virtual Enterprises Based on Structured Cognitive Computing, 2022, 2022, 1687-5273, 1, 10.1155/2022/4858434
    9. Jingjing Lu, Qiaoling Chen, Zhidong Teng, Tingting Zheng, Dynamical analysis of an age‐structured SIRE epidemic model with two routes of infection in environment, 2022, 148, 0022-2526, 461, 10.1111/sapm.12447
    10. Quentin Richard, Marc Choisy, Thierry Lefèvre, Ramsès Djidjou-Demasse, Human-vector malaria transmission model structured by age, time since infection and waning immunity, 2022, 63, 14681218, 103393, 10.1016/j.nonrwa.2021.103393
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1630) PDF downloads(105) Cited by(1)

Figures and Tables

Figures(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog