From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 29/CH9/EX9.10.15/exa9_10_15.sce | 32 ++++++++++++++++++++++++++++++++
 1 file changed, 32 insertions(+)
 create mode 100755 29/CH9/EX9.10.15/exa9_10_15.sce

(limited to '29/CH9/EX9.10.15/exa9_10_15.sce')

diff --git a/29/CH9/EX9.10.15/exa9_10_15.sce b/29/CH9/EX9.10.15/exa9_10_15.sce
new file mode 100755
index 000000000..34616873e
--- /dev/null
+++ b/29/CH9/EX9.10.15/exa9_10_15.sce
@@ -0,0 +1,32 @@
+//caption:obtain_time_response_using_diagonalization_process
+//example 9.10.15
+//page 403
+syms t;
+s=%s;
+poly(0,"l");
+A=[0 1;-6 -5]
+B=[0;1]
+C=[6 1]
+x0=[1 0]'
+[r c]=size(A);//size of matrix A
+I=eye(r,c);
+p=l*I-A;//l*I-A where I is identity matrix
+q=det(p)//determinant of lI-A
+//roots of q are -1 and -3
+l1=-2;
+l2=-3;
+//for determining vandermonde's matrix
+P=[1 1;l1 l2];
+P1=inv(P);
+k1=P1*A*P;
+z0=P1*x0;
+k2=P1*B
+Z=inv(s*I-k1)*z0+(inv(s*I-k1)*k2)*(1/s);
+X=P*Z
+X(2,1)=-5/(s^2+5*s+6)
+x=[0;0]
+x(1,1)=ilaplace(X(1,1),s,t);
+x(2,1)=ilaplace(X(2,1),s,t);
+y=C*x
+y=simple(y)
+disp(y,"output_equation_is,y(t)=")
\ No newline at end of file
-- 
cgit