A = magic(8)
A=A+A';
eig(A)
lam = poweriter(A);
plot(lam,'-o')
plot(520-lam,'-o')
shg
semilogy(abs(520-lam),'-o')
shg
[x,y]=ginput(2)
diff(log(y)) / diff(x)
exp(ans)
eig
eig(A)
(ans(end-1)/ans(end))^2
eig(A)
lam=inviter(A,160);
plot(lam,'-o')
shg
lam=inviter(A,160);
plot(lam,'-o')
shg
e=eig(A);
semilogy(abs(lam-e(7)),'-o')
shg
lam=inviter(A,164);
semilogy(abs(lam-e(7)),'-ok')
shg