Next Article in Journal
Length-Fuzzy Subalgebras in BCK/BCI-Algebras
Next Article in Special Issue
Existence of Solution, Filippov’s Theorem and Compactness of the Set of Solutions for a Third-Order Differential Inclusion with Three- Point Boundary Conditions
Previous Article in Journal
The Collapse of Ecosystem Engineer Populations
Previous Article in Special Issue
A Numerical Solution of Fractional Lienard’s Equation by Using the Residual Power Series Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Dynamics of Certain Mix Monotone Difference Equation

by
Senada Kalabušić
1,
Mehmed Nurkanović
2 and
Zehra Nurkanović
2,*
1
Department of Mathematics, University of Sarajevo, 71000 Sarajevo, Bosnia and Herzegovina
2
Department of Mathematics, University of Tuzla, 75000 Tuzla, Bosnia and Herzegovina
*
Author to whom correspondence should be addressed.
Mathematics 2018, 6(1), 10; https://doi.org/10.3390/math6010010
Submission received: 7 December 2017 / Revised: 3 January 2018 / Accepted: 7 January 2018 / Published: 12 January 2018
(This article belongs to the Special Issue Advances in Differential and Difference Equations with Applications)

Abstract

:
We investigate global dynamics of the following second order rational difference equation x n + 1 = x n x n 1 + α x n + β x n 1 a x n x n 1 + b x n 1 , where the parameters α , β , a , b are positive real numbers and initial conditions x 1 and x 0 are arbitrary positive real numbers. The map associated to the right-hand side of this equation is always decreasing in the second variable and can be either increasing or decreasing in the first variable depending on the corresponding parametric space. In most cases, we prove that local asymptotic stability of the unique equilibrium point implies global asymptotic stability.

1. Introduction and Preliminaries

