In this paper, we consider the explicit solution of the following system of nonlinear rational difference equations: x n + 1 = x n - 1 / x n - 1 + r , y n + 1 = x n - 1 y n / x n - 1 y n + r , with initial conditions x - 1 , x 0 and y 0 , which are arbitrary positive real numbers. By doing this, we encounter the hypergeometric function. We also investigate global dynamics of this system. The global dynamics of this system consists of two kind of bifurcations