CH8 exercises problem10

A=matrix(2,2,[1,1,1,1]) print A A.eigenvalues() A.eigenvectors_right() B=matrix([[1, 1], [1, -1]]) P=matrix([1/B.row(j).norm()*B.row(j) for j in range(0,2)]) print B 

[1 1]
[1 1]
[ 1  1]
[ 1 -1]

[1 1]
[1 1]
[ 1  1]
[ 1 -1]