In this paper, we investigate the local and global character of the equilibrium point of the following second order rational difference equation
x n + 1 = x n x n 1 + α x n + β x n 1 a x n x n 1 + b x n 1 , n = 0 , 1 ,
where the parameters α , β , a , b are nonnegative real numbers and initial conditions x 1 and x 0 are arbitrary nonnegative real numbers, such that a x n x n 1 + b x n 1 > 0 , n = 0 , 1 , .
Equation (1) is the special case of a general second order quadratic rational difference equation of the form
x n + 1 = A x n 2 + B x n x n 1 + C x n 1 2 + D x n + E x n 1 + F a x n 2 + b x n x n 1 + c x n 1 2 + d x n + e x n 1 + f , n = 0 , 1 ,
with nonnegative parameters and nonegative initial conditions such that A + B + C > 0 , a + b + c + d + e + f > 0 and a x n 2 + b x n x n 1 + c x n 1 2 + d x n + e x n 1 + f > 0 , n = 0 , 1 , . Several global asymptotic results for some special cases of Equation (2) were obtained in [1,2,3,4,5,6,7,8,9,10,11].
One interesting special case by (2) is the following rational difference equation studied in [12]:
x n + 1 = D x n + F d x n + e x n 1 , n = 0 , 1 ,
which represents discretization of the differential equation model in biochemical networks, see [13]. Notice that Equation (3) is an example of a rational difference equation, such that associated map is always strictly decreasing with respect to the second variable, and changes its monotonicity with respect to the first variable, i.e., can be increasing or decreasing depending on corresponding parametric space. Also, we see that Equation (3) is the special case of the linear rational difference equation
x n + 1 = D x n + E x n 1 + F d x n + e x n 1 + f , n = 0 , 1 ,
(which was investigated in detail in [12]) with well known and complicated dynamics, such as Lynes’ equation (see [14]).
There are not many papers that study in detail dynamics of the second order rational difference equations with quadratic terms such that associated map changes its monotonicity with respect to its variables. However, in [15] the behavior of the following rational difference equation has been investigated in great detail
x n + 1 = α + γ x n 1 B x n + D x n x n 1 + x n 1 , n = 0 , 1 ,
In both equations, (3) and (5), Theorems 1 and 2 were used in order to obtain the convergence results. In most cases of this paper we use the same results. However, in order to investigate the behaviour of the following four subsequences x 4 k k = 1 , x 4 k + 1 k = 0 , x 4 k + 2 k = 0 , x 4 k + 3 k = 0 that appears under the condition b α a β b = β b , where β > b a we cannot use this method, because the associated map in this case does not have the same monotonicity with respect to its variables in invariant interval. More precisely, the corresponding map changes its monotonicity in invariant interval with respect to the first variable. Instead of that, we use the brute-force method to show that each subsequence converges to the unique equilibrium point.
In the case when associated map of Equation (1) changes its monotonicity from “decreasing-decreasing” into “increasing-decreasing”, the problem of determining invariant interval appears. In all cases, we determine invariant interval and prove that the positive equilibrium of Equation (1), which is always locally asymptotically stable, is globally asymptotically stable for all values of the parameters, except in the case when α > α D (see Theorems 7–10).
The problem of determining invariant intervals in the case when the associated map changes its monotonicity with respect to its variables has been considered in [16,17]. Also, see [18,19,20].
Now, we state several well-known results.
Theorem 1.
[14] [Theorem 2.22] Let a , b be an interval, and suppose that f : a , b × a , b a , b is a continuous function. Consider the difference equation
x n + 1 = f x n , x n 1 , n = 0 , 1 ,
Assume that f satisfies the following two conditions:
a. f x , y is nondecreasing in x a , b for each y a , b , and f x , y is nonincreasing in y a , b for each x a , b ;
b. All solutions m , M a , b × a , b of the system
m = f m , M , M = f M , m
satisfy m = M .
Then, (6) has a unique equilibrium x ¯ a , b and every solution of (6) converges to x ¯ .
Theorem 2.
[12] [Theorem 1.4.7] Let a , b be an interval, and suppose that f : a , b × a , b a , b is a continuous function satisfying the following properties:
a. f x , y is nonincreasing in each of its arguments;
b. If m , M a , b × a , b is a solution of the system
m = f M , M , M = f m , m ,
then m = M .
Then, (6) has a unique equilibrium x ¯ a , b and every solution of (6) converges to x ¯ .
Theorem 3.
[21] [Theorem 1.4] Let f be the function from (6) with
1. 
f C 0 , × 0 , , 0 , ;
2. 
f u , v is nonincreasing in u and v respectively;
3. 
x f x , x is nondecreasing in x;
4. 
Equation (6) has a unique positive equilibrium x ¯ .
Then, every positive solution x n n = 1 of Equation (6) which is bounded from above and from below by positive constants converges to x ¯ .
Theorem 4.
[12] [Theorem 1.7.2] Assume that f C 0 , × 0 , , 0 , and that f x , y is decreasing in both arguments.
Let x ¯ be a positive equilibrium of Equation (6). Then, every oscillatory solution of Equation (6) has semicycles of length at most two.
Theorem 5.
[12] [Theorem 1.7.4] Assume that f C 0 , × 0 , , 0 , is such that: f x , y is increasing in x for each fixed y, and f x , y is decreasing in y for each fixed x.
Let x ¯ be a positive equilibrium of Equation (6). Then, except for the first semicycle, every oscillatory solution of Equation (6) has semicycles of length at least two.

2. Linearized Stability

