Let $f:[0,1]^2 \rightarrow \mathbb{R}^2,$ where $f_1(x,y) = g(y)-x$ and $f_2(x,y) = g(x)-y.$ Here $g(\cdot)$ is a strictly decreasing polynomial function such that $g(0)=1$ and $g(1)=0.$ I am interested in analyzing the asymptotic behavior of the following system of differential equations:
\begin{equation} \dot{x} = f_1(x,y) \ \text{and} \ \dot{y} = f_2(x,y). \end{equation}
I want to show that the above system's limiting behavior will be to one of the stationary states i.e., rule out limit cycles and other complicated behavior. Since $\frac{\partial f_1}{\partial x} + \frac{\partial f_2}{\partial y} = -2 \neq 0,$ can I use the Bendixson–Dulac theorem to rule out limit cycles and conclude that the above system converges to one of the stationary states?
The domain in the statement of Bendixson–Dulac's theorem is usually an open set. However my domain of interest $[0,1]^2$ is not open and so am confused. I would really appreciate help on this.