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