Next Article in Journal
A Computational Approach to a Model for HIV and the Immune System Interaction
Previous Article in Journal
Extended Gevrey Regularity via Weight Matrices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of a Patchy Predator–Prey Model with Fear Effect

Center for Applied Mathematics of Fujian Province, School of Mathematics and Statistics, Fuzhou University, Fuzhou 350108, China
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 577; https://doi.org/10.3390/axioms11100577
Submission received: 30 June 2022 / Revised: 10 October 2022 / Accepted: 12 October 2022 / Published: 21 October 2022

Abstract

:
In this paper, a predator–prey model with fear effect and dispersal is proposed. Assume that only the prey migrates at a constant rate between patches and the migration of prey on each patch is faster than the time scale of local predator–prey interaction. Using two time scales, an aggregation system of total prey density for two patches is constructed. Mathematical analysis shows that there may exist a trivial, a boundary and a unique positive equilibrium point. Under certain conditions, the corresponding unique equilibrium point is global asymptotically stable. The impact of the fear effect on the system is also investigated, i.e., the predator density decreases when the amount of fear effect increases. Moreover, dispersal has a great impact on the persistence of the predator and the prey. Numerical experiments are also presented to verify the feasibility of our conclusion.

1. Introduction

Predator-prey system [1] is one of the most basic relationships among populations, which plays an important role in maintaining the sustainable survival of populations, protecting biodiversity and preventing biological invasion. The predator–prey system with spatial heterogeneity has more complex dynamic properties and rich dynamic behaviors. Environmental heterogeneity can provide spatial shelter for the prey population. The population interaction model in a patchy environment is an effective method to study environmental heterogeneity [2,3,4,5]. Thus, the dynamic properties of the patchy predator–prey system have been widely studied [6,7,8,9]. Additionally, two time scales are often considered: a fast one in which migrations between patches occur, and a slow one in which interactions, reproduction, and mortality take place. In [10,11,12,13,14], by use of perturbation theory, the full system is aggregated into a reduced system for the total prey and predator populations. For example, Liu [14] proposes a two-patch Leslie–Gower predator–prey model with prey refuge. It is found that the unique positive equilibrium point of the system is always globally asymptotically stable. The above shows that the prey refuge does not change the stability of the system, though it has an effect on the density of the predators and the prey.
At present, not only the direct killing of the prey by the predator is considered, but also the physiological and psychological effects of the predator on the prey are investigated. For example, Wang et al. [15] firstly establish a predator-prey system with fear effect and the dynamic behavior of the system has become complicated under different functional response functions. Consequently, research in this direction has attracted great attention. One can see [16,17,18,19,20,21,22] for details. However, to the best knowledge of the authors, no scholar has studied the patchy model with fear effect under two different time scales. Actually, once the prey is suffering from the fear effect, it often leads to migration between patches. In other words, it is more reasonable to discuss the predator–prey system with fear effect and diffusion. Motivated by the above, in this paper, we assume that the prey in patch 1(2) will migrate to patch 2 due to the fear effect, and the prey in patch 2(1) will not be killed by predators. However, the prey in patch 2 will die exponentially due to lack of other sources of food. In addition, we suppose that migration between patches is much more rapid than the local interaction in the system. In this paper, we shall investigate the following two-patch predator–prey model with fear effect:
d x 1 d τ = α x 2 β x 1 + ε x 1 1 + k y ( r 1 b 1 x 1 ) ε a 1 x 1 y , d x 2 d τ = β x 1 α x 2 + ε ( d x 2 ) , d y d τ = ε ( r 2 + a 2 x 1 ) y ,
where x 1 , y is the density of the prey and the predator in patch 1 at time t , respectively. x 2 represents the density of the prey in patch 2. r 1 , b 1 , a 1 , d , r 2 , a 2 , α , β are all positive constants. Here, r 1 represents the intrinsic growth rate of prey in patch 1. d is the death rate of prey in patch 2. b 1 is the interspecific competition rate. a 1 represents the consumption rate of the predator in patch 1 and a 2 is the conversion ratio into fitness. r 2 is the death rate of the predator. β ( α ) represents the dispersal rate from patch 1(2) to patch 2(1). ε is an extremely small dimensionless positive number and represents the ratio of fast time scale to slow time scale. τ is a fast time variable.
As is known to all, based on perturbation techniques and the application of center manifold theorem, the aggregation method [10,11,12,13,14,23,24,25] reduces a system with a large number of variables involving different time scales into an aggregated system with few global variables. For the complete system (1), an aggregation method will be applied to obtain a reduced system. The reduced form presents a two-dimensional system of ordinary differential equation, which governs the total prey density and the predator density at slow time scale.
In this paper, we choose the aggregated variables as x = x 1 + x 2 , y = y . Adding the first two equations of system (1), we obtain
d x d τ = ε x 1 1 + k y ( r 1 b 1 x 1 ) d x 2 a 1 x 1 y , d y d τ = ε r 2 + a 2 x 1 y .
Let ε = 0 in system (1); we obtain that the solution of the fast part is x 1 * = α x α + β , x 2 * = β x α + β , and then substitute ( x 1 * , x 2 * ) into system (2). Thus, we can get the following aggregation system:
d x d t = α r 1 x ( α + β ) ( 1 + k y ) α 2 b 1 x 2 ( α + β ) 2 ( 1 + k y ) β d x α + β α a 1 x y α + β , d y d t = ( r 2 + α a 2 α + β x ) y ,
where t = ε τ .
As was pointed out in [10,11,12,13,14,24,25], the complete model is topologically equivalent to the reduced model, and it gives a good approximation solution of the complete model by using the aggregated model. Thus, in the following sections, we will consider the dynamic behaviors of the aggregated model (3).
In order to simplify system (3), we take the following transformations
τ = l t , x ¯ = g x , y ¯ = h y ,
where
l = r 2 , g = a 2 m r 2 , h = a 1 m r 2 m = α α + β .
We still reserve t, x, y to express τ , x ¯ , y ¯ , respectively. Then, we can obtain the system as follows.
d x d t = e x a + y c x 2 a + y f x x y , d y d t = y ( x 1 ) ,
where
a = a 1 m k r 2 , e = r 1 a 1 m 2 k r 2 2 , c = b 1 a 1 m 2 k r 2 a 2 , f = ( 1 m ) d r 2 .
It is not difficult to verify that for system (4), under non-negative initial values, i.e., x ( 0 ) 0 , y ( 0 ) 0 , the solution ( x ( t ) , y ( t ) ) is non-negative. Moreover, let w ( t ) = x ( t ) + y ( t ) . Suppose that 0 < σ < max { 1 , f } . Then, we have d w ( t ) d t + σ w ( t ) x a + y ( e c x ) e 2 4 a c . Thus, the solution ( x ( t ) , y ( t ) ) is ultimately bounded.
The rest of the paper is arranged as follows. In Section 2, the existence and stability of the equilibria are established. In Section 3, we show the global stability of the equilibria. In Section 4, we analyze the impact of fear effect and dispersal on predator density and give some numerical simulation results to verify these analytical results. We summarize the obtained theoretical results and biological application in Section 5.

