It looks like the system is
Compute the eigenvalues of the coefficient matrix.
For
, the corresponding eigenvector is
such that
Notice that the first row is 1 + 2i times the second row, so
Let
; then
, so that
The eigenvector corresponding to
is the complex conjugate of
.
So, the characteristic solution to the homogeneous system is
The characteristic solution contains
and
, both of which are linearly independent to
and
. So for the nonhomogeneous part, we consider the ansatz particular solution
Differentiating this and substituting into the ODE system gives
Then the general solution to the system is