A = matrix([[0,9,-12],[-9,0,20],[12,-20,0]])

print A


A.characteristic_polynomial()

A.eigenvalues()

# one of the eigenvalues is real, two are complex conjugate pairs and their absolute = 1 => orthogonal matrix

SageDays/Kreyszig-7.3-5 (last edited 2010-08-11 12:58:06 by bhanukiran)