y = np.mean(samples[:,1]) * np.array([[0,-1j],[1j,0]])