def elim(row, minus_row, scale, n=3): m = np.eye(n) m[row, minus_row] = -scale return m