code: def major_flux(self, fraction=0.9):
          r
          (paths, pathfluxes) = self.pathways(fraction=fraction)
          return self._pathways_to_flux(paths, pathfluxes, n=self.nstates)
docs: r"""Returns the main pathway part of the net flux comprising
          at most the requested fraction of the full flux.