Next Article in Journal
Polar Bear Optimization Algorithm: Meta-Heuristic with Fast Population Movement and Dynamic Birth and Death Mechanism
Next Article in Special Issue
Evaporation and Antievaporation Instabilities
Previous Article in Journal
Towards a New Proposal for the Time Delay in Gravitational Lensing
Previous Article in Special Issue
On Brane Solutions with Intersection Rules Related to Lie Algebras
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvability of the Class of Two-Dimensional Product-Type Systems of Difference Equations of Delay-Type (1, 3, 1, 1)

1
Mathematical Institute of the Serbian Academy of Sciences, Knez Mihailova 36/III, 11000 Beograd, Serbia
2
Operator Theory and Applications Research Group, Department of Mathematics, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
Symmetry 2017, 9(10), 200; https://doi.org/10.3390/sym9100200
Submission received: 1 September 2017 / Revised: 21 September 2017 / Accepted: 21 September 2017 / Published: 25 September 2017
(This article belongs to the Special Issue Symmetry: Feature Papers 2017)

Abstract

:
This paper essentially presents the last and important steps in the study of (practical) solvability of two-dimensional product-type systems of difference equations of the following form z n = α z n k a w n l b , w n = β w n m c z n s d , n N 0 , where k , l , m , s N , a , b , c , d Z , and where α , β and the initial values are complex numbers. It is devoted to the most complex case which has not been considered so far (the case k = l = s = 1 and m = 3 ). Closed form formulas for solutions to the system are found in all possible cases. The structure of the solutions to the system is considered in detail. The following five cases: (1) b = 0 ; (2) c = 0 ; (3) d = 0 ; (4) a c 0 ; (5) a = 0 , b c d 0 , are considered separately. Some of the situations appear for the first time in the literature.

1. Introduction

There has been a growing interest in difference equations and systems of difference equations (see, for example, [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43]). Among several subfields of recent interest we mention here two, on whose intersection is the study in this paper. The first one is the classical subfield of finding solutions to the equations and systems in closed form. Books [6,12,13,14,15,16] contain some old information in the subfield. Some recent ones can be found, for example, in the following papers: [5,27,31,32,33,34,35,36,37,38,39,40,41,42,43] (see also numerous references therein). For some related results and applications of solvable difference equations and systems, see also [1,2,3,4]. The recent interest has been considerably motivated by the fact that some new interesting classes of nonlinear difference equations and systems have been solved by transforming them to known solvable ones. One of our transformations has had some impact in the recent interest. For more information, see, for example, [5,27,34], and the related references therein. Beside the line of investigation, there have been some other several ones which also use some related ideas. The reader can consult the representative paper [34] and find many other related ones in its list of references. Generally speaking, above mentioned lines of investigations use the method of transformation in solving the equations and systems therein. In many of these papers obtained formulas for the solutions to the equations and systems studied therein are used in describing their long-term behavior (for example, in [5,39]).
The second one is the subfield on concrete systems of difference equations. Some of the papers which have had some impact on the growing interest in the subfield are [19,20,21] by Papaschinopoulos and Schinas. One of their main ideas is to consider symmetric systems of difference equations obtained from the following scalar one
x n = f ( x n k , x n l ) , n N 0 ,
where k , l N , k l , for concrete values of function f, that is, to study some concrete systems of difference equations of one of the following forms
x n = f ( x n k , y n l ) , y n = f ( y n k , x n l ) , x n = f ( y n k , x n l ) , y n = f ( x n k , y n l ) , x n = f ( y n k , y n l ) , y n = f ( x n k , x n l ) , n N 0 .
This, among other things, has motivated us to study the solvability of some concrete systems of the form, such as the ones in [5,39,40] (see also the related references therein). The idea naturally evolved into the investigation of more general symmetric systems of difference equations. For example, a few symmetric systems with three variables were studied in [23], while in [22] were studied, among others, the invariants of the following system of nonlinear difference equations
x n + 1 = a + b ( y n + x n 1 + + y n k + 1 + x n k ) y n k 1 y n + 1 = a + b ( x n + y n 1 + + x n k + 1 + y n k ) x n k 1
where k is an odd number and a , b ( 0 , ) , which means that more complex symmetric systems were studied therein. In fact, the study of the invariants of some equations and systems can be regarded as a kind of the study of their solvability, so that paper [22], as well as [21] essentially belong to both subfields.
It is highly expected that the methods and ideas used in the study of symmetric systems of difference equations can produce the same or related results for the systems which are not symmetric, but are close to them. For example, the following max-type system of difference equations
x n + 1 = max { a n y n , b n } y n x n 1 , y n + 1 = max { c n x n , d n } x n y n 1 , n N 0 ,
where a n , b n , c n , d n are sequences of positive numbers and x 1 , x 0 , y 1 , y 0 are positive numbers, is such one, and was studied in [25] (for another max-type system see also [28]). Such system are called close-to-symmetric systems of difference equations and are frequently studied (see, for example, [8,17,24,29,31,33] and the related references therein). In paper [11] was initiated study of cyclic systems of difference equations, which naturally evolved into the study of some close to cyclic systems of difference equations (see, for example, [7,18]).
Studying positive solutions to some classes of equations and systems, such as the ones of the special cases of the following equation
x n = α + x n k p x n l q , n N 0 ,
where α , p , q ( 0 , ) , k , l N , k l (see, for example, [30], as well as [9,10,26] and the related references therein), as well as of the corresponding two-dimensional symmetric systems of difference equations, such as
x n = α + y n k p x n l q , y n = α + x n k p y n l q , n N 0 ,
we came across some product-type equations and systems.
The solvability of the product-type equations and systems in this case is something which should be known to any expert. Namely, if all the initial values are positive a simple inductive argument shows that all other terms are also positive, so that by using the logarithm the equation/system is transformed to a linear one with constant coefficients, which is one of the most known solvable equation/system. If the domain is changed, then several problems occur which prevent using the standard method for solving the equations and systems for the case of positive solutions. Our study of the systems on the complex plane was set off in [39]. It turned out that finding solutions to some related systems is not so simple problem. The next product-type system was solved in [40], but without detailed analysis of the structure of its solutions. The forms of the systems studied in [39,40] are similar, which suggested an investigation of the extensions which include both of them. On the other hand, the occurrence of some multipliers in some cases of the one-dimensional equation in [34] has suggested an investigation of the related systems with additional constant multipliers, which has been done for the first time in [31].
Somewhat later, a detailed analysis has shown that complete lists of formulas can be given for some of concrete systems of difference equations of the following form
z n = α z n k a w n l b , w n = β w n m c z n s d , n N ,
with “small" values of delays k , l , m and s, which means that they are solvable. Since we have studied so far a number of the systems of type (1), to facilitate classification and terminology, from now on we will say that the system is of delay-type ( k , m , l , s ) .
The corresponding lists of formulas for solutions are given first for the systems in [33,43], unlike the systems in above mentioned paper [40] and for the system in [42]. For some systems such as the ones in [35,41] the solutions were obtained more easily, so the analysis was simpler and it was of a different character. Some technical problems in dealing with the systems of the sort in (1) lead us to devising another method for solving them in [32], which has been recently also used in [37]. Recently, we have done, for the first time a detailed analysis of the structure of solutions to a class of product-type systems with an associated polynomial (to the system) of the fourth order in [38], and quite recently in [36].
The main goal of the whole project is to classify solvable product-type systems of difference equations of the form in (1) and present their solutions in closed form in terms of the involved parameters and initial values. Here we continue the project. This paper is a natural continuation of our research in [31,32,33,35,36,37,38,39,40,41,42,43], and essentially presents the last and important steps in the finishing of the project.
Our task here is to show the solvability of system (1) of delay-type ( 1 , 3 , 1 , 1 ) (the case k = l = s = 1 and m = 3 ), that is, of
z n + 1 = α z n a w n b , w n + 1 = β w n 2 c z n d , n N 0 ,
where a , b , c , d Z , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . The case when some of the quantities α , β , z 0 , w 2 , w 1 , w 0 is zero we do not take into consideration because in the case are obtained solutions which are either not defined or trivial, so of not special interest.
The following five cases are considered separately in this paper: (1) b = 0 ; (2) c = 0 ; (3) d = 0 ; (4) a c 0 ; (5) a = 0 , b c d 0 . We would like to point out that the fifth case is not covered with the fourth one, since the condition a = 0 changes the order of an associated polynomial appearing in the study. If k , l Z , then the notation k , l ¯ , denotes the set of all j Z such that k j l , whereas we regard that j = l l 1 c j = 0 for each l Z .

2. Auxiliary Results

The following three lemmas are useful tools in our investigation and have been already used in some of our previous papers devoted to the project on product-type systems. The first one is a consequence of the Langrage formula applied to the functions f s ( t ) = t s , s N (see, for example, [13] or [44], as well as [40] for a proof based on complex analysis).
Lemma 1.
Let
p k ( t ) = c k j = 1 k ( t t j ) ,
c k 0 and t i t j , i j . Then
j = 1 k t j s p k ( t j ) = 0 , 0 s k 2 ,
and
j = 1 k t j k 1 p k ( t j ) = 1 c k .
Further, we need several closed form formulas for some sums which can be found in numerous books (see, for example, [13] or [16]). For a general method for calculating this type of sums consult our recent paper [33] where a recurrent formula for this type of sums is presented, and by using it the sums can be calculated.
Lemma 2.
Let
s n ( i ) ( z ) = j = 1 n j i z j 1 , n N ,
where i N 0 and z C .
Then
s n ( 0 ) ( z ) = 1 z n 1 z ,
s n ( 1 ) ( z ) = 1 ( n + 1 ) z n + n z n + 1 ( 1 z ) 2 ,
s n ( 2 ) ( z ) = 1 + z ( n + 1 ) 2 z n + ( 2 n 2 + 2 n 1 ) z n + 1 n 2 z n + 2 ( 1 z ) 3 ,
s n ( 3 ) ( z ) = n 3 z n ( z 1 ) 3 3 n 2 z n ( z 1 ) 2 + 3 n z n ( z 2 1 ) ( z n 1 ) ( z 2 + 4 z + 1 ) ( 1 z ) 4 ,
for every z C \ { 1 } and n N .
The following lemma describes the nature/type of the zeros of an arbitrary fourth order polynomial equation. The results in the lemma are certainly folklore and were essentially obtained, for example, in [45] (the lemma formulates the results appearing therein in a unified way, although the notation and some quantities are different).
Lemma 3.
Let
P 4 ( t ) = t 4 + b t 3 + c t 2 + d t + e ,
where b , c , d , e are real numbers,
Δ 0 = c 2 3 b d + 12 e , Δ 1 = 2 c 3 9 b c d + 27 b 2 e + 27 d 2 72 c e , Δ = 1 27 ( 4 Δ 0 3 Δ 1 2 ) , P = 8 c 3 b 2 , Q = b 3 + 8 d 4 b c , D = 64 e 16 c 2 + 16 b 2 c 16 b d 3 b 4 .
Then, the following statements hold.
(a) 
If Δ < 0 , then two zeros of P 4 are real and different, and two are complex conjugate.
(b) 
If Δ > 0 , then all the zeros of P 4 are real or none is. More precisely,
1
if P < 0 and D < 0 , then all four zeros of P 4 are real and different;
2
if P > 0 or D > 0 , then there are two pairs of complex conjugate zeros of P 4 .
(c) 
If Δ = 0 , then and only then P 4 has a multiple zero. The following cases can occur:
1
if P < 0 , D < 0 and Δ 0 0 , then two zeros of P 4 are real and equal and two are real and simple;
2
if D > 0 or ( P > 0 and ( D 0 or Q 0 )), then two zeros of P 4 are real and equal and two are complex conjugate;
3
if Δ 0 = 0 and D 0 , there is a triple zero of P 4 and one simple, all real;
4
if D = 0 , then
4.1
if P < 0 there are two double real zeros of P 4 ;
4.2
if P > 0 and Q = 0 there are two double complex conjugate zeros of P 4 ;
4.3
if Δ 0 = 0 , then all four zeros of P 4 are real and equal to b / 4 .

3. Main Results

In this section we state and prove our main results and by using them and some further analysis we get several corollaries. Before this we give a list of first several members of the sequences z n and w n , whose values will be used in the proofs of some of the results.
We have
z 1 = α z 0 a w 0 b , w 1 = β w 2 c z 0 d , z 2 = α ( α z 0 a w 0 b ) a ( β w 2 c z 0 d ) b = α 1 + a β b z 0 a 2 + b d w 2 b c w 0 a b , w 2 = β w 1 c z 1 d = β w 1 c ( α z 0 a w 0 b ) d = α d β z 0 a d w 1 c w 0 b d z 3 = α z 2 a w 2 b = α ( α 1 + a β b z 0 a 2 + b d w 2 b c w 0 a b ) a ( α d β z 0 a d w 1 c w 0 b d ) b = α 1 + a + a 2 + b d β b + a b z 0 a 3 + 2 a b d w 2 a b c w 1 b c w 0 b ( a 2 + b d ) .
As we have already mentioned we will consider the following five cases separately: (1) b = 0 ; (2) c = 0 ; (3) d = 0 ; (4) a c 0 ; (5) a = 0 , b c d 0 . Hence, we will prove five results on the solvability and by further analysis we will get several consequences from them.
Theorem 1.
Assume that a , c , d Z , b = 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then system (2) is solvable in closed form.
Proof. 
Since b = 0 , we have
z n + 1 = α z n a , w n + 1 = β w n 2 c z n d , n N 0 ,
from which it follows that
z n = α j = 0 n 1 a j z 0 a n ,
for n N .
From (10) we have that
z n = α 1 a n 1 a z 0 a n , n N ,
when a 1 , and
z n = α n z 0 , n N ,
when a = 1 .
Using (10) in the second equation in (9), is obtained
w n = α d j = 0 n 2 a j β z 0 d a n 1 w n 3 c , n 2 ,
which is equivalent to
w 3 m + i = α d j = 0 3 m + i 2 a j β z 0 d a 3 m + i 1 w 3 ( m 1 ) + i c ,
for m N and i = 1 , 0 , 1 .
From (13) and by induction it is proved that
w 3 m + i = α d l = 0 k 1 c l j = 0 3 ( m l ) + i 2 a j β l = 0 k 1 c l z 0 d l = 0 k 1 c l a 3 ( m l ) + i 1 w 3 ( m k ) + i c k ,
for m k and i = 1 , 0 , 1 , from which for k = m is obtained
w 3 m + i = α d l = 0 m 1 c l j = 0 3 ( m l ) + i 2 a j β l = 0 m 1 c l z 0 d l = 0 m 1 c l a 3 ( m l ) + i 1 w i c m ,
for m N and i = 1 , 0 , 1 .
From (8) and (14), we have
w 3 m 1 = α d l = 0 m 1 c l j = 0 3 ( m l 1 ) a j β l = 0 m 1 c l z 0 d l = 0 m 1 c l a 3 ( m l ) 2 w 1 c m ,
w 3 m = α d l = 0 m 1 c l j = 0 3 ( m l ) 2 a j β l = 0 m 1 c l z 0 d l = 0 m 1 c l a 3 ( m l ) 1 w 0 c m ,
w 3 m + 1 = α d l = 0 m 1 c l j = 0 3 ( m l ) 1 a j β l = 0 m 1 c l z 0 d l = 0 m 1 c l a 3 ( m l ) ( β w 2 c z 0 d ) c m = α d l = 0 m 1 c l j = 0 3 ( m l ) 1 a j β l = 0 m c l z 0 d l = 0 m c l a 3 ( m l ) w 2 c m + 1 ,
for m N .
Now we use the formulas in (15)–(17) in five subcases separately.
Case a 1 c a 3 . From (15)–(17) and by Lemma 2, we obtain
w 3 m 1 = α d l = 0 m 1 c l 1 a 3 ( m l ) 2 1 a β 1 c m 1 c z 0 a d a 3 m c m a 3 c w 1 c m = α d 1 a 1 c m 1 c a a 3 m c m a 3 c β 1 c m 1 c z 0 a d a 3 m c m a 3 c w 1 c m = α d ( a 3 c + ( 1 a ) ( a + a 2 + c ) c m + ( c 1 ) a 3 m + 1 ) ( 1 a ) ( 1 c ) ( a 3 c ) β 1 c m 1 c z 0 a d a 3 m c m a 3 c w 1 c m ,
w 3 m = α d l = 0 m 1 c l 1 a 3 ( m l ) 1 1 a β 1 c m 1 c z 0 d a 2 a 3 m c m a 3 c w 0 c m = α d 1 a 1 c m 1 c a 2 a 3 m c m a 3 c β 1 c m 1 c z 0 d a 2 a 3 m c m a 3 c w 0 c m = α d ( a 3 c + ( 1 a ) ( c + c a + a 2 ) c m + ( c 1 ) a 3 m + 2 ) ( 1 a ) ( 1 c ) ( a 3 c ) β 1 c m 1 c z 0 d a 2 a 3 m c m a 3 c w 0 c m ,
w 3 m + 1 = α d l = 0 m 1 c l 1 a 3 ( m l ) 1 a β 1 c m + 1 1 c z 0 d a 3 m + 3 c m + 1 a 3 c w 2 c m + 1 = α d 1 a 1 c m 1 c a 3 a 3 m c m a 3 c β 1 c m + 1 1 c z 0 d a 3 m + 3 c m + 1 a 3 c w 2 c m + 1 = α d ( a 3 c + ( 1 a 3 ) c m + 1 + ( c 1 ) a 3 m + 3 ) ( 1 a ) ( 1 c ) ( a 3 c ) β 1 c m + 1 1 c z 0 d a 3 m + 3 c m + 1 a 3 c w 2 c m + 1 ,
for m N .
Case a 1 c = a 3 . From (15)–(17) and by Lemma 2, we obtain
w 3 m 1 = α d l = 0 m 1 a 3 l j = 0 3 ( m l 1 ) a j β l = 0 m 1 a 3 l z 0 d l = 0 m 1 a 3 l a 3 ( m l ) 2 w 1 a 3 m = α d l = 0 m 1 a 3 l 1 a 3 ( m l ) 2 1 a β 1 a 3 m 1 a 3 z 0 d m a 3 m 2 w 1 a 3 m = α d 1 a 1 a 3 m 1 a 3 m a 3 m 2 β 1 a 3 m 1 a 3 z 0 d m a 3 m 2 w 1 a 3 m = α d ( m a 3 m + 1 a 3 m m a 3 m 2 + 1 ) ( 1 a ) ( 1 a 3 ) β 1 a 3 m 1 a 3 z 0 d m a 3 m 2 w 1 a 3 m ,
w 3 m = α d l = 0 m 1 a 3 l j = 0 3 ( m l ) 2 a j β l = 0 m 1 a 3 l z 0 d l = 0 m 1 a 3 l a 3 ( m l ) 1 w 0 a 3 m = α d l = 0 m 1 a 3 l 1 a 3 ( m l ) 1 1 a β 1 a 3 m 1 a 3 z 0 d m a 3 m 1 w 0 a 3 m = α d 1 a 1 a 3 m 1 a 3 m a 3 m 1 β 1 a 3 m 1 a 3 z 0 m d a 3 m 1 w 0 a 3 m = α d ( m a 3 m + 2 a 3 m m a 3 m 1 + 1 ) ( 1 a ) ( 1 a 3 ) β 1 a 3 m 1 a 3 z 0 d m a 3 m 1 w 0 a 3 m ,
w 3 m + 1 = α d l = 0 m 1 a 3 l j = 0 3 ( m l ) 1 a j β l = 0 m a 3 l z 0 d l = 0 m a 3 l a 3 ( m l ) w 2 a 3 m + 3 = α d l = 0 m 1 a 3 l 1 a 3 ( m l ) 1 a β 1 a 3 m + 3 1 a 3 z 0 d ( m + 1 ) a 3 m w 2 a 3 m + 3 = α d 1 a 1 a 3 m 1 a 3 m a 3 m β 1 a 3 m + 3 1 a 3 z 0 d ( m + 1 ) a 3 m w 2 a 3 m + 3 = α d ( m a 3 m + 3 ( m + 1 ) a 3 m + 1 ) ( 1 a ) ( 1 a 3 ) β 1 a 3 m + 3 1 a 3 z 0 d ( m + 1 ) a 3 m w 2 a 3 m + 3 ,
for m N .
Case a 1 = c . From (15)–(17) and by Lemma 2, we obtain
w 3 m 1 = α d l = 0 m 1 j = 0 3 ( m l 1 ) a j β l = 0 m 1 1 z 0 d l = 0 m 1 a 3 ( m l ) 2 w 1 = α d l = 0 m 1 1 a 3 ( m l ) 2 1 a β m z 0 a d a 3 m 1 a 3 1 w 1 = α d 1 a m a a 3 m 1 a 3 1 β m z 0 a d a 3 m 1 a 3 1 w 1 = α d ( a 3 m + 1 + m ( 1 a 3 ) a ) ( 1 a ) ( 1 a 3 ) β m z 0 a d a 3 m 1 a 3 1 w 1 ,
w 3 m = α d l = 0 m 1 j = 0 3 ( m l ) 2 a j β l = 0 m 1 1 z 0 d l = 0 m 1 a 3 ( m l ) 1 w 0 = α d l = 0 m 1 1 a 3 ( m l ) 1 1 a β m z 0 d a 2 a 3 m 1 a 3 1 w 0 = α d 1 a m a 2 a 3 m 1 a 3 1 β m z 0 d a 2 a 3 m 1 a 3 1 w 0 = α d ( a 3 m + 2 + m ( 1 a 3 ) a 2 ) ( 1 a ) ( 1 a 3 ) β m z 0 d a 2 a 3 m 1 a 3 1 w 0 ,
w 3 m + 1 = α d l = 0 m 1 j = 0 3 ( m l ) 1 a j β l = 0 m 1 z 0 d l = 0 m a 3 ( m l ) w 2 = α d l = 0 m 1 1 a 3 ( m l ) 1 a β m + 1 z 0 d a 3 m + 3 1 a 3 1 w 2 = α d 1 a m a 3 a 3 m 1 a 3 1 β m + 1 z 0 d a 3 m + 3 1 a 3 1 w 2 = α d ( a 3 m + 3 ( m + 1 ) a 3 + m ) ( 1 a ) ( 1 a 3 ) β m + 1 z 0 d a 3 m + 3 1 a 3 1 w 2 ,
for m N .
Case a = 1 c . From (15)–(17) and by Lemma 2, we obtain
w 3 m 1 = α d l = 0 m 1 c l j = 0 3 ( m l 1 ) 1 β l = 0 m 1 c l z 0 d l = 0 m 1 c l w 1 c m = α d l = 0 m 1 ( 3 ( m l 1 ) + 1 ) c l β 1 c m 1 c z 0 d 1 c m 1 c w 1 c m = α d ( 3 m 2 ) 1 c m 1 c 3 c 1 m c m 1 + ( m 1 ) c m ( 1 c ) 2 β 1 c m 1 c z 0 d 1 c m 1 c w 1 c m = α d ( 3 m 2 ( 3 m + 1 ) c + 2 c m + c m + 1 ) ( 1 c ) 2 β 1 c m 1 c z 0 d 1 c m 1 c w 1 c m ,
w 3 m = α d l = 0 m 1 c l j = 0 3 ( m l ) 2 1 β l = 0 m 1 c l z 0 d l = 0 m 1 c l w 0 c m = α d l = 0 m 1 ( 3 ( m l ) 1 ) c l β 1 c m 1 c z 0 d 1 c m 1 c w 0 c m = α d ( 3 m 1 ) 1 c m 1 c 3 c 1 m c m 1 + ( m 1 ) c m ( 1 c ) 2 β 1 c m 1 c z 0 d 1 c m 1 c w 0 c m = α d ( 3 m 1 ( 3 m + 2 ) c + c m + 2 c m + 1 ) ( 1 c ) 2 β 1 c m 1 c z 0 d 1 c m 1 c w 0 c m ,
w 3 m + 1 = α d l = 0 m 1 c l j = 0 3 ( m l ) 1 1 β l = 0 m c l z 0 d l = 0 m c l w 2 c m + 1 = α d l = 0 m 1 3 ( m l ) c l β 1 c m + 1 1 c z 0 d 1 c m + 1 1 c w 2 c m + 1 = α d 3 m 1 c m 1 c 3 c 1 m c m 1 + ( m 1 ) c m ( 1 c ) 2 β 1 c m + 1 1 c z 0 d 1 c m + 1 1 c w 2 c m + 1 = α 3 d ( m ( m + 1 ) c + c m + 1 ) ( 1 c ) 2 β 1 c m + 1 1 c z 0 d 1 c m + 1 1 c w 2 c m + 1 ,
for m N .
Case a = c = 1 . From (15)–(17), we obtain
w 3 m 1 = α d l = 0 m 1 j = 0 3 ( m l 1 ) 1 β l = 0 m 1 1 z 0 d l = 0 m 1 1 w 1 = α d l = 0 m 1 ( 3 ( m l 1 ) + 1 ) β m z 0 d m w 1 = α d m ( 3 m 1 ) 2 β m z 0 d m w 1 ,
w 3 m = α d l = 0 m 1 j = 0 3 ( m l ) 2 1 β l = 0 m 1 1 z 0 d l = 0 m 1 1 w 0 = α d l = 0 m 1 ( 3 ( m l ) 1 ) β m z 0 d m w 0 = α d m ( 3 m + 1 ) 2 β m z 0 d m w 0 ,
w 3 m + 1 = α d l = 0 m 1 j = 0 3 ( m l ) 1 1 β m + 1 z 0 d ( m + 1 ) w 2 = α 3 d l = 0 m 1 ( m l ) β m + 1 z 0 d ( m + 1 ) w 2 = α 3 d m ( m + 1 ) 2 β m + 1 z 0 d ( m + 1 ) w 2 ,
for m N .
From (11), (12), (18)–(32) the theorem follows. ☐
The following corollary follows from Theorem 1.
Corollary 1.
Assume that a , c , d Z , b = 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then the following statements are true.
(a) 
If a 1 c a 3 , then the general solution to system (2) is given by formulas (11), (18)–(20).
(b) 
If a 1 c = a 3 , then the general solution to system (2) is given by formulas (11), (21)–(23).
(c) 
If a 1 = c , then the general solution to system (2) is given by formulas (11), (24)–(26).
(d) 
If a = 1 c , then the general solution to system (2) is given by formulas (12), (27)–(29).
(e) 
If a = c = 1 , then the general solution to system (2) is given by formulas (12), (30)–(32).
Theorem 2.
Assume that a , b , c Z , d = 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then system (2) is solvable in closed form.
Proof. 
In this case we have
z n + 1 = α z n a w n b , w n + 1 = β w n 2 c , n N 0 .
Hence
w 3 n + i = β w 3 ( n 1 ) + i c ,
for n N and i = 2 , 0 ¯ , and consequently
w 3 n + i = β j = 0 n 1 c j w i c n ,
for n N and i = 2 , 0 ¯ .
When c 1 relation (35) implies
w 3 n + i = β 1 c n 1 c w i c n ,
for n N and i = 2 , 0 ¯ , and
w 3 n + i = β n w i ,
for n N and i = 2 , 0 ¯ , when c = 1 .
Further, we have
z n = α ( α z n 2 a w n 2 b ) a w n 1 b = α 1 + a z n 2 a 2 w n 2 b a w n 1 b = α 1 + a ( α z n 3 a w n 3 b ) a 2 w n 2 b a w n 1 b = α 1 + a + a 2 w n 3 b a 2 w n 2 b a w n 1 b z n 3 a 3 ,
for n 3 , which can be written as
z 3 n = α 1 + a + a 2 w 3 n 3 b a 2 w 3 n 2 b a w 3 n 1 b z 3 ( n 1 ) a 3
z 3 n + 1 = α 1 + a + a 2 w 3 n 2 b a 2 w 3 n 1 b a w 3 n b z 3 ( n 1 ) + 1 a 3
z 3 n + 2 = α 1 + a + a 2 w 3 n 1 b a 2 w 3 n b a w 3 n + 1 b z 3 ( n 1 ) + 2 a 3 ,
for n N .
Employing the method of induction we get
z 3 n = z 0 a 3 n j = 0 n 1 α 1 + a + a 2 w 3 j b a 2 w 3 j + 1 b a w 3 j + 2 b a 3 ( n j 1 )
z 3 n + 1 = z 1 a 3 n j = 0 n 1 α 1 + a + a 2 w 3 j + 1 b a 2 w 3 j + 2 b a w 3 j + 3 b a 3 ( n j 1 )
z 3 n + 2 = z 2 a 3 n j = 0 n 1 α 1 + a + a 2 w 3 j + 2 b a 2 w 3 j + 3 b a w 3 j + 4 b a 3 ( n j 1 ) ,
for n N 0 .
Using formula (35), along with the first equation in (33) with n = 0 , 1 , into (41)–(43), we get
z 3 n = z 0 a 3 n j = 0 n 1 α 1 + a + a 2 β b a 2 l = 0 j 1 c l w 0 b a 2 c j β b a l = 0 j c l w 2 b a c j + 1 β b l = 0 j c l w 1 b c j + 1 a 3 ( n j 1 ) = α j = 0 3 n 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j 1 c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l × z 0 a 3 n w 2 b a c j = 0 n 1 c j a 3 ( n j 1 ) w 1 b c j = 0 n 1 c j a 3 ( n j 1 ) w 0 b a 2 j = 0 n 1 c j a 3 ( n j 1 ) ,
z 3 n + 1 = ( α z 0 a w 0 b ) a 3 n j = 0 n 1 α 1 + a + a 2 β b a 2 l = 0 j c l w 2 b a 2 c j + 1 β b a l = 0 j c l w 1 b a c j + 1 β b l = 0 j c l w 0 b c j + 1 a 3 ( n j 1 ) = α j = 0 3 n a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l × z 0 a 3 n + 1 w 2 b a 2 c j = 0 n 1 c j a 3 ( n j 1 ) w 1 b a c j = 0 n 1 c j a 3 ( n j 1 ) w 0 b j = 0 n c j a 3 ( n j ) ,
z 3 n + 2 = ( α 1 + a β b z 0 a 2 w 2 b c w 0 a b ) a 3 n × j = 0 n 1 α 1 + a + a 2 β b a 2 l = 0 j c l w 1 b a 2 c j + 1 β b a l = 0 j c l w 0 b a c j + 1 β b l = 0 j + 1 c l w 2 b c j + 2 a 3 ( n j 1 ) = α j = 0 3 n + 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n a 3 ( n j ) l = 0 j c l × z 0 a 3 n + 2 w 2 b c j = 0 n c j a 3 ( n j ) w 1 b a 2 c j = 0 n 1 c j a 3 ( n j 1 ) w 0 a b j = 0 n c j a 3 ( n j ) ,
for n N 0 .
As in the proof of the previous theorem, there are five subcases, depending on the values of parameters a and c, for which we get closed form formulas for solutions to system (2).
Case a 1 c a 3 . From (44)–(46) and by Lemma 2, we have
z 3 n = α j = 0 3 n 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j 1 c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l × z 0 a 3 n w 2 b a c j = 0 n 1 c j a 3 ( n j 1 ) w 1 b c j = 0 n 1 c j a 3 ( n j 1 ) w 0 b a 2 j = 0 n 1 c j a 3 ( n j 1 ) = α 1 a 3 n 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 c j 1 c + b a j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c + b j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c × z 0 a 3 n w 2 b a c a 3 n c n a 3 c w 1 b c a 3 n c n a 3 c w 0 b a 2 a 3 n c n a 3 c = α 1 a 3 n 1 a β b ( a 2 + a + 1 ) 1 c 1 a 3 n 1 a 3 b ( a 2 + a c + c ) 1 c a 3 n c n a 3 c z 0 a 3 n w 2 b a c a 3 n c n a 3 c w 1 b c a 3 n c n a 3 c w 0 b a 2 a 3 n c n a 3 c = α 1 a 3 n 1 a β b ( a 3 c + ( c 1 ) a 3 n + 2 ( a 2 + a c + c ) ( a 1 ) c n ) ( 1 c ) ( a 3 c ) ( 1 a ) z 0 a 3 n w 2 b a c a 3 n c n a 3 c w 1 b c a 3 n c n a 3 c w 0 b a 2 a 3 n c n a 3 c ,
z 3 n + 1 = α j = 0 3 n a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l × z 0 a 3 n + 1 w 2 b a 2 c j = 0 n 1 c j a 3 ( n j 1 ) w 1 b a c j = 0 n 1 c j a 3 ( n j 1 ) w 0 b j = 0 n c j a 3 ( n j ) = α 1 a 3 n + 1 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c + b a j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c + b j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c × z 0 a 3 n + 1 w 2 b a 2 c a 3 n c n a 3 c w 1 b a c a 3 n c n a 3 c w 0 b a 3 n + 3 c n + 1 a 3 c = α 1 a 3 n + 1 1 a β b ( a 2 + a + 1 ) 1 c 1 a 3 n 1 a 3 b c ( a 2 + a + 1 ) 1 c a 3 n c n a 3 c z 0 a 3 n + 1 w 2 b a 2 c a 3 n c n a 3 c w 1 b a c a 3 n c n a 3 c w 0 b a 3 n + 3 c n + 1 a 3 c = α 1 a 3 n + 1 1 a β b ( a 3 c + ( c 1 ) a 3 n + 3 + ( 1 a 3 ) c n + 1 ) ( 1 c ) ( 1 a ) ( a 3 c ) z 0 a 3 n + 1 w 2 b a 2 c a 3 n c n a 3 c w 1 b a c a 3 n c n a 3 c w 0 b a 3 n + 3 c n + 1 a 3 c ,
z 3 n + 2 = α j = 0 3 n + 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j c l + b j = 0 n a 3 ( n j ) l = 0 j c l × z 0 a 3 n + 2 w 2 b c j = 0 n c j a 3 ( n j ) w 1 b a 2 c j = 0 n 1 c j a 3 ( n j 1 ) w 0 a b j = 0 n c j a 3 ( n j ) = α 1 a 3 n + 2 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c + b a j = 0 n 1 a 3 ( n j 1 ) 1 c j + 1 1 c + b j = 0 n a 3 ( n j ) 1 c j + 1 1 c × z 0 a 3 n + 2 w 2 b c a 3 n + 3 c n + 1 a 3 c w 1 b a 2 c a 3 n c n a 3 c w 0 a b a 3 n + 3 c n + 1 a 3 c = α 1 a 3 n + 2 1 a β b ( a 2 + a ) 1 c 1 a 3 n 1 a 3 + b ( 1 a 3 n + 3 ) ( 1 c ) ( 1 a 3 ) b c ( a 2 + a ) 1 c a 3 n c n a 3 c b c 1 c a 3 n + 3 c n + 1 a 3 c × z 0 a 3 n + 2 w 2 b c a 3 n + 3 c n + 1 a 3 c w 1 b a 2 c a 3 n c n a 3 c w 0 a b a 3 n + 3 c n + 1 a 3 c = α 1 a 3 n + 2 1 a β b ( a 3 c + ( c 1 ) a 3 n + 4 + ( 1 a ) ( a 2 + a + c ) c n + 1 ) ( 1 c ) ( 1 a ) ( a 3 c ) z 0 a 3 n + 2 × w 2 b c a 3 n + 3 c n + 1 a 3 c w 1 b a 2 c a 3 n c n a 3 c w 0 a b a 3 n + 3 c n + 1 a 3 c ,
for n N 0 .
Case a 1 c = a 3 . From (44)–(46) and by Lemma 2, we have
z 3 n = α j = 0 3 n 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j 1 a 3 l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l × z 0 a 3 n w 2 b a 4 j = 0 n 1 a 3 j a 3 ( n j 1 ) w 1 b a 3 j = 0 n 1 a 3 j a 3 ( n j 1 ) w 0 b a 2 j = 0 n 1 a 3 j a 3 ( n j 1 ) = α 1 a 3 n 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j 1 a 3 + b a j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 + b j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 × z 0 a 3 n w 2 b n a 3 n + 1 w 1 b n a 3 n w 0 b n a 3 n 1 = α 1 a 3 n 1 a β b a 2 1 a 3 1 a 3 n 1 a 3 n a 3 n 3 + b a 1 a 3 1 a 3 n 1 a 3 n a 3 n + b 1 a 3 1 a 3 n 1 a 3 n a 3 n × z 0 a 3 n w 2 b n a 3 n + 1 w 1 b n a 3 n w 0 b n a 3 n 1 = α 1 a 3 n 1 a β b n a 3 n + 2 a 3 n n a 3 n 1 + 1 ( 1 a ) ( 1 a 3 ) z 0 a 3 n w 2 b n a 3 n + 1 w 1 b n a 3 n w 0 b n a 3 n 1 ,
z 3 n + 1 = α j = 0 3 n a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l + b j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l × z 0 a 3 n + 1 w 2 b a 5 j = 0 n 1 a 3 j a 3 ( n j 1 ) w 1 b a 4 j = 0 n 1 a 3 j a 3 ( n j 1 ) w 0 b j = 0 n a 3 j a 3 ( n j ) = α 1 a 3 n + 1 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 + b a j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 + b j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 × z 0 a 3 n + 1 w 2 b n a 3 n + 2 w 1 b n a 3 n + 1 w 0 b ( n + 1 ) a 3 n = α 1 a 3 n + 1 1 a β b a 2 1 a 3 1 a 3 n 1 a 3 n a 3 n + b a 1 a 3 1 a 3 n 1 a 3 n a 3 n + b 1 a 3 1 a 3 n 1 a 3 n a 3 n × z 0 a 3 n + 1 w 2 b n a 3 n + 2 w 1 b n a 3 n + 1 w 0 b ( n + 1 ) a 3 n = α 1 a 3 n + 1 1 a β b ( n a 3 n + 3 ( n + 1 ) a 3 n + 1 ) ( 1 a ) ( 1 a 3 ) z 0 a 3 n + 1 w 2 b n a 3 n + 2 w 1 b n a 3 n + 1 w 0 b ( n + 1 ) a 3 n ,
z 3 n + 2 = α j = 0 3 n + 1 a j β b a 2 j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l + b a j = 0 n 1 a 3 ( n j 1 ) l = 0 j a 3 l + b j = 0 n a 3 ( n j ) l = 0 j a 3 l × z 0 a 3 n + 2 w 2 b a 3 j = 0 n a 3 j a 3 ( n j ) w 1 b a 5 j = 0 n 1 a 3 j a 3 ( n j 1 ) w 0 a b j = 0 n a 3 j a 3 ( n j ) = α 1 a 3 n + 2 1 a β b a 2 j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 + b a j = 0 n 1 a 3 ( n j 1 ) 1 a 3 j + 3 1 a 3 + b j = 0 n a 3 ( n j ) 1 a 3 j + 3 1 a 3 × z 0 a 3 n + 2 w 2 b ( n + 1 ) a 3 n + 3 w 1 b n a 3 n + 2 w 0 b ( n + 1 ) a 3 n + 1 = α 1 a 3 n + 2 1 a β b a 2 1 a 3 1 a 3 n 1 a 3 n a 3 n + b a 1 a 3 1 a 3 n 1 a 3 n a 3 n + b 1 a 3 1 a 3 n + 3 1 a 3 ( n + 1 ) a 3 n + 3 × z 0 a 3 n + 2 w 2 b ( n + 1 ) a 3 n + 3 w 1 b n a 3 n + 2 w 0 b ( n + 1 ) a 3 n + 1 = α 1 a 3 n + 2 1 a β b ( ( n + 1 ) a 3 n + 4 a 3 n + 3 ( n + 1 ) a 3 n + 1 + 1 ) ( 1 a ) ( 1 a 3 ) z 0 a 3 n + 2 w 2 b ( n + 1 ) a 3 n + 3 w 1 b n a 3 n + 2 w 0 b ( n + 1 ) a 3 n + 1 ,
for n N 0 .
Case If a 1 = c . From (44)–(46) and by Lemma 2, we have
z 3 n = α j = 0 3 n 1 a j β b a 2 j = 0 n 1 j a 3 ( n j 1 ) + b a j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) + b j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) × z 0 a 3 n w 2 b a j = 0 n 1 a 3 ( n j 1 ) w 1 b j = 0 n 1 a 3 ( n j 1 ) w 0 b a 2 j = 0 n 1 a 3 ( n j 1 ) = α 1 a 3 n 1 a β b a 2 a 3 n n a 3 + n 1 ( a 3 1 ) 2 + b a a 3 n + 3 ( n + 1 ) a 3 + n ( a 3 1 ) 2 + b a 3 n + 3 ( n + 1 ) a 3 + n ( a 3 1 ) 2 × z 0 a 3 n w 2 b a 1 a 3 n 1 a 3 w 1 b 1 a 3 n 1 a 3 w 0 b a 2 1 a 3 n 1 a 3 = α 1 a 3 n 1 a β b ( a 2 + a + 1 ) a 3 n + 2 n a 5 ( n + 1 ) a 4 ( n + 1 ) a 3 + ( n 1 ) a 2 + n a + n ( a 3 1 ) 2 × z 0 a 3 n w 2 b a 1 a 3 n 1 a 3 w 1 b 1 a 3 n 1 a 3 w 0 b a 2 1 a 3 n 1 a 3 = α 1 a 3 n 1 a β b a 3 n + 2 n a 3 + n a 2 ( a 1 ) ( a 3 1 ) z 0 a 3 n w 2 b a 1 a 3 n 1 a 3 w 1 b 1 a 3 n 1 a 3 w 0 b a 2 1 a 3 n 1 a 3 ,
z 3 n + 1 = α j = 0 3 n a j β b a 2 j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) + b a j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) + b j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) × z 0 a 3 n + 1 w 2 b a 2 j = 0 n 1 a 3 ( n j 1 ) w 1 b a j = 0 n 1 a 3 ( n j 1 ) w 0 b j = 0 n a 3 ( n j ) = α 1 a 3 n + 1 1 a β b a 3 n + 3 ( n + 1 ) a 3 + n ( a 1 ) ( a 3 1 ) z 0 a 3 n + 1 w 2 b a 2 1 a 3 n 1 a 3 w 1 b a 1 a 3 n 1 a 3 w 0 b 1 a 3 n + 3 1 a 3 ,
z 3 n + 2 = α j = 0 3 n + 1 a j β b a 2 j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) + b a j = 0 n 1 ( j + 1 ) a 3 ( n j 1 ) + b j = 0 n ( j + 1 ) a 3 ( n j ) × z 0 a 3 n + 2 w 2 b j = 0 n a 3 ( n j ) w 1 b a 2 j = 0 n 1 a 3 ( n j 1 ) w 0 a b j = 0 n a 3 ( n j ) = α 1 a 3 n + 2 1 a β b a 2 a 3 n + 3 ( n + 1 ) a 3 + n ( a 3 1 ) 2 + b a a 3 n + 3 ( n + 1 ) a 3 + n ( a 3 1 ) 2 + b a 3 n + 6 ( n + 2 ) a 3 + n + 1 ( a 3 1 ) 2 × z 0 a 3 n + 2 w 2 b 1 a 3 n + 3 1 a 3 w 1 b a 2 1 a 3 n 1 a 3 w 0 a b 1 a 3 n + 3 1 a 3 = α 1 a 3 n + 2 1 a β b a 3 n + 4 + ( 1 a 3 ) n ( a 3 + a 1 ) ( a 1 ) ( a 3 1 ) z 0 a 3 n + 2 w 2 b 1 a 3 n + 3 1 a 3 w 1 b a 2 1 a 3 n 1 a 3 w 0 a b 1 a 3 n + 3 1 a 3 ,
for n N 0 .
Case If a = 1 c . From (44)–(46) and by Lemma 2, we have
z 3 n = α 3 n β b j = 0 n 1 l = 0 j 1 c l + b j = 0 n 1 l = 0 j c l + b j = 0 n 1 l = 0 j c l z 0 w 2 b c j = 0 n 1 c j w 1 b c j = 0 n 1 c j w 0 b j = 0 n 1 c j = α 3 n β b j = 0 n 1 1 c j 1 c + b j = 0 n 1 1 c j + 1 1 c + b j = 0 n 1 1 c j + 1 1 c z 0 w 2 b c 1 c n 1 c w 1 b c 1 c n 1 c w 0 b 1 c n 1 c = α 3 n β b ( 2 c n + 1 + c n 3 ( c 1 ) n 2 c 1 ) ( 1 c ) 2 z 0 w 2 b c 1 c n 1 c w 1 b c 1 c n 1 c w 0 b 1 c n 1 c ,
z 3 n + 1 = α 3 n + 1 β b j = 0 n 1 l = 0 j c l + b j = 0 n 1 l = 0 j c l + b j = 0 n 1 l = 0 j c l z 0 w 2 b c j = 0 n 1 c j w 1 b c j = 0 n 1 c j w 0 b j = 0 n c j = α 3 n + 1 β 3 b j = 0 n 1 1 c j + 1 1 c z 0 w 2 b c 1 c n 1 c w 1 b c 1 c n 1 c w 0 b 1 c n + 1 1 c = α 3 n + 1 β 3 b c n + 1 ( c 1 ) n c ( 1 c ) 2 z 0 w 2 b c 1 c n 1 c w 1 b c 1 c n 1 c w 0 b 1 c n + 1 1 c ,
z 3 n + 2 = α 3 n + 2 β b j = 0 n 1 l = 0 j c l + b j = 0 n 1 l = 0 j c l + b j = 0 n l = 0 j c l × z 0 w 2 b c j = 0 n c j w 1 b c j = 0 n 1 c j w 0 b j = 0 n c j = α 3 n + 2 β b j = 0 n 1 1 c j + 1 1 c + b j = 0 n 1 1 c j + 1 1 c + b j = 0 n 1 c j + 1 1 c × z 0 w 2 b c 1 c n + 1 1 c w 1 b c 1 c n 1 c w 0 b 1 c n + 1 1 c = α 3 n + 2 β b ( c n + 2 + 2 c n + 1 3 ( c 1 ) n + 1 4 c ) ( 1 c ) 2 z 0 w 2 b c 1 c n + 1 1 c w 1 b c 1 c n 1 c w 0 b 1 c n + 1 1 c ,
for n N 0 .
Case If a = c = 1 . From (44)–(46), we have
z 3 n = α 3 n β b j = 0 n 1 j + b j = 0 n 1 ( j + 1 ) + b j = 0 n 1 ( j + 1 ) z 0 w 2 b n w 1 b n w 0 b n = α 3 n β b n ( 3 n + 1 ) 2 z 0 w 2 b n w 1 b n w 0 b n ,
z 3 n + 1 = α 3 n + 1 β b j = 0 n 1 ( j + 1 ) + b j = 0 n 1 ( j + 1 ) + b j = 0 n 1 ( j + 1 ) z 0 w 2 b n w 1 b n w 0 b ( n + 1 ) = α 3 n + 1 β 3 b n ( n + 1 ) 2 z 0 w 2 b n w 1 b n w 0 b ( n + 1 ) ,
z 3 n + 2 = α 3 n + 2 β b j = 0 n 1 ( j + 1 ) + b j = 0 n 1 ( j + 1 ) + b j = 0 n ( j + 1 ) z 0 w 2 b ( n + 1 ) w 1 b n w 0 b ( n + 1 ) = α 3 n + 2 β b ( n + 1 ) ( 3 n + 2 ) 2 z 0 w 2 b ( n + 1 ) w 1 b n w 0 b ( n + 1 ) ,
for n N 0 .
From all above presented formulas we see that system of difference equations (2) is solvable in this case, as claimed. ☐
The statements in the following corollary are direct consequences of the formulas presented in the proof of Theorem 2.
Corollary 2.
Assume that a , b , c Z , d = 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then the following statements are true.
(a) 
If a 1 c a 3 , then the general solution to system (2) is given by (36), (47)–(49).
(b) 
If a 1 c = a 3 , then the general solution to system (2) is given by (36), (50)–(52).
(c) 
If a 1 = c , then the general solution to system (2) is given by (37), (53)–(55).
(d) 
If a = 1 c , then the general solution to system (2) is given by (36), (56)–(58).
(e) 
If a = c = 1 , then the general solution to system (2) is given by (37), (59)–(61).
The following result concerns the case c = 0 . Note that under the condition system (2) of difference equations is
z n + 1 = α z n a w n b , w n + 1 = β z n d , n N 0 ,
which is system (3.15) in [33].
From (62) it follows that
z n + 1 = α β b z n a z n 1 b d ,
for n N .
How this product-type difference equation can be solved was explained in the proof of Theorem 3.3 in [33]. From the closed-form formulas obtained in the proof of Theorem 3.3 therein it follows directly that the following result holds. Hence, to avoid repeating, the proof of the theorem is omitted here.
Theorem 3.
Assume that a , b , d Z , c = 0 , α , β , z 0 , w 0 C \ { 0 } . Then the following statements are true.
(a) 
If b d 0 , a 2 + 4 b d 0 and a + b d 1 , then the general solution to system (2) is given by
z n = α ( λ 2 1 ) λ 1 n + 1 ( λ 1 1 ) λ 2 n + 1 + λ 1 λ 2 ( λ 1 λ 2 ) ( λ 1 1 ) ( λ 2 1 ) β b ( λ 2 1 ) λ 1 n ( λ 1 1 ) λ 2 n + λ 1 λ 2 ( λ 1 λ 2 ) ( λ 1 1 ) ( λ 2 1 ) z 0 λ 1 n + 1 λ 2 n + 1 λ 1 λ 2 w 0 b λ 1 n λ 2 n λ 1 λ 2 , w n = α d ( λ 2 1 ) λ 1 n ( λ 1 1 ) λ 2 n + λ 1 λ 2 ( λ 1 λ 2 ) ( λ 1 1 ) ( λ 2 1 ) β 1 + b d ( λ 2 1 ) λ 1 n 1 ( λ 1 1 ) λ 2 n 1 + λ 1 λ 2 ( λ 1 λ 2 ) ( λ 1 1 ) ( λ 2 1 ) z 0 d λ 1 n λ 2 n λ 1 λ 2 w 0 b d λ 1 n 1 λ 2 n 1 λ 1 λ 2 ,
for n N , where
λ 1 , 2 = a ± a 2 + 4 b d 2 .
(b) 
If b d 0 , a 2 + 4 b d = 0 and a + b d 1 , then the general solution to system (2) is given by
z n = α 1 ( n + 1 ) λ 1 n + n λ 1 n + 1 ( 1 λ 1 ) 2 β b 1 n λ 1 n 1 + ( n 1 ) λ 1 n ( 1 λ 1 ) 2 z 0 ( n + 1 ) λ 1 n w 0 b n λ 1 n 1 w n = α d 1 n λ 1 n 1 + ( n 1 ) λ 1 n ( 1 λ 1 ) 2 β 1 2 λ 1 + ( n 1 ) λ 1 n ( n 2 ) λ 1 n + 1 ( 1 λ 1 ) 2 z 0 d n λ 1 n 1 w 0 ( 1 n ) λ 1 n
for n N , where
λ 1 = a 2 .
(c) 
If b d 0 , a 2 + 4 b d 0 and a + b d = 1 , then the general solution to system (2) is given by
z n = α ( b d ) n + 1 + ( n + 1 ) b d + n ( 1 + b d ) 2 β b ( b d ) n + n b d + n 1 ( 1 + b d ) 2 z 0 1 ( b d ) n + 1 1 + b d w 0 b 1 ( b d ) n 1 + b d , w n = α d ( b d ) n + n b d + n 1 ( 1 + b d ) 2 β 1 + b d ( b d ) n 1 + ( n 1 ) b d + n 2 ( 1 + b d ) 2 z 0 d 1 ( b d ) n 1 + b d w 0 b d 1 ( b d ) n 1 1 + b d ,
for n N .
(d) 
If b d 0 , a 2 + 4 b d = 0 and a + b d = 1 , then the general solution to system (2) is given by
z n = α n ( n + 1 ) 2 β b ( n 1 ) n 2 z 0 n + 1 w 0 b n w n = α d ( n 1 ) n 2 β ( 3 n ) n 2 z 0 d n w 0 1 n
for n N .
(e) 
If b d = 0 and a 1 , then the general solution to system (2) is given by
z n = α 1 a n 1 a β b 1 a n 1 1 a z 0 a n w 0 b a n 1 w n = α d 1 a n 1 1 a β z 0 d a n 1 ,
for n N .
(f) 
If b d = 0 and a = 1 , then the general solution to system (2) is given by
z n = α n β b ( n 1 ) z 0 w 0 b w n = α d ( n 1 ) β z 0 d ,
for n N .
Theorem 4.
Assume that a , b , c , d Z , a c 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then system (2) is solvable in closed form.
Proof. 
To deal with the case we use and modify our method previously used, for example, in [40,43]. Since α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } , we have z n w n 0 for n N 0 . Hence,
w n b = z n + 1 α z n a , n N 0 ,
w n + 1 b = β b w n 2 b c z n b d , n N 0 ,
from which we get
z n + 2 = α 1 c β b z n + 1 a z n b d z n 1 c z n 2 a c ,
for n 2 .
Let δ = α 1 c β b ,
a 1 = a , b 1 = b d , c 1 = c , d 1 = a c , y 1 = 1 ,
then
z n + 2 = δ y 1 z n + 1 a 1 z n b 1 z n 1 c 1 z n 2 d 1 , n 2 .
From (67) is obtained
z n + 2 = δ y 1 ( δ z n a 1 z n 1 b 1 z n 2 c 1 z n 3 d 1 ) a 1 z n b 1 z n 1 c 1 z n 2 d 1 , = δ y 1 + a 1 z n a 1 a 1 + b 1 z n 1 b 1 a 1 + c 1 z n 2 c 1 a 1 + d 1 z n 3 d 1 a 1 = δ y 2 z n a 2 z n 1 b 2 z n 2 c 2 z n 3 d 2 ,
for n 3 , where
a 2 : = a 1 a 1 + b 1 , b 2 : = b 1 a 1 + c 1 , c 2 : = c 1 a 1 + d 1 , d 2 : = d 1 a 1 , y 2 : = y 1 + a 1 .
Suppose
z n + 2 = δ y k z n + 2 k a k z n + 1 k b k z n k c k z n k 1 d k ,
for a k 2 and n k + 1 , and
a k = a 1 a k 1 + b k 1 , b k = b 1 a k 1 + c k 1 , c k = c 1 a k 1 + d k 1 , d k = d 1 a k 1 ,
y k = y k 1 + a k 1 .
Using (67) in (70), we get
z n + 2 = δ y k ( δ z n + 1 k a 1 z n k b 1 z n k 1 c 1 z n k 2 d 1 ) a k z n + 1 k b k z n k c k z n k 1 d k = δ y k + a k z n + 1 k a 1 a k + b k z n k b 1 a k + c k z n k 1 c 1 a k + d k z n k 2 d 1 a k = δ y k + 1 z n + 1 k a k + 1 z n k b k + 1 z n k 1 c k + 1 z n k 2 d k + 1 ,
for n k + 2 , where
a k + 1 : = a 1 a k + b k , b k + 1 : = b 1 a k + c k , c k + 1 : = c 1 a k + d k , d k + 1 : = d 1 a k ,
y k + 1 : = y k + a k .
Hence, by the induction we get that (70)–(72) hold.
From (70) with k = n 1 and (8), it follows that
z n + 2 = δ y n 1 z 3 a n 1 z 2 b n 1 z 1 c n 1 z 0 d n 1 = ( α 1 c β b ) y n 1 ( α 1 + a + a 2 + b d β b + a b z 0 a 3 + 2 a b d w 2 a b c w 1 b c w 0 b ( a 2 + b d ) ) a n 1 × ( α 1 + a β b z 0 a 2 + b d w 2 b c w 0 a b ) b n 1 ( α z 0 a w 0 b ) c n 1 z 0 d n 1 = α ( 1 c ) y n 1 + ( 1 + a + a 2 + b d ) a n 1 + ( 1 + a ) b n 1 + c n 1 β b y n 1 + b ( 1 + a ) a n 1 + b b n 1 × z 0 ( a 3 + 2 a b d ) a n 1 + ( a 2 + b d ) b n 1 + a c n 1 + d n 1 w 2 a b c a n 1 + b c b n 1 w 1 b c a n 1 × w 0 b ( a 2 + b d ) a n 1 + a b b n 1 + b c n 1 = α y n + 2 c y n 1 β b y n + 1 z 0 a n + 2 c a n 1 w 2 b c a n w 1 b c a n 1 w 0 b a n + 1 ,
for n N 0 .
From (71) we have
a k = a 1 a k 1 + b 1 a k 2 + c 1 a k 3 + d 1 a k 4 , k 5 ,
and that b k , c k and d k also satisfy the equation, and using (74) and (75) for k = 0 , 1 , 2 , 3 is obtained
a 3 = 0 , a 2 = 0 , a 1 = 0 , a 0 = 1 ,
y 3 = y 2 = y 1 = y 0 = 0 , y 1 = 1 ,
and
y k = j = 0 k 1 a j ,
(see, for example, [33] for more details).
The solvability of (77) is a classical thing. Hence, by finding closed form formula for a k , employing it in (80), then using Lemma 2, is calculated y k . These two formulas and (76) give a closed-form formula for solution to (65).
Now note that
z n d = w n + 1 β w n 2 c , n N 0 ,
z n + 1 d = α d z n a d w n b d , n N 0 ,
from which we get
w n + 2 = α d β 1 a w n + 1 a w n b d w n 1 c w n 2 a c , n N 0 .
As above is get
w n + 2 = η y k w n + 2 k a k w n + 1 k b k w n k c k w n k 1 d k ,
for every k , n N , such that n k 1 , where η = α d β 1 a , ( a k ) k N , ( b k ) k N , ( c k ) k N and ( d k ) k N are defined by (66) and (71), while ( y k ) k N is defined by (72) and (79).
From (84) with k = n + 1 and (8), we have
w n + 2 = η y n + 1 w 1 a n + 1 w 0 b n + 1 w 1 c n + 1 w 2 d n + 1 = ( α d β 1 a ) y n + 1 ( β w 2 c z 0 d ) a n + 1 w 0 b n + 1 w 1 c n + 1 w 2 d n + 1 = α d y n + 1 β ( 1 a ) y n + 1 + a n + 1 z 0 d a n + 1 w 2 c a n + 1 + d n + 1 w 1 c n + 1 w 0 b n + 1 = α d y n + 1 β y n + 2 a y n + 1 z 0 d a n + 1 w 2 c ( a n + 1 a a n ) w 1 c ( a n a a n 1 ) w 0 a n + 2 a a n + 1 ,
for n N 0 .
Recall that closed form formulas for a k and y k can be found. Applying them into (85) we show the solvability of equation (83). It can be checked that (76) and (85) present a solution to (2), from which the theorem follows. ☐
Corollary 3.
Consider system (2) with a , b , c , d Z , a c 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then the general solution to (2) is given by (76) and (85), where ( a k ) k N is defined by (77) and (78), while ( y k ) k N is defined by (79) and (80)

