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



X = matrix(CDF, [[-2*I],[1 + I]])
Xc = X.conjugate()
Xct = Xc.transpose()
Xct*A*X

SageDays/Kreyszig-7.4-21 (last edited 2010-08-11 18:35:10 by bhanukiran)