This paper investigates the dynamics of non-autonomous cooperative systems of difference equations with asymptotically constant coefficients. We are mainly interested in global attractivity results for such systems and the application of such results to evolutionary population cooperation models. We use two methods to extend the global attractivity results for autonomous cooperative systems to related non-autonomous cooperative systems which appear in recent problems in evolutionary dynamics.

This paper investigates the rate of convergence of a certain mixed monotone rational second-order difference equation with quadratic terms. More precisely we give the precise rate of convergence for all attractors of the difference equation $x_{n+1}=\frac{Ax_{n}^{2}+Ex_{n-1}}{x_{n}^{2}+f}$, where all parameters are positive and initial conditions are non-negative.The mentioned methods are illustrated in several characteristic examples. 2020 Mathematics Subject Classification. 39A10, 39A20, 65L20.

We investigate the period-two trichotomies of solutions of the equation $$x_{n+1} = f(x_{n}, x_{n-1},x_{n-2}), \quad n=0, 1, \ldots $$ where the function $f$ satisfies certain monotonicity conditions. We give fairly general conditions for period-two trichotomies to occur and illustrate the results with numerous examples.   1991 Mathematics Subject Classification. 39A10, 39A11

We investigate the global character of the difference equation of the form $$ x_{n+1} = f(x_n, x_{n-1},\ldots, x_{n-k+1}), \quadn=0,1, \ldots $$ with several equilibrium points, where $f$ is increasing in all its variables. We show that a considerable number of well known difference equations can be embeded into this equation through the iteration process. We also show that a negative feedback condition can be used to determine a part of the basin of attraction of different equilibrium points, and that the boundaries of the basins of attractions of different locally asymptotically stable equilibrium points are in fact the global stable manifolds of neighboring saddle or non-hyperbolic equilibrium points.   2000 Mathematics Subject Classification. 39A10, 39A11

This paper investigates an autonomous discrete-time glycolytic oscillator model with a unique positive equilibrium point which exhibits chaos in the sense of Li–Yorke in a certain region of the parameters. We use Marotto’s theorem to prove the existence of chaos by finding a snap-back repeller. The illustration of the results is presented by using numerical simulations.

This paper investigates the dynamics of non-autonomous competitive systems of difference equations with asymptotically constant coefficients. We are mainly interested in global attractivity results for such systems and the application of such results to the evolutionary population of competition models of two species.

<jats:p>We present the bifurcation results for the difference equation <jats:inline-formula> <math xmlns="http://www.w3.org/1998/Math/MathML" id="M1"> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mi>n</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mrow> <msubsup> <mi>x</mi> <mi>n</mi> <mn>2</mn> </msubsup> </mrow> <mo>/</mo> <mfenced open="(" close=")" separators="|"> <mrow> <mi>a</mi> <msubsup> <mi>x</mi> <mi>n</mi> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>x</mi> <mrow> <mi>n</mi> <mo>−</mo> <mn>1</mn> </mrow> <mn>2</mn> </msubsup> <mo>+</mo> <mi>f</mi> </mrow> </mfenced> </math> </jats:inline-formula> where <jats:inline-formula> <math xmlns="http://www.w3.org/1998/Math/MathML" id="M2"> <mi>a</mi> </math> </jats:inline-formula> and <jats:inline-formula> <math xmlns="http://www.w3.org/1998/Math/MathML" id="M3"> <mi>f</mi> </math> </jats:inline-formula> are positive numbers and the initial conditions <jats:inline-formula> <math xmlns="http://www.w3.org/1998/Math/MathML" id="M4"> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msub> </math> </jats:inline-formula> and <jats:inline-formula> <math xmlns="http://www.w3.org/1998/Math/MathML" id="M5"> <msub> <mrow> <mi>x</mi> </mrow> <mrow> <mn>0</mn> </mrow> </msub> </math> </jats:inline-formula> are nonnegative numbers. This difference equation is one of the perturbations of the sigmoid Beverton–Holt difference equation, which is a major mathematical model in population dynamics. We will show that this difference equation exhibits transcritical and Neimark–Sacker bifurcations but not flip (period-doubling) bifurcation since this difference equation cannot have period-two solutions. Furthermore, we give the asymptotic approximation of the invariant manifolds, stable, unstable, and center manifolds of the equilibrium solutions. We give the necessary and sufficient conditions for global asymptotic stability of the zero equilibrium as well as sufficient conditions for global asymptotic stability of the positive equilibrium.</jats:p>

