Dynamics Analysis of a Discrete-Time Commensalism Model with Additive Allee for the Host Species

: We propose and study a class of discrete-time commensalism systems with additive Allee effects on the host species. First, the single species with additive Allee effects is analyzed for existence and stability, then the existence of ﬁxed points of discrete systems is given, and the local stability of ﬁxed points is given by characteristic root analysis. Second, we used the center manifold theorem and bifurcation theory to study the bifurcation of a codimension of one of the system at non-hyperbolic ﬁxed points, including ﬂip, transcritical, pitchfork, and fold bifurcations. Furthermore, this paper used the hybrid chaos method to control the chaos that occurs in the ﬂip bifurcation of the system. Finally, the analysis conclusions were veriﬁed by numerical simulations. Compared with the continuous system, the similarities are that both species’ densities decrease with increasing Allee values under the weak Allee effect and that the host species hastens extinction under the strong Allee effect. Further, when the birth rate of the beneﬁted species is low and the time is large enough, the beneﬁted species will be locally asymptotically stabilized. Thus, our new ﬁnding is that both strong and weak Allee effects contribute to the stability of the beneﬁted species under certain conditions.


Introduction
Interactions between different species play a crucial role in shaping and maintaining balance within ecosystems.Interspecies interactions hold the potential to deliver mutual advantages, inflict harm, or yield no discernible impact on each of the participating species.These interactions are categorized by paired outcomes, symbolized as a positive, negative, or neutral result [1].Commensal relationships, where one species gains a benefit while the other is unaffected, are often cited in the ecological literature, yet they receive surprisingly limited attention in research.In ecosystems, both commensalism and mutualism models involve positive interactions between populations.There is no harm in commensalism, in which one side is beneficial, while the other side is not beneficial and harmless.
A Belgian zoologist named Pierre-Joseph van Beneden was the pioneer in situating commensalism within a biological context.He contrasted commensalism with two other types of interactions: parasitism and mutualism [1].In mutualism, species are interdependent and both benefit, while in commensalism, one species gains advantages without strong dependency or notable harm to the other.Despite its milder impact, commensalism adds complexity to ecological relationships and community organization.According to [2,3], depending on the population density, in some cases, the mutualism relationship shifts to a resemblance closer to commensalism, e.g., the relation among the carnivorous plant Roridula dentata and hemipteran Pameridae marlothii.The larvae of midges and mosquitoes that live on pitcher plants exhibit a commensalism relationship.Midges, as upstream resource consumers, produce large quantities of bacteria and particles that contribute to the growth of mosquitoes.However, the experiments in [4] showed that mosquitoes do not affect midges.There is a symbiotic relationship between the two.Hari Prasad [5] pointed out that small green epiphytes that grow on other plants prepare themselves for survival by occupying the space of other plants to absorb the water and minerals that the roots obtain from the air.In turn, plants are unaffected.Squirrels burrow intooak trees for living space and food storage.Yet, the oak is neither beneficial nor harmful.To avoid it enemies, the clownfish will choose to hide in the tentacles of the sea anemone, and the tentacles are not affected by the clownfish, and so on.
From the standpoint of mathematical modeling, the investigation of the commensalism system was pioneered by Sun and Wei [6] in 2003.They were the first to explore this area using a two-species model.Over the past few years, an increasing number of researchers have directed their focus towards utilizing mathematical modeling methods to explore the commensalism model among populations (see [7][8][9][10][11][12][13][14][15][16][17] and the references therein).Afterwards, Han et al. [18] proposed a continuous time Lotka-Volterra-type commensalism system as follows: where x(0) > 0, y(0) > 0 and t ∈ [0, ∞).x(t) and y(t) denote the population densities of the commensal and host, respectively, at time t; b 1 and b 2 denote the intrinsic growth rates of the two species, respectively; a 11 and a 22 denote the intraspecific competition coefficient of the commensal and host; a 12 > 0 denotes the coefficient of influence that the host produces in the commensal.We should notice that b 1 , a 11 , a 12 , b 2 , a 22 are positive constants.The authors' study showed that Model (1) has a globally stable coexistence fixed point.Next, feedback control is considered by adding control variables.The results showed that the position of the coexistence fixed point is altered without affecting its global steady-state properties.
After that, many scholars began to consider the effects of various factors based on the system (1) and analyze the dynamic behaviors of the corresponding systems, for example periodic solutions [19,20], harvesting [21][22][23][24][25], the Allee effect [26][27][28], and delay [29,30].In [19], the authors investigated a modified commensal symbiosis model with Michaelis-Menten-type harvesting and observed positive periodic solutions.In [20], Chen et al. presented another discrete commensal symbiosis model with a Hassell-Varley-type functional response.Their study revealed a positive periodic solution of the system, which is both locally and globally stable.Considering a Lotka-Volterra-type commensal model subject to the Allee effect on the unaffected species, which occurs at a low population density, Xinyu Guan [27] observed that the Allee effect has no influence on the final densities of the unaffected species and the benefited species.The study involved numerical simulations that demonstrated a considerable elongation in the time required for the system subjected to an Allee effect to attain its stable steady-state solution, implying that the Allee effect has an unstable effect on the system.In [28], the authors explored a discrete-time commensal symbiosis model and observed period-doubling bifurcation phenomenon around the interior equilibrium state.
According to previous studies by W.C. Allee, population density has a positive impact on reproduction in a wide range of terrestrial and aquatic species.It could potentially contribute to the prolonged survival of the species in unfavorable environments and provide increased protection against harmful substances [31,32].When the population density drops below a specific threshold or when the population is too diminutive, the task of finding a mate can become more challenging.Similarly, in species where group activities such as fending off predators or searching for sustenance are impeded due to low population densities, these activities become less proficient.This positive correlation between population size and individual fitness is termed the Allee effect.In the past few years, populations from various ecosystems have become low or even at risk of becoming endangered due to escalating human predation, heightened occurrences of natural disasters, and rising global temperatures.In such situations, the populations may have trouble mating, foraging, and fending off natural enemies.Previous mathematical studies have shown that Allees add more complexity to the system's dynamic behavior [33,34].In 2022, He et al. [11] proposed the following continuous-time commensalism system incorporating an additive Allee to the host: where all of the parameters are positive constants and initial values are x(0) > 0, y(0) > 0, and t ∈ [0, ∞), the same as the previous definition in [11].He et al. [11] gave the existence and local stability of equilibrium points, then proved that, under suitable conditions, E r ( r b , 0) and E * 1 (x * 1 , y 1 ) are globally asymptotically stable, respectively.There are saddlenode bifurcation at E * 3 (x * 3 , y 3 ), E 3 (0, y 3 ) and transcritical bifurcation at E 0 (0, 0), E r ( r b , 0), respectively.
In recent years, discrete-time systems have also been one of the hotspots of scholars' research.The reason is that populations with non-overlapping generations are better described by difference equations, which are also easier to simulate numerically.Discrete-time systems also have more topological classifications and bifurcations at fixed points than corresponding continuous-time systems, with a codimension of one, such as flip, transcritical, pitchfork, fold, and Neimark-Sacker bifurcations, and a codimension of two, such as 1:1 and 1:2, and fold-flip bifurcations, e.g., [35][36][37][38][39][40][41].To the best of the authors' knowledge, limited research has been conducted on discrete models incorporating additive Allee effects, highlighting an ongoing need for such investigations.In this work, we explored the qualitative behavior of the continuous-time commensal model (2) by converting it into a discrete-time model.
Transformation is the first step in reducing the parameters, and it uses x = d b x, y = d e ȳ, t = τ d .Then, the bars are removed, and System (2) can be transformed as where α = r d , β = c e , M = me d 2 , and A = ae d .The initial value is x(0) > 0, y(0) > 0 and τ ∈ [0, ∞).By using the piecewise constant parameter method [42] to discretize the system (3), we have where 0 ≤ n ≤ τ < n + 1 and [τ] is the round-down function.First, the simple rounding function operation can be used to determine that the right-hand side of the ordinary differential Equation ( 4) is a constant value, when other parameters are fixed.Secondly, the integral of the interval n to τ is performed on both sides of the equation at the same time, and when τ tends to n + 1, the following is acquired: for n = {0, 1, 2, ...}, so we obtain the following a discrete-time commensalism mode with the additive Allee effect for the host species: where x(n): = x n and y(n): = y n , n = {0, 1, 2, ...}.The discrete-time commensalism system is now defined by mapping Then, the discrete-time single-species mode with the additive Allee is expressed as where M y+A represents the additive Allee.The Allee effect in the map (5) and the map ( 6) is weak when 0 < M < A. The Allee effect in the map (5) and the map (6) is strong when M > A.
In this paper, Section 2 provides the existence and local stability of the map (6) at fixed points.Then, Section 3 provides the existence and local stability of the map (5) at fixed points.Section 3.3 gives the bifurcation of different types of a codimension of one under certain conditionals including flip, fold, transcritical, and pitchfork bifurcations.Moreover, Section 3.4 provides a chaos-control system for the occurrence of chaos due to bifurcation.Finally, the correctness of the numerical simulations verifies the conclusions.

