diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /75/CH9/EX9.5/ex_5.sce | |
download | Scilab-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-x | 75/CH9/EX9.5/ex_5.sce | 25 |
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) |