A = matrix(CDF, [[1, I, 4],[-I, 3, 0],[4, 0, 2]])
Ac = A.conjugate() #complex conjugate
Act = Ac.transpose() #conjugate transpose
print A
print Act
# Act = A => A is Hermitian
X = matrix(CDF, [[1],[I],[-I]])
Xc = X.conjugate()
Xct = Xc.transpose()
Xct*A*X

SageDays/Kreyszig-7.4-24 (last edited 2010-08-11 18:38:56 by bhanukiran)