The Existence of Fixed Point
The fixed points of the map (6) satisfy the equations as follows: Obviously, 0 is a trivial fixed point of the map (6).Further, the positive fixed point of the map (6) needs to solve 1 − y − M y+A = 0, and we can obtain the following equation with the same solution through simplification: Let (M) be the discriminant function of G(y), then, according to the discriminant formula of the quadratic function, we have Let M * be the zero point of (M), then Thus, if 0 < M < M * , (M) > 0, then (7) has two roots of and (M) = 0, then (7) has a double-root denoted by y 3 = 1−A 2 ; when M > M * , (M) < 0, then (7) has no root.Then, we obtain the following result.
Theorem 1.The map always has a fixed point 0, regardless of how the parameters are altered in the map (6).The parameter requirements for the existence of positive fixed points are then stated: (1) (Case 1 : 0 < M < A, i.e., weak Allee) Existence of a positive fixed point y 2 in the map (6).
(ii) Map (6) has no other fixed point if A ≥ 1.
(3) (Case Firstly, it is easy to see that the map (6) at a fixed point y satisfies Then, by a simple computation, we have So, when 0 < M < A, d f dy (0) > 1, the trivial fixed point 0 is a source; when M > A, d f dy (0) ∈ (0, 1), the trivial fixed point 0 is a sink that is locally asymptotically stable; when M = A, d f dy (0) = 1, the trivial fixed point 0 is non-hyperbolic.Secondly, according to Theorem 1, if the positive fixed point y i , i = 1, 2, 3 exists, we can obtain and then, Thus, it follows that y i is a rigorously decreasing function of d f dy (y i ).From Theorem 1, when M = M * and A ∈ (0, 1), there is a positive fixed point y 3 .By a direct computation, we have Furthermore, it is possible to verify y 1 < y 3 < y 2 .Namely, d f dy (y 1 ) > 1 > d f dy (y 2 ).Hence, the positive fixed point y 1 is always the source; the positive fixed point y 3 is always non-hyperbolic.Now, in order to compare So, by the above analysis, the positive fixed point y 2 is a sink.Thus, we have the following result.
Theorem 2. When M < A, i.e., a weak Allee, or (M = A and 0 < A < 1), the map (6) has a unique positive fixed point y 2 that is globally asymptotically stable.
Firstly, let f (y) = y exp 1 − y − M y+A .Then, let us construct a Lyapunov function as Then, find the first derivative of V(y).We can obtain dy < 0, then V(y) is strictly monotonically decreasing with respect to y; if y > y 2 , d(V(y)) dy > 0, then V(y) is an increasing function with respect to y.Thus, V(y) has a minimum value V(y 2 ) = 0. Obviously, V(y) is positive definite for any y > 0 and for y = y 2 .
V(y) derives the solution of the map (6), and we have: To prove that the map ( 6) is globally stable is to prove that ∆V(y) is negative definite for all y > 0 and for y = y 2 .To test this condition, we require that ∆V(y) has a unique global maximum at y = y 2 , which is difficult to prove using pure analysis.For the specific parameter, this can often be performed computationally; first, we consider whether y = y 2 is the local maximum of ∆V(y), and we calculate: Hence, if M < A or M = A and 0 < A < 1, there is a unique positive fixed point y 2 , and ∆V(y) has global maximum at y = y 2 .Then, Theorem 2 is proven.

Parameter Conditions
Existence Stability

Existence and Local Stability of Fixed Points
The fixed points of the map (5) satisfy the equations as follows: Obviously, E 0 (0, 0) and E 1 (α, 0) are two boundary fixed points of the map (5).Further, the other fixed points of the map ( 5) are summed up in Table 2. Let E 2i (0, y i ), i = 1, 2, 3 be the boundary fixed point of the extinction of the commensal species and E * 3i (x * i , y i ), i = 1, 2, 3 be the fixed points of the coexistence of the two populations, where x * i = α + βy i .
Table 2. Boundary and positive fixed point of the map (5)'s existence.

Parameter Conditions Existence
no other fixed points Remark 1.To prepare for the prospect of bifurcations, we offer the following observations: ( At any fixed point E(x, y), the Jacobian matrix of the map ( 5) is given as where Let the eigenvalues of the matrix J(E) be λ 1 and λ 2 .
We classify fixed points topologically and examine their local stability by [43].
Note that, when M = 0, and for the case M > 0, Thus, the following conclusion can be drawn.

Theorem 4.
Assume that M = 0, and we have the following description: (1) E 0 (0, 0) is always a source since two eigenvalues are e α > 1 and e > 1.
By Theorem 4(iii), the two eigenvalues of J(E 1 ) are λ 1 = −1, λ 2 = e; hence, E 1 is non-hyperbolic.The conditions that satisfy Theorem 4(iii) for the set Ω E f are written as follows: Choose the bifurcation parameter α, and set it to α = 2 + ρ.A suitably tiny perturbation term is ρ.According to a quick calculation, the center manifold of the map (5) in Ω E f is y = 0.In this instance, the map (5) may be expressed as As a result, at E 1 , a flip bifurcation is formed.
Theorem 5. Assume that M > 0; in this case, the two eigenvalues of the map (5) at fixed point E 0 (0, 0) are e α and e 1− A M , then E 0 (0, 0) is: Theorem 6.The two eigenvalues of the map (5) at the fixed point E 1 (α, 0) are 1 − α and e 1− A M ; consequently, E 1 (α, 0) is: (1) A sink if and only if 0 < α < 2 and M > A; (2) A source if and only if α > 2 and 0 < M < A; (3) A saddle if and only if 0 < α < 2 and 0 < M < A or α > 2 and M > A; (4) Non-hyperbolic if M = A or α = 2. Proof.At the boundary fixed point E 1 (α, 0), whose eigenvalues are 1 − α < 1 and e Then, the result follows.Now, we give the Jacobian matrix of the map (5) at the boundary fixed points E 2i 0, y i (i = 1, 2, 3) as follows: where Clearly, the two eigenvalues of J(E 2i ) are e α+βy i > 1 and J i .Note that J i was analyzed in Section 2.2.Then, we have J 1 > 1, J 2 ∈ (−1, 1) and J 3 = 1.Thus, we obtain the following conclusions: Next, the Jacobian matrix at positive points E * 3i x * i , y i (i = 1, 2, 3) is obtained as where x * i = α + βy i , whose eigenvalues are 1 − x * i < 1, and J i was analyzed above.Then, we obtain This is equivalent to the following two situations.
Lemma 1. Suppose that the sequences x(n) satisfy x(n) > 0 and Lemma 2. Suppose that the sequences x(n) satisfy x(n) > 0 and where a and b are positive constants.Then, where M was given by Lemma 2.
Secondly, from Theorem 2, we obtain that, if 0 < M < A or M = A, 0 < A < 1 holds, y(n) is any positive solution of the map (6), then Now, we consider a system as follows: where x 1 (n) = α + βy 2 is any positive solution of System (9).We obtain the following result.
holds, the positive fixed point E * 32 (x * 2 , y 2 ) of the map ( 5) is globally attractive.In other words, where x 1 (n) is any positive solution of System (9).

Bifurcation Analysis of a Codimension of One
According to the prior analysis of the map (5) at fixed points, in this section, the center manifold theorem is used to reduce dimensionality, and then, bifurcation theory is employed to investigate the bifurcations of a codimension of one at non-hyperbolic fixed points by [44,45].
(ii) From Theorem 6(4), when (α, β, M, A) Taking α as a bifurcation parameter and assuming ε that is a sufficiently small perturbation parameter, namely ε 1, then, the perturbations corresponding to System (5) are mapped as follows: Then, we shift E 1 in the map ( 12) to (0,0) using a transformation u here, Taking the map (23) at origin (0, 0), the Jacobian matrix is Invertible matrices T are constructed as follows: the map (23) becomes where According to the center manifold theory, (X, Y) = (0, 0) is stable for ε = 0 and can be determined as follows: g(X, ε) must satisfy the following relation: Substituting ( 26) into (27) and comparing the coefficients of ε 2 , Xε and X 2 in ( 27), we obtain Consequently, the map restricted to the center manifold W c (0, 0, 0) is expressed as From [44], it can be seen that the conditions for flip bifurcation to occur are: Thus, by [44], the system (5) undergoes a flip bifurcation if (α, β, M, A) ∈ Ω E1 FB holds at E 1 (α, 0).
Proof.(a) From Theorem 8 (2)(iii), we can easily obtain that, when (α, β, M, A) Then, to obtain the flip bifurcation, we regard β as the bifurcation parameter and shift the fixed point and is a sufficiently small new variable.We then have here, and The linearization matrix for the map (28) is Find an invertible matrix T, for instance T = a 010 a 010 −1 − a 100 λ 2 − a 100 and make the transformation: the map (18) yields where where According to the center manifold theory, a one-parameter simplified system of equations restricted to the manifold determines the stability of (X, Y) = (0, 0) in the neighborhood of δ = 0, given by W c (0, 0, 0) = (X, Y, δ) ∈ R 3 : Y = g(X, δ), g(0, 0) = 0, Dg(0, 0) = 0 with X and δ sufficiently small.Moreover, let From ( 31) and ( 32), we obtain Hence, the map restricted to the center manifold becomes where We give the following two transversal conditions in [44]: Therefore, a flip bifurcation occurs when (α, β, M, A) ∈ Ω E32 FB at the fixed point 1, and it is a new variable.We transform the fixed point and then, the map (5) becomes where The coefficient matrix is given by 33) into the normal form with matrix translation as follows: where is an invertible matrix.From ( 33) and (34), we have where and , 3  .

Chaos Control
From Theorem 12(a), it can be seen that a flip bifurcation of the map ( 5) is shown at the coexistence fixed point E * 32 = (x * 2 , y 2 ).This will cause chaos in the system, so we use a hybrid chaos-control method in [46], which is a combination of parameter perturbation and feedback control, to delay or get rid of chaos.First, the control system corresponding to the map ( 5) is given by where 0 < ρ < 1 and ρ is an external control parameter.Moreover, the calculation shows that System (38) and the map (5) have the same fixed points.The linearized Jacobian matrix of the map (38) at fixed point Obviously, the two eigenvalues of the Jacobian matrix According to the previous analysis, J 2 < 1, then −1 + M (y 2 +A) 2 < 0. We have λ 2 = 1 + ρy 2 −1 + M (y 2 +A) 2 < 1.Therefore, with the help of Definition 1(1), the controlled system (38)

Numerical Simulations
In this section, the analysis results obtained above are further studied through numerical simulations.

Parameter
Figure 8 is a diagram of the transcritical bifurcation on the M − y plane in the local range of the fixed point E 1 (3, 0) when M ∈ [0, 0.331], and the other parameters are given in Case (3).From Figure 8, it can be found that, when M < 0.15, E 1 (3, 0) is unstable, but when M > 0.15, there is a stable fixed point Case (4): Further, based on Case (3), we change the parameter A and let A = 1.At this time, by Theorem 11(ii), the map (5) generates a pitchfork bifurcation at E 1 (3, 0).Namely, it is a pitchfork bifurcation point.
Figure 9 is a diagram of the pitchfork bifurcation on the M − y plane in the local area of the fixed point E 1 .Select the initial value of x 0 = 0.5, y 0 = 0.3, M ∈ [0, 1.5], α = 3, and β = 1.It can be seen from Figure 9 that the fixed point E 1 is unstable, when M < 1; the stable fixed point E 1 (3, 0) appears when M > 1.Additionally, there are two stable fixed points   (5).From Figure 10a, the coexistence equilibria E * 32 are stable when β < 2, and when β = 2, the map (5) is unstable, resulting in a stable two-period solution; when β > 3.065, it will further generate four-dimensional solutions and constantly flip to generate chaotic sets.
Figure 10b shows the maximum Lyapunov exponents of Figure 10a, and β ranges in [0, 6].When the exponent value is less than 0, this indicates that the map ( 5) is stable in the small field of E * 32 (2, 0.5).In other words, when the exponent value is greater than 0, the map ( 5) is chaotic in the small field of E * 32 (2, 0.5).From Figure 10b, when β varies in [0, 2], the index is negative, except for β = 2, which is 0. When β belongs to (2, 3.345], the index is only 0 at β = 3.065 and β = 3.345, and in other cases, the index is negative.When β is greater than 3.345, most of the indices are positive and very few are negative numbers, which indicates that the map (5) produces chaotic behavior near the fixed point E * 32 (2, 0.5).We also give local bifurcation diagrams in Figure 11.Further, when we choose β = 1, 2.5, 3.2, 3.345, 4.5, 5, the phase diagrams are depicted in Figure 12.As β increases, it can be seen from Figure 12 that there will be one point, two points, four points, eight points, etc., gradually becoming chaotic.
Case (6): To study the impact of the Allee effect on the dynamic behavior of the map (5), we used the matlab 2017a drawing software to give the bifurcation diagram on the M − x plane and M − y.We selected four sets of parameters.From Figure 13a,c, it is shown that, with the increase of M, it helps the stability of the commensal populations.Further, we also found that a strong Allee effect reduces the species density of the commensal.Figure 13b,d are graphs showing the impact of the Allee effect M on the population density of the host species.It can be seen from Figure 13b,d that, if M increases, the host population is accelerated to extinction.It can be seen from Figure 13e,f that, as M increases, the map (5) moves from chaos to a two-period solution, but the map (5) is unstable.Case (8): Set the parameters α = 2.5, β = 1, M = 0.1, A = 4, and the initial value x 0 = 0.5, y 0 = 0.3.We can simply calculate y 2 ≈ 0.96, x * 2 ≈ 3.48, and ρx * 2 < 2 is equivalent to ρ < 0.57.Moreover, ρy 2 (1 − M (y 2 +A) 2 ) < 2 implies that ρ < 2.04.So, let us take a smaller ρ.From Theorem 13, when 0 < ρ < 0.57, it is locally asymptotically stable.Figure 15 shows that when the external control parameters change in the range [0,0.57], the map (5) can be stabilized again.In addition, we set the parameter α = 2.5, β = 1, A = 4, then the stability region diagram is given of the control system (28) in Figure 16.

Summary
This paper investigated a discrete-time commensalism mode with an additive Allee effect for the host species.First, we gave the existence and local stability of fixed points in the map ( 6), and we found that, under the weak Allee effect, the system has a unique positive fixed point of global attraction.Under the strong Allee effect, the stability of the system is related to the attraction domain, which has similar properties to the corresponding continuous system.Secondly, we studied the existence and stability of the solutions of the map (5), then gave the possible bifurcation of the map (5) by the center manifold theorem and bifurcation theory.We found that the map (5) has new bifurcation phenomena, such as flip and pitchfork bifurcations, which could not be found in the corresponding continuoustime system.Wei et al. [47] presented a commensal model with additive Allees in the first species and discussed the occurrence of saddle-node bifurcation and transcritical bifurcation.Examining a Lotka-Volterra commensal model featuring an Allee effect within the first species, Lin [48] revealed that heightened Allee effects lead to an escalation in the final population density of the species.In [13], Chen examined a comparable two-species commensal symbiosis model and noted that the Allee effect introduces instability to the system; nevertheless, this effect can be managed and controlled.From our investigation, it was revealed that the Allee increases the system complexity through multiple bifurcations.Recently, the studies by He et al. [11] revealed that, in the case of the weak Allee effect, the additive Allee effect negatively correlated with the final population density of both species.For a strong Allee effect, the additive Allee effect played a significant role in the extinction of the second species.Our study in this domain showed more-complex dynamics in comparison to the continuous-time model (2), which were presented through bifurcation analysis.Through numerical simulations of our proposed model, in the case where the first population had a lower birth rate, we observed that the Allee effect enhanced the stability of the commensal species, while hastening the extinction of the host species (see Figure 13a-d).These intricate and multifaceted dynamic behaviors mark our model as a truly novel contribution.Moreover, the map (5) may also have a fold-flip bifurcation of a codimension of two at fixed point E * 33 (x * 3 , y 3 ), which is something we need to study in the future.When M = 0, the map (5) has at most four positive fixed points, and only flip bifurcation can occur.When M = 0, the system may have six positive fixed points, and they may undergo flip, fold, transcritical, pitchfork, and fold-flip bifurcations, which are left for future investigation.
one eigenvalue of the Jacobian matrix J(E 1 ) is 1 − α and the other is 1.First, we choose M as the bifurcation parameter.Set µ = M − A, where |µ| 1, and it is a new variable.Next, we translate System (5) to the originwith u = x − α, v = y − 0 and perform the Taylor expansion, leading to 

Table 1
gives the above analysis.
(5) boundary fixed point E 21 exists in the map(5), and the eigenvalues of J(E 21 ) are e α+βy 1 > 1 and J 1 > 1.So, the fixed point E 21 is always a source.(2)When the boundary fixed point E 22 exists, the eigenvalues of J(E 22 ) are e α+βy 2 > 1 and |J 2 | < 1.So, the fixed point E 22 is always a saddle.(3)WhenM = M * and 0 < A < 1, the boundary fixed point E 23 exists in the map(5), and the eigenvalues of J(E 23 ) are e α+βy 3 > 1 and J 3 = 1.Consequently, E 23 is always non-hyperbolic.