a,k = var('a,k')

A = matrix([[a,k,k],[k,a,k],[k,k,a]])

print A

A.characteristic_polynomial()

A.eigenvalues()

# all the eigenvalues are real => the matrix is symmetric, if a, k are real.

SageDays/Kreyszig-7.3-6 (last edited 2010-08-11 13:01:30 by bhanukiran)