code: def all_timescales_(self):
          
  
          us, lvs, rvs = self._get_eigensystem()
          # make sure to leave off equilibrium distribution
          timescales = -1 / us[:,1:]
          return timescales
docs: Implied relaxation timescales each sample in the ensemble
  
          Returns
          -------
          timescales : array-like, shape = (n_samples, n_timescales,)
              The longest implied relaxation timescales of the each sample in
              the ensemble.