x = np.real(np.trace(np.dot(rho,np.array([[0,1],[1,0]]))))