In this section, we prove that Equation (1) has a unique equilibrium point which is always locally asymptotically stable.
The equilibrium point x ¯ of Equation (1) satisfies
x ¯ = x ¯ 2 + α x ¯ + β x ¯ a x ¯ 2 + b x ¯ .
Equation (7) has the unique positive solution
x ¯ = 1 b + D 2 a ,
where
D = ( b 1 ) 2 + 4 a ( α + β ) .
Theorem 6.
The unique positive Equilibrium (8) of Equation (1) is always locally asymptotically stable.
Proof. 
The real function f ( u , v ) associated to Equation (1) is given by
f ( u , v ) = u v + α u + β v a u v + b v .
Derivatives with respect to u and v evaluated at the equilibrium point (8) are respectively
p = b α + ( b a β ) x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) , and q = α x ¯ + α + β .
We have that the function f ( u , v ) is always decreasing with respect to the second variable and can be either decreasing or increasing with respect to the first variable, depending on the sign of nominator, that is depending on corresponding parametric space. Now, we check the conditions of Theorem 1.1.1, see [12]. The condition | p | < 1 q < 2 becomes
b α + ( b a β ) x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) < 1 + α x ¯ + α + β < 2 .
The second inequality is equivalent to α < x ¯ + α + β , which is always true. The first inequality becomes
1 + α x ¯ + α + β < b α + ( b a β ) x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) < 1 + α x ¯ + α + β .
Now, we have
1 + α x ¯ + α + β < b α + ( b a β ) x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) 3 b α + b β + a x ¯ 2 + 2 b x ¯ + 2 a α x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) < 0 ,
which is always true.
Also, we have
b α + ( b a β ) x ¯ ( x ¯ + α + β ) ( b + a x ¯ ) < 1 + α x ¯ + α + β b α + b β + a x ¯ 2 + 2 a α x ¯ + 2 a x ¯ β ( x ¯ + α + β ) ( b + a x ¯ ) < 0 ,
which is always true. □
Lemma 1.
Every solution of Equation (1) is bounded.
Proof. 
From Equation (1), we have
x n + 1 = x n x n 1 + α x n + β x n 1 a x n x n 1 + b x n 1 x n x n 1 + β x n 1 a x n x n 1 + b x n 1 min 1 a , β b , x n + 1 = x n x n 1 + α x n + β x n 1 a x n x n 1 + b x n 1 = α x n a x n x n 1 + b x n 1 + x n x n 1 + β x n 1 a x n x n 1 + b x n 1 , x n + 1 α x n a x n x n 1 + max 1 a , β b , x n + 1 α a x n 1 + max 1 a , β b α + β b , 1 a < β b α b β a + 1 a , 1 a > β b α + 1 a , 1 a = β b
which completes the proof. □
Lemma 2.
Equation (1) does not posses a minimal period-two solution.
Proof. 
Period-two solution ( x , y ) , x > 0 , y > 0 , x y , satisfies the following system of algebraic equations
x f ( y , x ) = 0 y f ( x , y ) = 0 .
By replacing the function f , we have
a x 2 y + b x 2 β x x y α y = 0 a x y 2 + b y 2 α x x y β y = 0 .
By subtracting these two equations, we obtain
x y a x y + b x + y + α β = 0 .
If α β , then there does not exist a minimal period-two solution. However, if α < β , we have
a x y + b x + y = β α .
Similarly, by adding, from (10) we obtain
x + y a x y + b x + y α β = 2 x y 1 + b ,
from which
a x y + b x + y α β > 0 .
By using (11), we have
a x y + b x + y α β = β α α β = 2 α < 0 ,
which is a contradiction with (12). ☐

3. Global Attractivity Results

