code: def steady_state(self):
          
          v, w = linalg.eigs(self.P.transpose(), which='LR')
          ssf = np.real(w[:, v.argmax()])
          ssf /= ssf.sum()
          return ssf
docs: Steady-state frequencies of each [k-1]-mer