MATLAB example (you can cut and paste this into MATLAB):
A=[
2 0 0
-1 3 1
-1 1 3]
disp('all eigenvalues')
lambdas=eig(A)
lambda1= lambdas(1);
D1=A-lambda1*eye(3)
rref(D1)
disp('eigenvector for first eigenvalue')
v1=[-ans(1:2,3)
eye(1)]
disp('check eigenvector equation, i.e. norm of Av1-lambda1 v1')
norm(A*v1-lambda1*v1)
disp('now for second eigenvalue')
lambda2= lambdas(2);
D2=A-lambda2*eye(3)
rref(D2)
disp('two linearly indep. eigenvectors for second eigenvalue')
v2=[-ans(1,2:3)
eye(2)]
disp('check eigenvector equation, i.e. norm of Av2-lambda1 v2')
norm(A*v2-lambda2*v2)