code: def pathways(self, fraction=1.0, maxiter=1000):
          r
          return tptapi.pathways(self.net_flux, self.A, self.B,
                                 fraction=fraction, maxiter=maxiter)
docs: r"""Decompose flux network into dominant reaction paths.
  
          Parameters
          ----------
          fraction : float, optional
              Fraction of total flux to assemble in pathway decomposition
          maxiter : int, optional
              Maximum number of pathways for decomposition
  
          Returns
          -------
          paths : list
              List of dominant reaction pathways
          capacities: list
              List of capacities corresponding to each reactions pathway in paths
  
          References
          ----------
          .. [1] P. Metzner, C. Schuette and E. Vanden-Eijnden.
              Transition Path Theory for Markov Jump Processes.
              Multiscale Model Simul 7: 1192-1219 (2009)