summaryrefslogtreecommitdiff
path: root/75/CH9/EX9.5/ex_5.sce
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /75/CH9/EX9.5/ex_5.sce
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '75/CH9/EX9.5/ex_5.sce')
-rwxr-xr-x75/CH9/EX9.5/ex_5.sce25
1 files changed, 25 insertions, 0 deletions
diff --git a/75/CH9/EX9.5/ex_5.sce b/75/CH9/EX9.5/ex_5.sce
new file mode 100755
index 000000000..6982ab9db
--- /dev/null
+++ b/75/CH9/EX9.5/ex_5.sce
@@ -0,0 +1,25 @@
+ // PG 599
+
+ // e = 0.001
+ // From earlier example :
+ // eigen values of matrix A are 1 and 2. So,..
+
+ // inv(P)*A*P = [1 0;0 2]
+
+A = [101 -90;110 -98]
+B = [-1 -1;0 0]
+ // From the above equation, we get:
+
+P = [9/sqrt(181) -10/sqrt(221);10/sqrt(181) -11/sqrt(221)]
+inv(P)
+K = norm(P)*norm(inv(P)) // K is condition number
+u1 = P(:,1)
+u2 = P(:,2)
+Q = inv(P)
+R = Q'
+w1 = R(:,1)
+w2 = R(:,2)
+s1 = 1/norm(w1,2)
+norm(B)
+
+// abs(lam1(e) - lam1) <= sqrt(2)*e/0.005 + O(e^2) = 283*e + O(e^2)