2. Existence and Local Stability of Equilibrium

In this section, we shall discuss the existence and stability of the equilibrium for system (4). The equilibrium of system (4) is given by the following equation:
e x a + y c x 2 a + y f x x y = 0 , y ( x 1 ) = 0 .
By simple computation, there always exist a trivial equilibrium E 0 ( 0 , 0 ) and a boundary equilibrium E 1 ( e a f c , 0 ) . To obtain the positive equilibrium for the system, we consider the following equation:
y 2 + ( a + f ) y + a f + c e = 0 .
Let the discriminant of Equation (5) denoted by Δ = ( a + f ) 2 4 ( a f + c e ) . If e > a f + c , then Δ > 0 and system (4) has a unique positive equilibrium E 2 ( x * , y * ) . Here, x * = 1 and y * = ( a + f ) + Δ 2 .
Next, we will investigate the stability of the above equilibrium. The Jacobian matrix of system (4) at E ( x , y ) is calculated as
J E = e 2 c x a + y f y x c x e ( a + y ) 2 1 , y 1 + x . .
The dynamics of system (4) are given as follows.
Theorem 1.
(1) 
System (4) always has a trivial equilibrium E 0 ( 0 , 0 ) .
(a) 
If e < a f , E 0 ( 0 , 0 ) is a hyperbolic stable node;
(b) 
If e = a f , E 0 ( 0 , 0 ) is an attracting saddle-node;
(c) 
If e > a f , E 0 ( 0 , 0 ) is a hyperbolic saddle.
(2) 
If e > a f , system (4) has a boundary equilibrium E 1 ( x 1 , 0 ) . Here, x 1 = e a f c .
(a) 
If a f < e < a f + c , E 1 ( x 1 , 0 ) is a hyperbolic stable node;
(b) 
If e = a f + c , E 1 ( x 1 , 0 ) is an attracting saddle-node;
(c) 
If e > a f + c , E 1 ( x 1 , 0 ) is a hyperbolic saddle.
(3) 
When e > a f + c , system (4) has a unique positive equilibrium E 2 ( x * , y * ) , which is always asymptotically stable.
Proof. 
(1) The Jacobin matrix at E 0 ( 0 , 0 ) is
J E 0 = e a f 0 0 1 .
By direct calculation, we can obtain two eigenvalues of J E 0 , i.e., λ 1 = e a f and λ 2 = 1 < 0 . If e > a f , E 0 ( 0 , 0 ) is a hyperbolic saddle. If e < a f , we have λ 1 < 0 , then E 0 ( 0 , 0 ) is a hyperbolic stable node. If e = a f , then λ 1 = 0 . In order to investigate the stability property of the equilibrium E 0 , we introduce the new time variable τ = a a + y t , and then we have
d x d τ = 1 + f a x y + c a x 2 + 1 a x y 2 , d y d τ = y x y 1 a x y 2 + 1 a y 2 ,
the coefficient of x 2 is c a > 0 . Considering the new time variable τ and using Theorem 7.1 in [26], we can conclude that E 0 is an attracting saddle-node.
(2) The Jacobian matrix at E 1 ( x 1 , 0 ) is
J E 1 = e 2 c x 1 a f x 1 ( c x 1 e a 2 1 ) 0 1 + x 1 .
By simple calculation, we get two eigenvalues of J E 1 , i.e., λ 1 = a f e a < 0 , λ 2 = 1 + x 1 . It follows that if a f < e < a f + c , λ 2 < 0 , and so E 1 is a hyperbolic stable node. If e > a f + c , then E 1 is a saddle. When e = a f + c , it follows that λ 2 = 0 . To analyze the stability of E 1 ( x 1 , 0 ) , we make the transformation ( X , Y ) = ( x x 1 , y ) and expand the power series around the origin. Then, system (4) becomes
d X d t = c a X f + a a Y c a X 2 a 2 + a f c a 2 Y X + f a 2 Y 2 + P 1 ( X , Y ) , d Y d t = X Y .
Here, P 1 ( X , Y ) denote the power series with term U i V j satisfying i + j 3 .
Let U = X + f + a c Y , V = Y . Introducing a new time variable τ by τ = c a t and rewriting τ as t, we have
d U d t = U + 1 a U 2 a f ( f a + 1 ) f 2 c a 2 c V 2 a 2 + a f + c a c U V + P 2 ( U , V ) , d V d t = U V f + a c V 2 ,
where P 2 ( U , V ) denote the power series with term U i V j satisfying i + j 3 .
Hence, by Theorem 7.1 in Chapter 2 in [26], the coefficient of V 2 is f + a c < 0 , then the equilibrium E 1 is an attracting saddle node.
(3) The Jacobian matrix at E 2 is
J E 2 = e 2 c a + y * f y * c e ( a + y * ) 2 1 y * 0 .
Under the existence condition of E 2 , i.e., a f + c < e , it is easy to obtain that D e t ( E 2 ) = y * ( c e ( a + y * ) 2 1 ) > 0 and T r ( E 2 ) = e 2 c a + y * f y * = c a + y * < 0 . Thus, E 2 ( x * , y * ) is locally asymptotically stable. □
The corresponding illustration of Theorem 1 is shown in Figure 1 by matlab software.