In this section, we prove several global attractivity results in the corresponding parametric space.
We notice that the sign of the partial derivative with respect to the first variable at the equilibrium point depends on the sign of the b a β .
b a β 0 , In this case, the function f is increasing in the first variable and decreasing in the second variable.
Lemma 3.
If α β or α > β , b 1 or β < α α D = 4 a β 2 b 1 2 + β , b > 1 , then the system of algebraic equations
m = f m , M , M = f M , m
has a unique solution m , M = x ¯ , x ¯ .
Proof. 
System (13) is of the form
m = M β + m α + M m M b + a m , M = M α + m β + M m m b + M a ,
that is
M β + m α + M m = m M b + a m ,
M α + m β + M m = M m b + M a .
By subtracting Equations (14) and (15), we obtain
m M a m M α + β = 0 .
If α β , then system (13) has a unique solution m , M = x ¯ , x ¯ .
Suppose that α > β and M m . Then
M = α β a m .
By adding Equations (14) and (15), we have
m + M a m M α β = 2 1 b m M ,
from which, by using (16), we have
β m + M = b 1 m M .
Now, we see that Equation (17) has no positive solutions m , M if b 1 . This means that system (13) has a unique solution m , M = x ¯ , x ¯ in this case.
Assume that b > 1 . By substituting (16) into (17) we obtain quadratic equation
a β m 2 b 1 α β m + β α β = 0 ,
with solutions
m ± = b 1 α β ± α β b 1 2 α β 4 a β 2 2 a β .
It is easy to see that m ± R when α < 4 a β 2 b 1 2 + β = α D , which means that system (13) has a unique solution m , M = x ¯ , x ¯ in this case. Similarly, if α = α D , then m = M = 2 β b 1 = x ¯ .
Finally, if α > α D , then system (13) has more than one positive solution. □
Theorem 7.
Assume that one of the following conditions hold:
(1) 
1 b a β , α > 0 , a > 0 ,
(2) 
b > 1 , b a β , α α D , a > 0 .
Then, the unique Equilibrium (8) of Equation (1) is globally asymptotically stable.
Proof. 
In this case (see the proof of Lemma 1) the invariant interval (and an attracting interval) of Equation (1) is
1 a , β b + α .
Since f ( u , v ) u = b α + ( b a β ) v v ( b + a u ) 2 > 0 , then f is increasing in the first variable and decreasing in the second variable and we can apply Theorem 1. Also, we know that the equilibrium x ¯ is locally asymptotically stable, and consequently the proof will be completed by using Lemma 3 and Theorem 1. □
For some numerical values of parameters we give visual evidence for Theorem 7. (See Figure 1a).
b a β < 0 , Lemma 1 implies that
1 a x n + 1 α + β b .
Now, we have
f ( u , v ) u = b α + ( b a β ) v v ( b + a u ) 2 = 0 , v = v 0 = b α a β b , > 0 , v < v 0 = b α a β b , < 0 , v > v 0 = b α a β b .
Since
f b α a β b , b α a β b = β b ,
we have to consider the following three cases:
(i)
v 0 = b α a β b > β b ,
(ii)
v 0 = b α a β b = β b ,
(iii)
v 0 = b α a β b < β b .
For the case of v 0 = b α a β b > β b we have the following result about global behavior of solutions of Equation (1).
Theorem 8.
Assume that one of the following conditions hold:
(1) 
b 1 , β > b a α > α 0 , a > 0 ,
(2) 
b > 1 , β > b a , α 0 < α α D , a > 0 ,
where α 0 = β a β b b 2 . Then, the unique Equilibrium (8) of Equation (1) is globally asymptotically stable.
Proof. 
Since v 0 = b α a β b > β b , we have that b α a β b β b = b β a β 2 + b 2 α b a β b > 0 , i.e., b β a β 2 + b 2 α > 0 , from which α > α 0 = β a β b b 2 , β > b a .
In this case we have
f : β b , b α a β b 2 β b , b α a β b ,
which implies that the interval β b , b α a β b is an invariant interval. Indeed, since the function f ( u , v ) is continuous, then f ( u , v ) attains its extreme points at the end of closed interval or at the stationary point. Straightforward calculations show that all values
f b α a β b , b α a β b = β b , f β b , b α a β b = β b ,
f b α a β b , β b = b 2 α 2 + b α β b β 2 + a β 3 b β a α + a β b , f β b , β b = β + b α + b β b 2 + a β
are in β b , b α a β b .
Set
g ( x ) = a x 2 + ( b 1 ) x α β .
We know that g ( x ¯ ) = 0 . On the other hand, we have that
g β b = b β a β 2 + b 2 α b 2 < 0 , g b α a β b = a α + a β b b β a β 2 + b 2 α a β b 2 > 0 ,
which implies x ¯ β b , b α a β b .
Also, we know that f ( u , v ) u = b α + ( b a β ) v v ( b + a u ) 2 > 0 , which implies that f is increasing in the first variable and decreasing in the second variable. Since Equation (1) has the unique equilibrium point in invariant interval β b , b α a β b , we can apply Theorem 1. Also, we know that the equilibrium x ¯ is locally asymptotically stable, and consequently the proof will be completed by using Lemma 3 and Theorem 1. □
For some numerical values of parameters we give visual evidence for Theorem 8. (See Figure 1b).
For the case of v 0 = b α a β b < β b the following result holds.
Theorem 9.
Assume that v 0 = b α a β b < β b . Then, the invariant interval of Equation (1) is
b α a β b , β b ,
that is
f : b α a β b , β b 2 b α a β b , β b
and the unique Equilibrium (8) of Equation (1) is globally asymptotically stable.
Proof. 
Since v 0 = b α a β b < β b , then b α a β b β b = b β a β 2 + b 2 α b a β b < 0 , that is b β a β 2 + b 2 α < 0 , from which α < α 0 = β a β b b 2 , β > b a .
First, we prove that the invariant interval is given by b α a β b , β b . Since the function f ( u , v ) is continuous, then this function attains its extreme points at the end of closed interval or at the stationary point. Straightforward calculations show that
f b α a β b , b α a β b = β b , f β b , b α a β b = β b ,
f b α a β b , β b = b 2 α 2 + b α β b β 2 + a β 3 b β a α + a β b , f β b , β b = β + b α + b β b 2 + a β
are in b α a β b , β b .
Now, we prove that the equilibrium point is in interval b α a β b , β b . Set
g ( x ) = a x 2 + ( b 1 ) x α β .
We know that g ( x ¯ ) = 0 . On the other hand, we have
g β b = b β a β 2 + b 2 α b 2 > 0 , g b α a β b = ( a α + a β b ) b β a β 2 + b 2 α ( a β b ) 2 < 0 ,
which shows that x ¯ b α a β b , β b .
Since the function f ( u , v ) is decreasing in both variables, and Equation (1) has the unique equilibrium point in invariant interval b α a β b , β b , we can apply Theorem 2. System of algebraic equations
f ( M , M ) = m , f ( m , m ) = M
becomes
M + α + β b + a M = m , m + α + β b + a m = M .
It is easy to see that this system has a unique solution m = M , which completes the proof. □
For some numerical values of parameters we give visual evidence for Theorem 9. (See Figure 2a).
Remark 1.
Notice that we can prove Theorem 9 by using Theorem 3. In this case, we have
1 a x n + 1 α + β b , n = 0 , 1 , ,
which means that every solution x n n = 1 of Equation (1) is bounded from above and from below by positive constants.
Since
d d x x f x , x = a x 2 + 2 b x + b α + b β b + a x 2 > 0 ,
and f clearly satisfies the conditions 1. and 2. of Theorem 3. Now, by using Theorem 3, we have that every solution x n n = 1 of Equation (1) converges to x ¯ .
Now, consider the case v 0 = b α a β b = β b .
Lemma 4.
Assume that b α a β b = β b , that is α = α 0 = β a β b b 2 , where β > b a . Then, Equation (1) does not posses a minimal period-four solution.
Proof. 
Suppose the opposite, i.e., Equaiton (1) has a minimal period-four solution: x , y , z , t , x , y , z , t , , that is the following equalities hold
z = x y + β a β b b 2 y + β x a x y + b x , t = z y + β a β b b 2 z + β y a z y + b y , x = t z + β a β b b 2 t + β z a t z + b z , y = x t + β a β b b 2 x + β t a x t + b t .
By eliminating z and t we obtain
b x β Φ x , y , a , b , β b 2 x y + β + a β b y β Ω x , y , a , b , β = 0 , b y β Ψ x , y , a , b , β b + a x β a β b b 2 x y + a β b y β + b 2 x β + b 2 y 2 + y x b 4 y + β + b β + a y β = 0 ,
where
Φ x , y , a , b , β = A a β b 3 + B a β b 2 + C a β b + D , Ψ x , y , a , b , β = E a β b 2 + F a β b + G , Ω x , y , a , b , β = H + J
A = β 2 y 2 , B = β y b 1 + a β b y 2 + b + a β x y + b + a β x β , C = x b 3 b + a β + 3 a b β + a b 2 β y 3 + b + a x + a β + b 2 + b 3 + 2 a b β β y 2 + 2 y + β x a β 2 , D = y x 2 b 5 y + β a y + a β + a b y + 2 a b β + b 3 + a b 2 y + a 2 y β , E = y β x + β , F = b a x + b y 2 β + x + β y + β b x , G = b 3 x y y + β + b β + a y β b + a x H = β a a β b a β b y β + b 2 y + β x + 1 + b y 2 J = y x b 4 y + β a + a y b 2 + b + a β + 2 a b β + b 3 .
Since a β b > 0 , we have that Φ x , y , a , b , β > 0 and Ψ x , y , a , b , β > 0 . Therefore, system (18) has a unique solution of the form
x = β b , y = β b ,
z = x y + β a β b b 2 y + β x a x y + b x = β b 2 + β a β b b 2 β b + β β b a β b 2 + b β b = β b , t = z y + β a β b b 2 z + β y a z y + b y = β b 2 + β a β b b 2 β b + β β b a β b 2 + b β b = β b ,
which means that Equaiton (1) has no a minimal period-four solution. □
Theorem 10.
Assume that b α a β b = β b , that is α = α 0 = β a β b b 2 , where β > b a . Then, the unique Equilibrium (8) of Equaiton (1) is globally asymptotically stable. Furthermore, every solution oscillates about the equilibrium point x ¯ with semicycles of length two.
Proof. 
Notice that
x n + 1 β b = x n x n 1 + β a β b b 2 x n + β x n 1 a x n x n 1 + b x n 1 β b = b x n a β b x n 1 b 2 b + a x n β b x n 1 ,
which implies that the length of the semicycle is two. By using Eq.(1), we have
x n + 4 = x n + 3 x n + 2 + β a β b b 2 x n + 3 + β x n + 2 a x n + 3 x n + 2 + b x n + 2 = x n + 3 x n + 1 x n + β a β b b 2 x n + 1 + β x n a x n + 1 x n + b x n + β a β b b 2 x n + 3 + β x n + 1 x n + β a β b b 2 x n + 1 + β x n a x n + 1 x n + b x n a x n + 3 + b x n + 1 x n + β a β b b 2 x n + 1 + β x n a x n + 1 x n + b x n .
After straightforward calculations, we obtain
x n x n + 4 = β + b x n a β b β + x n + 3 x n + 1 + b x n β + x n + 1 b + a x n + 3 b + a x n + 3 a β 2 x n + 1 + b 2 x n x n + 1 + b 2 β x n b β x n + 1 .
From (19) we have
x n x n + 4 x n β b = b a β b β + x n + 3 x n + 1 + b x n β + x n + 1 b + a x n + 3 b + a x n + 3 a β b β x n + 1 + b 2 x n x n + 1 + b 2 β x n > 0 .
On the other hand,
b a β b β + x n + 3 x n + 1 + b x n β + x n + 1 b + a x n + 3 b + a x n + 3 a β b β x n + 1 + b 2 x n x n + 1 + b 2 β x n < 1
if and only if
x n + 3 x n + 1 a β b 2 < 0 ,
which is true.
Therefore, for a β b > 0 , we have that
0 < x n x n + 4 x n β b < 1 ,
and we see that x n and x n + 4 are always on the same side of the equilibrium point x ¯ . Namely, if x n < β b , then
x n < x n + 4 < β b , n N ,
and so, if x n > β b , then
x n > x n + 4 > β b , n N .
Therefore, every sequence x 4 k k = 1 , x 4 k + 1 k = 0 , x 4 k + 2 k = 0 , x 4 k + 3 k = 0 is monotone and bounded. This implies that each of the sequences is convergent. Since, by Lemma 4, Equation (1) has no minimal period-four solutions, we have that
lim k x 4 k = lim k x 4 k + 1 = lim k x 4 k + 2 = lim k x 4 k + 3 = x ¯ ,
which implies that the equilibrium x ¯ is an attractor. It means that x ¯ is globally asymptotically stable. □
For some numerical values of parameters we give visual evidence for Theorem 10. (See Figure 2b)
Remark 2.
Notice that in the case when these four subsequences exist, there is an invariant interval of this form 1 a , a β 2 b 2 , but we can not use any of Theorems 1 and 2 because the map associated with Equation (1) changes its monotonicity with respect to the first variable in this invariant interval.
Remark 3.
Based on our numerical simulations and Theorems 4 and 5, we believe that every solution of Equaiton (1) oscillates about the equilibrium x ¯ with semicycles of length two. (See Figure 1 and Figure 2).
Also, based on our numerical simulations, we give the following conjecture.
Conjecture 1.
The unique positive Equilibrium (8) of Equation (1) is always globally asymptotically stable.

