diff options
Diffstat (limited to '75/CH9/EX9.2/ex_2.sce')
-rwxr-xr-x | 75/CH9/EX9.2/ex_2.sce | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/75/CH9/EX9.2/ex_2.sce b/75/CH9/EX9.2/ex_2.sce new file mode 100755 index 000000000..2d28eadd5 --- /dev/null +++ b/75/CH9/EX9.2/ex_2.sce @@ -0,0 +1,25 @@ + // PG 591 + +n = 4 +A = [4 1 0 0;1 4 1 0;0 1 4 1;0 0 1 4] +lam = spec(A) + +// Since A is symmetric, all eigen values are real. +// The radii are all 1 or 2. +// The centers of all the circles are 4. +// All eigen values must all lie in the interval [2,6] +// Since the eigen values of inv(A) are the reciprocals of those of A, +// 1/6 <= mu <= 1/2 + +// Let inv(A) = B + +B=inv(A); +norm(B,2) +n +i = 1:n; +j = 1:n; + + // for j~i + // r = sum(abs(B(i,j))) + +// norm(B,2) = r(B) <= o.5
\ No newline at end of file |