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.11/ex_11.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.11/ex_11.sce')
-rwxr-xr-x | 75/CH9/EX9.11/ex_11.sce | 27 |
1 files changed, 27 insertions, 0 deletions
diff --git a/75/CH9/EX9.11/ex_11.sce b/75/CH9/EX9.11/ex_11.sce new file mode 100755 index 000000000..84d26d659 --- /dev/null +++ b/75/CH9/EX9.11/ex_11.sce @@ -0,0 +1,27 @@ + // PG (613) + +A = [4 1 1;1 4 1;1 1 4] +w1 = [0.985599 0.119573 0.119573]' +P1 = eye() - 2*w1*w1' +A2 = P1*A +w2 = [0 0.996393 0.0848572]' +P2 = eye() - 2*w2*w2' +R = P2*A2 +Q = P1*P2 +Q*R + +// A = Q * R + +abs(det(A)) +abs(det(Q)*det(R)) + +// |det(A)| = |det(Q)*det(R)| = |det(R)| = 54 (approx) + +lam = spec(A)' +lam1 = lam(1,1) +lam2 = lam(1,2) +lam3 = lam(1,3) +lam1 * lam2 * lam3 + +// Product of eigen values also comes out to be 54 + |