Di = D[i, np.concatenate((np.r_[0:i], np.r_[i + 1:n]))]