Boundary condition for np.linalg.solve?

Question

I have the following problem:

I have a system of 9 linear coupled equations, with 9 variables I want to solve for. I wrote it in Matrixform, so I can solve it like