A=matrix([[1,2,3],[0,4,5],[0,0,6]])
print A
[1 2 3]
[0 4 5]
[0 0 6]
A.eigenvalues()
[6, 4, 1]
A.jordan_form()
[6|0|0]
[-+-+-]
[0|4|0]
[-+-+-]
[0|0|1]
B=matrix([[1,1],[-1,-1]])
print B
[ 1  1]
[-1 -1]
B.eigenvalues()
[0, 0]
B.jordan_form()
[0 1]
[0 0]

Kreyszig-Appendix B-4-U (last edited 2010-12-18 07:56:55 by 10)