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