1. close all;
2. clear all;
3.
4. i = 0;
5. max = 6;
6.
7. A = [1 2 3
8. 2 4 5
9. 3 5 -1];
10.
11. v = [0
12. 0
13. 1];
14.
15. lambda = v' * A * v;
16.
17. while i < max
18. i++;
19. w = (A - lambda * eye(rows(A))) \ v;
20. v = w / norm(w);
21. lambda = v' * A * v;
22. endwhile
23.
24. save Rayleigh.txt;