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

SageDays/Kreyszig-7.4-25 (last edited 2010-08-11 18:42:01 by bhanukiran)