There are at least a couple of options:

1. Use the eigenvalue/eigenvector diagonalization approach.

2. Eliminate one variable to get a second-order DE.

You can check those out at Chris's DE tutorial. Is a method specified for this problem, or can you use any method?