Detailed Form of Solutions Given in (76) and (85)

Equation (77) is not only theoretically but also practically solvable. The reason for this is that its characteristic polynomial
p 4 ( λ ) = λ 4 a λ 3 b d λ 2 c λ + a c ,
for the case a c 0 , is of the forth degree, thus, solvable by radicals.
Note that the equation p 4 ( λ ) = 0 can be written as follows [44].
λ 2 a 2 λ + s 2 2 a 2 4 + s + b d λ 2 a s 2 c λ + s 2 4 a c = 0 ,
Now, choose s so that the expression in the second bracket in (87) is a perfect square. Thus, it must be
( a s 2 c ) 2 = ( a 2 + 4 s + 4 b d ) ( s 2 4 a c ) ,
that is,
s 3 + b d s 2 3 a c s a 3 c c 2 4 a b c d = 0 .
Hence, (87) can be written as
λ 2 a 2 λ + s 2 2 a 2 + 4 s + 4 b d 2 λ a s 2 c 2 a 2 + 4 s + 4 b d 2 = 0 ,
which is equivalent to
λ 2 a 2 + a 2 + 4 s + 4 b d 2 λ + s 2 + a s 2 c 2 a 2 + 4 s + 4 b d = 0 ,
λ 2 a 2 a 2 + 4 s + 4 b d 2 λ + s 2 a s 2 c 2 a 2 + 4 s + 4 b d = 0 .
Using the change of variables s = t b d 3 in (88), it follows that
t 3 b 2 d 2 3 + 3 a c t + 2 b 3 d 3 27 a 3 c c 2 3 a b c d = 0 .
Let
p = b 2 d 2 3 + 3 a c and q = 27 a 3 c + 27 c 2 + 81 a b c d 2 b 3 d 3 27 .
As usual, a solution to (92) is found in the form t = u + v . Putting it into (92) and requesting u v = p / 3 , is get u 3 + v 3 = q and u 3 v 3 = p 3 / 27 . Hence, u 3 and v 3 are solutions to the equation z 2 + q z p 3 / 27 = 0 , so, they must be ( q ± q 2 + 4 p 3 / 27 ) / 2 .
Hence
t = q 2 q 2 4 + p 3 27 3 + q 2 + q 2 4 + p 3 27 3 ,
where any of the three possible values of the right-hand side can be chosen. If p = Δ 0 / 3 and q = Δ 1 / 27 , then it can be written as
t = 1 3 2 3 Δ 1 Δ 1 2 4 Δ 0 3 3 + Δ 1 + Δ 1 2 4 Δ 0 3 3 .
For such chosen s, that is, t, (90) and (91) can be solved and by some calculation it is obtained that the zeros of p 4 are
λ 1 = a 4 + 1 2 a 2 4 + 2 b d 3 + t + 1 2 a 2 2 + 4 b d 3 t Q 4 a 2 4 + 2 b d 3 + t ,
λ 2 = a 4 + 1 2 a 2 4 + 2 b d 3 + t 1 2 a 2 2 + 4 b d 3 t Q 4 a 2 4 + 2 b d 3 + t ,
λ 3 = a 4 1 2 a 2 4 + 2 b d 3 + t + 1 2 a 2 2 + 4 b d 3 t + Q 4 a 2 4 + 2 b d 3 + t ,
λ 4 = a 4 1 2 a 2 4 + 2 b d 3 + t 1 2 a 2 2 + 4 b d 3 t + Q 4 a 2 4 + 2 b d 3 + t ,
where
Q : = a 3 8 c 4 a b d .
Recall, that the nature of these λ j ’s depends on the sign of the discriminant
Δ : = 1 27 ( 4 Δ 0 3 Δ 1 2 ) ,
where
Δ 0 : = b 2 d 2 + 9 a c ,
Δ 1 : = 2 b 3 d 3 + 81 a c b d + 27 a 3 c + 27 c 2 ,
and the signs of
P : = 8 b d 3 a 2
and
D : = 48 a c 16 b 2 d 2 16 a 2 b d 3 a 4 .
Zeros of p 4 are different and none of them is 1. If a, b, c and d are chosen such that
Δ 0 = b 2 d 2 + 9 a c < 0 ,
then it will be Δ < 0 , from which by Lemma 3 we have in this case that p 4 has four different zeros. Moreover, since Δ < 0 two zeros are real and two are complex-conjugate.
Zeros of p 4 are different and one of them is 1. Polynomial p 4 has a zero equal to 1 if p 4 ( 1 ) = 1 a b d c + a c = 0 , that is, if
( a 1 ) ( c 1 ) = b d ,
so that
p 4 ( λ ) = λ 4 a λ 3 ( a 1 ) ( c 1 ) λ 2 c λ + a c .
Thus, if we choose a and c such that
p 4 ( 1 ) = 4 3 a 2 ( a 1 ) ( c 1 ) c 0 ,
that is, ( 2 a 1 ) ( 2 c + 1 ) 3 , then p 4 will be such a polynomial if Δ 0 . For example, if a = 3 and c = 2 , then b d = 2 0 , Δ 0 , which means that polynomial p 4 has all zeros mutually different and exactly one of them is equal to 1
p 4 ( λ ) = λ 4 3 λ 3 2 λ 2 2 λ + 6 = ( λ 1 ) ( λ 3 2 λ 2 4 λ 6 ) .
Since in these two cases λ j λ i , i j , then the general solution to (77) is
a n = γ 1 λ 1 n + γ 2 λ 2 n + γ 3 λ 3 n + γ 4 λ 4 n , n N ,
where γ i , i = 1 , 4 ¯ , are arbitrary constants.
Lemma 1 implies
j = 1 4 λ j l p 4 ( λ j ) = 0 for l = 0 , 2 ¯ , and j = 1 4 λ j 3 p 4 ( λ j ) = 1 .
From initial conditions (78) and (109), it is obtained
a n = j = 1 4 λ j n + 3 p 4 ( λ j ) = λ 1 n + 3 ( λ 1 λ 2 ) ( λ 1 λ 3 ) ( λ 1 λ 4 ) + λ 2 n + 3 ( λ 2 λ 1 ) ( λ 2 λ 3 ) ( λ 2 λ 4 ) + λ 3 n + 3 ( λ 3 λ 1 ) ( λ 3 λ 2 ) ( λ 3 λ 4 ) + λ 4 n + 3 ( λ 4 λ 1 ) ( λ 4 λ 2 ) ( λ 4 λ 3 ) ,
for n 3 ([40]).
Combining (80) and (110), we get
y n = j = 0 n 1 i = 1 4 λ i j + 3 p 4 ( λ i ) = i = 1 4 λ i 3 ( λ i n 1 ) p 4 ( λ i ) ( λ i 1 ) , n N ,
when λ i 1 , i = 1 , 4 ¯ , and
y n = n 2 a + c 2 a c + i = 2 4 λ i 3 ( λ i n 1 ) p 4 ( λ i ) ( λ i 1 ) , n N ,
when one of the zeros is 1, say λ 1 , is equal to 1.
Note that if one of the zeros is equal to 1, then we have
p 4 ( λ ) = λ 4 a λ 3 ( a 1 ) ( c 1 ) λ 2 c λ + a c = ( λ 1 ) ( λ 3 + ( 1 a ) λ 2 + ( 1 a ) c λ a c ) .
By using the change of variables λ = t + a 1 3 in (113) we get that the equation
λ 3 + ( 1 a ) λ 2 + ( 1 a ) c λ a c = 0 ,
is transformed to the following one
t 3 + p ˜ t + q ˜ = 0 ,
with
p ˜ = ( a 1 ) ( a + 3 c 1 ) 3 and q ˜ = 2 ( a 1 ) 3 27 + c ( a 1 ) 2 3 + a c ,
whose solutions are given by
t j = ε j 1 q ˜ 2 q ˜ 2 4 + p ˜ 3 27 3 + ε j 1 q ˜ 2 + q ˜ 2 4 + p ˜ 3 27 3 , j = 1 , 3 ¯ ,
where ε 3 = 1 and ε 1 .
Hence,
λ j = a 1 3 + ε j 2 q ˜ 2 q 2 4 + p ˜ 3 27 3 + ε j 2 q ˜ 2 + q ˜ 2 4 + p ˜ 3 27 3 ,
for j = 2 , 4 ¯ , are the other three zeros of the equation p 4 ( λ ) = 0 , in this case.
A simple calculation along with (109) shows that (111) holds also for n = j , j = 0 , 3 ¯ .
The previous analysis along with Corollary 3 implies the following corollary.
Corollary 4.
Assume that a , b , c , d Z , a c 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } and Δ 0 . Then the following statements are true.
(a) 
If ( a 1 ) ( c 1 ) b d , then the general solution to (2) is given by (76) and (85), where ( a n ) n 3 is given by (110), ( y n ) n 3 is given by (111), while λ j ’s, j = 1 , 4 ¯ , are given by (95)–(98).
(b) 
If ( a 1 ) ( c 1 ) = b d and ( 2 a 1 ) ( 2 c + 1 ) 3 , then the general solution to (2) is given by (76) and (85), where ( a n ) n 3 is given by (110) with λ 1 = 1 , ( y n ) n 3 is given by (112), λ 1 = 1 , while λ j ’s, j = 2 , 4 ¯ , are given by (116) and (115).
p 4 has only one double zero which is equal to 1. Polynomial (86) has a double zero equal to 1 if (105) holds and
( 2 a 1 ) ( 2 c + 1 ) = 3 .
From (117) we have that it must be a = 2 and c = 0 , or a = 1 and c = 1 , or a = 0 and c = 2 , or a = 1 and c = 1 . If a = 0 or c = 0 , then a c = 0 , which contradicts to the assumption a c 0 .
If a = c = 1 , then b d = 0 , from which it follows that
p 4 ( λ ) = λ 4 λ 3 λ + 1 = ( λ 1 ) 2 ( λ 2 + λ + 1 ) ,
and consequently
λ 3 , 4 = 1 ± i 3 2 ,
are the non-unit zeros of polynomial p 4 in the case.
If a = c = 1 , then b d = 4 , from which it follows that
p 4 ( λ ) = λ 4 + λ 3 4 λ 2 + λ + 1 = ( λ 1 ) 2 ( λ 2 + 3 λ + 1 ) ,
and consequently
λ 3 , 4 = 3 ± 5 2 ,
are the non-unit zeros of polynomial p 4 in the case.
From this, we have proved in passing, that there are no a , c Z \ { 0 } , such that 1 is a triple zero of p 4 or that it has two pairs of double zeros one of which is 1.
In these two cases we have (see, for example, [38])
a n = n ( 1 λ 3 ) ( 1 λ 4 ) + 3 λ 3 λ 4 2 λ 3 2 λ 4 + 1 ( 1 λ 3 ) 2 ( 1 λ 4 ) 2 + λ 3 n + 3 ( λ 3 1 ) 2 ( λ 3 λ 4 ) + λ 4 n + 3 ( λ 4 1 ) 2 ( λ 4 λ 3 ) ,
and
y n = j = 0 n 1 j ( 1 λ 3 ) ( 1 λ 4 ) + 3 λ 3 λ 4 2 λ 3 2 λ 4 + 1 ( 1 λ 3 ) 2 ( 1 λ 4 ) 2 + λ 3 j + 3 ( λ 3 1 ) 2 ( λ 3 λ 4 ) + λ 4 j + 3 ( λ 4 1 ) 2 ( λ 4 λ 3 ) = ( n 1 ) n 2 ( 1 λ 3 ) ( 1 λ 4 ) + n ( 3 λ 3 λ 4 2 λ 3 2 λ 4 + 1 ) ( 1 λ 3 ) 2 ( 1 λ 4 ) 2 + λ 3 3 ( λ 3 n 1 ) ( λ 3 1 ) 3 ( λ 3 λ 4 ) + λ 4 3 ( λ 4 n 1 ) ( λ 4 1 ) 3 ( λ 4 λ 3 ) .
p 4 has a double zero different from 1. Let b d = 0 , then
p 4 ( λ ) = λ 4 a λ 3 c λ + a c = ( λ a ) ( λ 3 c ) .
If we take c = a 3 , then it is obtained
p 4 ( λ ) = ( λ a ) 2 ( λ 2 + a λ + a 2 ) ,
which for a Z \ { 0 , 1 } is a polynomial with a double zero different from 1 and two non-real complex-conjugate zeros.
Since, in the case λ 1 = λ 2 , λ i λ j , 2 i , j 4 , we have
a n = ( γ 1 + γ 2 n ) λ 2 n + γ 3 λ 3 n + γ 4 λ 4 n , n N ,
where γ i , i = 1 , 4 ¯ , are arbitrary constants, and the solution satisfying (78) can be obtained, for example, by letting λ 1 λ 2 in (110) [38]
a n = λ 2 n + 2 ( ( n + 3 ) ( λ 2 λ 3 ) ( λ 2 λ 4 ) λ 2 ( 2 λ 2 λ 3 λ 4 ) ) ( λ 2 λ 3 ) 2 ( λ 2 λ 4 ) 2 + λ 3 n + 3 ( λ 3 λ 2 ) 2 ( λ 3 λ 4 ) + λ 4 n + 3 ( λ 4 λ 2 ) 2 ( λ 4 λ 3 ) .
From (80), (123) and by Lemma 2, we get
y n = j = 0 n 1 ( λ 2 j + 2 ( ( j + 3 ) ( λ 2 λ 3 ) ( λ 2 λ 4 ) λ 2 ( 2 λ 2 λ 3 λ 4 ) ) ( λ 2 λ 3 ) 2 ( λ 2 λ 4 ) 2 + λ 3 j + 3 ( λ 3 λ 2 ) 2 ( λ 3 λ 4 ) + λ 4 j + 3 ( λ 4 λ 2 ) 2 ( λ 4 λ 3 ) ) = λ 2 3 n λ 2 n + 2 + ( n 1 ) λ 2 n + 3 ( λ 2 λ 3 ) ( λ 2 λ 4 ) ( 1 λ 2 ) 2 + ( λ 2 4 2 λ 2 3 λ 3 2 λ 2 3 λ 4 + 3 λ 2 2 λ 3 λ 4 ) ( λ 2 n 1 ) ( λ 2 λ 3 ) 2 ( λ 2 λ 4 ) 2 ( λ 2 1 ) + λ 3 3 ( λ 3 n 1 ) ( λ 3 λ 2 ) 2 ( λ 3 λ 4 ) ( λ 3 1 ) + λ 4 3 ( λ 4 n 1 ) ( λ 4 λ 2 ) 2 ( λ 4 λ 3 ) ( λ 4 1 ) .
From the previous analysis and Corollary 3 we obtain the following result.
Corollary 5.
Assume that a , b , c , d Z , a c 0 and α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then the following statements are true.
(a) 
If only one of the zeros of p 4 is double and different from 1, then the general solution to (2) is given by (76) and (85), where ( a n ) n 3 is given by (123), while ( y n ) n 3 is given by (124).
(b) 
If 1 is a unique double zero of polynomial p 4 , say λ 1 = λ 2 = 1 , then the general solution to (2) is given by (76) and (85), where ( a n ) n 3 is given by (120), ( y n ) n 3 is given by (121), while λ 3 , 4 are given by (118) if a = c = 1 or by (119) if a = c = 1 .
Two pairs of different double zeros. In this case it must be D = 0 which implies that
16 t 2 + 16 a 2 t + 3 a 4 48 a c = 0 ,
where t = b d . On the other hand, it must be Δ = 0 , which is equivalent to 4 Δ 0 3 = Δ 1 2 , that is,
( 2 b 3 d 3 + 81 a c b d + 27 a 3 c + 27 c 2 ) 2 = 4 ( b 2 d 2 + 9 a c ) 3 ,
from which it follows that
c 16 a t 4 + 4 ( a 3 + c ) t 3 207 a 2 c t 2 162 a c ( a 3 + c ) t 27 c ( a 3 c ) 2 ) = 0 .
The problem of the existence of a joint zero of the polynomials in (125) and (126) for some integers a , b, c and d, such that a c 0 , seems quite technical, so we leave it to the reader.
Solutions to (77) in this case are
a n = ( γ 1 + γ 2 n ) λ 1 n + ( γ 3 + γ 4 n ) λ 3 n , n N ,
where γ i , i = 1 , 4 ¯ are constants. The solution to (127) such that (78) holds is [38]
a n = λ 2 n + 2 ( n ( λ 2 λ 4 ) 2 + λ 2 2 4 λ 2 λ 4 + 3 λ 4 2 ) ( λ 2 λ 4 ) 4 + λ 4 n + 2 ( n ( λ 4 λ 2 ) 2 + λ 4 2 4 λ 2 λ 4 + 3 λ 2 2 ) ( λ 4 λ 2 ) 4 .
From (80), (128) and Lemma 2, we get
y n = j = 0 n 1 λ 2 j + 2 ( j ( λ 2 λ 4 ) 2 + λ 2 2 4 λ 2 λ 4 + 3 λ 4 2 ) ( λ 2 λ 4 ) 4 + λ 4 j + 2 ( j ( λ 4 λ 2 ) 2 + λ 4 2 4 λ 2 λ 4 + 3 λ 2 2 ) ( λ 4 λ 2 ) 4 = λ 2 3 n λ 2 n + 2 + ( n 1 ) λ 2 n + 3 ( λ 2 λ 4 ) 2 ( 1 λ 2 ) 2 + ( λ 2 4 4 λ 2 3 λ 4 + 3 λ 2 2 λ 4 2 ) ( λ 2 n 1 ) ( λ 2 λ 4 ) 4 ( λ 2 1 ) + λ 4 3 n λ 4 n + 2 + ( n 1 ) λ 4 n + 3 ( λ 4 λ 2 ) 2 ( 1 λ 4 ) 2 + ( λ 4 4 4 λ 2 λ 4 3 + 3 λ 2 2 λ 4 2 ) ( λ 4 n 1 ) ( λ 4 λ 2 ) 4 ( λ 4 1 ) .
Corollary 6.
Assume that a , b , c , d Z , a c 0 and α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then the following statements are true.
(a) 
If polynomial p 4 has two pairs of double zeros both different from 1, then the general solution to (2) is given by (76) and (85), where ( a n ) n 3 is given by (128), while ( y n ) n 3 is given by (129).
(b) 
The polynomial in (86) cannot have two pairs of double zeros such that one of them is equal to 1.
Triple zero case. In this case it must be Δ = Δ 0 = 0 , or equivalently, Δ 0 = Δ 1 = 0 . Hence,
a c = b 2 d 2 / 9 and 2 b 3 d 3 81 a c b d 27 a 3 c 27 c 2 = 0 ,
from which it follows that
11 ( b d ) 3 + 3 a 2 ( b d ) 2 ( b d ) 4 3 a 2 = 0 .
Since b d = 0 implies a c = 0 , which contradicts to the assumption a c 0 , from (130) it follows that
s 2 33 s 9 = 0 ,
where s = b d / a 2 . Hence, it must be
b d / a 2 = 33 + 1125 2 or b d / a 2 = 33 1125 2 ,
which is not possible since b d / a 2 is a rational number, whereas ( 33 ± 1125 ) / 2 are both irrational numbers. Hence, p 4 cannot have a triple, and consequently cannot have a quadruple zero.
Corollary 7.
Assume that a , b , c , d Z , a c 0 and α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then polynomial (86) cannot have a triple zero.
Theorem 5.
Assume that b , c , d Z , a = 0 , b c d 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } . Then system (2) is solvable in closed form.
Proof. 
In this case (2) is
z n + 1 = α w n b , w n + 1 = β w n 2 c z n d , n N 0 .
Using the first equation in (132) into the second one it follows that
w n + 1 = α d β w n 1 b d w n 2 c , n N .
Let η = α d β ,
β 1 = b d , γ 1 = c , δ 1 = 0 , x 1 = 1 ,
then, we have
w n + 1 = η w n 1 β 1 w n 2 γ 1 w n 3 δ 1 , n N .
Hence
w n + 1 = η x 1 ( η w n 3 β 1 w n 4 γ 1 ) β 1 w n 2 γ 1 w n 3 δ 1 = η x 1 + β 1 w n 2 γ 1 w n 3 β 1 β 1 + δ 1 w n 4 γ 1 β 1 = η x 2 w n 2 β 2 w n 3 γ 2 w n 4 δ 2 ,
for n 3 , where
x 2 : = x 1 + β 1 , β 2 : = γ 1 , γ 2 : = β 1 β 1 + δ 1 , δ 2 : = γ 1 β 1 .
Suppose that
w n + 1 = η x k w n k β k w n k 1 γ k w n k 2 δ k ,
for a k N \ { 1 } and every n k + 1 , and
β k = γ k 1 , γ k = β 1 β k 1 + δ k 1 , δ k = γ 1 β k 1 ,
x k = x k 1 + β k 1 .
Using (135) in (136), we have
w n + 1 = η x k ( η w n k 2 β 1 w n k 3 γ 1 ) β k w n k 1 γ k w n k 2 δ k = η x k + β k w n k 1 γ k w n k 2 β 1 β k + δ k w n k 3 γ 1 β k = η x k + 1 w n k 1 β k + 1 w n k 2 γ k + 1 w n k 3 δ k + 1 ,
for a k 2 and n k + 2 , and where
β k + 1 : = γ k , γ k + 1 : = β 1 β k + δ k , δ k + 1 : = γ 1 β k ,
x k + 1 : = x k + β k .
From this we see that hypotheses (136)–(138) are true.
Let k = n 1 , then from (136)–(138), it follows that
w n + 1 = η x n 1 w 1 β n 1 w 0 γ n 1 w 1 δ n 1 = ( α d β ) x n 1 ( β w 2 c z 0 d ) β n 1 w 0 γ n 1 w 1 δ n 1 = α d x n 1 β x n 1 + β n 1 z 0 d β n 1 w 2 c β n 1 w 1 δ n 1 w 0 γ n 1 = α d x n 1 β x n z 0 d β n 1 w 2 c β n 1 w 1 c β n 2 w 0 β n ,
for n 3 .
From (137) it follows that
β k = β 1 β k 2 + γ 1 β k 3 , k 4 ,
and that γ k and δ k also satisfy (140).
From (140) and since γ 1 = c 0 , we have that
β k 3 = β k β 1 β k 2 γ 1 ,
From (138) and (141) with k = 3 , 2 , 1 and some calculation is obtained
β 3 = 0 , β 2 = 0 , β 1 = 1 , β 0 = 0 ,
x 3 = x 2 = x 1 = 0 , x 0 = x 1 = 1 ,
and
x k = 1 + j = 1 k 1 β j .
Since equation (140) is solvable, we can calculate β k , from which along with (144) and Lemma 2, y k is calculated. These facts along with (139) gives a closed form formula for (133).
Using (139) in the first equation in (132) we get
z n + 1 = α 1 + b d x n 2 β b x n 1 z 0 b d β n 2 w 2 b c β n 2 w 1 b c β n 3 w 0 b β n 1 , n N 0 .
It is not difficult to see that (139) and (145) are solutions to (2) in the case. ☐
Theorem 5 solves theoretically system (2) when a = 0 and b c d 0 . Now we will practically solve it in terms of the parameters and initial values. The following polynomial
p 3 ( λ ) = λ 3 b d λ c ,
is the characteristic one associated to Equation (140), and its solutions are
λ j = 1 3 2 3 ε j Δ ^ 1 Δ ^ 1 2 4 Δ ^ 0 3 3 + ε ¯ j Δ ^ 1 + Δ ^ 1 2 4 Δ ^ 0 3 3 , j = 0 , 2 ¯ ,
where
Δ ^ 0 = 3 b d = : 3 p and Δ ^ 1 = 27 c = : 27 q ,
and ε 3 = 1 , ε 1 .
Zeros of p 3 are different and none of them is 1. In this case it must be Δ ^ 1 2 4 Δ ^ 0 3 , which can written as
27 c 2 4 ( b d ) 3 .
Hence, in the case all the zeros of p 3 are different. If additionally p 3 ( 1 ) 1 , that is, c + b d 1 , then none of the zeros is 1. For example, if c = b d = k N , then we have such a situation.
Zeros of p 3 are different and one of them is 1. Polynomial p 3 has a zero equal to 1 if c + b d = 1 . Then
p 3 ( λ ) = λ 3 + ( c 1 ) λ c = ( λ 1 ) ( λ 2 + λ + c ) .
Hence
λ 1 = 1 , λ 2 , 3 = 1 ± 1 4 c 2 ,
are the zeros of the equation p 3 ( λ ) = 0 , in this case. Since p 3 ( 1 ) = 2 + c , it follows that p 3 can have 1 as a double zero only if c = 2 .
The general solution to (140) in these two cases is
β n = α 1 λ 1 n + α 2 λ 2 n + α 3 λ 3 n , n N ,
for some constants α i , i = 1 , 3 ¯ , which due to γ 1 = c 0 can be prolonged for every non-positive index.
From Lemma 1 with Q ( t ) = p 3 ( t ) = j = 1 3 ( t λ j ) , we have
j = 1 3 λ j l p 3 ( λ j ) = 0 , for l = 0 , 1 , and j = 1 3 λ j 2 p 3 ( λ j ) = 1 .
From (142), (150) and (151), we get
β n = λ 1 n + 3 ( λ 1 λ 2 ) ( λ 1 λ 3 ) + λ 2 n + 3 ( λ 2 λ 1 ) ( λ 2 λ 3 ) + λ 3 n + 3 ( λ 3 λ 1 ) ( λ 3 λ 2 ) ,
for n 2 .
From (144) and (152), it follows that
x n = 1 + k = 1 n 1 j = 1 3 λ j k + 3 p 3 ( λ j ) ,
for n N .
If λ j 1 , j = 1 , 3 ¯ , then from formula (153), it follows that
x n = 1 + λ 1 4 ( λ 1 n 1 1 ) ( λ 1 λ 2 ) ( λ 1 λ 3 ) ( λ 1 1 ) + λ 2 4 ( λ 2 n 1 1 ) ( λ 2 λ 1 ) ( λ 2 λ 3 ) ( λ 2 1 ) + λ 3 4 ( λ 3 n 1 1 ) ( λ 3 λ 1 ) ( λ 3 λ 2 ) ( λ 3 1 ) ,
for n N , moreover, (154) holds for every n 2 .
If one of the zeros is 1, say λ 3 , then 1 λ 1 λ 2 1 , and we have
x n = 1 + λ 1 4 ( λ 1 n 1 1 ) ( λ 1 λ 2 ) ( λ 1 1 ) 2 + λ 2 4 ( λ 2 n 1 1 ) ( λ 2 λ 1 ) ( λ 2 1 ) 2 + n 1 ( λ 1 1 ) ( λ 2 1 ) .
for n N . Moreover, due to (151), (155) holds for every n 2 .
Corollary 8.
Assume that a , b , c , d Z , a = 0 , b c d 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } and Δ ^ 1 2 4 Δ ^ 0 3 . Then the following statements are true.
(a) 
If c + b d 1 , then the general solution to (2) is given by (139) and (145), where ( β n ) n 2 is given by (152), ( x n ) n 2 is given by (154), while λ j ’s, j = 1 , 3 ¯ are given by (147) and (148).
(b) 
If c + b d = 1 , then p 3 has a unique zero equal to 1, say λ 3 , and the general solution to (2) is given by formulas (139) and (145), where ( β n ) n 2 is given by (152) with λ 3 = 1 , ( x n ) n 2 is given by (155), while λ j ’s, j = 1 , 3 ¯ are given by (149).
One of the zeros is double. In this case it must be Δ ^ 1 2 = 4 Δ ^ 0 3 , that is, ( b d ) 3 = 27 c 2 / 4 . Assume that m is a double zero of p 3 , then it must be
m 3 b d m c = 0 and 3 m 2 b d = 0 ,
from which it follows that
p 3 ( λ ) = λ 3 3 m 2 λ + 2 m 3 = ( λ m ) 2 ( λ + 2 m ) .
Since b d 0 , we have m 0 . From this and since c 0 , from (156) we see that p 3 cannot have a triple zero. It also cannot have a unique zero equal to 1, since otherwise we would have 2 m = 1 , from which it would be b d = 3 / 4 Z , which would be a contradiction. Note also that the polynomial can have 1 as a double zero when b d = 3 and c = 2 .
If λ 1 λ 2 = λ 3 , then the general solution to (140) has the following form
β n = α ^ 1 λ 1 n + ( α ^ 2 + α ^ 3 n ) λ 2 n , n N ,
where α ^ i , i = 1 , 3 ¯ are constants. Since, in our case condition (142) must be satisfied, the solution ( β n ) n 2 to (140) can be found by letting λ 3 λ 2 in (152), so that
β n = λ 1 n + 3 + 2 λ 2 3 λ 1 + n ( λ 2 λ 1 ) λ 2 n + 2 ( λ 2 λ 1 ) 2 ,
for n 2 .
From (144) and (158), we have
x n = 1 + j = 1 n 1 β j = 1 + j = 1 n 1 λ 1 j + 3 + 2 λ 2 3 λ 1 + j ( λ 2 λ 1 ) λ 2 j + 2 ( λ 2 λ 1 ) 2 ,
for every n N .
From (159) and Lemma 2, we get
x n = 1 + λ 1 4 ( λ 1 n 1 1 ) ( λ 2 λ 1 ) 2 ( λ 1 1 ) + ( 2 λ 2 3 λ 1 ) λ 2 3 ( λ 2 n 1 1 ) ( λ 2 λ 1 ) 2 ( λ 2 1 ) + λ 2 3 ( 1 n λ 2 n 1 + ( n 1 ) λ 2 n ) ( λ 2 λ 1 ) ( λ 2 1 ) 2 ,
for n N (in fact, (160) hold also for every n 2 ).
If we assume that λ 1 1 and λ 2 = λ 3 = 1 , then from (159) it follows that
x n = 1 + λ 1 4 ( λ 1 n 1 1 ) ( λ 1 1 ) 3 + ( 2 3 λ 1 ) ( n 1 ) ( λ 1 1 ) 2 + ( n 1 ) n 2 ( 1 λ 1 ) .
A direct calculation shows that (161) holds also for every n 2 .
Corollary 9.
Assume that b , c , d Z , a = 0 , b c d 0 , α , β , z 0 , w 2 , w 1 , w 0 C \ { 0 } and Δ ^ 1 2 = 4 Δ ^ 0 3 . Then the following statements are true.
(a) 
If c + b d 1 , then the general solution to (2) is given by (139) and (145), where ( β n ) n 2 is given by (158), while ( x n ) n 2 is given by (160).
(b) 
If b d = 3 and c = 2 , then 1 is a double zero of p 3 , say, λ 2 = λ 3 = 1 , then the general solution to system (2) is given by (139) and (145), where ( β n ) n 2 is given by (158) with λ 2 = 1 , ( x n ) n 2 is given by (161), while λ 3 = 2 .
(c) 
It is not possible that 1 is a simple zero of p 3 .
Triple zero case. Since in this case p 3 must have the form in (156), we see that the only possibility that this polynomial has a triple zero is if m = 0 , which is impossible due to the condition c 0 .

