3 files changed, 44 insertions, 0 deletions
diff --git a/46/CH5/EX5.1/Example5_1.sce b/46/CH5/EX5.1/Example5_1.sce
new file mode 100755
index 000000000..5f5368ba0
--- /dev/null
+++ b/46/CH5/EX5.1/Example5_1.sce
@@ -0,0 +1,9 @@
+//Example 5.1
+tau=0.1;//min
+xs=90;//degrees
+A=10;//degrees
+Y_inf=10;//degrees
+Y_t=8;//degrees
+//Substituting into Eq.(5.12) the appropriate values of Y_t,A,and tau gives
+t=-0.1*logm(1-(Y_t/A));//min
+disp('min',t,'time=')
+\ No newline at end of file
diff --git a/46/CH5/EX5.2/Example5_2.sce b/46/CH5/EX5.2/Example5_2.sce
new file mode 100755
index 000000000..475584914
--- /dev/null
+++ b/46/CH5/EX5.2/Example5_2.sce
@@ -0,0 +1,35 @@
+//Example 5.2
+clear
+clc
+tau=0.1;//min
+xs=100;//Fahrenheit
+ys=100;//Fahrenheit
+A=2;//Fahrenheit
+f=10/%pi;//cycles/min
+w=2*%pi*f;//rad/min
+//From Eq.(5.25), the amplitude of the response and the phase angle are calculated;thus
+disp('Fahrenheit',A/sqrt((tau*w)^2+1),'A/sqrt((tau*w)^2+1)=')
+phi=atan(-w*tau);//radians
+phi=phi*180/%pi;//degrees
+disp('degrees',phi,'phase lag=')
+t=0:0.01:1;
+//From Eq.(5.19), the input of the thermometer is  therefore
+disp("X(t)=2*sin(20*t)");
+//or
+xt=xs+2*sin(20*t);
+//The response of the thermometer is therefore
+disp("Y(t)=0.8944*sin(20*t-63.4349)")
+//or
+yt=ys+0.8944*sin(20*t-63.4349);
+Lag=phi/(360*f);//min
+Lag=abs(Lag);//min
+disp('min',Lag,'Lag=')
+clf;
+plot(t,yt)
+plot(t,xt)
+xlabel('time')
+ylabel('x(t) , y(t)')
+title('x(t) , y(t) Vs time')
+xgrid
+
+
diff --git a/46/CH5/EX5.2/figure5_2.jpg b/46/CH5/EX5.2/figure5_2.jpg
new file mode 100755
index 000000000..a8b463c04
--- /dev/null
+++ b/46/CH5/EX5.2/figure5_2.jpg