1 files changed, 24 insertions, 0 deletions
diff --git a/257/CH13/EX13.19/Example_13_19.sce b/257/CH13/EX13.19/Example_13_19.sce
new file mode 100644
index 000000000..6785b98c2
--- /dev/null
+++ b/257/CH13/EX13.19/Example_13_19.sce
@@ -0,0 +1,24 @@
+A=[0 1; -2 -3]
+U=[1/s;1/(s+2)]
+B=[2 1; 0 1]
+phi=inv(s*eye(2,2)-A)
+X=[0;0]
+a1=ilaplace(phi(1,1),s,t)
+a2=ilaplace(phi(1,2),s,t)
+a3=ilaplace(phi(2,1),s,t)
+a4=ilaplace(phi(2,2),s,t)
+
+S=[a1 a2;a3 a4]
+disp(S,"%e^(A*t) = ")
+ZIR=S*X
+disp(ZIR," ZIR = ")
+
+k=phi*B*U
+b1=ilaplace(k(1,1),s,t)
+b3=ilaplace(k(2,1),s,t)
+
+ZSR=[b1;b3]
+disp(ZSR,"ZSR = ")
+
+X=ZIR+ZSR;
+