Ch-10-Prob-8-New-박종우

1840 days ago by matrix

A=matrix(5,5,[1,3,-6,4,3,0,1,2,8,2,0,0,2,0,-8,0,0,0,1,1,0,0,0,0,4]) print A.eigenvectors_right() 
       
[(4, [
(1, -2/15, -12/25, 1/25, 3/25)
], 1), (2, [
(0, 1, 1/2, 0, 0)
], 1), (1, [
(1, 0, 0, 0, 0)
], 3)]
[(4, [
(1, -2/15, -12/25, 1/25, 3/25)
], 1), (2, [
(0, 1, 1/2, 0, 0)
], 1), (1, [
(1, 0, 0, 0, 0)
], 3)]
M=matrix(5,5,[1,1,1,0,1,0,1/3,5/18,1,-2/15,0,0,0,1/2,-12/25,0,0,1/24,0,1/25,0,0,0,0,3/25]) J=matrix(5,5,[1,1,0,0,0,0,1,1,0,0,0,0,1,0,0,0,0,0,2,0,0,0,0,0,4]) print A*M==M*J 
       
True
True