Author Contributions

All three authors have significant contribution to this paper and the final form of this paper is approved by all three authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bektešević, J.; Hadžiabdić, V.; Kalabušić, S.; Mehuljić, M. Global Asymptotic Behavior of Some Quadratic Rational Second-Order Difference Equations. Int. J. Differ. Equ. 2017, 12, 169–183. [Google Scholar]
  2. Garić-Demirović, M.; Kulenović, M.R.S.; Nurkanović, M. Global Dynamics of Certain Homogeneous Second-Order Quadratic Fractional Difference Equations. Sci. World J. 2013, 2013, 10. [Google Scholar] [CrossRef] [PubMed]
  3. Garić-Demirović, M.; Kulenović, M.R.S.; Nurkanović, M. Basins of Attraction of Certain Homogeneous Second Order Quadratic Fractional Difference Equation. J. Concr. Appl. Math. 2015, 13, 35–50. [Google Scholar]
  4. Garić-Demirović, M.; Nurkanović, M.; Nurkanović, Z. Stability, Periodicity and Neimark–Sacker Bifurcation of Certain Homogeneous Fractional Difference Equations. Int. J. Differ. Equ. 2017, 12, 27–53. [Google Scholar]
  5. Hrustić, S.J.; Kulenović, M.R.S.; Nurkanović, M. Global Dynamics and Bifurcations of Certain Seond Order Rational Difference Equation with Quadratic Terms. Qual. Theory Dyn. Syst. 2016, 15, 283. [Google Scholar] [CrossRef]
  6. Jašarević-Hrustić, S.; Nurkanović, Z.; Kulenović, M.R.S.; Pilav, E. Birkhoff normal forms, KAM theory and symmetries for certain second order rational difference equation with quadratic term. Int. J. Differ. Equ. 2015, 10, 181–199. [Google Scholar]
  7. Kent, C.M.; Sedaghat, H. Global attractivity in a quadratic-linear rational difference equation with delay. J. Differ. Equ. Appl. 2009, 15, 913–925. [Google Scholar] [CrossRef]
  8. Kent, C.M.; Sedaghat, H. Global attractivity in a rational delay difference equation with quadratic terms. J. Differ. Equ. Appl. 2011, 17, 457–466. [Google Scholar] [CrossRef]
  9. Kulenović, M.R.S.; Moranjkić, S.; Nurkanović, Z. Naimark-Sacker Bifurcation of Second Order Rational Difference Equation with Quadratic Terms. J. Nonlinear Sci. Appl. 2017, 10, 3477–3489. [Google Scholar] [CrossRef]
  10. Moranjkić, S.; Nurkanović, Z. Local and Global Dynamics of Certain Second-Order Rational Difference Equations Containing Quadratic Terms. Adv. Dyn. Syst. Appl. 2017, 12, 123–157. [Google Scholar]
  11. Sedaghat, H. Global behaviours of rational difference equations of orders two and three with quadratic terms. J. Differ. Equ. Appl. 2009, 15, 215–224. [Google Scholar] [CrossRef]
  12. Kulenović, M.R.S.; Ladas, G. Dynamics of Second Order Rational Difference Equations with Open Problems and Conjectures; Chapman and Hall/CRC: Boca Raton, FL, USA; London, UK, 2001. [Google Scholar]
  13. Enciso, E.G.; Sontag, E.D. Global attractivity, I/O monotone small-gain theorems, and biological delay systems. Discrete Contin. Dyn. Syst. 2006, 14, 549–578. [Google Scholar]
  14. Kulenović, M.R.S.; Merino, O. Discrete Dynamical Systems and Difference Equations with Mathematica; Chapman and Hall/CRC Press: Boca Raton, FL, USA, 2002. [Google Scholar]
  15. Kostrov, Y.; Kudlak, Z. On a Second-Order Rational Difference Equation with a Quadratic Term. Int. J. Differ. Equ. 2016, 11, 179–202. [Google Scholar]
  16. Kulenović, M.R.S.; Nurkanović, M. Asymptotic behavior of a two dimensional linear fractional system of difference equations. Radovi Matematički 2002, 11, 59–78. [Google Scholar]
  17. Kulenović, M.R.S.; Nurkanović, M. Asymptotic Behavior of a Linear Fractional System of Difference Equations. J. Inequal. Appl. 2005, 127–143. [Google Scholar] [CrossRef]
  18. Burgić, D.; Kulenović, M.R.S.; Nurkanović, M. Global Dynamics of a Rational System of Difference Equations in the plane. Commun. Appl. Nonlinear Anal. 2008, 15, 71–84. [Google Scholar]
  19. Moranjkić, S.; Nurkanović, Z. Basins of attraction of certain rational anti-competitive system of difference equations in the plane. Adv. Differ. Equ. 2012, 2012, 153. [Google Scholar] [CrossRef]
  20. Smith, H.L. Non-monotone systems decomposable into monotone systems with negative feedback. J. Math. Biol. 2006, 53, 747–758. [Google Scholar] [CrossRef] [PubMed]
  21. Drymonis, E.; Ladas, G. On the global character of the rational system x n + 1 = α A 1 + B 1 x n + y n , y n + 1 = α 2 + β 2 x n A 2 + B 2 x n + C 2 y n . Sarajevo J. Math. 2012, 8, 293–309. [Google Scholar] [CrossRef]
