Given the following matrix $$A = \begin{bmatrix} 8 & 1 & 0\\ 1 & 4 & \epsilon\\ 0 & \epsilon & 1\\ \end{bmatrix}, |\epsilon| < 1.$$ Gerschgorin's theorem states that each of the $\lambda_i$ eigenvalues of A will be placed inside a circular disk with center $a_{i,i}$ and radius $\sum^n_{j = 1, j\neq i} |a_{i,j}|$. As a direct consequence, the lowest eigenvalue $\lambda_3$ of A is such that $|\lambda_3 -1| < \epsilon$.
However I'm now asked to prove the tighter bound $|\lambda_3-1| < \epsilon^2$ using diagonal similarity transformations. My initial idea was to use the first step of the QR algorithm for Hessemberg diagonalization as usual to obtain a similar matrix $A_1 = Q^*AQ$ (with $Q$ computed by Householder transformations) preserving the eigenvalues and hopefully the tighter bound appears using Gerschgorin's theorem once again. I've made some numerical tests using Octave for some particular $\epsilon$ and it seems to be indeed the case that $A_1$ satisfies the bound. The problem is that the exact computation of $Q$ and $A_1$ have been absolutely painful, specially when dealing with the unknown symbol $\epsilon$. I also tried to use a shift of $\rho = 1$ with no better results. So now I'm suspicious there is a simpler similarity transformation that could lead to the desired result. Any suggestions would be highly appreciated.
Subtract off the identity matrix so the upper left term is now $7.$ The characteristic polynomial is now $$ x^3 - 10 x^2 + (20 -\epsilon^2)x + 7 \epsilon^2 $$
We wish to show that this has a root near zero, between $- \epsilon^2$ and $\epsilon^2.$ I guess we take $\epsilon \neq 0$ and treat that as a separate case.
The value of the characteristic polynomial when $x = \epsilon^2$ is $$ \epsilon^2 \left( 27 - 11 \epsilon^2 + \epsilon^4 \right) $$ which is positive, as $\epsilon^2 < 1.$
The value of the characteristic polynomial when $x = - \epsilon^2$ is $$ -\epsilon^2 \left( 13 +9 \epsilon^2 + \epsilon^4 \right) $$ which is negative, as $\epsilon^2 < 1.$
Te characteristic polynomial is continuous in $x,$ thus it has a root with $$ - \epsilon^2 < x < \epsilon^2 $$
The original matrix has an eigenvalue $x$ with $$ 1- \epsilon^2 < x < 1 +\epsilon^2 $$
ADDED: As Robert points out, we can simply evaluate the shifted characteristic polynomial $ x^3 - 10 x^2 + (20 -\epsilon^2)x + 7 \epsilon^2 $ at $x = t \epsilon^2,$ with $t$ real. When $t=0$ we get (positive) $7 \epsilon^2.$ There is some cancellation available if we take $t = -\frac{7}{20},$ and the polynomial comes out as $$ \epsilon^2 \left( - \frac{350}{400} \epsilon^2 - \frac{343}{8000} \epsilon^4\right) $$ which is negative. Thus the shifted eigenvalue is between $-\frac{7}{20} \epsilon^2$ and $0,$ the original eigenvalue between $1-\frac{7}{20} \epsilon^2$ and $1.$