docs: def remove_constraint(self, name):
          """Remove a constraint from the problem"""
          index = self._get_constraint_index(name)
          # Remove from matrix
          self._A = np.delete(self.A, index, 0)
          # Remove from upper_bounds
          self.upper_bounds = np.delete(self.upper_bounds, index)
          # Remove from constraint list
          del self._constraints[name]
          self._update_constraint_indices()
          self._reset_solution()