In this paper, we represent that the following three-dimensional system of difference equations
x
n
+
1
=
α
y
n
+
a
y
n
y
n
−
β
z
n
−
1
,
y
n
+
1
=
β
z
n
+
b
z
n
z
n
−
γ
x
n
−
1
,
z
n
+
1
=
γ
x
n
+
c
x
n
x
n
−
α
y
n
−
1
,
n
∈
ℕ
0
,
$$\matrix{{{x_{n + 1}} = \alpha {y_n} + {{a{y_n}} \over {{y_n} - \beta {z_{n - 1}}}},\quad {y_{n + 1}} = \beta {z_n} + {{b{z_n}} \over {{z_n} - \gamma {x_{n - 1}}}},\quad {z_{n + 1}} = \gamma {x_n} + {{c{x_n}} \over {{x_n} - \alpha {y_{n - 1}}}},\qquad n \in {{\mathbb N}_0},} \cr} $$
where the parameters a, b, c, α, β, γ and the initial values x
−i
, y
−i
, z
−i
, i ∈ {0, 1}, are real numbers, can be solved in closed form by using transformation. We analyzed the solutions in 10 different cases depending on whether the parameters are zero or nonzero. It is noteworthy to depict that the solutions of some particular cases of this system are presented in terms of generalized Fibonacci numbers. Note that our results considerably extend and improve some recent results in the literature.