<s>[INST] Solve system of first order differential equations I have to solve differential systems like this:
$$
\left\{ 
\begin{array}{c}
x' = 3x - y + z \\
y' = x + 5y - z \\
z' = x - y + 3z
\end{array}
\right. 
$$
Until now I computed the eigenvalues $k = \{2,4,5\}$ by solving the equation resulted from this determinant of this matrix:
$$
\begin{pmatrix}
3 && -1 && 1 \\
1 && 5 && -1 \\
1 && -1 && 3
\end{pmatrix} - kI_3 = 
\begin{pmatrix}
3-k && -1 && 1 \\
1 && 5-k && -1 \\
1 && -1 && 3-k
\end{pmatrix}
$$
I don't know what to do next.
NOTE: This is an example but from it I want to learn the method to solve any system of this kind. I learn better from particular examples than directly from generalization.
 [/INST] To find the eigenvalues of: $$A =\begin{bmatrix}3 & -1 & 1\\1 & 5 & -1\\1 & -1 &3\end{bmatrix},$$  
we set up $|A - \lambda I| = 0$ and solve the characteristic polynomial, so we have:
$$|A -\lambda I| = \begin{bmatrix} 3-\lambda & -1 & 1\\1 & 5-\lambda & -1\\1 & -1 &3 -\lambda \end{bmatrix} = 0$$  
From this, we get the characteristic polynomial as: $$-\lambda^3+11 \lambda^2-38 \lambda+40 = -(\lambda-5) (\lambda-4) (\lambda-2)= 0$$
This gives us three eigenvalues: $λ_1 = 2, λ_2 = 4$  and $λ_3 = 5$
To find the eigenvectors, we set up and solve $[A - \lambda_i I]v_i = 0$, so lets do $\lambda_1 = 2$. We have:
$$[A - \lambda_i I]v_i = \begin{bmatrix} 3-2 & -1 & 1\\1 & 5-2 & -1\\1 & -1 &3 -2 \end{bmatrix}v_1 = \begin{bmatrix} 1 & -1 & 1\\1 & 3 & -1\\1 & -1 & 1 \end{bmatrix}v_1  = 0$$
This leads to the row-reduced-echelon form:
$$\begin{bmatrix} 1 & 0 & \dfrac{1}{2}\\0 & 1 & -\dfrac{1}{2}\\0 & 0 & 0 \end{bmatrix}v_1  = 0$$
This gives us a solution of:
$b = \dfrac{1}{2}c \rightarrow c = 2 \rightarrow b = 1$, and $a = -\dfrac{1}{2} c \rightarrow a = -1$.
Thus the eigenvector for the eigenvalue $\lambda_1 = 2$ is $v_1 = (-1, 1, 2)$.
If we repeat this process two more times for the other two eigenvalues, we end up with the eigenvalue/eigenvector pairs:
$$\lambda_1 = 2, ~v_1 = (-1, 1, 2)$$  
$$\lambda_2 = 4, ~v_2 = (1, 0, 1)$$ 
$$\lambda_3 = 5, ~v_3 = (1, -1, 1)$$  
Note that different approaches are needed when we get into repeated eigenvalues, but that is for another problem. 
Try to derive those last two eigenvectors and report back if they do not work out.
To write out the solution for the three equation, we would use a linear combination of the the eigenvalues, eigenvectors and some unknown constants, as:
$$X(t) = \begin{bmatrix} x(t) \\ y(t) \\ z(t) \end{bmatrix} = c_1e^{2t} \begin{pmatrix} -1 \\ 1 \\ 2 \end{pmatrix} + c_2 e^{4t} \begin{pmatrix} 1 \\ 0 \\ 1 \end{pmatrix} + c_3 e^{5t} \begin{pmatrix} 1 \\ -1 \\ 1 \end{pmatrix}$$
</s>