import numpy as np

from scipy import linalg

A $=$ np$.$array$([[1,2],[3,4]])$

v$0=$ np$.$array$([[7],[-2]])$

v $=$ v$0$

d $=1$

tol$=1$e$-4$

i$=0$;

while d $>$ tol:

i $=$ i$+1$

v$0=$ v

z $=$ A@v

v $=$ z$/($abs$($z$[0])+$abs$($z$[1]))$

d $=$ abs$($v$[0]-$v$0[0])+$abs$($v$[1]-$v$0[1])$

print$("$iterations $=\%$d$"\%$ i$)$

print$("$distance $=\%$e$"\%$ d$)$

lam $=$ float$($v$.$T@A@v$)/$float$($v$.$T@v$)$

print$("$lambda $=\%\mathrm{.}8$f$"\%$ lam $)$

print$("$eigenvector $=",$v$)$

print$("$Spectrum of A $=",$ linalg.eigvals$($A$))$