sigma(:,:,k) = cov(X(idx==k,:));