A=matrix(5, 5, [4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4])
print A.charpoly(). factor()
print A.eigenvalues()
E=identity_matrix(5)
print (A-0*E). rank()
print ((A-0*E)^2). rank()
print ((A-0*E)^3). rank()
print ((A-0*E)^4). rank()
|
|
(x - 20) * x^4
[20, 0, 0, 0, 0]
1
1
1
1
(x - 20) * x^4
[20, 0, 0, 0, 0]
1
1
1
1
|