Figure 1. Numerical simulation of stability for: (a) b a β 0 ; (b) b a β < 0 and v 0 > β b .
Figure 1. Numerical simulation of stability for: (a) b a β 0 ; (b) b a β < 0 and v 0 > β b .
Mathematics 06 00010 g001
Figure 2. Numerical simulation of stability for b a β < 0 (a) v 0 < β b ; (b) v 0 = β b .
Figure 2. Numerical simulation of stability for b a β < 0 (a) v 0 < β b ; (b) v 0 = β b .
Mathematics 06 00010 g002

Share and Cite

MDPI and ACS Style

Kalabušić, S.; Nurkanović, M.; Nurkanović, Z. Global Dynamics of Certain Mix Monotone Difference Equation. Mathematics 2018, 6, 10. https://doi.org/10.3390/math6010010

AMA Style

Kalabušić S, Nurkanović M, Nurkanović Z. Global Dynamics of Certain Mix Monotone Difference Equation. Mathematics. 2018; 6(1):10. https://doi.org/10.3390/math6010010

Chicago/Turabian Style

Kalabušić, Senada, Mehmed Nurkanović, and Zehra Nurkanović. 2018. "Global Dynamics of Certain Mix Monotone Difference Equation" Mathematics 6, no. 1: 10. https://doi.org/10.3390/math6010010

APA Style

Kalabušić, S., Nurkanović, M., & Nurkanović, Z. (2018). Global Dynamics of Certain Mix Monotone Difference Equation. Mathematics, 6(1), 10. https://doi.org/10.3390/math6010010

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

Article Metrics

Back to TopTop