3. Global Stability of Equilibrium

In this section, for system (4), we will investigate the global stability of equilibria.
Theorem 2.
If e < a f , E 0 ( 0 , 0 ) is globally asymptotically stable.
Proof. 
When e < a f , system (4) has a unique stable equilibrium E 0 ( 0 , 0 ) . The solution of system (4) is bounded, so the solution cannot extend to the infinity. Additionally, there exists no limit cycle, since system (4) has no interior equilibrium. Thus, E 0 ( 0 , 0 ) is globally asymptotically stable. □
Theorem 3.
If a f < e < a f + c , E 1 ( x 1 , 0 ) is globally asymptotically stable.
The proof of Theorem 3 is similar to that of Theorem 2; thus, we omit it for brevity.
Theorem 4.
If e > a f + c , E 2 ( x * , y * ) is globally asymptotically stable.
Proof. 
Here, we consider the Dulac function B ( x , y ) = 1 x y . Then,
( B P ) x + ( B Q ) y = c ( a + y ) y < 0 ,
where
P = e x a + y c x 2 a + y f x y x , Q = y + x y .
According to the Bendixson–Dulac discriminant, system (4) has no limit cycle in the first quadrant. Thus, E 2 ( x * , y * ) is globally asymptotically stable. The proof of Theorem 4 is finished. □
Notice that in Theorems 2–4, the conditions e < a f , a f < e < a f + c and e > a f + c are equivalent to r 1 < β α d , β α d < r 1 < β α d + b 1 r 2 a 2 and r 1 > β α d + b 1 r 2 a 2 , respectively. Thus the intrinsic growth rate of the prey in patch 1 plays an important role in the sustainable development of the species. In detail, when r 1 is not large, both the prey and the predator are led to extinction. Meanwhile, intermediate intrinsic growth rate r 1 can be favorable to the survival of the prey. However, an extremely large intrinsic growth rate could result in the coexistence of both the prey and the predator.

