Processing math: 21%
Home Complex dynamics of a nonlinear discrete predator-prey system with Allee effect
Article Open Access

Complex dynamics of a nonlinear discrete predator-prey system with Allee effect

  • Jing Wang and Ceyu Lei EMAIL logo
Published/Copyright: June 20, 2024

Abstract

The transition between strong and weak Allee effects in prey provides a simple regime shift in ecology. In this article, we study a discrete predator-prey system with Holling type II functional response and Allee effect. First, the number of fixed points of the system, local stability, and global stability is discussed. The population changes of predator and prey under strong or weak Allee effects are proved using the nullclines and direction field, respectively. Second, using the bifurcation theory, the bifurcation conditions for the system to undergo transcritical bifurcation and Neimark-Sacker bifurcation at the equilibrium point are obtained. Finally, the dynamic behavior of the system is analyzed by numerical simulation of bifurcation diagram, phase diagram, and maximum Lyapunov exponent diagram. The results show that the system will produce complex dynamic phenomena such as periodic state, quasi-periodic state, and chaos.

MSC 2010: 35K57; 37N25; 39A30; 92D25

1 Introduction

Population refers to a collection of the same kind of individuals living in a certain area at the same time. In nature, each population belongs to a certain level, and often, there are both predators at the upper level and preys at the lower level. As a result, in order to grasp the link between populations at all levels, assess and predict population persistence and extinction, a population model to explain the population system’s evolution law must be established.

The general population model can be divided into two categories: differential equation model and difference equation model. In general, when the population number is relatively large or generations overlap, it can be described by differential equation. The difference equation is suitable for describing the population with long life cycle, few population, and non-overlapping generations. For example, herbs often flower in their first year and then die, roots and all, after setting seed, and a small fraction ( <1% ) of the 22,000 species of teleost fish are semelparous and die soon after spawning (see [1] for more examples). In 1976, May [2] showed that although the first-order difference equation model is simple, it can also show a series of surprising dynamic behaviors, such as from stable point to a bifurcation hierarchy of stable cycles, and finally produce chaos. However, in a continuous-time model, a minimum of three species are needed for exhibiting chaos [3]. It can be seen that the discrete model will show more interesting dynamic behavior. Therefore, the discrete dynamic system model has attracted the attention of many scholars [49]. AlSharawi et al. [10] considered the influence of vigilance of prey on dynamics of a discrete-time predator-prey system. They studied the stability, persistence, flip bifurcation, and Neimark-Sacker bifurcation of the discrete system. The results show that with the increase of prey vigilance, the density of predator population continuously decreases and high vigilance will have a detrimental role for the prey population. Streipert et al. [11] deduced a discrete predator-prey model through the first principles in economics. They extended standard phase plane analysis by introducing the next iterate root-curve associated with the nontrivial prey nullcline. The stability of the system is proved by combining the curve with the nullclines and direction field. Finally, it is proved that the system will have a transcritical bifurcation at the boundary equilibrium point and a Neimark-Sacker bifurcation at the internal equilibrium point.

The Allee effect, a reduction of the per capita growth rate of a population of biological species at densities smaller than a critical value, was first introduced by Allee in 1931 [12]. Allee effects are mainly classified into two ways: strong and weak Allee effects. There is a critical value for the strong Allee effect. When the population density is less than this critical value, the population shows a negative growth trend, and then, the population will become extinct. When the population density is higher than this critical value, the population growth is positive, and then, it can develop. The weak Allee effect does not have a critical threshold. Its characteristic is that the individual growth rate at low density is lower than that at high density and always keeps a positive value, so the two populations can coexist permanently. Empirical evidence of Allee effects has been observed in many natural species, for example, plants [13,14], insects [15], marine invertebrates [16], and birds and mammals [17], etc. Therefore, many researchers have studied the bifurcation and stability analysis for discrete-time predator-prey system with the Allee effect [7,1821]. However, in these articles, they did not consider the impact of strong and weak Allee effects on the local stability of the equilibrium point of the predation system, nor did they prove the global stability of the equilibrium point. So, in this article, we consider a discrete-time predator-prey model with Holling type II functional response and Allee effect in preys, which is given by

