A = Matrix([[1,4],[-4,1]])

print A

A.characteristic_polynomial()

A.eigenvalues()

# eigenvalues are imaginary => the matrix is not symmetric

# eigenvalues are not purely imaginary => the matrix is not skew-symmetric

# eigenvalues are complex conjugates with modulii = 1 => the matrix is not orthogonal

SageDays/Kreyszig-7.3-3 (last edited 2010-08-11 12:27:02 by bhanukiran)