1 files changed, 26 insertions, 0 deletions
diff --git a/29/CH12/EX12.56/exa12_56.sce b/29/CH12/EX12.56/exa12_56.sce
new file mode 100755
index 000000000..e7c9302a6
--- /dev/null
+++ b/29/CH12/EX12.56/exa12_56.sce
@@ -0,0 +1,26 @@
+//caption:obtain_state_transistion_matrix
+//example 12_56
+//page 586
+s=%s;
+syms t
+A=[0 1;0 -3]
+[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
+    if(i==2 & j==1)
+      else
+    ip(i,j)=ilaplace(r(i,j),s,t);
+    j=j+1;
+    end
+  end
+  i=i+1;
+end
+r(2,1)=0
+disp(ip,"state transistion matrix,ip(t)=");
+\ No newline at end of file