Research article Special Issues

Bifurcation structure of nonconstant positive steady states for a diffusive predator-prey model

  • In this paper, we make a detailed descriptions for the local and global bifurcation structure of nonconstant positive steady states of a modified Holling-Tanner predator-prey system under homogeneous Neumann boundary condition. We first give the stability of constant steady state solution to the model, and show that the system exhibits Turing instability. Second, we establish the local structure of the steady states bifurcating from double eigenvalues by the techniques of space decomposition and implicit function theorem. It is shown that under certain conditions, the local bifurcation can be extended to the global bifurcation.

    Citation: Dongfu Tong, Yongli Cai, Bingxian Wang, Weiming Wang. Bifurcation structure of nonconstant positive steady states for a diffusive predator-prey model[J]. Mathematical Biosciences and Engineering, 2019, 16(5): 3988-4006. doi: 10.3934/mbe.2019197

    Related Papers:

    [1] Yasir Nadeem Anjam, Mehmet Yavuz, Mati ur Rahman, Amna Batool . Analysis of a fractional pollution model in a system of three interconnecting lakes. AIMS Biophysics, 2023, 10(2): 220-240. doi: 10.3934/biophy.2023014
    [2] Mohammed Alabedalhadi, Mohammed Shqair, Ibrahim Saleh . Analysis and analytical simulation for a biophysical fractional diffusive cancer model with virotherapy using the Caputo operator. AIMS Biophysics, 2023, 10(4): 503-522. doi: 10.3934/biophy.2023028
    [3] Shaimaa A. M. Abdelmohsen, D. Sh. Mohamed, Haifa A. Alyousef, M. R. Gorji, Amr M. S. Mahdy . Mathematical modeling for solving fractional model cancer bosom malignant growth. AIMS Biophysics, 2023, 10(3): 263-280. doi: 10.3934/biophy.2023018
    [4] Stephen Haller, Rabin Gerrah, Sandra Rugonyi . Towards virtual surgery planning: the modified Blalock-Taussig Shunt. AIMS Biophysics, 2020, 7(3): 169-188. doi: 10.3934/biophy.2020014
    [5] Omar El Deeb, Joseph El Khoury Edde . COVID19 vaccines as boosters or first doses: simulating scenarios to minimize infections and deaths. AIMS Biophysics, 2024, 11(2): 239-254. doi: 10.3934/biophy.2024014
    [6] Etimad Alattar, Eqbal Radwan, Khitam Elwasife . Improvement in growth of plants under the effect of magnetized water. AIMS Biophysics, 2022, 9(4): 346-387. doi: 10.3934/biophy.2022029
    [7] Mary Jane Beilby, Sabah Al Khazaaly . Re-modeling Chara action potential: II. The action potential form under salinity stress. AIMS Biophysics, 2017, 4(2): 298-315. doi: 10.3934/biophy.2017.2.298
    [8] Mehmet Yavuz, Kübra Akyüz, Naime Büşra Bayraktar, Feyza Nur Özdemir . Hepatitis-B disease modelling of fractional order and parameter calibration using real data from the USA. AIMS Biophysics, 2024, 11(3): 378-402. doi: 10.3934/biophy.2024021
    [9] Natarajan Mala, Arumugam Vinodkumar, Jehad Alzabut . Passivity analysis for Markovian jumping neutral type neural networks with leakage and mode-dependent delay. AIMS Biophysics, 2023, 10(2): 184-204. doi: 10.3934/biophy.2023012
    [10] Jean Berthier, Kenneth A. Brakke, David Gosselin, Maxime Huet, Erwin Berthier . Metastable capillary filaments in rectangular cross-section open microchannels. AIMS Biophysics, 2014, 1(1): 31-48. doi: 10.3934/biophy.2014.1.31
  • In this paper, we make a detailed descriptions for the local and global bifurcation structure of nonconstant positive steady states of a modified Holling-Tanner predator-prey system under homogeneous Neumann boundary condition. We first give the stability of constant steady state solution to the model, and show that the system exhibits Turing instability. Second, we establish the local structure of the steady states bifurcating from double eigenvalues by the techniques of space decomposition and implicit function theorem. It is shown that under certain conditions, the local bifurcation can be extended to the global bifurcation.


    In the middle of the 20th century, the disease of Poliomyelitis (polio) was one of the contagious diseases which has been transmitted from the western regions of the Asian continent. Jonas Salk, in the year 1952 developed a vaccination process for this outbreak. Poliovirus became the cause of this outbreak among most societies. It affects the nerves of the human body and may cause whole-body paralysis within a short duration. The associated disease mainly affects children in the first three years of age but it may also affect anyone at any age. The poliovirus is transmitted to the children's bodies through the mouth, eventually reaching the intestines of the digestive system. The poliovirus can be transferred from one person to the other. The starting signs of the associated infection are fever, headache, fatigue, vomiting, and pain produced in the limbs. There is still no effective treatment for the virus, but its spread can be controlled through the process of vaccination and immunizations which commonly save many children's lives.

    Mostly, assumptions regarding dynamic problems of daily life are necessary to manage these problems; therefore, in every mathematical model the infections are only transmitted to the susceptible or healthy class due to the infective population. Several diseases like polio are contagious during the time of incubation. So one can conclude that the interactions between the healthy population classes exposed classes have their importance in the spread of diseases. Thus, scholars have discussed the relationship between the susceptible and the exposed classes. Many scientists have attempted to construct mathematical models of Polio. The authors [1] developed a mathematical system spread of polio after the re-introduction of poliovirus. Garfinkel and Sarewitz et al. [2] briefly studied the conditions that would lead to polio eradication. Furthermore, to prove efficacy of polio vaccines in the control of polio infection, most scholars have constructed various models for the spread of polio [3][5]. The authors of [6] specifically considered the SIR problem with pulse vaccination and have shown that pulse vaccination can lead to eradication under specific conditions related to the magnitude of vaccination proportion. In some of the articles like [7], the authors dealt with a two-dimensional SIS polio system with the addition of vaccination terms for bifurcation in the opposite direction. The authors of [8] developed a polio mathematical system to investigate its dynamical behavior in the presence of vaccinations strategy with nonlinear incidence, which also deals with the optimality of the model that converges the vaccine quantity to the threshold value. Agarwal and Bhadauria [9] developed a polio epidemic model with vaccine quantity to investigate the effect of vaccination once a strategy is applied to the susceptible and exposable populations. In this study, the polio model was adapted from [9], and it has been divided into five subgroups, namely susceptible persons S, exposed persons E, acutely infected persons I, chronically infected persons V and recovered persons R:

    ((1.1))

    with the following initial conditions

    The parameters used in model (1.1) are described as follows: A is the rate of immigration in the society, µ is the rate of death that occurs naturally, β is the probability of disease transmission because of the infectious population, is the rate of chance of infection because of the exposed population, r is the rate of reduction in the transmission of infection by exposed class, ν is the rate of individuals moving from the susceptible to vaccinated class, ν1 is the rate at which exposed individuals are vaccinated, b is the rate at which exposed individuals move to the infective class and α is the rate of death due to the disease.

    The utilization of fractional calculus in the fractional-order model, which combines differentiation and integration, offers a more effective approach to comprehending real-world problems as compared to classical derivatives [10][16]. The concept of fractional derivatives, originally introduced by Riemann-Liouville based on the power law, has recently faced scrutiny from researchers due to the presentation of experimental results questioning the adequacy of a single fractional operator, such as the Caputo, Caputo-Fabrizio and Atangana-Baleanu operators, as a tool to describe complex phenomena in science and engineering [17][21]. These findings suggest that a broader range of fractional operators may be necessary to accurately capture the behavior of diverse systems [22][26]. The authors of [27] developed a model by using a new definition of the constant proportional Caputo operator, which describes the generalized memory effects. A new recursive algorithm has been constructed to solve certain initial value problems involving a fractional differential equation [28]. This new approach Adomian decomposition method (ADM) is based on the application of the ADM, and it involves combining the decomposition with a recurrence formula and utilizing the solutions of the generalized Abel equation. According to the researchers, a new fractional derivative can be defined by utilizing the exponential kernel [29], [30]. The utilization of non-singular kernel fractional derivatives in the modeling of epidemics offers valuable insights into the dynamics of infectious disease outbreaks, particularly when applied to trigonometric and exponential functions [31][33]. Initialization can pose challenges when working with non-singular kernel fractional derivatives, especially in cases involving non-singular kernels. Regardless of their type, equations in the following form have been observed to present similar challenges, as noted by the authors of [34]:

    xaF(x,τ)g(τ)dτ=y(x),axb.

    The condition F(a,a)0 introduces distinct limitations on differential equations involving non-singular kernels, implying that if F(x,τ) and y(x) are continuous and F(a,a)0, then y(a)0. This condition gives rise to peculiar behaviors in these equations, necessitating careful consideration during the development of models in various research fields. As mentioned in [34], this issue remains unsolved. However, the problem associated with the operator introduced in [17] has been solved by various researchers [35], [36]. To overcome the challenges associated with non-singular operators, the authors of [37] proposed a modification to the operator that involves utilizing a Mittag-Leffler kernel. They demonstrated that the resulting fractional differential equations, based on this modified operator, are easier to initialize than those based on non-singular kernels. Additionally, they showcased that the modified operator, known as the modified Atangana-Baleanu-Caputo (mABC) derivative, is capable of solving multiple fractional differential problems, which is not possible for the Atangana-Baleanu-Caputo (ABC) derivative. The mABC derivative possesses an integrable singularity at the origin. In [38], the authors developed a unique numerical method for the mABC derivative, utilizing finite differences; the method streamlines the initialization of the corresponding fractional differential equations. Researchers have also developed modified fractional difference operators by employing the mABC derivative and Mittag-Leffler kernels [39].

    We propose a memory-affected mABC model to model the spread of polio. Although the epidemic model described by Eq (1.1) utilizes classical derivatives and should be considered, our memory-affected mABC model provides additional insights and enhances the accuracy of hepatitis B spread predictions. To convert the ordinary system of equations given by Eq (1.1) into the mABC operator, we make the following transformation:

    ((1.3))

    with the following initial conditions

    The rest of the paper is arranged as follows. In Section 2, we recall some basic definitions and statements from the literature. In Section 3, we give some formulas and show the non-zero solution to the homogeneous fractional initial value problem. With the use of the Laplace Adomian decomposition method, we find the approximate solution for the mABC derivative in Section 4. In the same section, we also establish the stability and uniqueness results for the considered system. The approximate series solution is graphically presented in Section 5. We conclude our work in Section 6.

    Here, we recall basic results from the literature on fractional calculus.

    Definition 2.1. [37] Let f(t)L1(0,T) be a function; then, the mABC derivation is presented as follows:

    ((2.1))

    where Eθ is known as the Mittag-Leffler function for one parameter while Eθ,θ denotes the Mittag-Leffler for two parameters; from the above definition one may prove that $^\mathrm{mABC}\mathbf{D}_t^\vartheta\mathrm{f}=0.

    Definition 2.2. [37] Let f(t)L1(0,T) be a function, then, the integral of mABC operator is given as

    ((2.2))

    Lemma 2.2.1. For fL1(0,) and the order ϑ(0,1), we obtain

    ((2.3))

    The Laplace transform for the mABC is defined as

    ((2.4))

    We state that there is a non-zero solution to the homogeneous fractional initial value problem. To achieve that, we use the following formulas:

    ((3.1))
    ((3.2))

    Lemma 3.0.1. [37] Suppose that the fractional initial value problem is as follows:

    ((3.3))

    where 0<ϑ<1.

    (1) For Ω=M(ϑ)1ϑ, the solution is as follows:

    ((3.4))

    (2). For ΩM(ϑ)1ϑ, the solution is as follows:

    ((3.5))

    where, τϑ=Ω(1ϑ)M(ϑ).

    Proof. (1) Given that

    t0(ty)ϑ1Eϑ,ϑ(µϑ(1y)ϑ)yϑdy=Γ(1ϑ)Eϑ(µtϑ),

    for 0<t, we have

    ((3.7))
    ((3.8))
    ((3.9))
    ((3.10))
    ((3.11))

    which shows that the proof is complete.

    Applying Eq.(3.1) and Eq.(3.2) for t>0, we have

    Thus, the proof is complete.

    Lemma 3.0.2. [37] Assume the following fractional differential equations

    where 0<ϑ<1 and ΩM(ϑ)1ϑ; the solution of the above initial value system is given by

    ((3.12))

    where

    ((3.13))

    and zϑ=M(ϑ+Ω(1ϑ)).

    Proof. Utilizing Eq (3.1) and (3.2) we may verify that

    ((3.14))

    by Eq (2.4), we obtain

    ((3.15))

    Further calculation then yields the following:

    ((3.16))

    By substituting Eq (3.14) into Eq. (3.16), the following results are obtained:

    ((3.17))
    ((3.18))

    This finishes the proof.      �

    Remark 1. Suppose that GiC[0,T], with i=1,2,3,4,; then, we have

    ((3.19))

    subject to the following conditions:

    ((3.20))

    Therefore, ˜S=S0,˜E=E0,˜I=I0 and ˜V=V0; hence, the soultion given by Eq (3.12) presents the condition of continuity. On the basis of the above conditions the solution exists for the proposed model.

    Applying the Laplace transform to the mABC, as given in [37] to the considered model (1.3) gives

    ((4.1))

    Proceeding we get

    ((4.2))

    By application of the initial conditions we get

    ((4.3))

    Decomposing each quantity as yields

    We write the non-linear term as follows:

    where

    Gi=1Γ(i+1)didpi[k=0pkSjVj]|p=0.

    Plugging all of the above values into Eq (4.3) we get

    ((4.4))

    after calculating the above Eq. (4.4) we get

    ((4.5))
    ((4.6))
    ((4.7))
    ((4.8))

    Considering the first three terms and by applying the inverse Laplace transform we get

    ((4.9))
    ((4.10))
    ((4.11))
    ((4.12))

    By putting the values S1, E1, I1, V1 into the second term of each quantity, we get

    ((4.13))
    ((4.14))
    ((4.15))
    ((4.16))

    Lastly, each quantity can be written as follows:

    S(t)=S0+S1(t)+S2(t)\hdots
    E(t)=E0+E1(t)+E2(t)\hdots
    I(t)=I0+I1(t)+I2(t)\hdots
    V(t)=V0+V1(t)+V2(t)\hdots

    Theorem 1. Suppose that a Banach space is denoted by (B,|.|) with a mapping T:BB which satisfies

    TxTyΘTxTy+πxy

    for every x,yB; also, 0Θand0π<1. Then, T is Picard T-stable.

    Theorem 2. Suppose T to be a self map defined as given below:

    ((4.17))

    Thus, the iteration is T-stable in L1(x,y), if we obtained the following results:

    (1a(N1+N4)w1(θ)(η2+µ)w2(θ))<1,(1+a(N1+N4)w1(θ)(α+µ)w3(θ))<1,(1+αw4(θ)(µ+a1+a2)w5(θ))<1,(1+a2w5(θ)(η1+a3+µ)w6(θ))<1,

    Proof. We need to prove that T has a fixed point; for this we use (q,p)N×N

    ((4.19))

    Taking the first equation of Eq (4.19) and calculating the norm on both sides, we have

    ((4.20))

    with the help of the triangular inequality, and by solving Eq (4.20), we obtain

    ((4.21))

    By establishing the relation as follows:

    ((4.22))

    replacing the above relation in Eq (4.21), we get the following relation

    ((4.23))

    Furthermore, the convergent sequence Vq is bounded.

    Additionally, one may obtain different constants N1, N2, N3 and N4 for every t such that

    ((4.24))

    Therefore, considering Eq (4.23) and Eq (4.24), we obtain

    ((4.25))

    here, w1 and w2 are the functions of L1[(1ϑ)(sϑ+μϑ)M(ϑ)sϑL]. Considering the same procedure for the remaining equations, we have

    ((4.26))

    Hence, T has a fixed point. Applying Eq (4.25) and (4.26), we assume that

    ψ=(0,0,0,0),

    Ψ={(1a(N1+N4)w1(θ)(η2+µ)w2(θ)),(1+a(N1+N4)w1(θ)(α+µ)w3(θ)),(1+αw4(θ)(µ+a1+a2)w5(θ)),(1+a2w5(θ)(η1+a3+µ)w6(θ)),

    Therefore, the conditions of Theorem 1 are satisfied; hence, the proof is complete.

    Theorem 3. The aforementioned method gives a unique solution for the considered model.   �

    The graphical representations for the proposed model under the modified operator have been established by using the data taken from [9]. Four studied compartments of the model were tested on different fractional orders by using the obtained numerical scheme. The initial values and corresponding parameter descriptions are given in Table 1.

    Table 1.  Parameters and their numerical values for model (1.3).
    Notation Value Source Notation Value Source
    S0 500 [9] E0 200 [9]
    I0 150 [9] E0 600 [9]
    A 1000 [9] β 0.002, 0.003, 0.0016, 0.05 [9]
    r 0.5 [9] µ 0.5 [9]
    ν 0.6 [9] ν1 0.001 [9]
    b 0.9 [9] α 0.6 [9]

     | Show Table
    DownLoad: CSV
    Figure 1.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.002.

    In Figure 1(a–d), the dynamics of all four compartments are shown for different fractional orders under the modified fractional operator which has extra terms as compared to the ABC operator. The susceptible population gains stability after increasing to some peak value and the stability is best achieved for relatively small fractional orders instead of the higher order. The exposed class declines and then increases to the equilibrium point to achieve stability. During this period, the exposed population is transferred to the infected population, whose density is increasing and then stabilizes, or decreases as a result of individuals recovering from the infection after the vaccination is applied to the population. The vaccination process is terminated once the infection is controlled.

    Figure 2.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.003.

    In Figure 2(a–d), the dynamics of all four compartments are shown for different arbitrary orders under the modified fractional operator with extra informative terms. The susceptible class gains stability after declining to some small value and the stability is achieved faster with small fractional orders. The exposed class here declines quickly and then increases to the equilibrium point to achieve the convergence. During this period, the susceptible population is transferred to the exposed class, and then to the infected class whose density increases and then stabilizes, or decreases as a result of recovering from the infection after the vaccination is injected into the total population. The vaccination process is reduced as the transmission is increased and the infection is controlled after more populations.

    Figure 3.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.0016.

    In Figure 3(a–d), the rate of transmission is decreased which affects the dynamics of all four compartments as shown for different arbitrary orders. The susceptible class increases and gains stability after reaching some high value. The exposed class here declines much more quickly and then converges to its equilibrium point. In this case, after a small increase, the infection is controlled quickly. Here the vaccination rate also increases, i.e. as the vaccination rate increases the infection may be easier to control.

    Figure 4.  Plots for all four agents on six different fractional orders ϑ lying between 0 and 1; β = 0.0002.

    In Figure 4(a–d), the rate of transmission is again decreased which affects the dynamics of all the four compartments, as shown for different arbitrary orders. In this case, the dynamics are about the same as observed for the previous dynamics but the dynamics can be controlled more quickly. Such dynamics provide the total density of each compartment in the form of a continuous spectrum lying between 0 and 1.

    This study developed the analysis of the impact of vaccination on the spread of polio in the human population by using a four-compartment system with the generalized novel mABC fractional operator. The dynamical behaviors of all compartments have been successfully examined for different fractional orders to demonstrate the validity of an extra degree of freedom in the selection of the derivative order. The fixed-point theory shows the uniqueness of the solution in the generalized novel modified format. The stability analysis of the model has been conducted by T-Picard-type stability techniques. The numerical solution for the model has been achieved by using the Laplace transform, along with the Adomian technique in the format of the decomposition process. Such a study involved the analysis of the complex geometry in the dynamical system. The study was performed by using various fractional orders and iterations with different transmission rate; also, each curve has been plotted for six different fractional orders and compared with the integer-order result. The study also shows that the impact of vaccination is more significant on the total population. In terms of stability, relatively smaller fractional orders yield high accuracy. In all of the numerical simulations, decreasing the transmission rate led to better control of the polio infection.

    The authors declare that they have not used artificial intelligence tools in the creation of this article.



    [1] J. Jang, W.-M. Ni and M. Tang, Global bifurcation and structure of Turing patterns in the 1-D Lengyel-Epstein model, J. Dyn. Diff. Equa., 2 (2004), 297–320.
    [2] A. M. Turing, The chemical basis of morphogenesis, Philos. T. R. SOC. B, 237 (1852), 37–72.
    [3] N. F. Britton, Essential Mathematical Biology, Springer, New York, 2003.
    [4] L. A. Segel and J. L. Jackson, Dissipative structure: an explanation and an ecological example, J. Theor. Biol., 37 (1972), 545–559.
    [5] D. Alonso, F. Bartumeus and J. Catalan, Mutual interference between predators can give rise to Turing spatial patterns, Ecology, 83 (2002), 28–34.
    [6] W.-M. Wang, L. Zhang, H. Wang, et al., Pattern formation of a predator-prey system with Ivlev-type functional response, Ecol. Model., 221 (2008), 131–140.
    [7] C. Neuhauser, Mathematical challenges in spatial ecology, Notices AMS, 48 (2001), 1304–1314.
    [8] A. Okubo and S. A. Levin, Diffusion and Ecological Problems: Modern Perspectives, Springer, New York, 2001.
    [9] A. B. Medvinsky, S. V. Petrovskii, I. A. Tikhonova, et al., Spatiotemporal complexity of plankton and fish dynamics. SIAM Rev., 44 (2002), 311–370.
    [10] J. D. Murray, Mathematical biology. II: Spatial models and biomedical applications, Springer, New York, 2003.
    [11] P. Y. H. Pang and M. Wang, Qualitative analysis of a ratio-dependent predator-prey system with diffusion, Proc. R. Soc. Edinb., 133 (2003), 919–942.
    [12] K. Kuto and Y. Yamada. Multiple coexistence states for a prey-predator system with crossdi ffusion, J. Differ. Equations, 197 (2004), 315–348.
    [13] K. Kuto, Stability of steady-state solutions to a prey-predator system with cross-diffusion, J. Differ. Equations, 197 (2004), 293–314.
    [14] X. Zeng and Z. Liu, Non-constant positive steady states of a prey-predator system with cross-diffusions, J. Math. Anal. Appl., 332 (2007), 989–1009.
    [15] R. Peng and J. Shi, Non-existence of non-constant positive steady states of two Holling type-II predator-prey systems: Strong interaction case, J. Differ. Equations, 247 (2009), 866–886.
    [16] F. Yi, J. Wei and J. Shi, Bifurcation and spatiotemporal patterns in a homogeneous diffusive predator-prey system, J. Differ. Equations, 246 (2009), 1944–1977.
    [17] H. Shi,W.-T. Li and G. Lin. Positive steady states of a diffusive predator-prey system with modified Holling-Tanner functional response, Nonl. Anal. Real, 11 (2010), 3711–3721.
    [18] Y. Cai, M. Banerjee, Y. Kang, et al., Spatiotemporal complexity in a predator-prey model with weak Allee effects, Math. Biosci. Eng., 11 (2014), 1247–1274.
    [19] S. Li, J. Wu and Y. Dong, Turing patterns in a reaction-diffusion model with the Degn-Harrison reaction scheme, J. Differ. Equations, 259 (2015), 1990–2029.
    [20] H. Shi and S. Ruan. Spatial, temporal and spatiotemporal patterns of diffusive predator-prey models with mutual interference, IMA J. Appl. Math., 80 (2015), 1534–1568.
    [21] Y. Cai and W.-M. Wang, Fish-hook bifurcation branch in a spatial heterogeneous epidemic model with cross-diffusion, Nonl. Anal. Real, 30 (2016), 99–125.
    [22] T. Kuniya and J.Wang, Global dynamics of an SIR epidemic model with nonlocal diffusion, Nonl. Anal. Real, 43 (2018), 262–282.
    [23] J.Wang, J.Wang and T. Kuniya, Analysis of an age-structured multi-group heroin epidemic model.Appl. Math. Comp., 347 (2019), 78–100.
    [24] Y. Cai, Z. Ding, B. Yang, et al., Transmission dynamics of Zika virus with spatial structure–A case study in Rio de Janeiro, Brazil. Phys. A, 514 (2019), 729–740.
    [25] Y. Cai, K. Wang and W.M. Wang, Global transmission dynamics of a Zika virus model, Appl. Math. Lett., 92 (2019), 190–195.
    [26] Y. Cai, X. Lian, Z. Peng, et al., Spatiotemporal transmission dynamics for influenza disease in a heterogenous environment, Nonl. Anal. Real., 46 (2019), 178–194.
    [27] Y. Cai, Z. Gui, X. Zhang, et al., Bifurcations and pattern formation in a predator-prey model,. Inter. J. Bifur. Chaos, 28 (2018), 1850140.
    [28] H. Zhang, Y. Cai, S. Fu, et al., Impact of the fear effect in a prey-predator model incorporating a prey refuge, Appl. Math. Comp., 356 (2019), 328–337.
    [29] X. Cao, Y. Song and T. Zhang. Hopf bifurcation and delay-induced Turing instability in a diffusive Iac Operon model, Inter. J. Bifur. Chaos, 26 (2016), 1650167.
    [30] J. Jiang, Y. Song and P. Yu, Delay-induced Triple-Zero bifurcation in a delayed Leslie-type predator-prey model with additive Allee effect, Inter. J. Bifur. Chaos, 26 (2016), 1650117.
    [31] Y. Song, H. Jiang, Q. Liu, et al., Spatiotemporal dynamics of the diffusive Mussel-Algae model near Turing-Hopf bifurcation. SIAM J. Appl. Dyn. Sys., 16 (2017), 2030–2062.
    [32] Y. Song and X. Tang, Stability, steady-state bifurcations, and Turing patterns in a predator-prey model with herd behavior and prey-taxis, Stud. Appl. Math., 139 (2017), 371–404.
    [33] S. Wu and Y. Song, Stability and spatiotemporal dynamics in a diffusive predator Cprey model with nonlocal prey competition, Nonl. Anal. Real, 48 (2019), 12–39.
    [34] S.-B. Hsu and T.-W. Huang, Global stability for a class of predator-prey systems, SIAM J. Appl. Math., 55 (1995), 763–783.
    [35] Y. Lou and W.-M. Ni. Diffusion, self-diffusion and cross-diffusion, J. Differ. Equations, 131 (1996), 79–131.
    [36] W.-M. Ni and M. Tang, Turing patterns in the Lengyel-Epstein system for the CIMA reaction, T. Am. Math. Soc., 357 (2005), 3953–3969.
    [37] P. H. Rabinowitz, Some global results for nonlinear eigenvalue problems, J. Func. Anal., 7 (1971), 487–513.
    [38] Y. Nishiura, Global structure of bifurcating solutions of some reaction-diffusion systems, SIAM J. Math. Anal., 13 (1982), 555–593.
    [39] I. Takagi, Point-condensation for a reaction-diffusion system, J. Differ. Equations, 61 (1986), 208–249.
    [40] A.Chertock, A. Kurganov, X. Wang, et al., On a chemotaxis model with saturated chemotactic flux, Kinet. Relat. Models, 5 (2012), 51–95.
  • This article has been cited by:

    1. Elham Raeisi, Mehmet Yavuz, Mohammadreza Khosravifarsani, Yasin Fadaei, Mathematical modeling of interactions between colon cancer and immune system with a deep learning algorithm, 2024, 139, 2190-5444, 10.1140/epjp/s13360-024-05111-4
    2. Elhoussine Azroul, Sara Bouda, Mathematical modeling of poliomyelitis virus with vaccination and post‐paralytic syndrome dynamics using Caputo and ABC fractional derivatives, 2024, 0170-4214, 10.1002/mma.10406
    3. Mohamed A. Abdoon, Abdulrahman B. M. Alzahrani, Comparative Analysis of Influenza Modeling Using Novel Fractional Operators with Real Data, 2024, 16, 2073-8994, 1126, 10.3390/sym16091126
    4. Elkhateeb S. Aly, Mohammed A. Almalahi, Khaled A. Aldwoah, Kamal Shah, Criteria of existence and stability of an n-coupled system of generalized Sturm-Liouville equations with a modified ABC fractional derivative and an application to the SEIR influenza epidemic model, 2024, 9, 2473-6988, 14228, 10.3934/math.2024691
    5. Lahcen Boulaasair, Hassane Bouzahir, Mehmet Yavuz, Global mathematical analysis of a patchy epidemic model, 2024, 14, 2146-5703, 365, 10.11121/ijocta.1558
    6. Mehmet Yavuz, Feyza Nur Özdemir, Kübra Akyüz, Naime Büşra Bayraktar, The relationship between colon cancer and immune system: a fractional order modelling approach, 2025, 27, 1301-7985, 126, 10.25092/baunfbed.1498218
    7. Isah Abdullahi, D.G. Yakubu, M.Y. Adamu, Musa Ali, A.M. Kwami, Inclined magnetic fields heat transfer and thermal radiation on fractionalized EMHD Burgers’ fluid flow via bifurcated artery for tumor treatments, 2025, 26668181, 101093, 10.1016/j.padiff.2025.101093
    8. R. Rajaraman, Modeling Immobilized Enzyme Reactions: Nonlinear Kinetics With Fractional‐ and Integer‐Order Analysis, 2025, 0170-4214, 10.1002/mma.10791
    9. Ateq Alsaadi, Advancing water quality management: A synergistic approach using fractional differential equations and neural networks, 2025, 10, 2473-6988, 5332, 10.3934/math.2025246
  • Reader Comments
  • © 2019 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(5020) PDF downloads(569) Cited by(2)

Figures and Tables

Figures(1)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog