1 files changed, 26 insertions, 0 deletions
diff --git a/46/CH6/EX6.1/Example6_1.sce b/46/CH6/EX6.1/Example6_1.sce
new file mode 100755
index 000000000..cc04a958a
--- /dev/null
+++ b/46/CH6/EX6.1/Example6_1.sce
@@ -0,0 +1,26 @@
+//Example 6.1
+clc;
+syms s t;
+tau=1;//min
+R=1/9;//ft/cfm
+A=9;
+//from Equation 6.8
+g=R/(tau*s+1);
+disp(g,'H(s)/Q(s)=')
+//from Example 4.5 
+disp('Q(t)=90[u(t)-u(t-0.1)')
+//where u(t) is a unit step function, the laplace transform of it gives
+Qs=90*(1-exp(-0.1*s))/s
+disp(Qs,'Q(s)=')
+Hs=Qs*g;
+disp(Hs,'H(s)=')
+//taking first term for t<0.1, the second term goes equals to zero
+Ht=ilaplace('10*(1/(s*(s+1)))',s,t);//t<0.1
+disp(Ht,'H(t)=')
+disp('H(t)=10(1-yexp(-(-t-0.1)))')//t>0.1
+Ht=10*((1-exp(-t))-(1-exp(-(-t-0.1))));
+disp(Ht,'H(t)=')
+//from Eq.(5.16)
+Ht=R*A*exp(-(t/tau));//impulse 
+disp(Ht,'H(t)=')
+