solution = linalg.lstsq(a_matrix, rhs, rcond=-1)