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