: $\, {{\boldsymbol x}}_{k+1} \, = \, A{{\boldsymbol x}}_k, \, k \ge 0$. In most of practical cases it converges to a leading eigenvector ${{\boldsymbol v}}$ linearly with the rate $O \bigl( \bigl| \frac{\lambda_2}{\lambda_1}\bigr|^k\bigr)$, where $\lambda_1, \lambda_2$ are the first and the second largest by modulus eigenvalues of $A$ or $O\bigl( \frac{1}{k} \bigr)$, if the leading eigenvalue has non-trivial (i.e., of size bigger than one) Jordan blocks. In a rare case when there are several largest by modulus eigenvalues, then the “averaged” sequence $\frac1r\, \sum_{j=0}^{r-1} x_{k+j}$ converges to ${{\boldsymbol v}}$, where $r$ is the [*imprimitivity index*]{} (see [@G chapter 13]). In all cases we will say that the power method converges and, for the sake of simplicity, assume that there is only one largest by modulus eigenvalue, maybe multiple, i.e., that $r=1$.

\[th.40\] The selected leading eigenvector of a non-negative matrix $A$ is proportional to the limit of the power method applied to the matrix $A$ with the initial vector ${{\boldsymbol x}}_0 = {{\boldsymbol e}}$.

[Proof.]{} It may be assumed that $\rho(A) = 1$. The spectral radius of the matrix $A_{\varepsilon}$ strictly increases in $\varepsilon$, so $\rho(A_{\varepsilon}) =  1 + \delta$ for some $\delta > 0$. We have $$\bigl( A\, + \, \varepsilon \, {{\boldsymbol e}}\, {{\boldsymbol e}}^T \bigr)\, {{\boldsymbol v}}_{\varepsilon} \ = \ 
(1+\delta){{\boldsymbol v}}_{\varepsilon}\,.$$ Denoting by $S_{\varepsilon} = ({{\boldsymbol e}}, {{\boldsymbol v}}_{\varepsilon})$ the sum of entries of the vector ${{\boldsymbol v}}_{\varepsilon}$, we obtain $$A{{\boldsymbol v}}_{\varepsilon}\ + \ \varepsilon \, S_{\varepsilon}\, {{{\boldsymbol e}}}   \ = \ 
(1+\delta)\, {{\boldsymbol v}}_{\varepsilon}\,,$$ and hence $$\left( I \, - \, \frac{1}{(1+\delta)}\, A   \right)\,  {{\boldsymbol v}}_{\
: $ \, { { \boldsymbol x}}_{k+1 } \, = \, A{{\boldsymbol x}}_k, \, k \ge 0$. In most of practical cases it converges to a leading eigenvector   $ { { \boldsymbol v}}$ linearly with the pace $ oxygen \bigl (\bigl| \frac{\lambda_2}{\lambda_1}\bigr|^k\bigr)$, where $ \lambda_1, \lambda_2 $ are the beginning and the second largest by modulus eigenvalue of   $ A$ or $ O\bigl (\frac{1}{k } \bigr)$, if the leading eigenvalue have non - trivial (i.e., of size big than one) Jordan blocking. In a rare case when there embody several largest by modulus eigenvalue, then the “ averaged ” sequence $ \frac1r\, \sum_{j=0}^{r-1 } x_{k+j}$ converges to $ { { \boldsymbol v}}$, where $ r$ is the [ * imprimitivity index * ] { } (see   [ @G chapter 13 ]). In all cases we will read that the power method acting converges and, for the sake of simplicity, wear that there embody only one largest by modulus eigenvalue, maybe multiple, i.e., that $ r=1$. 

 \[th.40\ ] The choose leading eigenvector of a non - negative matrix $ A$ is proportional to the terminus ad quem of the power method acting applied to the matrix   $ A$ with the initial vector   $ { { \boldsymbol x}}_0 = { { \boldsymbol e}}$. 

 [ Proof. ] { } It may be assumed that $ \rho(A) = 1$. The spectral radius of the matrix $ A_{\varepsilon}$ strictly increases in $ \varepsilon$, therefore $ \rho(A_{\varepsilon }) =   1 + \delta$ for some $ \delta > 0$. We have $ $ \bigl (A\, + \, \varepsilon \, { { \boldsymbol e}}\, { { \boldsymbol e}}^T \bigr)\, { { \boldsymbol v}}_{\varepsilon } \ = \ 
 (1+\delta){{\boldsymbol v}}_{\varepsilon}\,.$$ Denoting by $ S_{\varepsilon } = ({ { \boldsymbol e } }, { { \boldsymbol v}}_{\varepsilon})$ the sum of entries of the vector   $ { { \boldsymbol v}}_{\varepsilon}$, we obtain $ $ A{{\boldsymbol v}}_{\varepsilon}\ + \ \varepsilon \, S_{\varepsilon}\, { { { \boldsymbol e } } }    \ = \ 
 (1+\delta)\, { { \boldsymbol v}}_{\varepsilon}\,,$$ and hence $ $ \left (I \, - \, \frac{1}{(1+\delta)}\, A    \right)\,   { { \boldsymbol v}}_{\
: $\, {{\blldsymbol x}}_{k+1} \, = \, A{{\boldsymnol x}}_k, \, k \ge 0$. In most oh practjcal casds it converges to a leading eugenvtbtor ${{\boldsymbol v}}$ linexrly with the ratw $O \uigl( \bigl| \frac{\lambda_2}{\lambda_1}\bigr|^i\nigr)$, chxre $\lambda_1, \lambca_2$ are the first and the sdclnd largest by modulus eigenvalues jf $A$ or $P\blgl( \frac{1}{k} \bigr)$, if uhe leaspnn eigenvalue has non-trivial (i.e., or size uigger than one) Jordan blocks. In a rare cwse ahen there are sevfral largesj by nodulus eigevvalues, thtn the “averafed” sequence $\frac1r\, \sum_{j=0}^{r-1} x_{k+j}$ covvergzs to ${{\boldsimgop v}}$, where $r$ ms the [*imprimitivibj index*]{} (see [@G cnapter 13]). In all caves we will say that the power method convgrges and, xox the sake of simplicutt, assome tvat gyerd ia pnmy one lacgest by mosulus eigencalue, maybe multiplt, i.q., that $r=1$.

\[th.40\] The seleceeq leading eigenvector of a non-negative katdix $A$ is proportional ti the limit of the poaer methoq applied to the matrix $A$ with the initial vector ${{\bmldsyjcol x}}_0 = {{\vopdsymbol e}}$.

[Proof.]{} It may be assumed that $\rho(A) = 1$. Fht siectral radius of the matrix $A_{\gatgpsilon}$ strictuy incxsaaes in $\varepsilon$, do $\rho(A_{\darepwilon}) =  1 + \delya$ for some $\delta > 0$. We have $$\bigl( A\, + \, \vagepsulon \, {{\boldsymbol e}}\, {{\boldsymbol e}}^T \nigr)\, {{\noldsymbol v}}_{\varepsilon} \ = \ 
(1+\dslta){{\boldsymhol v}}_{\varelrilon}\,.$$ Denoting bh $S_{\nare[silon} = ({{\boldsymbol e}}, {{\boldsrmbol v}}_{\vacepsikon})$ the sum of eneries of tje vegdor ${{\boldsymbol v}}_{\varfpsilln}$, we obtain $$A{{\boldsymbol v}}_{\varepsilon}\ + \ \vare'xilon \, S_{\vareksinon}\, {{{\boldsymyol e}}}   \ = \ 
(1+\delta)\, {{\bojdsymbol v}}_{\vareksilon}\,,$$ anb hencd $$\left( I \, - \, \frac{1}{(1+\dxlta)}\, A   \rigrt)\,  {{\boldsymbon v}}_{\
: $\, {{\boldsymbol x}}_{k+1} \, = \, \, \ge 0$. most of practical leading ${{\boldsymbol v}}$ linearly the rate $O \bigl| \frac{\lambda_2}{\lambda_1}\bigr|^k\bigr)$, where $\lambda_1, \lambda_2$ are first and the second largest by modulus eigenvalues of $A$ or $O\bigl( \frac{1}{k} if the leading eigenvalue has non-trivial (i.e., of size bigger than one) Jordan In rare when are several largest by modulus eigenvalues, then the “averaged” sequence $\frac1r\, \sum_{j=0}^{r-1} x_{k+j}$ converges to ${{\boldsymbol where $r$ is the [*imprimitivity index*]{} (see [@G 13]). In all cases will say that the power converges for the of assume there is only largest by modulus eigenvalue, maybe multiple, i.e., that $r=1$. \[th.40\] The selected leading eigenvector of a non-negative $A$ is the limit the method to the matrix the initial vector ${{\boldsymbol x}}_0 = It may be assumed that $\rho(A) = 1$. spectral radius the matrix $A_{\varepsilon}$ strictly increases in so $\rho(A_{\varepsilon}) = 1 + \delta$ for some > 0$. We have $$\bigl( A\, + \, \varepsilon \, {{\boldsymbol e}}\, {{\boldsymbol e}}^T \bigr)\, \ = \ (1+\delta){{\boldsymbol Denoting by $S_{\varepsilon} ({{\boldsymbol {{\boldsymbol the of entries the vector ${{\boldsymbol v}}_{\varepsilon}$, we obtain $$A{{\boldsymbol v}}_{\varepsilon}\ + \ \varepsilon S_{\varepsilon}\, {{{\boldsymbol e}}} \ = \ (1+\delta)\, {{\boldsymbol v}}_{\varepsilon}\,,$$ and I - \, \frac{1}{(1+\delta)}\, \right)\, {{\boldsymbol v}}_{\
: $\, {{\boldsymbol x}}_{k+1} \, = \, A{{\boldsymbol x}}_k, \, K \ge 0$. In most oF pracTicAl cAsEs it ConvErges to a leadinG EigeNvector ${{\boldsymbol v}}$ lineArly wItH The rATe $o \bigl( \Bigl| \fraC{\LaMBDa_2}{\lAmBdA_1}\biGr|^K\BiGr)$, wheRe $\lAmbda_1, \laMbda_2$ are the FirSt And the second LArGest by moduLus Eigenvalues oF $A$ oR $O\bigl( \FrAc{1}{k} \BIgr)$, if The LeadiNg eigeNValue hAs non-trivIaL (I.e., of siZE bigger THAn One) JOrdan blocks. In a rarE CaSE when there are sEveral LaRGeST By mOduLus eigenvaLuEs, theN The “averAGeD” SEQueNCe $\frac1r\, \sum_{j=0}^{r-1} x_{K+j}$ converges TO ${{\boLdsymbOl V}}$, whERe $r$ is tHe [*impRiMItiVity index*]{} (seE [@G chApter 13]). In alL cases WE will saY That the Power mEthOd cOnveRGeS aNd, fOr THe sAKe Of sIMplIcity, assUmE tHat thEre iS ONLY one LarGest By modUlus eigenvaluE, maYbe mULtiPle, i.e., That $r=1$.

\[Th.40\] ThE sElectEd leadIng eiGeNvector of a non-neGatiVe matrix $A$ Is pRoPorTiOnal tO The limIt oF thE power mEthod apPLieD tO THE mAtrix $A$ with the initiAl VECtOr ${{\boldsyMbol x}}_0 = {{\bOLdSyMBol e}}$.

[ProoF.]{} IT maY be aSSUmed tHat $\rHO(A) = 1$. the spectRal radIUs Of The matrIx $a_{\varepSiLon}$ StrIctly INcreAses in $\VarepsilOn$, so $\rHO(A_{\varepsilon}) =  1 + \deLTa$ for some $\deltA > 0$. we HAVe $$\BIgl( A\, + \, \VarEpsilon \, {{\boldSymbOL e}}\, {{\boLdsyMBoL e}}^T \BIgr)\, {{\boLdsymBoL V}}_{\vARepsilon} \ = \ 
(1+\delta){{\boldsyMbOl v}}_{\varEpsilOn}\,.$$ Denoting by $S_{\Varepsilon} = ({{\BOLDsymbol e}}, {{\BoldSYmBOl v}}_{\varepsilon})$ tHe sum Of entries oF The vectoR ${{\boldSymbol v}}_{\vArepsilon}$, WE Obtain $$A{{\bOldSymBol V}}_{\vaREPsIlon}\ + \ \varepsiloN \, s_{\VarePsIlon}\, {{{\bolDsyMbol e}}}   \ = \ 
(1+\deLta)\, {{\BolDsyMboL v}}_{\VarepsiloN}\,,$$ and hencE $$\lEfT( I \, - \, \FrAc{1}{(1+\dElta)}\, A   \RIght)\,  {{\boldSyMboL v}}_{\
: $\, {{\boldsymbol x}}_{k +1} \, = \ , A{{ \bo lds ym bolx}}_ k, \, k \ge 0$ .  Inmost of practical case s itco n verg e sto aleading ei g e nve ct or  ${ {\ b ol dsymb olv}}$ li nearly wit h t he  rate $O \bi g l(  \bigl| \f rac {\lambda_2}{ \la mbda_1 }\ big r |^k\b igr )$, w here $ \ lambda _1, \lamb da _ 2$ are the fir s t  a nd t he second largest by modulus eigenv aluesof $A $ or$O\ bigl( \fra c{ 1}{k} \bigr)$ ,  i f t hel eading eigenv alue has no n -tr ivial(i .e. ,  of si ze bi gg e r t han one) Jo rdan  blocks.In a r a re case when th ere ar e s eve rall ar ge stby mod u lu s e i gen values,th en  the“ave r a g e d” s equ ence  $\fr ac1r\, \sum_{ j=0 }^{r - 1}x_{k+ j}$ c onve rg es to  ${{\b oldsy mb ol v}}$, where$r$is the [* imp ri mit iv ity i n dex*]{ } ( see  [@G ch apter 1 3 ]).  I n a ll  cases we will say  t h a tthe powe r meth o dco n verges a nd , f or t h e  sake  ofs im plicity,  assum e  t ha t there  i s only  o nelar gestb y mo duluseigenval ue, m a ybe multiple,i .e., that $r= 1 $. \[ t h.40 \]The selecte d le a ding  eig e nv ect o r ofa non -n e ga t ive matrix $A$ is p ro portio nal t o the limit o f the powe r m ethod ap plie d  t o  the matrix $A $ wit h the init i al vecto r ${{ \boldsym bol x}}_0 =  {{\bold sym bol  e} }$. [P roof.]{} It m a y  beas sumed t hat  $\rho( A)= 1 $.The  s pectral r adius of  t he  m at rix  $A_{ \ varepsil on }$st ric tly i n crease s in$\va re ps i lon $, so $ \ rh o ( A_{\ va re psil on} )=  1+ \d e lta $ for s ome $\del ta>  0$.  W ehave $$ \bigl( A\, +\,  \varepsil on  \,  {{\bo l d symbol e }}\, {{\boldsymbol e}}^ T  \bigr) \,{{\bo ldsy mbol v}}_ {\v arepsi lon }  \ = \  
(1+\ delta ){ {\b o l dsymb o l  v }}_ {\ varepsilon } \ ,.$ $ Den ot ingby $S_{ \varepsilon} = ({{ \ bol dsymbol e}},{{\ bold s y mb olv }} _ {\v ar e psi l o n})$ the sum of  entries o ft he  vector ${ { \bo ld symbolv}}_{\v areps i lon}$,we obtain  $$A{{\bo ld symb o l  v} }_{\vareps ilon}\ +  \ \varep s ilon\ ,S_{\v are psilon }\ , { {{\bo ldsymb o l e }}}  \ = \  (1+\de lta)\ ,{{\bolds ymbol v}}_{\varepsilon} \,,$$and h enc e $$\left ( I \,- \, \fra c{1} {(1+\delta )}\ , A    \r igh t )\, {{\b o ld sym b ol v} }_{\ 
: $\,_{{\boldsymbol x}}_{k+1}_\, = \, A{{\boldsymbol_x}}_k, \,_k_\ge 0$._In_most of practical_cases it converges_to a leading eigenvector ${{\boldsymbol_v}}$ linearly with_the_rate $O \bigl( \bigl| \frac{\lambda_2}{\lambda_1}\bigr|^k\bigr)$, where $\lambda_1, \lambda_2$ are the first and the second_largest_by modulus_eigenvalues_of $A$_or $O\bigl( \frac{1}{k} \bigr)$, if_the leading eigenvalue has non-trivial_(i.e., of_size bigger than one) Jordan blocks. In a_rare_case when there_are several largest by modulus eigenvalues, then the “averaged”_sequence $\frac1r\, \sum_{j=0}^{r-1} x_{k+j}$ converges to_${{\boldsymbol v}}$, where_$r$_is_the [*imprimitivity index*]{} (see [@G_chapter 13]). In all cases we_will say that the power method_converges and, for the sake of simplicity,_assume that there is only one_largest by modulus eigenvalue, maybe_multiple, i.e.,_that $r=1$.

\[th.40\] The selected leading_eigenvector of a_non-negative matrix_$A$ is proportional_to the limit of the power_method applied to_the matrix $A$ with the initial vector ${{\boldsymbol_x}}_0_= {{\boldsymbol e}}$.

[Proof.]{}_It_may_be assumed_that $\rho(A) =_1$._The spectral_radius_of the matrix $A_{\varepsilon}$ strictly increases_in_$\varepsilon$, so $\rho(A_{\varepsilon}) =  1 +_\delta$ for some $\delta_>_0$. We have $$\bigl(_A\, + \, \varepsilon \,_{{\boldsymbol e}}\, {{\boldsymbol e}}^T \bigr)\, {{\boldsymbol_v}}_{\varepsilon} \_= \_
(1+\delta){{\boldsymbol v}}_{\varepsilon}\,.$$ Denoting by $S_{\varepsilon} = ({{\boldsymbol e}}, {{\boldsymbol v}}_{\varepsilon})$ the_sum of entries of the vector ${{\boldsymbol_v}}_{\varepsilon}$, we obtain $$A{{\boldsymbol_v}}_{\varepsilon}\ +_\_\varepsilon \, S_{\varepsilon}\,_{{{\boldsymbol_e}}} _ \ = \ 
(1+\delta)\, {{\boldsymbol v}}_{\varepsilon}\,,$$_and hence_$$\left( I \, - \, \frac{1}{(1+\delta)}\,_A   \right)\,__{{\boldsymbol v}}_{\