It is shown that locally asymptotically stable equilibria of planar cooperative or competitive maps have basin of attraction \begin{document}$ \mathcal{B} $\end{document} with relatively simple geometry: the boundary of each component of \begin{document}$ \mathcal{B} $\end{document} consists of the union of two unordered curves, and the components of \begin{document}$ \mathcal{B} $\end{document} are not comparable as sets. The boundary curves are Lipschitz if the map is of class \begin{document}$ C^1 $\end{document} . Further, if a periodic point is in \begin{document}$ \partial \mathcal{B} $\end{document} , then \begin{document}$ \partial\mathcal{B} $\end{document} is tangential to the line through the point with direction given by the eigenvector associated with the smaller characteristic value of the map at the point. Examples are given.

<jats:p>In this paper, certain dynamic scenarios for general competitive maps in the plane are presented and applied to some cases of second-order difference equation <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>f</mml:mi><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:mrow></mml:mfenced><mml:mo>,</mml:mo><mml:mtext> </mml:mtext><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>0,1</mml:mn><mml:mo>,</mml:mo><mml:mo>…</mml:mo></mml:math>, where <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"><mml:mi>f</mml:mi></mml:math> is decreasing in the variable <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub></mml:math> and increasing in the variable <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M4"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub></mml:math>. As a case study, we use the difference equation <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi><mml:mo>+</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>/</mml:mo><mml:mfenced open="(" close=")" separators="|"><mml:mrow><mml:mi>c</mml:mi><mml:msubsup><mml:mi>x</mml:mi><mml:mrow><mml:mi>n</mml:mi><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mn>2</mml:mn></mml:msubsup><mml:mo>+</mml:mo><mml:mi>d</mml:mi><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>n</mml:mi></mml:mrow></mml:msub><mml:mo>+</mml:mo><mml:mi>f</mml:mi></mml:mrow></mml:mfenced></mml:mrow></mml:mfenced><mml:mo>,</mml:mo><mml:mtext> </mml:mtext><mml:mi>n</mml:mi><mml:mo>=</mml:mo><mml:mn>0,1</mml:mn><mml:mo>,</mml:mo><mml:mo>…</mml:mo></mml:math>, where the initial conditions <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6"><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mo>−</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msub><mml:mo>≥</mml:mo><mml:mn>0</mml:mn></mml:math> and the parameters satisfy <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7"><mml:mi>c</mml:mi><mml:mo>,</mml:mo><mml:mi>d</mml:mi><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mo>></mml:mo><mml:mn>0</mml:mn></mml:math>. In this special case, we characterize completely the global dynamics of this equation by finding the basins of attraction of its equilibria and periodic solutions. We describe the global dynamics as a sequence of global transcritical or period-doubling bifurcations.</jats:p>

We investigate the global dynamics of the following rational difference equation of second order\begin{equation*}x_{n+1}=\frac{Ax_{n}^{2}+Ex_{n-1}}{x_{n}^{2}+f},\quad n=0,1,\ldots ,\end{equation*}where the parameters $A$ and $E$ are positive real numbers and the initial conditions $x_{-1}$ and $x_{0}$ are arbitrary non-negative real numbers such that $x_{-1}+x_{0}>0$. The transition function associated with the right-hand side of this equation is always increasing in the second variable and can be either increasing or decreasing in the first variable depending on the parametric values. The unique feature of this equation is that the second iterate of the map associated with this transition function changes from strongly competitive to strongly cooperative. Our main tool for studying the global dynamics of this equation is the theory of monotone maps while the local stability is determined by using center manifold theory in the case of the nonhyperbolic equilibrium point.

We investigate the period-doubling and Naimark–Sacker bifurcations of the equilibrium of the difference equation xn+1 = γxn−1 + δxn Cxn−1 + xn where the parameters γ, δ, C are positive numbers and the initial conditions x−1 and x0 are arbitrary nonnegative numbers such that x−1 + x0 > 0. AMS Subject Classifications: 39A10, 39A20, 39A60, 37B25.