Skip to main content

On a product-type system of difference equations of second order solvable in closed form

Abstract

It is shown that the following system of difference equations

$$z_{n+1}=\frac{z_{n}^{a}}{w_{n-1}^{b}},\qquad w_{n+1}=\frac{w_{n}^{c}}{z_{n-1}^{d}}, \quad n\in \mathbb {N}_{0}, $$

where \(a,b,c,d\in \mathbb {Z}\), \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\), is solvable in closed form.

1 Introduction

Recently there has been a great interest in studying nonlinear difference equations and systems not stemming from differential ones (see, e.g., [130]). The old area of solving difference equations and systems has re-attracted recent attention (see, e.g., [16, 12, 15, 1926, 2830]). Recent Stević’s idea of transforming complicated equations and systems into simpler solvable ones, used for the first time in explaining the solvability of the equation appearing in [6] (an extension of the original result can be found in [20], see also [22]), was employed in several other papers (see, e.g., [1, 2, 4, 12, 15, 19, 21, 2426, 29, 30] and the related references therein). Another area of some recent interest, essentially initiated by Papaschinopoulos and Schinas, is studying symmetric and close to symmetric systems of difference equations (see, e.g., [5, 810, 13, 14, 19, 21, 23, 25, 2729]).

Stević also essentially triggered a systematic study of non-rational concrete difference equations and systems, from one side those obtained by using the translation operator (see, e.g., [16] and also [11]) and from the other side those obtained by using max-type operators (see, e.g., [17, 18, 27]), see also the related references cited therein. We would like to point out that for the equations and systems in [1618, 27] only long-term behavior of their positive solutions are studied. For instance, the boundedness of positive solutions to the system

$$\begin{aligned} x_{n+1}=\max \biggl\{ a,\frac{y_{n}^{p}}{x_{n-1}^{q}} \biggr\} ,\qquad y_{n+1}=\max \biggl\{ a,\frac{x_{n}^{p}}{y_{n-1}^{q}} \biggr\} ,\quad n\in \mathbb {N}_{0}, \end{aligned}$$
(1)

with \(\min\{a, p, q\}>0\), was investigated in [27]. System (1) is obviously obtained from the next product-type one

$$\begin{aligned} x_{n+1}=\frac{y_{n}^{p}}{x_{n-1}^{q}},\qquad y_{n+1}=\frac{x_{n}^{p}}{y_{n-1}^{q}}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(2)

by acting with the max-type operator \(m_{a}(t)=\max\{a,t\}\) onto the right-hand sides of both equations in (2) (see also [17] and [18] for related scalar equations). Note that for the case of positive initial values, system (2) can be solved by taking the logarithm to the both sides of both equations therein, since this transforms the system to a linear second order system of difference equations with constant coefficients, which is solvable. Note that the method does not work if initial values are not positive. Let us also mention here that positive solutions to difference equations and systems are often studied since many real-life models produce such solutions (see, e.g., [7, 15, 31]). It is also interesting to note that there are max-type systems of difference equations which are solvable (see [23]). Finally, we want to note that the long-term behavior of solutions to product-type systems and those obtained from them by acting with some ‘reasonable good’ transformations are frequently closely related, which is another reason for studying these systems.

Hence, a natural problem is to investigate the solvability of product-type difference equations and systems with real and/or complex initial values. In [26], Stević and his collaborators started studying the problem with an approach different from the ones in [5, 21, 24, 25], but which can be regarded as a modification of some of the methods in [1618, 27]. They showed therein that the system

$$\begin{aligned} z_{n+1}=\frac{w_{n}^{a}}{z_{n-1}^{b}},\qquad w_{n+1}=\frac{z_{n}^{c}}{w_{n-1}^{d}}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(3)

where \(a,b,c,d\in \mathbb {Z}\) and \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\), is solvable in closed form and presented numerous applications of obtained formulas.

In this paper we continue our investigation by studying the solvability of the following system of difference equations:

$$\begin{aligned} z_{n+1}=\frac{z_{n}^{a}}{w_{n-1}^{b}},\qquad w_{n+1}=\frac{w_{n}^{c}}{z_{n-1}^{d}}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(4)

where \(a,b,c,d\in \mathbb {Z}\) and \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\).

Let us mention here that although systems (3) and (4) are similar in appearance, the methods used in dealing with them are quite different.

It is easy to see that the domain of undefinable solutions [24] to system (4) is the set

$${\mathcal{U}}=\bigl\{ (z_{-1}, z_{0}, w_{-1}, w_{0})\in \mathbb {C}^{4}: z_{-1}=0 \mbox{ or } z_{0}=0 \mbox{ or } w_{-1}=0 \mbox{ or } w_{0}=0\bigr\} . $$

Hence, from now on we will assume that our initial values belong to the set \(\mathbb {C}^{4}\setminus{\mathcal{U}}\).

A solution \((z_{n},w_{n})_{n\ge-1}\) of system (4) is called periodic (or eventually periodic) with period \(p\in \mathbb {N}\) if there is \(n_{0}\ge-1\) such that

$$(z_{n+p},w_{n+p})=(z_{n},w_{n}) \quad \mbox{for }n\geq n_{0}. $$

Period p is prime if there is no \(\hat{p}\in \mathbb {N}\), \(\hat{p}< p\) which is a period of the solution. For \(p=1\), the solution is called eventually constant (see, e.g., [32]). For some results on the topic, see, e.g., [7, 15] and the related references therein. If it is said that a solution of system (4) is periodic with period p, it will need not mean that it is prime.

A system of difference equations of the form

$$\begin{aligned} &z_{n}=f(z_{n-1},\ldots,z_{n-k},w_{n-1}, \ldots,w_{n-k}) \\ &w_{n}=g(z_{n-1},\ldots,z_{n-k},w_{n-1}, \ldots,w_{n-k}),\quad n\in \mathbb {N}_{0}, \end{aligned}$$

where \(k\in \mathbb {N}\), is said to be solvable in closed form if its general solution can be found in terms of initial values \(z_{-i}\), \(w_{-i}\), \(i=\overline {1,k}\), delay k and index n only.

2 Main result

The main result in this paper is proved in this section.

Theorem 1

Assume that \(a,b,c,d\in \mathbb {Z}\) and \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\setminus\{0\}\). Then system (4) is solvable in closed form.

