docs: def factorize_and_solve(self, A, b):
          """
          Factorizes A and solves Ax=b.
  
          Returns
          -------
          x : vector
          """
  
          self.factorize(A)
  
          return self.solve(b)