A = matrix(QQ,[[0.96, -0.28],[0.28,0.96]])

print A

A.determinant()


evals = A.eigenvalues()
print evals

evals[1].abs()

evals[0].abs()

# 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 orthogonal

SageDays/Kreyszig-7.3-1 (last edited 2010-08-11 12:22:56 by bhanukiran)