Author Contributions

The author has contributed solely to the writing of this paper. He read and approved the manuscript.

Conflicts of Interest

The author declares that he has no conflicts of interest.

References

  1. Andruch-Sobilo, A.; Migda, M. Further properties of the rational recursive sequence xn+1 = axn−1/(b+cxnxn−1). Opusc. Math. 2006, 26, 387–394. [Google Scholar]
  2. Andruch-Sobilo, A.; Migda, M. On the rational recursive sequence xn+1 = axn−1/(b+cxnxn−1). Tatra Mt. Math. Publ. 2009, 43, 1–9. [Google Scholar]
  3. Berezansky, L.; Braverman, E. On impulsive Beverton-Holt difference equations and their applications. J. Differ. Equ. Appl. 2004, 10, 851–868. [Google Scholar] [CrossRef]
  4. Berezansky, L.; Migda, M.; Schmeidel, E. Some stability conditions for scalar Volterra difference equations. Opusc. Math. 2016, 36, 459–470. [Google Scholar] [CrossRef]
  5. Berg, L.; Stević, S. On some systems of difference equations. Appl. Math. Comput. 2011, 218, 1713–1718. [Google Scholar] [CrossRef]
  6. Brand, L. Differential and Difference Equations; John Wiley & Sons, Inc.: New York, NY, USA, 1966. [Google Scholar]
  7. Fotiades, N.; Papaschinopoulos, G. Asymptotic behavior of the positive solutions of a system of k difference equations of exponential form. Dyn. Contin. Discret. Impuls. Syst. Ser. A 2012, 19, 585–597. [Google Scholar]
  8. Fotiades, N.; Papaschinopoulos, G. On a system of difference equations with maximum. Appl. Math. Comput. 2013, 221, 684–690. [Google Scholar] [CrossRef]
  9. Iričanin, B. Global stability of some classes of higher-order nonlinear difference equations. Appl. Math. Comput. 2010, 216, 1325–1328. [Google Scholar] [CrossRef]
  10. Iričanin, B. On a higher-order nonlinear difference equation. Abstr. Appl. Anal. 2010, 2010, 8. [Google Scholar] [CrossRef]
  11. Iričanin, B.; Stević, S. Some systems of nonlinear difference equations of higher order with periodic solutions. Dyn. Contin. Discret. Impuls. Syst. Ser. A Math. Anal. 2006, 13, 499–508. [Google Scholar]
  12. Jordan, C. Calculus of Finite Differences; Chelsea Publishing Company: New York, NY, USA, 1956. [Google Scholar]
  13. Krechmar, V.A. A Problem Book in Algebra; Mir Publishers: Moscow, Russia, 1974. [Google Scholar]
  14. Levy, H.; Lessman, F. Finite Difference Equations; Dover Publications: New York, NY, USA, 1992. [Google Scholar]
  15. Mitrinović, D.S.; Adamović, D.D. Sequences and Series; Naučna Knjiga: Beograd, Serbia, 1980. [Google Scholar]
  16. Mitrinović, D.S.; Kečkić, J.D. Methods for Calculating Finite Sums; Naučna Knjiga: Beograd, Serbia, 1984. [Google Scholar]
  17. Papaschinopoulos, G.; Fotiades, N.; Schinas, C.J. On a system of difference equations including negative exponential terms. J. Differ. Equ. Appl. 2014, 20, 717–732. [Google Scholar] [CrossRef]
  18. Papaschinopoulos, G.; Psarros, N.; Papadopoulos, K.B. On a system of m difference equations having exponential terms. Electron. J. Qual. Theory Differ. Equ. 2015, 2015, 13. [Google Scholar] [CrossRef]
  19. Papaschinopoulos, G.; Schinas, C.J. On a system of two nonlinear difference equations. J. Math. Anal. Appl. 1998, 219, 415–426. [Google Scholar] [CrossRef]
  20. Papaschinopoulos, G.; Schinas, C.J. On the behavior of the solutions of a system of two nonlinear difference equations. Commun. Appl. Nonlinear Anal. 1998, 5, 47–59. [Google Scholar]
  21. Papaschinopoulos, G.; Schinas, C.J. Invariants for systems of two nonlinear difference equations. Differ. Equ. Dyn. Syst. 1999, 7, 181–196. [Google Scholar] [CrossRef]
  22. Papaschinopoulos, G.; Schinas, C.J. Invariants and oscillation for systems of two nonlinear difference equations. Nonlinear Anal. TMA 2001, 46, 967–978. [Google Scholar] [CrossRef]
  23. Papaschinopoulos, G.; Schinas, C.J. Oscillation and asymptotic stability of two systems of difference equations of rational form. J. Differ. Equ. Appl. 2001, 7, 601–617. [Google Scholar] [CrossRef]
  24. Papaschinopoulos, G.; Schinas, C.J. On the dynamics of two exponential type systems of difference equations. Comput. Math. Appl. 2012, 64, 2326–2334. [Google Scholar] [CrossRef]
  25. Papaschinopoulos, G.; Schinas, C.J.; Hatzifilippidis, V. Global behavior of the solutions of a max-equation and of a system of two max-equations. J. Comput. Anal. Appl. 2003, 5, 237–254. [Google Scholar]
  26. Papaschinopoulos, G.; Schinas, C.J.; Stefanidou, G. On the nonautonomous difference equation x n + 1 = A n + ( x n 1 p / x n q ) . Appl. Math. Comput. 2011, 217, 5573–5580. [Google Scholar]
  27. Papaschinopoulos, G.; Stefanidou, G. Asymptotic behavior of the solutions of a class of rational difference equations. Int. J. Differ. Equ. 2010, 5, 233–249. [Google Scholar]
  28. Stefanidou, G.; Papaschinopoulos, G.; Schinas, C.J. On a system of max difference equations. Dyn. Contin. Discret. Impuls. Syst. Ser. A 2007, 14, 885–903. [Google Scholar]
  29. Stefanidou, G.; Papaschinopoulos, G.; Schinas, C.J. On a system of two exponential type difference equations. Commun. Appl. Nonlinear Anal. 2010, 17, 1–13. [Google Scholar]
  30. Stević, S. Boundedness character of a class of difference equations. Nonlinear Anal. TMA 2009, 70, 839–848. [Google Scholar] [CrossRef]
  31. Stević, S. First-order product-type systems of difference equations solvable in closed form. Electron. J. Differ. Equ. 2015, 2015, 1–11. [Google Scholar] [CrossRef]
  32. Stević, S. New solvable class of product-type systems of difference equations on the complex domain and a new method for proving the solvability. Electron. J. Qual. Theory Differ. Equ. 2016, 2016, 120. [Google Scholar] [CrossRef]
  33. Stević, S. Solvability of a product-type system of difference equations with six parameters. Adv. Nonlinear Anal. 2016. [Google Scholar] [CrossRef]
  34. Stević, S. Solvable subclasses of a class of nonlinear second-order difference equations. Adv. Nonlinear Anal. 2016, 5, 147–165. [Google Scholar] [CrossRef]
  35. Stević, S. Third order product-type systems of difference equations solvable in closed form. Electron. J. Differ. Equ. 2016, 2016, 1–11. [Google Scholar]
  36. Stević, S. Product-type system of difference equations with complex structure of solutions. Adv. Differ. Equ. 2017, 2017. [Google Scholar] [CrossRef]
  37. Stević, S. Solution to the solvability problem for a class of product-type systems of difference equations. Adv. Differ. Equ. 2017, 2017. [Google Scholar] [CrossRef]
  38. Stević, S. Solvable product-type system of difference equations whose associated polynomial is of the fourth order. Electron. J. Qual. Theory Differ. Equ. 2017, 2017. [Google Scholar] [CrossRef]
  39. Stević, S.; Alghamdi, M.A.; Alotaibi, A.; Elsayed, E.M. Solvable product-type system of difference equations of second order. Electron. J. Differ. Equ. 2015, 2015, 1–20. [Google Scholar] [CrossRef]
  40. Stević, S.; Iričanin, B.; Šmarda, Z. On a product-type system of difference equations of second order solvable in closed form. J. Inequal. Appl. 2015, 2015. [Google Scholar] [CrossRef]
  41. Stević, S.; Iričanin, B.; Šmarda, Z. Solvability of a close to symmetric system of difference equations. Electron. J. Differ. Equ. 2016, 2016, 1–13. [Google Scholar]
  42. Stević, S.; Iričanin, B.; Šmarda, Z. Two-dimensional product-type system of difference equations solvable in closed form. Adv. Differ. Equ. 2016, 2016. [Google Scholar] [CrossRef]
  43. Stević, S.; Ranković, D. On a practically solvable product-type system of difference equations of second order. Electron. J. Qual. Theory Differ. Equ. 2016, 2016. [Google Scholar] [CrossRef]
  44. Faddeyev, D.K. Lectures on Algebra; Nauka: Moscow, Russia, 1984. [Google Scholar]
  45. Rees, E.L. Graphical discussion of the roots of a quartic equation. Am. Math. Mon. 1922, 29, 51–55. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Stević, S. Solvability of the Class of Two-Dimensional Product-Type Systems of Difference Equations of Delay-Type (1, 3, 1, 1). Symmetry 2017, 9, 200. https://doi.org/10.3390/sym9100200

AMA Style

Stević S. Solvability of the Class of Two-Dimensional Product-Type Systems of Difference Equations of Delay-Type (1, 3, 1, 1). Symmetry. 2017; 9(10):200. https://doi.org/10.3390/sym9100200

Chicago/Turabian Style

Stević, Stevo. 2017. "Solvability of the Class of Two-Dimensional Product-Type Systems of Difference Equations of Delay-Type (1, 3, 1, 1)" Symmetry 9, no. 10: 200. https://doi.org/10.3390/sym9100200

APA Style

Stević, S. (2017). Solvability of the Class of Two-Dimensional Product-Type Systems of Difference Equations of Delay-Type (1, 3, 1, 1). Symmetry, 9(10), 200. https://doi.org/10.3390/sym9100200

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