This is easily solved by using methods from linear algebra. The system can be written as:

The eigenvalues that belong to the transformation matrix are -4 and -7, so you end up with the eigenvectors:

With inverse:

The new transformation matrix becomes:

And the system is transformed to:

Then just solve the two systems independently and transform the solution to the old coordinates by using: