First of all I guess that what you wrote is the Hamiltonian and not the Lagrangian of the system and $\dot{x}$ stays for $p_x$ and $\dot{y}$ stays for $p_y$.
You can decouple the problem redefining $$(X,Y)^t = R(x,y)^t$$ for a suitable $R\in O(2)$ diagonalizing the symmetric matrix in the potential part of your Hamiltonian. This way you see the final Hamiltonian is $$\left(\frac{p^2_X}{2m} + \lambda_+ X^2\right) + \left(\frac{p^2_Y}{2m} + \lambda_- Y^2\right)$$ where $\lambda_\pm$ are the eigenvalues of the above symmetric matrix.
In the considered case (up to the dimensional problem already stressed) you find that $ \lambda_+ \lambda_-<0$ (because the determinant of the symmetric matrix is negative nomatter the sign in front of $xy$).
So you have a pair of 1D non-mutually interacting particles, one subjected to a standard harmonic potential and the other subjected to a repulsive harmonic potential (like the one of centrifugal force for a constant angular speed).