Proof

Case \(b=0\). In this case system (4) becomes

$$\begin{aligned} z_{n+1}=z_{n}^{a},\qquad w_{n+1}= \frac{w_{n}^{c}}{z_{n-1}^{d}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(5)

From the first equation in (5) we easily obtain

$$\begin{aligned} z_{n}=z_{0}^{a^{n}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(6)

Employing (6) into the second equation in (5), we get

$$\begin{aligned} w_{n}=\frac{w_{n-1}^{c}}{z_{0}^{da^{n-2}}} \end{aligned}$$
(7)

for \(n\ge 2\).

Hence, by using (7), we have that

$$\begin{aligned} w_{n}&=\frac{1}{z_{0}^{da^{n-2}}} \biggl(\frac {w_{n-2}^{c}}{z_{0}^{da^{n-3}}} \biggr)^{c}=\frac {w_{n-2}^{c^{2}}}{z_{0}^{da^{n-2}+dca^{n-3}}} \\ &=\frac{1}{z_{0}^{da^{n-2}+dca^{n-3}}} \biggl(\frac {w_{n-3}^{c}}{z_{0}^{da^{n-4}}} \biggr)^{c^{2}}= \frac {w_{n-3}^{c^{3}}}{z_{0}^{da^{n-2}+dca^{n-3}+dc^{2}a^{n-4}}} \end{aligned}$$

for \(n\ge4\).

Assume that we have proved

$$\begin{aligned} w_{n}=\frac{w_{n-k}^{c^{k}}}{z_{0}^{da^{n-2}+dca^{n-3}+dc^{2}a^{n-4}+\cdots +dc^{k-1}a^{n-k-1}}} \end{aligned}$$
(8)

for \(n\ge k+1\).

Then, by using (7) with \(n\to n-k\) into (8), we get

$$\begin{aligned} w_{n}&=\frac{1}{z_{0}^{da^{n-2}+dca^{n-3}+dc^{2}a^{n-4}+\cdots +dc^{k-1}a^{n-k-1}}} \biggl(\frac{w_{n-k-1}^{c}}{z_{0}^{da^{n-k-2}}} \biggr)^{c^{k}} \\ &=\frac {w_{n-k-1}^{c^{k+1}}}{z_{0}^{da^{n-2}+dca^{n-3}+dc^{2}a^{n-4}+\cdots +dc^{k-1}a^{n-k-1}+dc^{k}a^{n-k-2}}} \end{aligned}$$
(9)

for \(n\ge k+2\).

From (7), (9) and the method of induction we see that (8) holds for every k such that \(1\le k\le n-1\).

By taking \(k=n-1\) into (8) we get

$$\begin{aligned} w_{n}&=\frac{w_{1}^{c^{n-1}}}{z_{0}^{da^{n-2}+dca^{n-3}+dc^{2}a^{n-4}+\cdots +dc^{n-2}}} \end{aligned}$$
(10)

for \(n\ge2\).

Now we have two subcases to consider.

Subcase \(a\ne c\). In this case from (10) we get

$$\begin{aligned} w_{n}=\frac{w_{1}^{c^{n-1}}}{z_{0}^{d\frac{a^{n-1}-c^{n-1}}{a-c}}},\quad n\ge2. \end{aligned}$$
(11)

Using the next relation

$$\begin{aligned} w_{1}=\frac{w_{0}^{c}}{z_{-1}^{d}} \end{aligned}$$
(12)

in (11) we get

$$\begin{aligned} w_{n}=\frac{w_{0}^{c^{n}}}{z_{0}^{d\frac {a^{n-1}-c^{n-1}}{a-c}}z_{-1}^{dc^{n-1}}},\quad n\in \mathbb {N}. \end{aligned}$$
(13)

Subcase \(a=c\). In this case from (10) we get

$$\begin{aligned} w_{n}=\frac{w_{1}^{a^{n-1}}}{z_{0}^{d(n-1)a^{n-2}}} \end{aligned}$$
(14)

for \(n\ge2\).

Using (12) with \(a=c\) into (14), we get

$$\begin{aligned} w_{n}=\frac{w_{0}^{a^{n}}}{z_{0}^{d(n-1)a^{n-2}}z_{-1}^{da^{n-1}}},\quad n\in \mathbb {N}. \end{aligned}$$
(15)

Case \(d=0\). In this case system (4) becomes

$$\begin{aligned} z_{n+1}=\frac{z_{n}^{a}}{w_{n-1}^{b}}, \qquad w_{n+1}=w_{n}^{c}, \quad n\in \mathbb {N}_{0}. \end{aligned}$$
(16)

From the second equation in (16) we have that

$$\begin{aligned} w_{n}=w_{0}^{c^{n}}, \quad n\in \mathbb {N}_{0}. \end{aligned}$$
(17)

Employing (17) into the first equation in (16), we get

$$\begin{aligned} z_{n}=\frac{z_{n-1}^{a}}{w_{0}^{bc^{n-2}}} \end{aligned}$$
(18)

for \(n\ge 2\).

Hence, by using (18), we have

$$\begin{aligned} z_{n}&=\frac{1}{w_{0}^{bc^{n-2}}} \biggl(\frac {z_{n-2}^{a}}{w_{0}^{bc^{n-3}}} \biggr)^{a}=\frac {z_{n-2}^{a^{2}}}{w_{0}^{bc^{n-2}+bac^{n-3}}} \\ &=\frac{1}{w_{0}^{bc^{n-2}+bac^{n-3}}} \biggl(\frac {z_{n-3}^{a}}{w_{0}^{bc^{n-4}}} \biggr)^{a^{2}}= \frac {z_{n-3}^{a^{3}}}{w_{0}^{bc^{n-2}+bac^{n-3}+ba^{2}c^{n-4}}} \end{aligned}$$

for \(n\ge4\).

Assume that we have proved

$$\begin{aligned} z_{n}&=\frac{z_{n-k}^{a^{k}}}{w_{0}^{bc^{n-2}+bac^{n-3}+ba^{2}c^{n-4}+\cdots +ba^{k-1}c^{n-k-1}}} \end{aligned}$$
(19)

for \(n\ge k+1\).

Then, by using (18) with \(n\to n-k\) into (19), we get

$$\begin{aligned} z_{n}&=\frac{1}{w_{0}^{bc^{n-2}+bac^{n-3}+ba^{2}c^{n-4}+\cdots +ba^{k-1}c^{n-k-1}}} \biggl(\frac{z_{n-k-1}^{a}}{w_{0}^{bc^{n-k-2}}} \biggr)^{a^{k}} \\ &=\frac {z_{n-k-1}^{a^{k+1}}}{w_{0}^{bc^{n-2}+bac^{n-3}+ba^{2}c^{n-4}+\cdots +ba^{k-1}c^{n-k-1}+ba^{k}c^{n-k-2}}} \end{aligned}$$
(20)

for \(n\ge k+2\).

From (18), (20) and the method of induction we see that (19) holds for every k such that \(1\le k\le n-1\).

By taking \(k=n-1\) into (19) we get

$$\begin{aligned} z_{n}=\frac{z_{1}^{a^{n-1}}}{w_{0}^{bc^{n-2}+bac^{n-3}+ba^{2}c^{n-4}+\cdots +bca^{n-3}+ba^{n-2}}} \end{aligned}$$
(21)

for \(n\ge2\).

Now we have two subcases to consider.

Subcase \(a\ne c\). In this case from (21) we get

$$\begin{aligned} z_{n}=\frac{z_{1}^{a^{n-1}}}{w_{0}^{b\frac{a^{n-1}-c^{n-1}}{a-c}}},\quad n\ge2. \end{aligned}$$
(22)

Using the next relation

$$\begin{aligned} z_{1}=\frac{z_{0}^{a}}{w_{-1}^{b}} \end{aligned}$$
(23)

in (22) we get

$$\begin{aligned} z_{n}=\frac{z_{0}^{a^{n}}}{w_{0}^{b\frac {a^{n-1}-c^{n-1}}{a-c}}w_{-1}^{ba^{n-1}}},\quad n\in \mathbb {N}. \end{aligned}$$
(24)

Subcase \(a=c\). In this case from (21) we get

$$\begin{aligned} z_{n}=\frac{z_{1}^{a^{n-1}}}{w_{0}^{b(n-1)a^{n-2}}} \end{aligned}$$
(25)

for \(n\ge2\).

Using (23) into (25) we get

$$\begin{aligned} z_{n}=\frac{z_{0}^{a^{n}}}{w_{0}^{b(n-1)a^{n-2}}w_{-1}^{ba^{n-1}}},\quad n\in \mathbb {N}. \end{aligned}$$
(26)

Case \(bd\ne0\). First note that from the first equation in (4), for every well-defined solution, we have that

$$\begin{aligned} w_{n-1}^{b}=\frac{z_{n}^{a}}{z_{n+1}},\quad n\in \mathbb {N}_{0}, \end{aligned}$$
(27)

while from the second one it follows that

$$\begin{aligned} w_{n+1}^{b}=\frac{w_{n}^{bc}}{z_{n-1}^{bd}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(28)

Using (27) into (28) we obtain

$$\frac{z_{n+2}^{a}}{z_{n+3}}=\frac {z_{n+1}^{ac}}{z_{n+2}^{c}z_{n-1}^{bd}},\quad n\in \mathbb {N}_{0}, $$

which can be written as

$$\begin{aligned} z_{n+3}=z_{n+2}^{a+c}z_{n+1}^{-ac}z_{n-1}^{bd}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(29)

which is a fourth order product-type difference equation.

Note also that

$$\begin{aligned} z_{1}=\frac{z_{0}^{a}}{w_{-1}^{b}}\quad\mbox{and}\quad z_{2}= \frac{z_{1}^{a}}{w_{0}^{b}}=\frac{z_{0}^{a^{2}}}{w_{0}^{b}w_{-1}^{ab}}. \end{aligned}$$
(30)

Let

$$\begin{aligned} a_{1}=a+c,\qquad b_{1}=-ac,\qquad c_{1}=0, \qquad d_{1}=bd. \end{aligned}$$
(31)

Then equation (29) can be written as

$$\begin{aligned} z_{n+3}=z_{n+2}^{a_{1}}z_{n+1}^{b_{1}}z_{n}^{c_{1}}z_{n-1}^{d_{1}}, \quad n\in \mathbb {N}_{0}. \end{aligned}$$
(32)

From (32) with \(n\to n-1\) we get

$$\begin{aligned} z_{n+2}=z_{n+1}^{a_{1}}z_{n}^{b_{1}}z_{n-1}^{c_{1}}z_{n-2}^{d_{1}}, \quad n\in \mathbb {N}. \end{aligned}$$
(33)

Employing (33) into (32) we get

$$\begin{aligned} z_{n+3}&=\bigl(z_{n+1}^{a_{1}}z_{n}^{b_{1}}z_{n-1}^{c_{1}} z_{n-2}^{d_{1}}\bigr)^{a_{1}}z_{n+1}^{b_{1}}z_{n}^{c_{1}}z_{n-1}^{d_{1}} \\ &=z_{n+1}^{a_{1}a_{1}+b_{1}}z_{n}^{a_{1}b_{1}+c_{1}}z_{n-1}^{a_{1}c_{1}+d_{1}} z_{n-2}^{a_{1}d_{1}} \\ &=z_{n+1}^{a_{2}}z_{n}^{b_{2}}z_{n-1}^{c_{2}}z_{n-2}^{d_{2}} \end{aligned}$$
(34)

for \(n\in \mathbb {N}\), where

$$\begin{aligned} a_{2}:=a_{1}a_{1}+b_{1},\qquad b_{2}:=a_{1}b_{1}+c_{1}, \qquad c_{2}:=a_{1}c_{1}+d_{1},\qquad d_{2}:=a_{1}d_{1}. \end{aligned}$$
(35)

From (32) with \(n\to n-2\) we get

$$\begin{aligned} z_{n+1}=z_{n}^{a_{1}}z_{n-1}^{b_{1}}z_{n-2}^{c_{1}}z_{n-3}^{d_{1}} \end{aligned}$$
(36)

for \(n\ge2\).

Employing (36) into (34) we get

$$\begin{aligned} z_{n+3}&=\bigl(z_{n}^{a_{1}}z_{n-1}^{b_{1}} z_{n-2}^{c_{1}}z_{n-3}^{d_{1}} \bigr)^{a_{2}}z_{n}^{b_{2}}z_{n-1}^{c_{2}}z_{n-2}^{d_{2}} \\ &=z_{n}^{a_{1}a_{2}+b_{2}}z_{n-1}^{b_{1}a_{2}+c_{2}} z_{n-2}^{c_{1}a_{2}+d_{2}}z_{n-3}^{d_{1}a_{2}} \\ &=z_{n}^{a_{3}}z_{n-1}^{b_{3}}z_{n-2}^{c_{3}}z_{n-3}^{d_{3}} \end{aligned}$$
(37)

for \(n\ge2\), where

$$\begin{aligned} a_{3}:=a_{1}a_{2}+b_{2},\qquad b_{3}:=b_{1}a_{2}+c_{2}, \qquad c_{3}:=c_{1}a_{2}+d_{2},\qquad d_{3}:=d_{1}a_{2}. \end{aligned}$$
(38)

Assume that for some \(2\le k \le n\), we have proved that

$$\begin{aligned} z_{n+3}=z_{n+3-k}^{a_{k}}z_{n+2-k}^{b_{k}}z_{n+1-k}^{c_{k}}z_{n-k}^{d_{k}} \end{aligned}$$
(39)

for \(n\ge k-1\), and that

$$ \begin{aligned} &a_{k}=a_{1}a_{k-1}+b_{k-1}, \qquad b_{k}=b_{1}a_{k-1}+c_{k-1}, \\ &c_{k}=c_{1}a_{k-1}+d_{k-1},\qquad d_{k}=d_{1}a_{k-1}. \end{aligned} $$
(40)

Then, by using the relation

$$z_{n+3-k}=z_{n+2-k}^{a_{1}}z_{n+1-k}^{b_{1}}z_{n-k}^{c_{1}}z_{n-1-k}^{d_{1}}, $$

for \(n\ge k\), into (39) we obtain

$$\begin{aligned} z_{n+3} &=\bigl(z_{n+2-k}^{a_{1}}z_{n+1-k}^{b_{1}}z_{n-k}^{c_{1}} z_{n-1-k}^{d_{1}}\bigr)^{a_{k}}z_{n+2-k}^{b_{k}}z_{n+1-k}^{c_{k}}z_{n-k}^{d_{k}} \\ &=z_{n+2-k}^{a_{1}a_{k}+b_{k}}z_{n+1-k}^{b_{1}a_{k}+c_{k}} z_{n-k}^{c_{1}a_{k}+d_{k}}z_{n-1-k}^{d_{1}a_{k}} \\ &=z_{n+2-k}^{a_{k+1}}z_{n+1-k}^{b_{k+1}}z_{n-k}^{c_{k+1}}z_{n-1-k}^{d_{k+1}} \end{aligned}$$
(41)

for \(n\ge k\), where

$$ \begin{aligned} &a_{k+1}:=a_{1}a_{k}+b_{k}, \qquad b_{k+1}:=b_{1}a_{k}+c_{k}, \\ &c_{k+1}:=c_{1}a_{k}+d_{k},\qquad d_{k+1}:=d_{1}a_{k}. \end{aligned} $$
(42)

This along with (34), (35) and the method of induction shows that (39) and (40) hold for every \(2\le k\le n+1\).

Hence, for \(k=n+1\), we have

$$\begin{aligned} z_{n+3}&=z_{2}^{a_{n+1}}z_{1}^{b_{n+1}}z_{0}^{c_{n+1}}z_{-1}^{d_{n+1}} \\ &= \biggl(\frac{z_{0}^{a^{2}}}{w_{0}^{b}w_{-1}^{ab}} \biggr)^{a_{n+1}} \biggl(\frac{z_{0}^{a}}{w_{-1}^{b}} \biggr)^{b_{n+1}}z_{0}^{c_{n+1}}z_{-1}^{d_{n+1}} \\ &=z_{0}^{a^{2}a_{n+1}+ab_{n+1}+c_{n+1}}z_{-1}^{d_{n+1}} w_{0}^{-ba_{n+1}}w_{-1}^{-aba_{n+1}-bb_{n+1}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(43)

From the recurrent relations (40) we easily obtain that the sequence \((a_{k})_{k\ge5}\) satisfies the difference equation

$$\begin{aligned} a_{k}=a_{1}a_{k-1}+b_{1}a_{k-2}+c_{1}a_{k-3}+d_{1}a_{k-4}. \end{aligned}$$
(44)

Since \(b_{k-1}=a_{k}-a_{1}a_{k-1}\) and equation (44) is linear, we have that the sequence \((b_{k})_{k\in \mathbb {N}}\) is also a solution to equation (44). From this, the linearity of equation (44) and since \(c_{k-1}=b_{k}-b_{1}a_{k-1}\), we have that the sequence \((c_{k})_{k\in \mathbb {N}}\) is also a solution to equation (44). Finally, since \(d_{k}=d_{1}a_{k-1}\), the linearity of equation (44) shows that \((d_{k})_{k\in \mathbb {N}}\) is also a solution to the equation.

Now, we show that these four sequences can be prolonged for some negative indices of use. This enables easier getting formulas for solutions to system (4).

From (42) with \(k=0\) we get

$$\begin{aligned} a_{1}&=a_{1}a_{0}+b_{0},\qquad b_{1}=b_{1}a_{0}+c_{0}, \qquad c_{1}=c_{1}a_{0}+d_{0},\qquad d_{1}=d_{1}a_{0}. \end{aligned}$$
(45)

Since \(bd=d_{1}\ne0\), from the last equation in (45) we get \(a_{0}=1\). Using this fact in the first three equalities in (45), we get \(b_{0}=c_{0}=d_{0}=0\).

From this and by (42) with \(k=-1\) we get

$$ \begin{aligned} &1=a_{0}=a_{1}a_{-1}+b_{-1}, \qquad 0=b_{0}=b_{1}a_{-1}+c_{-1}, \\ &0=c_{0}=c_{1}a_{-1}+d_{-1},\qquad 0=d_{0}=d_{1}a_{-1}. \end{aligned} $$
(46)

Since \(d_{1}\ne0\), from the last equation in (46) we get \(a_{-1}=0\). Using this fact in other three equalities in (46), we get \(b_{-1}=1\), \(c_{-1}=d_{-1}=0\).

From this and by (42) with \(k=-2\) we get

$$ \begin{aligned} &0=a_{-1}=a_{1}a_{-2}+b_{-2}, \qquad 1=b_{-1}=b_{1}a_{-2}+c_{-2}, \\ &0=c_{-1}=c_{1}a_{-2}+d_{-2},\qquad 0=d_{-1}=d_{1}a_{-2}. \end{aligned} $$
(47)

Since \(d_{1}\ne0\), from the last equation in (47) we get \(a_{-2}=0\). Using this fact in other three equalities in (47), we get \(b_{-2}=0\), \(c_{-2}=1\) and \(d_{-2}=0\).

From this and by (42) with \(k=-3\) we get

$$ \begin{aligned} &0=a_{-2}=a_{1}a_{-3}+b_{-3}, \qquad 0=b_{-2}=b_{1}a_{-3}+c_{-3}, \\ &1=c_{-2}=c_{1}a_{-3}+d_{-3},\qquad 0=d_{-2}=d_{1}a_{-3}. \end{aligned} $$
(48)

Since \(d_{1}\ne0\), from the last equation in (48) we get \(a_{-3}=0\). Using this fact in other three equalities in (48), we get \(b_{-3}=0\), \(c_{-3}=0\) and \(d_{-3}=1\).

Hence, sequences \((a_{k})_{k\ge-3}\), \((b_{k})_{k\ge-3}\), \((c_{k})_{k\ge-3}\) and \((d_{k})_{k\ge-3}\) are solutions to linear difference equation (44) satisfying the following initial conditions:

$$ \begin{aligned} &a_{-3}=0,\qquad a_{-2}=0, \qquad a_{-1}=0,\qquad a_{0}=1; \\ &b_{-3}=0, \qquad b_{-2}=0,\qquad b_{-1}=1, \qquad b_{0}=0; \\ &c_{-3}=0,\qquad c_{-2}=1, \qquad c_{-1}=0,\qquad c_{0}=0; \\ &d_{-3}=1,\qquad d_{-2}=0,\qquad d_{-1}=0, \qquad d_{0}=0, \end{aligned} $$
(49)

respectively.

Since difference equation (44) is solvable, it follows that closed form formulas for \((a_{k})_{k\ge-3}\), \((b_{k})_{k\ge-3}\), \((c_{k})_{k\ge-3}\) and \((d_{k})_{k\ge-3}\) can be found. From this fact and (43) we see that equation (29) is solvable too.

From the second equation in (4), for every well-defined solution, we have that

$$\begin{aligned} z_{n-1}^{d}=\frac{w_{n}^{c}}{w_{n+1}},\quad n\in \mathbb {N}_{0}, \end{aligned}$$
(50)

while from the first one it follows that

$$\begin{aligned} z_{n+1}^{d}=\frac{z_{n}^{ad}}{w_{n-1}^{bd}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(51)

Using (50) into (51) we obtain

$$\frac{w_{n+2}^{c}}{w_{n+3}}=\frac {w_{n+1}^{ac}}{w_{n+2}^{a}w_{n-1}^{bd}}, \quad n\in \mathbb {N}_{0}, $$

which can be written as

$$\begin{aligned} w_{n+3}=w_{n+2}^{a+c}w_{n+1}^{-ac}w_{n-1}^{bd}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(52)

which is nothing but difference equation (29). However, the sequence \((w_{n})_{n\ge-1}\) satisfies the following initial conditions:

$$\begin{aligned} w_{1}=\frac{w_{0}^{c}}{z_{-1}^{d}}\quad\mbox{and}\quad w_{2}= \frac{w_{1}^{c}}{z_{0}^{d}}=\frac{w_{0}^{c^{2}}}{z_{0}^{d}z_{-1}^{cd}}. \end{aligned}$$
(53)

Hence, the above presented procedure can be repeated, and it can be obtained that for \(1\le k \le n+1\),

$$\begin{aligned} w_{n+3}=w_{n+3-k}^{a_{k}}w_{n+2-k}^{b_{k}}w_{n+1-k}^{c_{k}}w_{n-k}^{d_{k}}, \quad n\in \mathbb {N}_{0}, \end{aligned}$$
(54)

where \((a_{k})_{k\in \mathbb {N}}\), \((b_{k})_{k\in \mathbb {N}}\), \((c_{k})_{k\in \mathbb {N}}\) and \((d_{k})_{k\in \mathbb {N}}\) satisfy recurrent relations (40) with initial conditions (31).

From (54) with \(k=n+1\) and by using (53) we get

$$\begin{aligned} w_{n+3}&=w_{2}^{a_{n+1}}w_{1}^{b_{n+1}}w_{0}^{c_{n+1}}w_{-1}^{d_{n+1}} \\ &= \biggl(\frac{w_{0}^{c^{2}}}{z_{0}^{d}z_{-1}^{cd}} \biggr)^{a_{n+1}} \biggl(\frac{w_{0}^{c}}{z_{-1}^{d}} \biggr)^{b_{n+1}}w_{0}^{c_{n+1}}w_{-1}^{d_{n+1}} \\ &=w_{0}^{c^{2}a_{n+1}+cb_{n+1}+c_{n+1}}w_{-1}^{d_{n+1}} z_{0}^{-da_{n+1}}z_{-1}^{-cda_{n+1}-db_{n+1}},\quad n\in \mathbb {N}_{0}. \end{aligned}$$
(55)

Also the sequences \((a_{k})_{k\in \mathbb {N}}\), \((b_{k})_{k\in \mathbb {N}}\), \((c_{k})_{k\in \mathbb {N}}\) and \((d_{k})_{k\in \mathbb {N}}\) satisfy the difference equation (44) with initial conditions in (49), respectively.

As above the solvability of equation (44) shows that closed form formulas for \((a_{k})_{k\ge-3}\), \((b_{k})_{k\ge-3}\), \((c_{k})_{k\ge-3}\) and \((d_{k})_{k\ge-3}\) can be found. This fact along with (55) implies that equation (52) is solvable too. A direct calculation shows that the sequences \((z_{n})_{n\ge-1}\) in (43) and \((w_{n})_{n\ge-1}\) in (55) are solutions to system (4) with initial values \(z_{-1}\), \(z_{0}\), that is, \(w_{-1}\), \(w_{0}\) respectively. Hence, system (4) is also solvable in this case, finishing the proof of the theorem. □

Remark 1

Note that difference equation (44) is not only theoretically but also practically solvable since the characteristic polynomial

$$\begin{aligned} p_{4}(\lambda )=\lambda ^{4}-a_{1}\lambda ^{3}-b_{1} \lambda ^{2}-c_{1}\lambda -d_{1} \end{aligned}$$
(56)

associated to the difference equation is of fourth order, which means that we can explicitly find its roots.

Remark 2

Since we are interested in those initial values \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\) which uniquely define solutions to system (4), to avoid multi-valued solutions to the system, we posed the condition \(a,b,c,d\in \mathbb {Z}\).

From the proof of Theorem 1 we obtain the following corollary.

Corollary 1

Consider system (4) with \(a,b,c,d\in \mathbb {Z}\). Assume that \(z_{-1}, z_{0}, w_{-1}, w_{0}\in \mathbb {C}\setminus\{0\}\). Then the following statements are true.

  1. (a)

    If \(b=0\) and \(a\ne c\), then the general solution to system (4) is given by (6) and (13).

  2. (b)

    If \(b=0\) and \(a=c\), then the general solution to system (4) is given by (6) and (15).

  3. (c)

    If \(d=0\) and \(a\ne c\), then the general solution to system (4) is given by (17) and (24).

  4. (d)

    If \(d=0\) and \(a=c\), then the general solution to system (4) is given by (17) and (26).

  5. (e)

    If \(bd\ne0\), then the general solution to system (4) is given by (43) and (55).

Let \(\lambda _{i}\), \(i=\overline {1,4}\), be the roots of the characteristic polynomial (56) of difference equation (44). If they satisfy the condition

$$\lambda _{i}\ne \lambda _{j}\quad\mbox{for }i\ne j, $$

then it is known that a general solution to equation (44) has the following form:

$$\begin{aligned} u_{n}=\alpha _{1}\lambda _{1}^{n}+\alpha _{2}\lambda _{2}^{n}+\alpha _{3} \lambda _{3}^{n}+\alpha _{4}\lambda _{4}^{n}, \quad n\in \mathbb {N}, \end{aligned}$$
(57)

where \(\alpha _{i}\), \(i=\overline {1,4}\), are arbitrary constants. Since for the case \(d_{1}\ne0\) the solution can be prolonged for nonpositive indices, then we may assume that formula (57) holds also for \(n\ge-3\) (or \(n\ge-4\) if necessary).

In order to find, in this case, a general solution to system (4) in closed form, we will need the following known lemma. We give a proof of it for the completeness and benefit of the reader.

Lemma 1

Assume that \(\lambda _{j}\), \(j=\overline {1,k}\), are pairwise different zeros of the polynomial

$$P(z)=\alpha _{k}z^{k}+\alpha _{k-1}z^{k-1}+ \cdots+\alpha _{1}z+\alpha _{0}. $$

Then

$$\sum_{j=1}^{k}\frac{\lambda _{j}^{l}}{P'(\lambda _{j})}=0 $$

for \(l=\overline {0,k-2}\), and

$$\sum_{j=1}^{k}\frac{\lambda _{j}^{k-1}}{P'(\lambda _{j})}= \frac{1}{\alpha _{k}}. $$

Proof

The functions

$$f_{l}(z)=\frac{z^{l}}{P(z)},\quad l\in \mathbb {N}, $$

are meromorphic on the Riemann sphere. Hence, by the residue theorem, we have that

$$\begin{aligned} \sum_{j=1}^{k}\operatorname{Res}_{z=\lambda _{j}}f_{l}(z) +\operatorname{Res}_{z=\infty}f_{l}(z)=0 \end{aligned}$$
(58)

for every \(l\in \mathbb {N}\).

Now note that the Laurent expansion of \(f_{l}\) at zero is

$$f_{l}(z)=\frac{z^{l}}{\alpha _{k}\prod_{j=1}^{k}(z-\lambda _{j})}=\frac {z^{l-k}}{\alpha _{k}\prod_{j=1}^{k}(1-\lambda _{j}/z)}=\frac{1}{\alpha _{k}}z^{l-k}+ \sum_{s=1}^{\infty}b_{-s}z^{l-k-s} $$

for some complex numbers \(b_{-s}\), \(s\in \mathbb {N}\).

On the other hand, since \(\lambda _{j}\), \(j=\overline {1,k}\), are simple poles of \(f_{l}\), we have that

$$\operatorname{Res}_{z=\lambda _{j}}f_{l}(z)=\frac{\lambda _{j}^{l}}{P'(\lambda _{j})}, \quad j=\overline {1,k}. $$

From this and since \(\operatorname{Res}_{z=\infty}f_{l}(z)\) is equal to the negative value of the coefficient at \(1/z\) in the Laurent expansion, it follows that \(\operatorname{Res}_{z=\infty}f_{l}(z)=0\) when \(l=\overline {0,k-2}\) and \(\operatorname{Res}_{z=\infty}f_{k-1}(z)=-1/\alpha _{k}\). Using these facts in (58) the lemma follows. □

If we apply Lemma 1 to polynomial \(p_{4}\) in (56), and since \(p_{4}(t)=\prod_{l=1}^{4}(t-\lambda _{j})\) (note that \(\alpha _{4}=1\)), we have

$$\sum_{j=1}^{4}\frac{\lambda _{j}^{l}}{p'_{4}(\lambda _{j})}=0 $$

for \(l=\overline {0,2}\), and

$$\sum_{j=1}^{4}\frac{\lambda _{j}^{3}}{p'_{4}(\lambda _{j})}=1. $$

From this, since from (49) we have \(a_{-3}=a_{-2}=a_{-1}=0\) and \(a_{0}=1\), and a general solution of (4) has the form in (57), we obtain

$$\begin{aligned} a_{n}={}&\sum_{j=1}^{4} \frac{\lambda _{j}^{n+3}}{p_{4}'(\lambda _{j})} \\ ={}&\frac{\lambda _{1}^{n+3}}{(\lambda _{1}-\lambda _{2})(\lambda _{1}-\lambda _{3})(\lambda _{1}-\lambda _{4})}+\frac {\lambda _{2}^{n+3}}{(\lambda _{2}-\lambda _{1})(\lambda _{2}-\lambda _{3})(\lambda _{2}-\lambda _{4})} \\ &{}+\frac{\lambda _{3}^{n+3}}{(\lambda _{3}-\lambda _{1})(\lambda _{3}-\lambda _{2})(\lambda _{3}-\lambda _{4})}+\frac {\lambda _{4}^{n+3}}{(\lambda _{4}-\lambda _{1})(\lambda _{4}-\lambda _{2})(\lambda _{4}-\lambda _{3})} \end{aligned}$$
(59)

for \(n\ge-3\).

On the other hand, from (40) we get

$$\begin{aligned} &b_{n}=a_{n+1}-a_{1}a_{n}, \end{aligned}$$
(60)
$$\begin{aligned} &c_{n}=b_{n+1}-b_{1}a_{n}, \end{aligned}$$
(61)
$$\begin{aligned} &d_{n}=d_{1}a_{n-1} \end{aligned}$$
(62)

for \(n\ge-3\).

By using (59) into (60) we get

$$\begin{aligned} b_{n}=\sum_{j=1}^{4} \frac{\lambda _{j}-a_{1}}{p_{4}'(\lambda_{j})}\lambda _{j}^{n+3} \end{aligned}$$
(63)

for \(n\ge-3\).

By using (59) and (63) into (61) we get

$$\begin{aligned} c_{n}=\sum_{j=1}^{4} \frac{(\lambda _{j}-a_{1})\lambda _{j}-b_{1}}{p_{4}'(\lambda_{j})}\lambda _{j}^{n+3} \end{aligned}$$
(64)

for \(n\ge-3\).

By using (59) into (62) we get

$$\begin{aligned} d_{n}=\sum_{j=1}^{4} \frac{d_{1}}{p_{4}'(\lambda_{j})}\lambda _{j}^{n+2} \end{aligned}$$
(65)

for \(n\ge-3\), where we have used the fact that (59) also holds for \(n=-4\) (in fact, we may assume that equality (59) holds for every \(n\ge-s\), for any fixed \(s\in \mathbb {N}\), since due to the assumption \(d_{1}\ne0\), any solution of equation (44) can be prolonged for any nonpositive value of index n).

By using (59), (63), (64) and (65) into (43) and (55), we get formulas for general solutions to system (4) in closed form.

Formulas obtained in this section can be used in describing the long-term behavior of solutions to system (4) in many cases. We will formulate and prove here only one result, just as an example. The formulations and proofs of other results, which are similar and whose proofs use standard techniques, we leave to the reader as some exercises.

Theorem 2

Assume that \(b=c=0\) and \(a,d\in \mathbb {Z}\). Then the following statements hold:

  1. (a)

    If \(a=1\), then every solution to system (4) is eventually constant.

  2. (b)

    If \(a=0\), then \(z_{n}=w_{n}=1\), \(n\ge3\).

  3. (c)

    If \(a=-1\), then every solution to system (4) is two-periodic.

  4. (d)

    If \(a>1\) and \(|z_{0}|<1\), then \(z_{n}\to0\) as \(n\to\infty\).

  5. (e)

    If \(a>1\) and \(|z_{0}|>1\), then \(|z_{n}|\to\infty\) as \(n\to\infty\).

  6. (f)

    If \(a>1\) and \(|z_{0}^{d}|<1\), then \(|w_{n}|\to\infty\) as \(n\to\infty\).

  7. (g)

    If \(a>1\) and \(|z_{0}^{d}|>1\), then \(w_{n}\to0\) as \(n\to\infty\).

  8. (h)

    If \(a<-1\) and \(|z_{0}|<1\), then \(z_{2n}\to0\) as \(n\to\infty\) and \(|z_{2n+1}|\to\infty\) as \(n\to\infty\).

  9. (i)

    If \(a<-1\) and \(|z_{0}|>1\), then \(z_{2n+1}\to0\) as \(n\to\infty\) and \(|z_{2n}|\to\infty\) as \(n\to\infty\).

  10. (j)

    If \(a<-1\) and \(|z_{0}^{d}|>1\), then \(w_{2n}\to0\) as \(n\to\infty\) and \(|w_{2n+1}|\to\infty\) as \(n\to\infty\).

  11. (k)

    If \(a<-1\) and \(|z_{0}^{d}|<1\), then \(w_{2n+1}\to0\) as \(n\to\infty\) and \(|w_{2n}|\to\infty\) as \(n\to\infty\).

Proof

(a) If we replace \(a=1\) and \(c=0\) in (6) and (7), we obtain \(z_{n}=z_{0}\) and \(w_{n+2}=1/z_{0}^{d}\), \(n\in \mathbb {N}_{0}\), from which the statement follows.

(b) By replacing \(a=0\) and \(c=0\) in (6) and (7), we get \(z_{n}=1\), \(n\in \mathbb {N}\) and \(w_{n}=1\), \(n\ge3\), from which the statement follows.

(c) By replacing \(a=-1\) and \(c=0\) in (6) and (7), we get \(z_{2n}=z_{0}\), \(z_{2n+1}=\frac{1}{z_{0}}\), \(w_{2n}=1/{z_{0}^{d}}\) and \(w_{2n+1}=z_{0}^{d}\), \(n\in \mathbb {N}\), from which the statement follows.

(d)-(k) From (6) and (7) with \(c=0\) we get

$$\begin{aligned} z_{n}=z_{0}^{a^{n}}, \qquad w_{n}= \frac{1}{z_{0}^{da^{n-2}}},\quad n\ge 2. \end{aligned}$$
(66)

Using the formulas in (66) all these statements easily follow. □

References

  1. Aloqeili, M: Dynamics of a kth order rational difference equation. Appl. Math. Comput. 181, 1328-1335 (2006)

    Article  MATH  MathSciNet  Google Scholar 

  2. Aloqeili, M: Dynamics of a rational difference equation. Appl. Math. Comput. 176, 768-774 (2006)

    Article  MATH  MathSciNet  Google Scholar 

  3. Andruch-Sobilo, A, Migda, M: On the rational recursive sequence \(x_{n+1}=ax_{n-1}/(b+cx_{n}x_{n-1})\). Tatra Mt. Math. Publ. 43, 1-9 (2009)

    MATH  MathSciNet  Google Scholar 

  4. Bajo, I, Liz, E: Global behaviour of a second-order nonlinear difference equation. J. Differ. Equ. Appl. 17, 1471-1486 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  5. Berg, L, Stević, S: On some systems of difference equations. Appl. Math. Comput. 218, 1713-1718 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  6. Cinar, C: On the positive solutions of difference equation. Appl. Math. Comput. 150(1), 21-24 (2004)

    Article  MATH  MathSciNet  Google Scholar 

  7. Grove, EA, Ladas, G: Periodicities in Nonlinear Difference Equations. Chapman & Hall/CRC Press, Boca Raton (2005)

    MATH  Google Scholar 

  8. Papaschinopoulos, G, Schinas, CJ: On a system of two nonlinear difference equations. J. Math. Anal. Appl. 219(2), 415-426 (1998)

    Article  MATH  MathSciNet  Google Scholar 

  9. Papaschinopoulos, G, Schinas, CJ: On the behavior of the solutions of a system of two nonlinear difference equations. Commun. Appl. Nonlinear Anal. 5(2), 47-59 (1998)

    MATH  MathSciNet  Google Scholar 

  10. Papaschinopoulos, G, Schinas, CJ: Invariants for systems of two nonlinear difference equations. Differ. Equ. Dyn. Syst. 7(2), 181-196 (1999)

    MATH  MathSciNet  Google Scholar 

  11. Papaschinopoulos, G, Schinas, CJ, Stefanidou, G: On the nonautonomous difference equation \(x_{n+1}=A_{n}+(x_{n-1}^{p}/x_{n}^{q})\). Appl. Math. Comput. 217, 5573-5580 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  12. Papaschinopoulos, G, Stefanidou, G: Asymptotic behavior of the solutions of a class of rational difference equations. Int. J. Difference Equ. 5(2), 233-249 (2010)

    MathSciNet  Google Scholar 

  13. Stefanidou, G, Papaschinopoulos, G, Schinas, C: On a system of max difference equations. Dyn. Contin. Discrete Impuls. Syst. 14(6), 885-903 (2007)

    MATH  MathSciNet  Google Scholar 

  14. Stefanidou, G, Papaschinopoulos, G, Schinas, CJ: On a system of two exponential type difference equations. Commun. Appl. Nonlinear Anal. 17(2), 1-13 (2010)

    MATH  MathSciNet  Google Scholar 

  15. Stević, S: A short proof of the Cushing-Henson conjecture. Discrete Dyn. Nat. Soc. 2006, Article ID 37264 (2006)

    Google Scholar 

  16. Stević, S: Boundedness character of a class of difference equations. Nonlinear Anal. TMA 70, 839-848 (2009)

    Article  MATH  Google Scholar 

  17. Stević, S: On a generalized max-type difference equation from automatic control theory. Nonlinear Anal. TMA 72, 1841-1849 (2010)

    Article  MATH  Google Scholar 

  18. Stević, S: On a nonlinear generalized max-type difference equation. J. Math. Anal. Appl. 376, 317-328 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  19. Stević, S: On a system of difference equations. Appl. Math. Comput. 218, 3372-3378 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  20. Stević, S: On the difference equation \(x_{n}=x_{n-2}/(b_{n}+c_{n}x_{n-1}x_{n-2})\). Appl. Math. Comput. 218, 4507-4513 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  21. Stević, S: On a third-order system of difference equations. Appl. Math. Comput. 218, 7649-7654 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  22. Stević, S: On the difference equation \(x_{n}=x_{n-k}/(b+cx_{n-1}\cdots x_{n-k})\). Appl. Math. Comput. 218, 6291-6296 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  23. Stević, S: Solutions of a max-type system of difference equations. Appl. Math. Comput. 218, 9825-9830 (2012)

    Article  MATH  MathSciNet  Google Scholar 

  24. Stević, S: Domains of undefinable solutions of some equations and systems of difference equations. Appl. Math. Comput. 219, 11206-11213 (2013)

    Article  MATH  MathSciNet  Google Scholar 

  25. Stević, S: Representation of solutions of bilinear difference equations in terms of generalized Fibonacci sequences. Electron. J. Qual. Theory Differ. Equ. 2014, Article No. 67 (2014)

    Article  Google Scholar 

  26. Stević, S, Alghamdi, MA, Alotaibi, A, Elsayed, EM: Solvable product-type system of difference equations of second order. Electron. J. Differ. Equ. 2015, Article No. 169 (2015)

    Article  Google Scholar 

  27. Stević, S, Alghamdi, MA, Alotaibi, A, Shahzad, N: Boundedness character of a max-type system of difference equations of second order. Electron. J. Qual. Theory Differ. Equ. 2014, Article No. 45 (2014)

    Article  Google Scholar 

  28. Tollu, DT, Yazlik, Y, Taskara, N: On fourteen solvable systems of difference equations. Appl. Math. Comput. 233, 310-319 (2014)

    Article  MathSciNet  Google Scholar 

  29. Stević, S, Diblik, J, Iričanin, B, Šmarda, Z: On a solvable system of rational difference equations. J. Differ. Equ. Appl. 20(5-6), 811-825 (2014)

    Article  MATH  Google Scholar 

  30. Stević, S, Diblik, J, Iričanin, B, Šmarda, Z: Solvability of nonlinear difference equations of fourth order. Electron. J. Differ. Equ. 2014, Article No. 264 (2014)

    Article  Google Scholar 

  31. Brauer, F, Castillo-Chavez, C: Mathematical Models in Population Biology and Epidemiology. Springer, Berlin (2012)

    Book  MATH  Google Scholar 

  32. Iričanin, B, Stević, S: Eventually constant solutions of a rational difference equation. Appl. Math. Comput. 215, 854-856 (2009)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

The work of Stevo Stević is supported by the Serbian Ministry of Education and Science projects III 41025 and III 44006. The work of Bratislav Iričanin is supported by the Serbian Ministry of Education and Science projects III 41025 and OI 171007. The work of Zdeňek Šmarda was realized in CEITEC - Central European Institute of Technology with research infrastructure supported by the project CZ.1.05/1.1.00/02.0068 financed from European Regional Development Fund. He was also supported by the project FEKT-S-14-2200 of Brno University of Technology.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Stevo Stević.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the manuscript.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

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, 327 (2015). https://doi.org/10.1186/s13660-015-0835-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13660-015-0835-9

MSC

Keywords