260/CH5/EX5.8/5_8.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19

//Eg-5.8
//pg-231

clear
clc

A=[1 2 4;3 1 2;4 2 5];
invA=inv(A);
z=[1;1;1];
for i=1:25
    a=invA*z;
    b=(sum(a.^2))^.5;
      z=a/b;
end

la_d=b;
minval=1/b;
disp("minimum eigen value of A")
disp(minval)