Eigenvalue, Eigenvector Analysis of a Square Matrix A

Example 2.

`Example2.  Eigenvalue, eigenvector analysis of `*A*`:`

A = matrix([[3, 0, -5], [1/5, -1, 0], [1, 1, -2]]), `  `*(A-I*lambda) = matrix([[3-lambda, 0, -5], [1/5, -1-lambda, 0], [1, 1, -2-lambda]])

det(A-I*lambda)*` = `*lambda*(lambda^2-2) = 0

evecsAval1 = 0

evecsAmul1 = 1

evecsAvec1 = vector([5, 1, 3])

lambda = 0

A-I*lambda = matrix([[3, 0, -5], [1/5, -1, 0], [1, 1, -2]])

matrix([[3, 0, -5], [1/5, -1, 0], [1, 1, -2]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[1, 0, -5/3], [0, 1, -1/3], [0, 0, 0]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[a-5/3*c], [b-1/3*c], [0]]) = matrix([[0], [0], [0]])

v[1] = matrix([[5], [1], [3]])

evecsAval2 = 2^(1/2)

evecsAmul2 = 1

evecsAvec2 = vector([5*2^(1/2)+5, 1, 1+2*2^(1/2)])

lambda = 2^(1/2)

A-I*lambda = matrix([[3-2^(1/2), 0, -5], [1/5, -1-2^(1/2), 0], [1, 1, -2-2^(1/2)]])

matrix([[3-2^(1/2), 0, -5], [1/5, -1-2^(1/2), 0], [1, 1, -2-2^(1/2)]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[1, 0, 5/(-3+2^(1/2))], [0, 1, (-1+2^(1/2))/(-3+2^(1/2))], [0, 0, 0]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[a+5/(-3+2^(1/2))*c], [b+(-1+2^(1/2))/(-3+2^(1/2))*c], [0]]) = matrix([[0], [0], [0]])

v[2] = matrix([[5*2^(1/2)+5], [1], [1+2*2^(1/2)]])

evecsAval3 = -2^(1/2)

evecsAmul3 = 1

evecsAvec3 = vector([-5*2^(1/2)+5, 1, 1-2*2^(1/2)])

lambda = -2^(1/2)

A-I*lambda = matrix([[3+2^(1/2), 0, -5], [1/5, -1+2^(1/2), 0], [1, 1, -2+2^(1/2)]])

matrix([[3+2^(1/2), 0, -5], [1/5, -1+2^(1/2), 0], [1, 1, -2+2^(1/2)]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[1, 0, -5/(3+2^(1/2))], [0, 1, (1+2^(1/2))/(3+2^(1/2))], [0, 0, 0]])*matrix([[a], [b], [c]]) = matrix([[0], [0], [0]])

matrix([[a-5/(3+2^(1/2))*c], [b+(1+2^(1/2))/(3+2^(1/2))*c], [0]]) = matrix([[0], [0], [0]])

v[3] = matrix([[-5*2^(1/2)+5], [1], [1-2*2^(1/2)]])

`Use (column) eigenvectors of `*A*` to construct P:`

v[1] = matrix([[5], [1], [3]]), v[2] = matrix([[5*2^(1/2)+5], [1], [1+2*2^(1/2)]]), v[3] = matrix([[-5*2^(1/2)+5], [1], [1-2*2^(1/2)]])

P = matrix([[5, 5*2^(1/2)+5, -5*2^(1/2)+5], [1, 1, 1], [3, 1+2*2^(1/2), 1-2*2^(1/2)]]), P^`-1` = matrix([[-1/5, 1/2, 1/2], [1/20*(1+2^(1/2))*2^(1/2), 1/8*(-2+2^(1/2))*2^(1/2), -1/4], [1/20*(-1+2^(1/2))...

P^`-1`*P*` = `*matrix([[-1/5, 1/2, 1/2], [1/20*(1+2^(1/2))*2^(1/2), 1/8*(-2+2^(1/2))*2^(1/2), -1/4], [1/20*(-1+2^(1/2))*2^(1/2), 1/8*(2+2^(1/2))*2^(1/2), -1/4]])*matrix([[5, 5*2^(1/2)+5, -5*2^(1/2)+5],...

P^`-1`*A*P*` = `*matrix([[-1/5, 1/2, 1/2], [1/20*(1+2^(1/2))*2^(1/2), 1/8*(-2+2^(1/2))*2^(1/2), -1/4], [1/20*(-1+2^(1/2))*2^(1/2), 1/8*(2+2^(1/2))*2^(1/2), -1/4]])*matrix([[3, 0, -5], [1/5, -1, 0], [1,...