diff options
Diffstat (limited to '260/CH5/EX5.7/5_7.sce')
-rw-r--r-- | 260/CH5/EX5.7/5_7.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/260/CH5/EX5.7/5_7.sce b/260/CH5/EX5.7/5_7.sce new file mode 100644 index 000000000..326e7e921 --- /dev/null +++ b/260/CH5/EX5.7/5_7.sce @@ -0,0 +1,27 @@ +//Eg-5.7
+//pg-230
+
+clear
+clc
+
+A=[0 1 1 1;1 0 1 1;1 1 0 1;1 1 1 0];
+ z=[1;0;0;0];
+ for i=1:20
+ printf('iteration number=%f\n',i);
+ a=A*z;
+ b=(sum(a.^2))^.5;
+ printf('dominant eigen value=%f\n',b);
+ z=a/b;
+ printf('z=%f\n',z);
+end
+
+disp("from iterations dominant eigen value converged to ")
+disp(b)
+disp("after 20 iterations")
+
+
+if round((A-b*eye(4,4))*z)==0 then
+ disp("dominant eigen value is positive")
+else
+ disp("dominant eigen value is negative")
+end
|