EN
We study k th order systems of two rational difference equations $$ x_n = \frac{{\alpha + \sum\nolimits_{i = 1}^k {\beta _i x_{n - 1} + } \sum\nolimits_{i = 1}^k {\gamma _i y_{n - 1} } }} {{A + \sum\nolimits_{j = 1}^k {B_j x_{n - j} + } \sum\nolimits_{j = 1}^k {C_j y_{n - j} } }}, y_n = \frac{{p + \sum\nolimits_{i = 1}^k {\delta _i x_{n - i} + } \sum\nolimits_{i = 1}^k {\varepsilon _i y_{n - i} } }} {{q + \sum\nolimits_{j = 1}^k {D_j x_{n - j} + } \sum\nolimits_{j = 1}^k {E_j y_{n - j} } }} n \in \mathbb{N} $$. In particular, we assume non-negative parameters and non-negative initial conditions, such that the denominators are nonzero. We develop several approaches which allow us to extend well known boundedness results on the k th order rational difference equation to the setting of systems in certain cases.