1 files changed, 29 insertions, 0 deletions
diff --git a/1967/CH13/EX13.9/13_9.sce b/1967/CH13/EX13.9/13_9.sce
new file mode 100755
index 000000000..3dea93238
--- /dev/null
+++ b/1967/CH13/EX13.9/13_9.sce
@@ -0,0 +1,29 @@
+clc
+//initialisation of variables
+clear
+a= -15.84 
+b= 22.84*10^-3
+c= -80.97*10^-7 
+T= 25 //C
+H1= -48.1 //kcal
+H2= -26.4
+dS= 53.09
+T1= 327 //C
+r1= 0.58
+r2= 1.1
+r3= 1.13
+//CALCULATIONS
+dH= (H1-H2)*1000-a*(273.2+T)-0.5*b*(273.2+T)^2-0.33*c*(273.2+T)^3
+dF= (H1-H2)*1000+(273.2+T)*dS
+I= (dF-dH+a*(273.2+T)*log(273.2+T)+0.5*b*(273.2+T)^2+0.166*c*(273.2+T)^3)/(273.2+T)
+dF1= (dH-a*(273+T1)*log(273+T1)-0.5*b*(273+T1)^2-0.166*c*(273+T1)^3)+I*(273+T1)
+Kf= 10^(-dF1/(4.576*(273+T1)))
+Jr= r1/(r2^2*r3)
+Kp= Kf/Jr
+//RESULTS
+printf ('heat of formation = %.f cal',dH)
+printf ('\n Entropy = %.f cal',dF)
+printf ('\n Inertia = %.f gm cm^2',I)
+printf ('\n Entropy = %.f cal',dF1)
+printf ('\n Kf = %.1e ',Kf)
+printf ('\n Kp = %.1e ',Kp)