Quiz: How does `scipy.linalg.solve(A, b)` solve the linear system A * x = b compared to computing `np.linalg.inv(A) @ b`?
Never invert matrices directly: scipy.linalg.solve uses efficient LU/Cholesky decomposition , avoiding instability and O(n³) cost.