14.1. THE POWER METHOD FOR EIGENVALUES 369

Now pick the eigenvalue λq which is closest to q. Then

|Ax− qx|2 =

n∑k=1

|ak|2 (λk − q)2 ≥ (λq − q)

2n∑

k=1

|ak|2 = (λq − q)2 |x|2

which implies 14.5. ■

Example 14.1.10 Consider the symmetric matrix A =

 1 2 3

2 2 1

3 1 4

 . Let x =(1, 1, 1)T.

How close is the Rayleigh quotient to some eigenvalue of A? Find the eigenvector and eigen-value to several decimal places.

Everything is real and so there is no need to worry about taking conjugates. Therefore,the Rayleigh quotient is

(1 1 1

) 1 2 3

2 2 1

3 1 4

 1

1

1

3

=19

3

According to the above theorem, there is some eigenvalue of this matrix λq such that

∣∣∣∣λq − 19

3

∣∣∣∣ ≤

∣∣∣∣∣∣∣ 1 2 3

2 2 1

3 1 4

 1

1

1

− 193

 1

1

1

∣∣∣∣∣∣∣

√3

=1√3

 − 13

− 43

53

=

√19 +

(43

)2+(53

)2√3

= 1. 247 2

Could you find this eigenvalue and associated eigenvector? Of course you could. This iswhat the shifted inverse power method is all about.

Solve  1 2 3

2 2 1

3 1 4

− 19

3

 1 0 0

0 1 0

0 0 1

 x

y

z

 =

 1

1

1

In other words solve  − 16

3 2 3

2 − 133 1

3 1 − 73

 x

y

z

 =

 1

1

1

and divide by the entry which is largest, 3. 870 7, to get

u2 =

 . 699 25

. 493 89

1.0



14.1. THE POWER METHOD FOR EIGENVALUES 369Now pick the eigenvalue \, which is closest to g. Then2 2 2 2 2 2) 12|Ax — gx|* = SO lax! Or — a = Og — a)” do lanl? = Ag — a |x|k=1 k=1which implies 14.5.1 2 3Example 14.1.10 Consider the symmetric matria A= | 2 2 1 |. Letx=(1,1, 1)" .3.1 4How close is the Rayleigh quotient to some eigenvalue of A? Find the eigenvector and eigen-value to several decimal places.Everything is real and so there is no need to worry about taking conjugates. Therefore,the Rayleigh quotient is123 1(1 11)/ 221 13°14 1 193 ~ 31 2 1 122 1 1 ae ,19 314 1 1 1 —3Ay-—| < = _432 2ot (3) +)= = 1.2472v3Could you find this eigenvalue and associated eigenvector? Of course you could. This iswhat the shifted inverse power method is all about.Solve1 2 8 1 0 0 1192 2 1 - 3 01 0 y = 13 14 0 1 z 1In other words solve~ 2 3 xv 113 _2 -z 1 y = 13 1 -# z 1and divide by the entry which is largest, 3.8707, to get. 699 25uz = | .493891.0