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

SageDays/Kreyszig-7.4-23 (last edited 2010-08-11 18:37:33 by bhanukiran)