def METHOD_NAME(self):
    output = num.zeros((self.n + 1,)).astype(np.float64)
    msg = "Output shape mismatch"
    with pytest.raises(ValueError, match=msg):
        num.linalg.solve(self.a, self.b, out=output)