1 files changed, 32 insertions, 0 deletions
diff --git a/965/CH4/EX4.13/13.sci b/965/CH4/EX4.13/13.sci
new file mode 100644
index 000000000..71609efaa
--- /dev/null
+++ b/965/CH4/EX4.13/13.sci
@@ -0,0 +1,32 @@
+clc;
+clear all;
+disp("spherical thermocouple junction")
+R=8/2000;//m
+hg=40;//W/m^2.C
+ha=10;//W/m^2.C
+k=40;// W/m.C
+c=420;// J/kg.C
+rho=8000;// kg/m^3
+
+As=4*3.1416*R*R;//m^2
+V=4*3.1416*(R^3)/3;//m^3
+
+
+tau1=rho*V*c/(hg*As);
+disp("s",tau1,"tau*  =")
+
+ta=300;// degree C
+ti=40;//degree C
+tau2=10;//s
+X=-1*tau2/tau1;
+t=ta+(ti-ta)*exp(X);
+disp("degree C",t,"temperature during heating t =")
+
+tau3=20;//s
+tau4=rho*R*c/(3*ha);//s
+disp("s",tau4,"tau* =")
+ti=t;// degree C
+X=-1*tau3/tau4;
+ta=30;//degree C
+t=ta+(ti-ta)*exp(X);//degree C
+disp("degree C",t,"temperature during cooling t =")