ix = np.r_[1:(n + 1) * (n + 2):2]