1. close all;
2. clear all;
3.
4. max = 18;
5. i = 0;
6.
7. A =[1 2 3
8. 2 4 5
9. 3 5 -1];
10.
11. while i < max
12. i++;
13. [Q, R] = qr(A);
14. A = R*Q;
15. endwhile
16.
17. B = diag(A);
18.
19. save qr.txt