1 files changed, 27 insertions, 0 deletions
diff --git a/29/CH9/EX9.10.11/exa9_10_11.sce b/29/CH9/EX9.10.11/exa9_10_11.sce
new file mode 100755
index 000000000..c99f11012
--- /dev/null
+++ b/29/CH9/EX9.10.11/exa9_10_11.sce
@@ -0,0 +1,27 @@
+//caption:obtain_time_response
+//example 9.10.11
+//page 394
+s=%s;
+syms t
+A=[0 1;-2 0]
+B=[1 -1]
+x0=[1 1]'
+[r c]=size(A);//size of matrix A
+p=s*eye(r,c)-A;//s*I-A where I is identity matrix
+q=det(p)//determinant of sI-A
+r=inv(p)//inverse of sI-A
+//for calculating state transistion matrix
+ip=[0 0;0 0]
+i=1;j=1;
+for i=1:2
+  for j=1:2
+    ip(i,j)=ilaplace(r(i,j),s,t);
+    j=j+1;
+  end
+  i=i+1;
+end
+disp(ip,"state transistion matrix,ip(t)=");
+x=ip*x0
+y=x(1,1)-x(2,1)
+y=simple(y) //output
+disp(y,"time response of the system,y(t)=");