4. The Impact of Fear Effect and Dispersal

4.1. The Impact of Fear Effect on the Predator Density

Next, we will discuss the impact of the fear effect on the predator density. Computing the derivation along the predator y * with respect to k, we have
d y * d k = 2 a Δ ( 2 ( a f ) a 4 ( c e ) ) 4 k Δ = f 2 ( Δ a ) 2 4 k Δ .
Under the existence condition of E 2 , i.e., e > a f + c , one can obtain Δ a > f . Thus, we can easily deduce d y * d k < 0 . Additionally,
lim k d y * d k = 0
which shows that y * is a strictly decreasing function of k, i.e., increasing the fear effect k may decrease the predator density.

4.2. The Impact of Dispersal on the Predator Density

The derivation of y * with respect to m along system (4) is as follows.
d y * d m = ( 1 m ) · a Δ + ( a f ) a 4 ( c e ) 2 Δ · m · ( 1 m ) + f m ( Δ + a f ) 2 Δ · m · ( 1 m ) .
Under the existence condition of E 2 , i.e., e > a f + c , one can get Δ > a + f . Additionally, notice that: a Δ + ( a f ) a 4 ( c e ) > f ( Δ + a f ) > 0 . Based on the above analysis, we can get d y * d m > 0 . In other words, y * increases as m increases.
Figure 2a,b show how fear effect and dispersal have a direct influence on the predator density, respectively. In other words, the predator density decreases when k increases or m decreases. In fact, it is not difficult to understand the above observation biologically. When the fear effect k increases or m decreases, the prey in patch 1 is more likely to move to refuge patch 2, which will lead to a decease in the total prey density. Consequently, the predation density decreases due to the reduced amount of food.

5. Conclusions

In this work, a patchy predator–prey system with a fear effect has been proposed and analyzed. Only the prey population migrates at a constant rate between patches, and there are two time scales, i.e., a fast one for migration of prey between patches, and a slow one corresponding to local predator–prey interaction. The aggregated model is studied analytically, and the threshold conditions for the existence and stability of various steady states are worked out in Theorems 1–4. In detail, the fear effect does not change the stability of the unique positive equilibrium of the system, which is in agreement with [15]. However, unlike [15], dispersal has a vital influence on the species’ survival. In fact, when m < d r 1 + d , i.e., β > r 1 d α , both the prey and the predator species are led to extinction. When d r 1 + d < m < d r 1 + d b 1 r 2 a 2 , i.e., r 1 d α < β < r 1 b 1 r 2 a 2 d α , the prey species is permanent, but the predator species is led to extinction. When m > d r 1 + d b 1 r 2 a 2 , i.e., β > r 1 b 1 r 2 a 2 d α , both the prey and the predator can persist. From the above, whether the species can survive or not depends on the dispersal. The results show that the dispersal and fear effect play an important role in the dynamic behaviours of the system, which is a good extension and supplements those in [14,27]. We would also like to point out here that several potential directions can be investigated. For example, in this paper, we only discuss two patches for the prey. It is also interesting to propose a model with dispersal among n > 2 patches. Moreover, similar to [28,29,30], once there is noise or the conformable derivative in the system, respectively, how about the corresponding dynamic behaviour? We leave this for future work.

