x = np.real(np.dot(psi.conj(),np.dot(np.array([[0,1],[1,0]]),psi)))