From 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Tue, 10 Oct 2017 12:27:19 +0530
Subject: initial commit / add all books

---
 1499/CH7/EX7.19/q19.sce | 43 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)
 create mode 100644 1499/CH7/EX7.19/q19.sce

(limited to '1499/CH7/EX7.19/q19.sce')

diff --git a/1499/CH7/EX7.19/q19.sce b/1499/CH7/EX7.19/q19.sce
new file mode 100644
index 000000000..75985e2a7
--- /dev/null
+++ b/1499/CH7/EX7.19/q19.sce
@@ -0,0 +1,43 @@
+syms m11 m12 m13 m21 m22 m23 m31 m32 m33 ^ t m
+s=%s;
+poly(0,"l");
+A=[0 1 0;3 0 2;-12 -7 -6]
+B=[1;0;2]
+C=[1 0 0]
+[r c]=size(A)
+I=eye(r,c);
+p=l*I-A;
+q=det(p); // determinant of li-p
+// roots of q are
+l1=-1;
+l2=-2;
+l3=-3;
+x1=[m11;m21;m31];
+q1=(l1*I-A)*1
+
+//on solving
+m11=1;
+m21=-1;
+m31=-1;
+x2=[m12;m22;m32];
+q2=(l2*I-A)*1
+//on solving;
+m12=2;
+m22=-4;
+m32=1;
+
+x3=[m13;m23;m33];
+q3=(l3*I-A)*1
+//on solving
+m13=1;
+m23=-3;
+m33=3;
+
+// modal matrix is
+M=[m11 m12 m13;m21 m22 m23;m31 m32 m33]
+
+
+q=inv(M)
+A1=real(q*A*M)
+B1=(q*B)
+C1=C*M
-- 
cgit