code: def capacity(F, path):
      r
      F = F.todok()
      L = len(path)
      currents = np.zeros(L - 1)
      for l in range(L - 1):
          i = path[l]
          j = path[l + 1]
          currents[l] = F[i, j]
  
      return currents.min()
docs: r"""Compute capacity (min. current) of path.
  
      Paramters
      ---------
      F : (M, M) scipy.sparse matrix
          The flux network (matrix of netflux values)
      path : list
          Reaction path
  
      Returns
      -------
      c : float
          Capacity (min. current of path)