Author Contributions

Writing—original draft preparation, T.L.; writing—review and editing, L.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grant (11601085) and the Natural Science Foundation of Fujian Province (2021J01614, 2021J01613).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the anonymous reviewers for their valuable comments, which greatly improved the final expression of the paper.

Conflicts of Interest

The authors declare that there are no conflict of interest.

References

  1. Xiao, Z.W.; Li, Z. Stability and bifurcation in a stage-structured predator-prey model with Allee effect and time delay. IAENG Int. J. Appl. Math. 2019, 49, 6–13. [Google Scholar]
  2. McNair, J.N. The effects of refuges on predator-prey interactions: A reconsideration. Theor. Popul. Biol. 1986, 29, 38–63. [Google Scholar] [CrossRef]
  3. Briggs, C.J.; Hoopes, M.F. Stabilizing effects in spatial parasitoid-host and predator-prey models: A review. Theor. Popul. Biol. 2004, 65, 299–315. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, F.D.; Chen, L.J.; Xie, X.D. On a Leslie-Gower predator-prey model incorporating a prey refuge. Nonlinear Anal. Real World Appl. 2009, 10, 2905–2908. [Google Scholar] [CrossRef]
  5. Lima, S.L. Nonlethal Effects in the Ecology of Predator-Prey Interactions. BioScience 1998, 48, 25–34. [Google Scholar] [CrossRef] [Green Version]
  6. Beretta, E.; Solimano, F.; Takeuchi, Y. Global stability and periodic orbits for two-patch predator-prey diffusion-delay models. Math. Biosci. 1987, 85, 153–183. [Google Scholar] [CrossRef]
  7. Kuang, Y.; Takeuchi, Y. Predator-prey dynamics in models of prey dispersal in two-patch environments. Math. Biosci. 1994, 120, 77–98. [Google Scholar] [CrossRef]
  8. Aly, S.; Farkas, M. Bifurcations in a predator–prey model in patchy environment with diffusion. Nonlinear Anal. Real World Appl. 2004, 5, 519–526. [Google Scholar] [CrossRef]
  9. Souna, F.; Djilali, S.; Lakmeche, A. Spatiotemporal behavior in a predator–prey model with herd behavior and cross-diffusion and fear effect. Eur. Phys. J. Plus 2021, 136, 1–21. [Google Scholar] [CrossRef]
  10. Michalski, J.; Poggiale, J.-C.; Arditi, R.; Auger, P. Macroscopic Dynamic Effects of Migrations in Patchy Predator-prey Systems. J. Theor. Biol. 1997, 185, 459–474. [Google Scholar] [CrossRef]
  11. El Abdllaoui, A.; Auger, P.; Kooi, B.W.; de la Parra, R.B.; Mchich, R. Effects of density-dependent migrations on stability of a two-patch predator–prey model. Math. Biosci. 2007, 210, 335–354. [Google Scholar] [CrossRef] [PubMed]
  12. Ma, Z.; Wang, S.; Li, W.; Li, Z. The effect of prey refuge in a patchy predator–prey system. Math. Biosci. 2013, 243, 126–130. [Google Scholar] [CrossRef] [PubMed]
  13. Rani, R.; Gakkhar, S.; Moussaoui, A. Dynamics of a fishery system in a patchy environment with nonlinear harvesting. Math. Methods Appl. Sci. 2019, 42, 7192–7209. [Google Scholar] [CrossRef]
  14. Liu, Z.G. The effect of prey refuge in a patchy Leslie-Gower predation system. Acta Ecol. Sin. 2018, 38, 2958–2964. [Google Scholar]
  15. Wang, X.; Zanette, L.; Zou, X. Modelling the fear effect in predator–prey interactions. J. Math. Biol. 2016, 73, 1179–1204. [Google Scholar] [CrossRef]
  16. Wang, X.; Zou, X. Modeling the Fear Effect in Predator–Prey Interactions with Adaptive Avoidance of Predators. Bull. Math. Biol. 2017, 79, 1325–1359. [Google Scholar] [CrossRef]
  17. Sasmal, S. Population dynamics with multiple Allee effects induced by fear factors-a mathematical study on prey-predator. Appl. Math. Model. 2018, 64, 1–14. [Google Scholar] [CrossRef]
  18. Xiao, Z.W.; Li, Z. Stability analysis of a mutual interference predator-prey model with the fear effect. J. Appl. Sci. Eng. 2019, 6, 205–211. [Google Scholar]
  19. Kumar, V.; Kumar, N. Stability and bifurcation analysis of Hassell-Varley prey-predator system with fear effect. Appl. Math. Comput. 2020, 22, 150–170. [Google Scholar] [CrossRef]
  20. Das, M.; Samanta, G.P. A prey-predator fractional order model with fear effect and group defense. Int. J. Dyn. Control 2020, 9, 334–349. [Google Scholar] [CrossRef]
  21. Liu, T.T.; Chen, L.J.; Chen, F.D.; Li, Z. Stability analysis of a Leslie-Gower model with strong Allee effect on prey and fear effect on predator. Int. J. Bifurc. Chaos 2022, 32, 2250082. [Google Scholar] [CrossRef]
  22. Biktashev, V.N.; Brindley, J.; Holden, A.V.; Tsyganov, M.A. Pursuit-evasion predator-prey waves in two spatial dimensions. Chaos: Interdiscip. J. Nonlinear Sci. 2004, 14, 988–994. [Google Scholar] [CrossRef] [Green Version]
  23. Bogoliubov, N.N.; Mitropolsky, Y.A. Asymptotic Methods in the Theory of Nonlinear Oscillations, Gordon and Breach; CRC Press: New York, NY, USA, 1961. [Google Scholar]
  24. Iwasa, Y.; Andreasen, V.; Levin, S. Aggregation in model ecosystems. I. Perfect aggregation. Ecol. Model. 1987, 37, 287–302. [Google Scholar] [CrossRef]
  25. Iwasa, Y.; Levin, S.; Andreasen, V. Aggregation in Model Ecosystems II. Approximate Aggregation. Math. Med. Biol. A J. IMA 1989, 6, 1–23. [Google Scholar] [CrossRef]
  26. Zhang, Z.F.; Ding, T.R.; Huang, W.Z.; Dong, Z.X. Qualitative Theory of Differential Equation; Science Press: Beijing, China, 1992; Volume 101. (In Chinese) [Google Scholar]
  27. Zhang, H.; Cai, Y.; Fu, S.; Wang, W. Impact of the fear effect in a prey-predator model incorporating a prey refuge. Appl. Math. Comput. 2019, 356, 328–337. [Google Scholar] [CrossRef]
  28. Waezizadeh, T.; Mehrpooya, A. A stochastic model for dynamics of two populations and its stability. In Proceedings of the 47th Annual Iranlan Mathematics Conference, Kharazmi University, Karaj, Iran, 28–31 August 2016. [Google Scholar]
  29. Gokila, C.; Sambath, M.; Balachandran, K.; Ma, Y.-K. Analysis of Stochastic Predator-Prey Model with Disease in the Prey and Holling Type II Functional Response. Adv. Math. Phys. 2020, 2020, 3632091. [Google Scholar] [CrossRef]
  30. Hosseini1, K.; Sadri1, K.; Mirzazadeh, M.; Salahshour, S.; Park, C.; Lee, J.R. The Guava Model Involving the Conformable Derivative and its Mathematical Analysis. Fractals 2022. [Google Scholar] [CrossRef]
Figure 1. The phase portraits of system (4).
Figure 1. The phase portraits of system (4).
Axioms 11 00577 g001
Figure 2. The total predator density y * with respect to k and m.
Figure 2. The total predator density y * with respect to k and m.
Axioms 11 00577 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, T.; Chen, L. Stability Analysis of a Patchy Predator–Prey Model with Fear Effect. Axioms 2022, 11, 577. https://doi.org/10.3390/axioms11100577

AMA Style

Liu T, Chen L. Stability Analysis of a Patchy Predator–Prey Model with Fear Effect. Axioms. 2022; 11(10):577. https://doi.org/10.3390/axioms11100577

Chicago/Turabian Style

Liu, Tingting, and Lijuan Chen. 2022. "Stability Analysis of a Patchy Predator–Prey Model with Fear Effect" Axioms 11, no. 10: 577. https://doi.org/10.3390/axioms11100577

APA Style

Liu, T., & Chen, L. (2022). Stability Analysis of a Patchy Predator–Prey Model with Fear Effect. Axioms, 11(10), 577. https://doi.org/10.3390/axioms11100577

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop