%
% Example 3 from Sec 5.1
%
A = [1 6; 5 2]
v = [6 -5]'
A * v
(-4)*v
[X,D] = eig(A)
rref([A + 4*eye(2) zeros(2,1)])
%
% Example 4 from Sec 5.1
%
A = [4 -1 6; 2 1 6; 2 -1 8]
v = 2 * [1 2 0]' + 3 * [-3 0 1]'
A * v
2 * v
[X,D] = eig(A)
rref([A - 2*eye(3) zeros(3,1)])