(1) {xn+1=xnexp[r(Kxn)(xnA)βyna+cxn],yn+1=ynexp[αβxna+cxnd],

where x and y represent the densities of prey and predator population, respectively, and f(x)=βxa+cx is the generalized Holling type II functional response. In System (1), r denotes the intrinsic growth rate of the prey population, K denotes the environmental carrying capacity of the prey, d is the predator’s mortality rate, 0<α<1 represents the conversion efficiency of intake of prey into new predators, β denotes the maximal predator capita consumption rate, a represents the half-saturation constant, and c denotes the handling time. The parameters r , K , β , c , and d are assumed to be positive. The parameter A represents the threshold of multiplicative Allee effect. When 0<A<K , it indicates the the strong Allee effect; when K<A<0 , it indicates the weak Allee effect. Note that A=0 is the transition between the weak and strong Allee effects.

In this article, we study a discrete predator-prey system with Holling type-II functional response and Allee effect. There are several highlights in our analysis: (a) the influence of Allee effect parameter A on the local stability of equilibrium point is analyzed; (b) the conditions of global stability of discrete System (1) at the positive equilibrium point are obtained; (c) the population changes of predator and prey under strong or weak Allee effect are proved using the nullclines and direction field, respectively; and (d) the transcritical bifurcation and Neimark-Sacker bifurcation of System (1) are completely and rigorously analyzed.

This article is organized as follows. In Section 2, we analyze the dynamics of System (1), including the local stability of the equilibrium points, the global stability of the positive equilibrium point, and the bifurcation analysis at the equilibrium point. In Section 3, we verify our analytical results through numerical simulations. In Section 4, this article is ended with a brief conclusion.

2 Model dynamics

Lemma 1

Solutions of System (1) with nonnegative initial conditions remain nonnegative. If x0=0 , then xn=0 for all n0 . If y0=0 , then yn=0 for all n0 . If x0>0 and y00 , then xn>0 for all n0 . If x00 and y0>0 , then yn>0 for all n0 .

Proof

It can be directly proved by the model structure.□

Lemma 2

  1. System (1) always has a trivial equilibrium point E0(0,0) .

  2. System (1) always has a positive semi-trivial equilibrium point EK(K,0) .

  3. System (1) always has a semi-trivial equilibrium point EA(A,0) .

  4. Under the strong Allee effect, if 0<A<adαβcd<K , then System (1) has a positive nontrivial equilibrium point E*(x*,y*) . Under the weak Allee effect, if 0<adαβcd<K , then System (1) has a positive nontrivial equilibrium point E*(x*,y*) , where x*=adαβcd and y*=rβ(Kx*)(x*A)(a+cx*) .

Proof

Direct computation.□

Lemma 3

Consider System (1) with initial conditions (x0,y0) . If x0=0 and y00 , then (xn,yn) converges to E0 .

Proof

Since xn=0 for all n0 if x0=0 . In this case, yn+1=1edny0 . This converges to zero for d>0 . This completes the proof.□

The linearized form of System (1) is then

(2) {xt+1=fxxt+fyyt,yt+1=gxxt+gyyt,

which has the Jacobian matrix

(3) J=[fxfygxgy](x,y),

where (x,y) denotes the fixed points of System (1).

To obtain the stable results of (2), we consider the algebraic equation

(4) λ2Tλ+D=0,

where

T=TrJ=fx+gyandD=DetJ=fxgyfygx.

It is well-known that it has two roots of the form

λ1,2=(T±T24D)/2.

Lemma 4

[22] The equilibrium point (x,y) is called

  1. sink (locally asymptotically stable) if λ1<1 and λ2<1 ;

  2. source (locally unstable) if λ1>1 and λ2>1 ;

  3. saddle if λ1>1 and λ2<1 (or λ1<1 and λ2>1 );

  4. non-hyperbolic if λ1=1 or λ2=1 .

Lemma 5

[23] Let E*(x*,y*) be the unique positive equilibrium point of System (1), then the following propositions hold:

  1. It is a sink if

    D<1andT<D+1.

  2. It is a source if

    D>1andT<D+1,orT>D+1.

  3. It is a saddle if

    0<T+D+1<2T.

  4. It is non-hyperbolic if

    T=D+1,orD=1andT2.

Jacobian matrix can be evaluated at E0(0,0) as

(5) J(E0)=(erAK00ed).

The eigenvalues of the Jacobian are λ1=ed and λ2=erAK at trivial equilibrium point E0(0,0) . The results regarding dynamical behaviors are listed in Table 1.

Table 1

Properties of origin equilibrium point E0(0,0)

Conditions Eigenvalues Properties
λ 1 = e d λ 2 = e r A K
d > 0 0 < A < K λ 1 < 1 λ 2 < 1 Sink
K < A < 0 λ 2 > 1 Saddle
A = 0 λ 2 = 1 Non-hyperbolic

From Table 1, we can obtain the following theorem.

Theorem 1

E 0 ( 0 , 0 ) is always locally asymptotically stable under strong Allee effect, while it is always unstable under weak Allee effect.

The Jacobian matrix computed at E K ( K , 0 ) is

(6) J ( E K ) = 1 + r A K r K 2 β K a + c K 0 exp α β K a + c K d .

The eigenvalues of the Jacobian are λ 1 = 1 + r A K r K 2 and λ 2 = exp [ α β K a + c K d ] . The properties of semi-trivial equilibrium point E K ( K , 0 ) are summarized in Table 2.

Table 2

Properties of semi-trivial equilibrium point E K ( K , 0 )

Conditions Eigenvalues Properties
λ 1 = 1 + r A K r K 2 λ 2 = exp [ α β K a + c K d ]
0 < r K ( K A ) < 2 ( α β c d ) K = a d λ 1 < 1 λ 2 = 1 Non-hyperbolic
( α β c d ) K < a d λ 2 < 1 Sink
( α β c d ) K > a d λ 2 > 1 Saddle
r K ( K A ) > 2 ( α β c d ) K = a d λ 1 > 1 λ 2 = 1 Non-hyperbolic
( α β c d ) K < a d λ 2 < 1 Saddle
( α β c d ) K > a d λ 2 > 1 Source
r K ( K A ) = 2 ( α β c d ) K = a d λ 1 = 1 λ 2 = 1 Non-hyperbolic
( α β c d ) K < a d λ 2 < 1 Non-hyperbolic
( α β c d ) K > a d λ 2 > 1 Non-hyperbolic

From Table 2, we can obtain the following theorem.

Theorem 2

When 0 < r K ( K A ) < 2 and ( α β c d ) K < a d are satisfied, the semi-trivial equilibrium point E K ( K , 0 ) is locally asymptotically stable.

Jacobian matrix can be evaluated at E A ( A , 0 ) as

(7) J ( E 2 ) = 1 + r A K r A 2 β A a + c A 0 exp α β A a + c A d .

The eigenvalues of the Jacobian are λ 1 = 1 + r A K r A 2 and λ 2 = exp [ α β A a + c A d ] at the semi-trivial equilibrium point E A ( A , 0 ) . The results regarding dynamical behaviors are listed in Table 3.

Table 3

Properties of semi-trivial equilibrium point E A ( A , 0 )

Conditions Eigenvalues Properties
λ 1 = 1 + r A K r A 2 λ 2 = exp [ α β A a + c A d ]
Strong Allee effect: 0 < A < K
0 < A < K α β A a + c A = d λ 1 > 1 λ 2 = 1 Non-hyperbolic
α β A a + c A < d λ 2 < 1 Saddle
α β A a + c A > d λ 2 > 1 Source
Weak Allee effect: K < A < 0
0 < r A ( A K ) < 2 α β A a + c A = d λ 1 < 1 λ 2 = 1 Non-hyperbolic
α β A a + c A < d λ 2 < 1 Sink
α β A a + c A > d λ 2 > 1 Saddle
r A ( A K ) > 2 α β A a + c A = d λ 1 > 1 λ 2 = 1 Non-hyperbolic
α β A a + c A < d λ 2 < 1 Saddle
α β A a + c A > d λ 2 > 1 Source

From Table 3, we can obtain the following theorem.

Theorem 3

The semi-trivial equilibrium point E A ( A , 0 ) is always unstable under the strong Allee effect. Under the weak Allee effect, when 0 < r A ( A K ) < 2 and α β A a + c A < d are satisfied, it is locally asymptotically stable.

J ( x , y ) evaluated at the positive equilibrium point E * ( x * , y * ) is

(8) J ( E * ) = f x * f y * g x * g y * = 1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 β x * a + c x * a α β y * ( a + c x * ) 2 1 .

Then characteristic equation of J ( E * ) is given by

(9) λ 2 Tr J λ + Det J = 0 ,

where

T = Tr J = 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 , D = Det J = 1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 .

In summary, we have the following theorem.

Theorem 4

System (1) at the positive equilibrium point E * ( x * , y * ) is local asymptotically stable when the conditions

4 + 2 r x * ( K + A ) 4 r x * 2 + 2 c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 > 0

and

r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 < 0

hold.

Proof

It can be proven by Lemma 5.□

Theorem 5

The positive equilibrium point E * ( x * , y * ) of System (1) is a global attractor if 1 + r ( K + A ) x * 2 r x * 2 + c β x * y * ( a + c x * ) 2 0 .

Proof

Let f ( x , y ) = x exp [ r ( K x ) ( x A ) β y a + c x ] , g ( x , y ) = y exp [ α β x a + c x d ] . Then, f ( x , y ) is non-decreasing in x if 1 + r ( K + A ) x * 2 r x * 2 + c β x * y * ( a + c x * ) 2 0 and non-increasing in y for all ( x , y ) . Moreover, g ( x , y ) is non-decreasing in both arguments x and y for all ( x , y ) . Let ( m 1 , M 1 , m 2 , M 2 ) be a positive solution of the system

m 1 = f ( m 1 , M 2 ) , M 1 = f ( M 1 , m 2 ) ,

m 2 = g ( m 1 , m 2 ) , M 2 = g ( M 1 , M 2 ) .

Then, we have

m 1 = m 1 exp r ( K m 1 ) ( m 1 A ) β M 2 a + c m 1 ,

M 1 = M 1 exp r ( K M 1 ) ( M 1 A ) β m 2 a + c M 1 ,

and

m 2 = m 2 exp α β m 1 a + c m 1 d , M 2 = M 2 exp α β M 1 a + c M 1 d .

Through calculation, we can obtain

m 1 = M 1 , m 2 = M 2 .

Hence, the unique positive equilibrium point E * ( x * , y * ) is a global attractor if 1 + r ( K + A ) x * 2 r x * 2 + c β x * y * ( a + c x * ) 2 0 .□

From Theorems 4 and 5, we can obtain the following conclusions.

Theorem 6

Assume that 1 + r ( K + A ) x * 2 r x * 2 + c β x * y * ( a + c x * ) 2 0 , D < 1 , and T < D + 1 . Then, the unique positive equilibrium point E * ( x * , y * ) of System (1) is globally asymptotically stable.

In fact, x n is increasing if the forward operator, x n = x n + 1 x n , is positive. For System (1), the forward operators are

(10) x n = x n + 1 x n = x n exp r ( K x n ) ( x n A ) β y n a + c x n 1

and

(11) y n = y n + 1 y n = y n exp α β x n a + c x n d 1 .

From (10), we see that x n > 0 is increasing if and only if y n < r β ( K x n ) ( x n A ) ( a + c x n ) , since x n + 1 > 0 for x n > 0 . Thus, if ( x n , y n ) is above the line

(12) y = l ( x ) r β ( K x ) ( x A ) ( a + c x ) ,

y n > l ( x n ) , then the sequence of iterates, x n is decreasing, and if y n < l ( x n ) , then x n is increasing. Similarly, from (11), it follows that y n is increasing as long as x n > a d α β c d and decreasing if x n < a d α β c d .

By Lemma 1, solutions remain in the first quadrant for all nonnegative initial conditions. We divide the first quadrant into regions based on the component-wise monotonicity obtained by solving x n + 1 = x n and y n + 1 = y n . The curves along which (1) satisfying x n = 0 and y n = 0 are given by the line

(13) y n = l ( x n ) = r β ( K x n ) ( x n A ) ( a + c x n ) and x n = a d α β c d ,

respectively. We refer to these curves as nullclines. These two curves divide the first quadrant into three or four regions R i ( i = 1 , 2 , 3 , 4 ) .

We define the regions

R 1 ( x n , y n ) ( 0 , + ) 2 : x n > a d α β c d and y n l ( x n ) , R 2 ( x n , y n ) ( 0 , + ) 2 : x n a d α β c d and y n > l ( x n ) , R 3 ( x n , y n ) ( 0 , + ) 2 : x n < a d α β c d and y n l ( x n ) , R 4 ( x n , y n ) ( 0 , + ) 2 : x n a d α β c d and y n < l ( x n ) .

In Figure 1, when 0 < A < K , for a d α β c d > K or 0 < a d α β c d < A , R 4 = , and R 3 = { ( x n , y n ) ( 0 , + ) 2 : A < x n < K and y n l ( x n ) } . It can be seen from Figure 1 that only when A < a d α β c d < K , System (1) has a unique positive equilibrium point. By observing the arrow direction in the image, we can find that when the number of prey population is small, the prey will be extinct and the predators will also be extinct. This is due to x n + 1 x n < 0 and y n + 1 y n < 0 , i.e., x n + 1 < x n and y n + 1 < y n . Thus, the trivial equilibrium E 0 is stable under the strong Allee effect, which is consistent with the conclusion of Theorem 1.

Figure 1 
               Phase diagram of System (1) when 
                     
                        
                        
                           0
                           <
                           A
                           <
                           K
                        
                        0\lt A\lt K
                     
                  . The red dashed line corresponds to the predator nullclines, and the blue curves correspond to the prey nullclines. A horizontal arrow pointing to the right (left) represents 
                     
                        
                        
                           
                              
                                 x
                              
                              
                                 n
                                 +
                                 1
                              
                           
                           ‒
                           
                              
                                 x
                              
                              
                                 n
                              
                           
                           >
                           0
                           
                              (
                              
                                 <
                                 0
                              
                              )
                           
                        
                        {x}_{n+1}&#x2012;{x}_{n}\gt 0\left(\lt 0)
                     
                  , and a vertical arrow pointing up (down) represents 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 n
                                 +
                                 1
                              
                           
                           ‒
                           
                              
                                 y
                              
                              
                                 n
                              
                           
                           >
                           0
                           
                              (
                              
                                 <
                                 0
                              
                              )
                           
                        
                        {y}_{n+1}&#x2012;{y}_{n}\gt 0\left(\lt 0)
                     
                  . Subfigure (a) is a schematic image if 
                     
                        
                        
                           K
                           <
                           
                              
                                 a
                                 d
                              
                              
                                 α
                                 β
                                 ‒
                                 c
                                 d
                              
                           
                        
                        K\lt \frac{ad}{\alpha \beta &#x2012;cd}
                     
                  , (b) is a schematic image if 
                     
                        
                        
                           0
                           <
                           
                              
                                 a
                                 d
                              
                              
                                 α
                                 β
                                 ‒
                                 c
                                 d
                              
                           
                           <
                           A
                        
                        0\lt \frac{ad}{\alpha \beta &#x2012;cd}\lt A
                     
                  , while (c) shows the case when 
                     
                        
                        
                           A
                           <
                           
                              
                                 a
                                 d
                              
                              
                                 α
                                 β
                                 ‒
                                 c
                                 d
                              
                           
                           <
                           K
                        
                        A\lt \frac{ad}{\alpha \beta &#x2012;cd}\lt K
                     
                  . Note that in (a) and (b), there is no region 
                     
                        
                        
                           
                              
                                 R
                              
                              
                                 4
                              
                           
                        
                        {R}_{4}
                     
                   in the first quadrant.
Figure 1

Phase diagram of System (1) when 0 < A < K . The red dashed line corresponds to the predator nullclines, and the blue curves correspond to the prey nullclines. A horizontal arrow pointing to the right (left) represents x n + 1 x n > 0 ( < 0 ) , and a vertical arrow pointing up (down) represents y n + 1 y n > 0 ( < 0 ) . Subfigure (a) is a schematic image if K < a d α β c d , (b) is a schematic image if 0 < a d α β c d < A , while (c) shows the case when A < a d α β c d < K . Note that in (a) and (b), there is no region R 4 in the first quadrant.

In Figure 2, when K < A < 0 , for a d α β c d > K , R 4 = and R 3 = { ( x n , y n ) ( 0 , + ) 2 : 0 < x n < K and y n l ( x n ) } . Similarly, we can find that when the prey population is small, under the influence of the weak Allee effect, the prey population will not be extinct, and the population will increase. However, the number of predator population will gradually decrease and become extinct. This is different from the strong Allee effect. And by observing Figures 1 and 2, it can be found that when x > K , the growth rate of the prey population shows a negative growth, which is consistent with the biological significance.

Figure 2 
               Phase diagram of System (1) when 
                     
                        
                        
                           ‒
                           K
                           <
                           A
                           <
                           0
                        
                        &#x2012;K\lt A\lt 0
                     
                  . The red dashed line corresponds to the predator nullclines, and the blue curves correspond to the prey nullclines. A horizontal arrow pointing to the right (left) represents 
                     
                        
                        
                           
                              
                                 x
                              
                              
                                 n
                                 +
                                 1
                              
                           
                           ‒
                           
                              
                                 x
                              
                              
                                 n
                              
                           
                           >
                           0
                           
                              (
                              
                                 <
                                 0
                              
                              )
                           
                        
                        {x}_{n+1}&#x2012;{x}_{n}\gt 0\left(\lt 0)
                     
                  , and a vertical arrow pointing up (down) represents 
                     
                        
                        
                           
                              
                                 y
                              
                              
                                 n
                                 +
                                 1
                              
                           
                           ‒
                           
                              
                                 y
                              
                              
                                 n
                              
                           
                           >
                           0
                           
                              (
                              
                                 <
                                 0
                              
                              )
                           
                        
                        {y}_{n+1}&#x2012;{y}_{n}\gt 0\left(\lt 0)
                     
                  . Subfigure (a) is a schematic image if 
                     
                        
                        
                           K
                           <
                           
                              
                                 a
                                 d
                              
                              
                                 α
                                 β
                                 ‒
                                 c
                                 d
                              
                           
                        
                        K\lt \frac{ad}{\alpha \beta &#x2012;cd}
                     
                  , while (b) shows the case when 
                     
                        
                        
                           0
                           <
                           
                              
                                 a
                                 d
                              
                              
                                 α
                                 β
                                 ‒
                                 c
                                 d
                              
                           
                           <
                           K
                        
                        0\lt \frac{ad}{\alpha \beta &#x2012;cd}\lt K
                     
                  . Note that in (a), there is no region 
                     
                        
                        
                           
                              
                                 R
                              
                              
                                 4
                              
                           
                        
                        {R}_{4}
                     
                   in the first quadrant.
Figure 2

Phase diagram of System (1) when K < A < 0 . The red dashed line corresponds to the predator nullclines, and the blue curves correspond to the prey nullclines. A horizontal arrow pointing to the right (left) represents x n + 1 x n > 0 ( < 0 ) , and a vertical arrow pointing up (down) represents y n + 1 y n > 0 ( < 0 ) . Subfigure (a) is a schematic image if K < a d α β c d , while (b) shows the case when 0 < a d α β c d < K . Note that in (a), there is no region R 4 in the first quadrant.

Theorem 7

If A = 0 , then E 0 = E A and System (1) undergoes a transcritical bifurcation at the trivial equilibrium point E 0 .

Proof

For A = 0 , the equilibria E 0 and E A coalesce. The Jacobian evaluated at E 0 given in (5) has eigenvalues λ 1 = e d and λ 2 = 1 , the trivial equilibrium point E 0 is non-hyperbolic. It can be concluded that the central manifold of the map x x exp [ r x ( K x ) ] . According to the literature [22], it manifests that System (1) exhibits a transcritical bifurcation at E 0 .□

Theorem 8

If a d α β c d = K , then E K = E * and System (1) undergoes a transcritical bifurcation at the semi-trivial equilibrium point E K .

Proof

For a d α β c d = K , the equilibria E K and E * coalesce. The Jacobian evaluated at E K given in (5) has eigenvalues λ 1 = 1 + r A K r A 2 and λ 2 = 1 , and the trivial equilibrium point E K is non-hyperbolic. Similarly, according to the literature [23], it manifests that System (1) exhibits a transcritical bifurcation at E K .□

Theorem 9

If 0 < a d α β c d = A , then E A = E * and System (1) undergoes a transcritical bifurcation at the semi-trivial equilibrium point E A .

Proof

For a d α β c d = A , the equilibria E A and E * coalesce. The Jacobian evaluated at E A given in (5) has eigenvalues λ 1 = 1 + r A K r K 2 and λ 2 = 1 , and the trivial equilibrium point E A is non-hyperbolic. Similarly, according to the literature [23], it manifests that System (1) exhibits a transcritical bifurcation at E A .□

Theorem 10

The interior equilibrium point E * loses its stability via Neimark-Sacker bifurcation if

(14) r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 = 0 .

Proof

Neimark-Sacker bifurcation occurs in the system when a pair of complex eigenvalues with unit modulus [24], i.e.,

Det ( J ) = 1 .

It is obtained

1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 = 1 .

This gives the condition for the Neimark-Sacker bifurcation.□

Now, we discuss the Neimark-Sacker bifurcation of the equilibrium point E * ( x * , y * ) . Here, we choose α as a bifurcation parameter. Neimark-Sacker bifurcation in a discrete system is the birth of a closed invariant curve from an equilibrium point. The bifurcation can be supercritical when closed invariant curve is stable and subcritical, when it is unstable.

The characteristic polynomial (8) of Jacobian matrix of linearized system of (1) about the positive equilibrium point ( x * , y * ) can be rewritten as:

(15) F ( λ ) = λ 2 p ( x * , y * ) λ + q ( x * , y * ) ,

where

p ( x * , y * ) = 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 , q ( x * , y * ) = 1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 .

Consider (15), then F ( λ ) = 0 has two complex conjugate roots with modulus one if the following conditions are satisfied:

1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α β 2 x * y * ( a + c x * ) 3 = 1

and

(16) T = 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 < 2 .

Let

Ω N S = { ( r , K , A , a , c , β , α , d ) R + 8 : ( 14 ) and ( 16 ) are satisfied } .

Then, the dynamic analysis of System (1) is analyzed when the parameters change in the small field of Ω N S . Select parameter ( r , K , A , a , c , β , α , d ) Ω N S , and consider the following system:

x n + 1 = x n exp r ( K x n ) ( x n A ) β y n a + c x n , y n + 1 = y n exp [ ( α + α * ) β x n a + c x n d ] ,

where α * 1 is a small perturbation parameter.

Let u = x x * , v = y y * , then we obtain

(17) u n + 1 = A 11 u + A 12 v + A 13 u 2 + A 14 u v + A 15 v 2 + A 16 u 3 + A 17 u 2 v + A 18 u v 2 + A 19 v 3 + O ( ( u + v ) 4 ) , v n + 1 = A 21 u + A 22 v + A 23 u 2 + A 24 u v + A 25 v 2 + A 26 u 3 + A 27 u 2 v + A 28 u v 2 + A 29 v 3 + O ( ( u + v ) 4 ) ,

where

A 11 = 1 + x * Φ 1 , A 12 = β x * a + c x * , Φ 1 = r ( K + A ) 2 r x * + c β y * ( a + c x * ) 2 , A 13 = 1 2 r ( K + A ) 2 r x * + c β y * ( a c x * ) 2 ( a + c x * ) 3 + 1 2 A 11 Φ 1 , A 14 = c β x * ( a + c x * ) 2 β A 11 a + c x * , A 15 = β 2 x * 2 ( a + c x * ) 2 , A 16 = 1 2 4 r + 2 c 3 β x * y * 4 a c 2 β y * ( a + c x * ) 3 + Φ 1 [ r ( K + A ) 4 r x * + c β y * ( a c x * ) ( a + c x * ) 3 ] A 11 ( 2 r + 2 c 2 β y * ( a + c x * ) 3 ) + 2 A 13 Φ 1 , A 17 = 1 2 β y * ( a c x * ) ( a + c x * ) 3 + c β x * Φ 1 ( a + c x * ) 2 + c β A 11 ( a + c x * ) 2 β A 13 a + c x * , A 18 = c β 2 x * 2 ( a + c x * ) 3 β A 14 2 ( a + c x * ) , A 19 = β 3 x * 6 ( a + c x * ) 3 , A 21 = a α β y * ( a + c x * ) 2 , A 22 = 1 , A 23 = a 2 α 2 β 2 y * 2 ( a + c x * ) 4 a c α β y * ( a + c x * ) 3 , A 24 = a α β ( a + c x * ) 2 , A 25 = A 28 = A 29 = 0 , A 26 = a c 2 α β y * ( a + c x * ) 4 + a α β A 23 3 ( a + c x * ) 2 2 a 2 c α 2 β 2 y * 3 ( a + c x * ) 5 , A 27 = a 2 α 2 β 2 2 ( a + c x * ) 4 a c α β ( a + c x * ) 3 .

The characteristic equation of System (17) at ( u , v ) = ( 0 , 0 ) is as follows:

λ 2 p ( α * ) λ + q ( α * ) = 0 ,

where

p ( α * ) = 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 , q ( α * ) = 1 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 + a α * β 2 x * y * ( a + c x * ) 3 .

Since parameters ( r , K , A , a , c , β , α , d ) Ω N S , the roots of the characteristic equation are

λ 1 , 2 = p ( α * ) 2 ± i 2 4 q ( α * ) p 2 ( α * ) ,

we have

λ 1 , 2 = q ( α * )

and

L = d λ 1 , 2 d α * α * = 0 0 .

In addition, it is required that α * = 0 , λ 1 , 2 j 1 ( j = 1 , 2 , 3 , 4 ) , which is equivalent to p ( 0 ) 2 , 0, 1, 2. Because ( r , K , A , a , c , β , α , d ) Ω N S , p ( 0 ) 2 , 2. We only require p ( 0 ) 0 , 1 , so that

(18) 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 0 and 2 + r x * ( K + A ) 2 r x * 2 + c β x * y * ( a + c x * ) 2 1 .

Therefore, eigenvalues λ 1 and λ 2 of the equilibrium point (0,0) of System (17) do not lay in the intersection of the unit circle with the coordinate axes when α * = 0 and the condition (18) holds.

Let η = p ( 0 ) 2 , ω = 4 q ( 0 ) p 2 ( 0 ) 2 , we use the following transformation:

u v = T W Z = A 12 0 η A 11 ω W Z ,

and System (17) becomes

(19) W t + 1 Z t + 1 = η ω ω η W t Z t + f ¯ ( W , Z ) g ¯ ( W , Z ) ,

where

f ¯ ( W , Z ) = A 13 A 12 u 2 + A 14 A 12 u v + A 15 A 12 v 2 + A 16 A 12 u 3 + A 17 A 12 u 2 v + A 18 A 12 u v 2 + A 19 A 12 v 3 + O ( ( W + Z ) 4 ) , g ¯ ( W , Z ) = A 13 ( η A 11 ) ω A 12 A 23 ω u 2 + A 14 ( η A 11 ) ω A 12 A 24 ω u v + A 15 ( η A 11 ) ω A 12 A 25 ω v 2 + A 16 ( η A 11 ) ω A 12 A 26 ω u 3 + A 17 ( η A 11 ) ω A 12 A 27 ω u 2 v + A 18 ( η A 11 ) ω A 12 A 28 ω u v 2 + A 19 ( η A 11 ) ω A 12 A 29 ω v 3 + O ( ( W + Z ) 4 ) , u = A 12 W , v = ( η A 11 ) W ω Z .

System (17) undergoes the Neimark-Sacker bifurcation if the following quantity is not zero:

(20) L = Re ( 1 2 λ 1 ) λ 2 2 1 λ 1 L 11 L 12 1 2 L 11 2 L 21 2 + Re ( λ 2 L 22 ) ,

where

L 11 = 1 4 [ ( f ¯ W W + f ¯ Z Z ) + i ( g ¯ W W + g ¯ Z Z ) ] , L 12 = 1 8 [ ( f ¯ W W f ¯ Z Z + 2 g ¯ W Z ) + i ( g ¯ W W g ¯ Z Z 2 f ˜ W Z ) ] , L 21 = 1 8 [ ( f ¯ W W f ¯ Z Z 2 g ¯ W Z ) + i ( g ¯ W W g ¯ Z Z + 2 f ¯ W Z ) ] , L 22 = 1 16 [ ( f ¯ W W W + f ¯ Z Z Z + g ¯ W W Z + g ¯ Z Z Z ) + i ( g ¯ W W W + g ¯ W Z Z f ¯ W W Z f ¯ Z Z Z ) ] .

If L 0 , the Neimark-Sacker bifurcation will occur in System (1), and the following theorem holds:

Theorem 11

System (1) undergoes a Neimark-Sacker bifurcation at the positive equilibrium point E * ( x * , y * ) if conditions in (18) are satisfied and L 0 in (20). Moreover, if L < 0 (resp., L > 0 ), an attracting (resp., repelling) invariant closed curve bifurcates from the steady state for α > α * (resp., α < α * ).

Remark 1

When the system has the Neimark-Sacker bifurcation, it will dispose of the expenditure invariant curve from the fixed point, which indicates that prey and predator can coexist, and the dynamic behavior can be periodic or quasi-periodic.

3 Numerical simulations

This section will show the bifurcation diagram, phase diagram, and maximum Lyapunov exponent diagram with the Allee effect model to verify the correctness of theoretical analysis.

3.1 Weak Allee effect

Assuming that the parameter is r = 1.7 , K = 1 , A = 0.1 , β = 0.7 , a = 1 , c = 1 , and d = 0.1 , α is the bifurcation parameter, and the initial value of the system is ( x 0 , y 0 ) = ( 0.8 , 0.5 ) . Lemma 2 shows that when α < 0.2857 , the coexistence equilibrium point does not exist; when α > 0.2857 , there is a unique coexistence equilibrium point. It can be seen from Theorem 8 and Figures 3(a) and (b) when α = 0.2857 , System (1) has a transcritical bifurcation at the boundary equilibrium point E K . And according to Theorem 10, when α = 0.3981 , System (1) will have a supercritical Neimark-Sacker bifurcation at the coexistence equilibrium point. The bifurcation diagram and the maximum Lyapunov exponent diagram are shown in Figure 3. Combined with Figures 3 and 4, when α < 0.3981 , the fixed point is stable. When α = 0.39 , the phase and evolution diagrams of predator and prey with time are given in Figures 4(b) and 5(a). It can be seen from Theorem 6 that the coexistence equilibrium points of System (1) are globally asymptotically stable. However, when α > 0.3981 , the coexistence equilibrium loses its stability, and a stable invariant loop appears. At this time, the periodic solution of System (1) appears (Figures 4(c) and (d) and 5(b)). When α increases, it can be seen from Figures 3(c), 4(e) and (f), and 5(c) that System (1) will have quasi-periodic solutions and chaos.

Figure 3 
                  Neimark-Sacker bifurcation diagram and maximum Lyapunov exponent diagram for System (1).
Figure 3

Neimark-Sacker bifurcation diagram and maximum Lyapunov exponent diagram for System (1).

Figure 4 
                  Phase diagram corresponding to Figure 3 when 
                        
                           
                           
                              α
                           
                           \alpha 
                        
                      takes different values: (a) 
                        
                           
                           
                              α
                              =
                              0.25
                           
                           \alpha =0.25
                        
                     , (b) 
                        
                           
                           
                              α
                              =
                              0.39
                           
                           \alpha =0.39
                        
                     , (c) 
                        
                           
                           
                              α
                              =
                              0.3981
                           
                           \alpha =0.3981
                        
                     , (d) 
                        
                           
                           
                              α
                              =
                              0.4
                           
                           \alpha =0.4
                        
                     , (e) 
                        
                           
                           
                              α
                              =
                              0.6
                           
                           \alpha =0.6
                        
                     , and (f) 
                        
                           
                           
                              α
                              =
                              0.9
                           
                           \alpha =0.9
                        
                     .
Figure 4

Phase diagram corresponding to Figure 3 when α takes different values: (a) α = 0.25 , (b) α = 0.39 , (c) α = 0.3981 , (d) α = 0.4 , (e) α = 0.6 , and (f) α = 0.9 .

Figure 5 
                  Evolution of predator and prey population with time: (a) 
                        
                           
                           
                              α
                              =
                              0.39
                           
                           \alpha =0.39
                        
                     , (b) 
                        
                           
                           
                              α
                              =
                              0.39
                           
                           \alpha =0.39
                        
                     , and (c) 
                        
                           
                           
                              α
                              =
                              0.6
                           
                           \alpha =0.6
                        
                     .
Figure 5

Evolution of predator and prey population with time: (a) α = 0.39 , (b) α = 0.39 , and (c) α = 0.6 .

3.2 Strong Allee effect

When A = 0.1 and other parameters remain unchanged, it can also be obtained that when α < 0.2857 , the coexistence equilibrium point does not exist; when α > 0.2857 , there is a unique coexistence equilibrium point, and when α = 0.2857 , System (1) has a transcritical bifurcation at the boundary equilibrium point E K . According to Theorem 10, when α = 0.3729 , System (1) has the Neimark-Sacker bifurcation at the coexistence equilibrium point. Figure 6(a) and (b) is the bifurcation graph of α on [0, 1], and Figure 6(c) is the maximum Lyapunov exponent graph corresponding to Figure 6(a). Figure 6(d–f) is a locally enlarged view of Figure 6(a–c). It can be seen from Figures 6 and 7 that when α < 0.3729 , the fixed point is stable; when α > 0.3729 , the fixed point loses its stability and a stable invariant loop appears. At this time, System (1) generates a periodic solution (Figures 7(d) and 8(b)). When α increases, System (1) produces quasi-periodic solutions and chaotic phenomena. However, when we continue to increase α , the population will become extinct (Figure 8(c)). Combining Theorem 6 and Figure 8(a), when 0.2875 < α < 0.3729 , the coexistence equilibrium point of System (1) exists and is globally asymptotically stable.

Figure 6 
                  Neimark-Sacker bifurcation diagram and maximum Lyapunov exponent diagram for System (1).
Figure 6

Neimark-Sacker bifurcation diagram and maximum Lyapunov exponent diagram for System (1).

Figure 7 
                  Phase diagram corresponding to Figure 6 when 
                        
                           
                           
                              α
                           
                           \alpha 
                        
                      takes different values: (a) 
                        
                           
                           
                              α
                              =
                              0.25
                           
                           \alpha =0.25
                        
                     , (b) 
                        
                           
                           
                              α
                              =
                              0.37
                           
                           \alpha =0.37
                        
                     , (c) 
                        
                           
                           
                              α
                              =
                              0.3729
                           
                           \alpha =0.3729
                        
                     , (d) 
                        
                           
                           
                              α
                              =
                              0.375
                           
                           \alpha =0.375
                        
                     , (e) 
                        
                           
                           
                              α
                              =
                              0.379
                           
                           \alpha =0.379
                        
                     , and (f) 
                        
                           
                           
                              α
                              =
                              0.4
                           
                           \alpha =0.4
                        
                     .
Figure 7

Phase diagram corresponding to Figure 6 when α takes different values: (a) α = 0.25 , (b) α = 0.37 , (c) α = 0.3729 , (d) α = 0.375 , (e) α = 0.379 , and (f) α = 0.4 .

Figure 8 
                  Evolution of predator and prey population with time. (a) 
                        
                           
                           
                              α
                              =
                              0.37
                           
                           \alpha =0.37
                        
                     , (b) 
                        
                           
                           
                              α
                              =
                              0.375
                           
                           \alpha =0.375
                        
                     , and (c) 
                        
                           
                           
                              α
                              =
                              0.4
                           
                           \alpha =0.4
                        
                     .
Figure 8

Evolution of predator and prey population with time. (a) α = 0.37 , (b) α = 0.375 , and (c) α = 0.4 .

It can be found from Figures 3 and 6 that under the weak Allee effect, the growth rate of the population at low density is always positive and the population will not become extinct. Under the strong Allee effect, the population will become extinct at low density. This is consistent with the conclusion obtained from Figures 1 and 2.

4 Conclusion

In this article, we study the stability and bifurcation of equilibrium points in a discrete predator-prey model with the Allee effect. The stability analysis shows that the model with strong Allee effect has two positive boundary equilibrium points, and the boundary equilibrium point corresponding to the Allee effect threshold is always unstable. Under the strong Allee effect, when 0 < A < a d α β c d < K , the discrete System (1) has a unique coexistence equilibrium, and when the population density is low, the predator and prey populations will be extinct. Under the weak Allee effect, when 0 < a d α β c d < K , the discrete System (1) has a unique coexistence equilibrium, and when the population density is low, the growth rate of the prey population always remains positive, and the population will not be extinct. The bifurcation analysis shows that when a d α β c d = K , the boundary equilibrium point E K will have a transcritical bifurcation, and when the coexistence equilibrium E * exists and loses stability, System (1) will have a Neimark-Sacker bifurcation. The numerical simulation reveals that when the energy conversion rate of predator α increases gradually, System (1) will produce periodic, quasi-periodic windows, and chaos.

  1. Funding information: This work was supported by Gansu Province University Teacher Innovation Fund (No. 2023A-178) and Graduate Research Support project of Northwest Normal University (No. 2023KYZZ-B051).

  2. Author contributions: JW: conceptualization, methodology, investigation, supervision, writing – review and editing, visualization. CL: conceptualization, methodology, software, formal analysis, investigation, writing – original draft.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

References

[1] M. Kot, Elements of Mathematical Ecology, Cambridge University Press, New York, 2001. 10.1017/CBO9780511608520Search in Google Scholar

[2] R. M. May, Simple mathematical models with very complicated dynamics, Nature 261 (1976), no. 5560, 459–467, DOI: https://doi.org/10.1038/261459a0. 10.1038/261459a0Search in Google Scholar PubMed

[3] A. Hastings and T. Powell, Chaos in three-species food chain, Ecology 72 (1991), no. 3, 896–903, DOI: https://doi.org/10.2307/1940591. 10.2307/1940591Search in Google Scholar

[4] A. Q. Khan and T. Khalique, Neimark-Sacker bifurcation and hybrid control in a discrete-time Lotka-Volterra model, Math. Methods Appl. Sci. 43 (2020), no. 9, 5887–5904, DOI: https://doi.org/10.1002/mma.6331. 10.1002/mma.6331Search in Google Scholar

[5] H. N. Agiza, E. M. Elabbasy, H. EL-Metwally, and A. A. Elsadany, Chaotic dynamics of a discrete prey-predator model with Holling type II, Nonlinear Anal. Real World Appl. 10 (2009), no. 1, 116–129, DOI: https://doi.org/10.1016/j.nonrwa.2007.08.029. 10.1016/j.nonrwa.2007.08.029Search in Google Scholar

[6] E. M. Elabbasy, A. A. Elsadany, and Y. Zhang, Bifurcation analysis and chaos in a discrete reduced Lorenz system, Appl. Math. Comput. 228 (2014), 184–194, DOI: https://doi.org/10.1016/j.amc.2013.11.088. 10.1016/j.amc.2013.11.088Search in Google Scholar

[7] L. Cheng and H. Cao, Bifurcation analysis of a discrete-time ratio-dependent predator-prey model with Allee effect, Commun. Nonlinear Sci. Numer. Simul. 38 (2016), 288–302, DOI: https://doi.org/10.1016/j.cnsns.2016.02.038. 10.1016/j.cnsns.2016.02.038Search in Google Scholar

[8] Q. Din, Complexity and chaos control in a discrete-time prey-predator model, Commun. Nonlinear Sci. Numer. Simul. 49 (2017), 113–134, DOI: https://doi.org/10.1016/j.cnsns.2017.01.025. 10.1016/j.cnsns.2017.01.025Search in Google Scholar

[9] C. Lei, X. Han, and W. Wang, Bifurcation analysis and chaos control of a discrete-time prey-predator model with fear factor, Math. Biosci. Eng. 19 (2022), no. 7, 6659–6679, DOI: https://doi.org/10.3934/mbe.2022313. 10.3934/mbe.2022313Search in Google Scholar PubMed

[10] Z. AlSharawi, N. Pal, and J. Chattopadhyay, The role of vigilance on a discrete-time predator-prey model, Discrete Contin. Dyn. Syst. Ser. B 27 (2022), no. 11, 6723–6744, DOI: https://doi.org/10.3934/dcdsb.2022017. 10.3934/dcdsb.2022017Search in Google Scholar

[11] S. H. Streipert, G. S. K. Wolkowicz, and M. Bohner, Derivation and analysis of a discrete predator-prey model, Bull. Math. Biol. 84 (2022), no. 7, 67, DOI: https://doi.org/10.1007/s11538-022-01016-4. 10.1007/s11538-022-01016-4Search in Google Scholar PubMed

[12] W. C. Allee, Animal Aggregations: A Study in General Sociology, University of Chicago Press, Chicago, 1931. 10.5962/bhl.title.7313Search in Google Scholar

[13] J.-B. Ferdy, F. Austerlitz, J. Moret, P.-H. Gouyon, and B. Godelle, Pollinator-induced density dependence in deceptive species, Oikos 87 (1999), no. 3, 549–560, DOI: https://doi.org/10.2307/3546819. 10.2307/3546819Search in Google Scholar

[14] M. J. Groom, Allee effects limit population viability of an annual plant, Am. Nat. 151 (1998), no. 6, 487–496, DOI: https://doi.org/10.1086/286135. 10.1086/286135Search in Google Scholar PubMed

[15] M. Kuussaari, I. Saccheri, M. Camara, and I. Hanski, Allee effect and population dynamics in the Glanville fritillary butterfly, Oikos 82 (1998), no. 2, 384–392, DOI: https://doi.org/10.2307/3546980. 10.2307/3546980Search in Google Scholar

[16] A. W. Stoner and M. Ray-Culp, Evidence for Allee effects in an over-harvested marine gastropod: density-dependent mating and egg production, Mar. Ecol. Prog. Ser. 202 (2000), 297–302, DOI: https://doi.org/10.3354/meps202297. 10.3354/meps202297Search in Google Scholar

[17] F. Courchamp, B. T. Grenfell, and T. H. Clutton-Brock, Impact of natural enemies on obligately cooperative breeders, Oikos 91 (2000), no. 2, 311–322, DOI: https://doi.org/10.1034/j.1600-0706.2000.910212.x. 10.1034/j.1600-0706.2000.910212.xSearch in Google Scholar

[18] S. Pal, S. K. Sasmal, and N. Pal, Chaos control in a discrete-time predator-prey model with weak Allee effect, Int. J. Biomath. 11 (2018), no. 7, 1850089, DOI: https://doi.org/10.1142/S1793524518500894. 10.1142/S1793524518500894Search in Google Scholar

[19] S. Işik, A study of stability and bifurcation analysis in discrete-time predator-prey system involving the Allee effect, Int. J. Biomath. 12 (2019), no. 1, 1950011, DOI: https://doi.org/10.1142/S1793524519500116. 10.1142/S1793524519500116Search in Google Scholar

[20] P. Chakraborty, S. Sarkar, and U. Ghosh, Stability and bifurcation analysis of a discrete prey-predator model with sigmoid functional response and Allee effect, Rend. Circ. Mat. Palermo (2) 70 (2021), no. 1, 253–273, DOI: https://doi.org/10.1007/s12215-020-00495-5. 10.1007/s12215-020-00495-5Search in Google Scholar

[21] L. Zhang and L. Zou, Bifurcations and control in a discrete predator-prey model with strong Allee effect, Int. J. Bifurcat. Chaos 28 (2018), no. 5, 1850062, DOI: https://doi.org/10.1142/S0218127418500621. 10.1142/S0218127418500621Search in Google Scholar

[22] L.-G. Yuan and Q.-G. Yang, Bifurcation, invariant curve and hybrid control in a discrete-time predator-prey system, Appl. Math. Model. 39 (2015), no. 8, 2345–2362, DOI: https://doi.org/10.1016/j.apm.2014.10.040. 10.1016/j.apm.2014.10.040Search in Google Scholar

[23] N. Britton, Essential Mathematical Biology, Springer, London, 2003. 10.1007/978-1-4471-0049-2Search in Google Scholar

[24] S. N. Elaydi, Discrete Chaos: With Applications in Science and Engineering, Chapman and Hall/CRC, London, 2007. 10.1201/9781420011043Search in Google Scholar

Received: 2024-02-05
Revised: 2024-04-11
Accepted: 2024-04-11
Published Online: 2024-06-20

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Special Issue on Contemporary Developments in Graph Topological Indices
  2. On the maximum atom-bond sum-connectivity index of graphs
  3. Upper bounds for the global cyclicity index
  4. Zagreb connection indices on polyomino chains and random polyomino chains
  5. On the multiplicative sum Zagreb index of molecular graphs
  6. The minimum matching energy of unicyclic graphs with fixed number of vertices of degree two
  7. Special Issue on Convex Analysis and Applications - Part I
  8. Weighted Hermite-Hadamard-type inequalities without any symmetry condition on the weight function
  9. Scattering threshold for the focusing energy-critical generalized Hartree equation
  10. (pq)-Compactness in spaces of holomorphic mappings
  11. Characterizations of minimal elements of upper support with applications in minimizing DC functions
  12. Some new Hermite-Hadamard-type inequalities for strongly h-convex functions on co-ordinates
  13. Global existence and extinction for a fast diffusion p-Laplace equation with logarithmic nonlinearity and special medium void
  14. Extension of Fejér's inequality to the class of sub-biharmonic functions
  15. On sup- and inf-attaining functionals
  16. Regularization method and a posteriori error estimates for the two membranes problem
  17. Rapid Communication
  18. Note on quasivarieties generated by finite pointed abelian groups
  19. Review Articles
  20. Amitsur's theorem, semicentral idempotents, and additively idempotent semirings
  21. A comprehensive review of the recent numerical methods for solving FPDEs
  22. On an Oberbeck-Boussinesq model relating to the motion of a viscous fluid subject to heating
  23. Pullback and uniform exponential attractors for non-autonomous Oregonator systems
  24. Regular Articles
  25. On certain functional equation related to derivations
  26. The product of a quartic and a sextic number cannot be octic
  27. Combined system of additive functional equations in Banach algebras
  28. Enhanced Young-type inequalities utilizing Kantorovich approach for semidefinite matrices
  29. Local and global solvability for the Boussinesq system in Besov spaces
  30. Construction of 4 x 4 symmetric stochastic matrices with given spectra
  31. A conjecture of Mallows and Sloane with the universal denominator of Hilbert series
  32. The uniqueness of expression for generalized quadratic matrices
  33. On the generalized exponential sums and their fourth power mean
  34. Infinitely many solutions for Schrödinger equations with Hardy potential and Berestycki-Lions conditions
  35. Computing the determinant of a signed graph
  36. Two results on the value distribution of meromorphic functions
  37. Zariski topology on the secondary-like spectrum of a module
  38. On deferred f-statistical convergence for double sequences
  39. About j-Noetherian rings
  40. Strong convergence for weighted sums of (α, β)-mixing random variables and application to simple linear EV regression model
  41. On the distribution of powered numbers
  42. Almost periodic dynamics for a delayed differential neoclassical growth model with discontinuous control strategy
  43. A new distributionally robust reward-risk model for portfolio optimization
  44. Asymptotic behavior of solutions of a viscoelastic Shear beam model with no rotary inertia: General and optimal decay results
  45. Silting modules over a class of Morita rings
  46. Non-oscillation of linear differential equations with coefficients containing powers of natural logarithm
  47. Mutually unbiased bases via complex projective trigonometry
  48. Hyers-Ulam stability of a nonlinear partial integro-differential equation of order three
  49. On second-order linear Stieltjes differential equations with non-constant coefficients
  50. Complex dynamics of a nonlinear discrete predator-prey system with Allee effect
  51. The fibering method approach for a Schrödinger-Poisson system with p-Laplacian in bounded domains
  52. On discrete inequalities for some classes of sequences
  53. Boundary value problems for integro-differential and singular higher-order differential equations
  54. Existence and properties of soliton solution for the quasilinear Schrödinger system
  55. Hermite-Hadamard-type inequalities for generalized trigonometrically and hyperbolic ρ-convex functions in two dimension
  56. Endpoint boundedness of toroidal pseudo-differential operators
  57. Matrix stretching
  58. A singular perturbation result for a class of periodic-parabolic BVPs
  59. On Laguerre-Sobolev matrix orthogonal polynomials
  60. Pullback attractors for fractional lattice systems with delays in weighted space
  61. Singularities of spherical surface in R4
  62. Variational approach to Kirchhoff-type second-order impulsive differential systems
  63. Convergence rate of the truncated Euler-Maruyama method for highly nonlinear neutral stochastic differential equations with time-dependent delay
  64. On the energy decay of a coupled nonlinear suspension bridge problem with nonlinear feedback
  65. The limit theorems on extreme order statistics and partial sums of i.i.d. random variables
  66. Hardy-type inequalities for a class of iterated operators and their application to Morrey-type spaces
  67. Solving multi-point problem for Volterra-Fredholm integro-differential equations using Dzhumabaev parameterization method
  68. Finite groups with gcd(χ(1), χc (1)) a prime
  69. Small values and functional laws of the iterated logarithm for operator fractional Brownian motion
  70. The hull-kernel topology on prime ideals in ordered semigroups
  71. ℐ-sn-metrizable spaces and the images of semi-metric spaces
  72. Strong laws for weighted sums of widely orthant dependent random variables and applications
  73. An extension of Schweitzer's inequality to Riemann-Liouville fractional integral
  74. Construction of a class of half-discrete Hilbert-type inequalities in the whole plane
  75. Analysis of two-grid method for second-order hyperbolic equation by expanded mixed finite element methods
  76. Note on stability estimation of stochastic difference equations
  77. Trigonometric integrals evaluated in terms of Riemann zeta and Dirichlet beta functions
  78. Purity and hybridness of two tensors on a real hypersurface in complex projective space
  79. Classification of positive solutions for a weighted integral system on the half-space
  80. A quasi-reversibility method for solving nonhomogeneous sideways heat equation
  81. Higher-order nonlocal multipoint q-integral boundary value problems for fractional q-difference equations with dual hybrid terms
  82. Noetherian rings of composite generalized power series
  83. On generalized shifts of the Mellin transform of the Riemann zeta-function
  84. Further results on enumeration of perfect matchings of Cartesian product graphs
  85. A new extended Mulholland's inequality involving one partial sum
  86. Power vector inequalities for operator pairs in Hilbert spaces and their applications
  87. On the common zeros of quasi-modular forms for Γ+0(N) of level N = 1, 2, 3
  88. One special kind of Kloosterman sum and its fourth-power mean
  89. The stability of high ring homomorphisms and derivations on fuzzy Banach algebras
  90. Integral mean estimates of Turán-type inequalities for the polar derivative of a polynomial with restricted zeros
  91. Commutators of multilinear fractional maximal operators with Lipschitz functions on Morrey spaces
  92. Vector optimization problems with weakened convex and weakened affine constraints in linear topological spaces
  93. The curvature entropy inequalities of convex bodies
  94. Brouwer's conjecture for the sum of the k largest Laplacian eigenvalues of some graphs
  95. High-order finite-difference ghost-point methods for elliptic problems in domains with curved boundaries
  96. Riemannian invariants for warped product submanifolds in Q ε m × R and their applications
  97. Generalized quadratic Gauss sums and their 2mth power mean
  98. Euler-α equations in a three-dimensional bounded domain with Dirichlet boundary conditions
  99. Enochs conjecture for cotorsion pairs over recollements of exact categories
  100. Zeros distribution and interlacing property for certain polynomial sequences
  101. Random attractors of Kirchhoff-type reaction–diffusion equations without uniqueness driven by nonlinear colored noise
  102. Study on solutions of the systems of complex product-type PDEs with more general forms in ℂ2
  103. Dynamics in a predator-prey model with predation-driven Allee effect and memory effect
  104. A note on orthogonal decomposition of 𝔰𝔩n over commutative rings
  105. On the δ-chromatic numbers of the Cartesian products of graphs
  106. Binomial convolution sum of divisor functions associated with Dirichlet character modulo 8
  107. Commutator of fractional integral with Lipschitz functions related to Schrödinger operator on local generalized mixed Morrey spaces
  108. System of degenerate parabolic p-Laplacian
  109. Stochastic stability and instability of rumor model
  110. Certain properties and characterizations of a novel family of bivariate 2D-q Hermite polynomials
  111. Stability of an additive-quadratic functional equation in modular spaces
  112. Monotonicity, convexity, and Maclaurin series expansion of Qi's normalized remainder of Maclaurin series expansion with relation to cosine
  113. On k-prime graphs
  114. On the existence of tripartite graphs and n-partite graphs
  115. Classifying pentavalent symmetric graphs of order 12pq
  116. Almost periodic functions on time scales and their properties
  117. Some results on uniqueness and higher order difference equations
  118. Coding of hypersurfaces in Euclidean spaces by a constant vector
  119. Cycle integrals and rational period functions for Γ0+(2) and Γ0+(3)
  120. Degrees of (L, M)-fuzzy bornologies
  121. A matrix approach to determine optimal predictors in a constrained linear mixed model
  122. On ideals of affine semigroups and affine semigroups with maximal embedding dimension
  123. Solutions of linear control systems on Lie groups
  124. A uniqueness result for the fractional Schrödinger-Poisson system with strong singularity
  125. On prime spaces of neutrosophic extended triplet groups
  126. On a generalized Krasnoselskii fixed point theorem
  127. On the relation between one-sided duoness and commutators
  128. Non-homogeneous BVPs for second-order symmetric Hamiltonian systems
  129. Erratum
  130. Erratum to “Infinitesimals via Cauchy sequences: Refining the classical equivalence”
  131. Corrigendum
  132. Corrigendum to “Matrix stretching”
  133. Corrigendum to “A comprehensive review of the recent numerical methods for solving FPDEs”
Downloaded on 17.6.2025 from https://www.degruyterbrill.com/document/doi/10.1515/math-2024-0013/html
Scroll to top button