Once the characteristic
equation hands you an eigenvalue \lambda, the matching
eigenvectors are the non-zero solutions of
(A - \lambda I)\vec{v} = \vec{0}.
Subtract \lambda from the diagonal of A,
then find the direction that gets sent to zero — a
linear system
you already know how to solve. Because A - \lambda I is singular by
construction, there's a whole line of solutions: every eigenvector comes with all its
scalar multiples, so we just report the direction.