Steady-state frequencies of each [k-1]-mer

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