12 files changed, 191 insertions, 0 deletions
diff --git a/1967/CH13/EX13.1/13_1.sce b/1967/CH13/EX13.1/13_1.sce
new file mode 100755
index 000000000..dd37883ce
--- /dev/null
+++ b/1967/CH13/EX13.1/13_1.sce
@@ -0,0 +1,12 @@
+clc
+//initialisation of variables
+clear
+R= 4.576 //cal deg^-1 mole^-1
+T= 700 //C
+Kp= 0.71
+p1= 1.5 //atm
+p2= 5 //atm
+//CALCULATIONS
+dF= -R*(273+T)*(log(Kp)-log((p1*p2)/(10*p2)))*0.77
+//RESULTS
+printf ('Free energy = %.f cal',dF-10)
diff --git a/1967/CH13/EX13.10/13_10.sce b/1967/CH13/EX13.10/13_10.sce
new file mode 100755
index 000000000..3c123d514
--- /dev/null
+++ b/1967/CH13/EX13.10/13_10.sce
@@ -0,0 +1,18 @@
+clc
+//initialisation of variables
+clear
+F1= 24.423//cal deg^-1
+F2= 21.031 //cal deg^-1
+F3= 37.172 //cal deg^-1
+H1= 2.024 //kcal
+H2= 1.035 //kcal
+H3= 2.365 //kcal
+H= -57.8 //kcal
+T= 25 //C
+//CALCULATIONS
+dF= F3-F1-F2
+dH= H3-H1-H2
+Hf= H-dH
+F= Hf-((273.2+T)*dF*10^-3) 
+//RESULTS
+printf ('Standard free energy = %.2f kcal',F)
diff --git a/1967/CH13/EX13.11/13_11.sce b/1967/CH13/EX13.11/13_11.sce
new file mode 100755
index 000000000..5da1c40e9
--- /dev/null
+++ b/1967/CH13/EX13.11/13_11.sce
@@ -0,0 +1,15 @@
+clc
+//initialisation of variables
+clear
+T= 1000 //C
+j= 1.5
+Q= 35840 //cal
+I= 743*10^-40 //g cm^2
+w= 214 //cm^-2
+Kf= 0.184
+//RESULTS
+printf ('Equilibrium constant = %.3f ',Kf)
+
+//NO SOLUTION IS GIVEN TO SOLVE Kf
+//INCOMPLETE SOLUTION IN THE TEXTBOOK
+
diff --git a/1967/CH13/EX13.12/13_12.sce b/1967/CH13/EX13.12/13_12.sce
new file mode 100755
index 000000000..941a9282a
--- /dev/null
+++ b/1967/CH13/EX13.12/13_12.sce
@@ -0,0 +1,18 @@
+clc
+//initialisation of variables
+clear
+dH= 83 //cal
+R= 1.98 //cal mole K^-1
+T= 25 //C
+M1= 128 //gms
+M2= 4 //gms
+M3= 2 //gms
+M4= 129 //gms
+I1= 4.31 //g cm^2
+I2= 0.920 //g cm^2
+I3= 0.459 //g cm^2
+I4= 8.555 // g cm^2
+//CALCULATIONS
+K= 10^((-dH/(2.303*R*(273.2+T)))+1.5*log10(M1^2*M2/(M3*M4^2))+log10(I1^2*I2/(I3*I4^2)))
+//RESULTS
+printf ('Equilibrium constant = %.2f ',K)
diff --git a/1967/CH13/EX13.2/13_2.sce b/1967/CH13/EX13.2/13_2.sce
new file mode 100755
index 000000000..94801f070
--- /dev/null
+++ b/1967/CH13/EX13.2/13_2.sce
@@ -0,0 +1,13 @@
+clc
+//initialisation of variables
+clear
+k1= 4600
+k2= -8.64
+k3= 1.86*10^-3
+k4= -0.12*10^-6
+k5= 12.07
+T= 600 //K
+//CALCULATIONS
+Kf= %e^(k1*(1/T)+k2*log10(T)+k3*T+k4*T^2+k5)
+//RESULTS
+printf ('Kf = %.3f ',Kf)
diff --git a/1967/CH13/EX13.3/13_3.sce b/1967/CH13/EX13.3/13_3.sce
new file mode 100755
index 000000000..81b834856
--- /dev/null
+++ b/1967/CH13/EX13.3/13_3.sce
@@ -0,0 +1,12 @@
+clc
+//initialisation of variables
+clear
+k= -8810 //cal
+k1= -7.46 //cal K^-1
+k2= 3.69*10^-3 //cal K^-2
+k3= -0.47*10^-6 //cak K^-3
+T= 298 //K
+//CALCULAATIONS
+dH= k+k1*T+k2*T^2+k3*T^3
+//RESULTS
+printf ('Standard heat of reaction = %.f cal',dH)
diff --git a/1967/CH13/EX13.4/13_4.sce b/1967/CH13/EX13.4/13_4.sce
new file mode 100755
index 000000000..a17842e3a
--- /dev/null
+++ b/1967/CH13/EX13.4/13_4.sce
@@ -0,0 +1,14 @@
+clc
+//initialisation of variables
+clear
+k1= -9130 //cal
+k2= 7.46 //cal K^-1
+k3= -3.69*10^-3 //K^-2
+k4= 0.235*10^-6 //K^-3
+k5= -12.07
+T= 298 //K
+R= 1.987 //cal deg^-1 mole^-1
+//CALCULATIONS
+dF= k1+k2*T*log(T)+k3*T^2+k4*T^3+k5*R*T
+//RESULTS
+printf ('Free energy = %.f cal',dF)
diff --git a/1967/CH13/EX13.5/13_5.sce b/1967/CH13/EX13.5/13_5.sce
new file mode 100755
index 000000000..fbfbf7522
--- /dev/null
+++ b/1967/CH13/EX13.5/13_5.sce
@@ -0,0 +1,11 @@
+clc
+//initialisation of variables
+clear
+T= 25 //C
+dF1= 61.44 //kcal
+dF= 54.65 //kcal
+R= 4.576 //cal deg^-1 mole^-1
+//CALCULATIONS
+Kf= 10^(-(dF1-dF)*10^3/(R*(273.2+T)))
+//RESULTS
+printf ('Kf at this temperature = %.2e ',Kf)
diff --git a/1967/CH13/EX13.6/13_6.sce b/1967/CH13/EX13.6/13_6.sce
new file mode 100755
index 000000000..99b949248
--- /dev/null
+++ b/1967/CH13/EX13.6/13_6.sce
@@ -0,0 +1,13 @@
+clc
+//initialisation of variables
+clear
+R= 4.576 //cal mole^-1 K^-1
+T= 25 //C
+p1= 122 //mm
+F1= -5.88 //kcal
+F2= -33 //kcal
+//CALCULATIONS
+dF= R*(273.2+T)*log10(p1/760)
+F= F2+F1+(dF/1000)
+//RESULTS
+printf ('Standard free energy change = %.f kcal',F)
diff --git a/1967/CH13/EX13.7/13_7.sce b/1967/CH13/EX13.7/13_7.sce
new file mode 100755
index 000000000..eeb5af902
--- /dev/null
+++ b/1967/CH13/EX13.7/13_7.sce
@@ -0,0 +1,15 @@
+clc
+//initialisation of variables
+clear
+r= 3.38*10^-4 //volt deg^-1
+F= 23070 //cal volt^-1 deg^-1
+Sagcl= 23 //E.U.mole^-1
+Shg= 18.5 //E.U.mole^-1
+Sag= 10.2 //E.U.mole^-1
+//CALCULATIONS
+dS= F*r
+shgcl= 2*-(dS-Sagcl-Shg+Sag)
+//RESULTS
+printf ('dS  = %.1f E.U.cal deg^-1',dS)
+printf ('\n molar entropy  = %.f E.U.mole^-1',shgcl)
+
diff --git a/1967/CH13/EX13.8/13_8.sce b/1967/CH13/EX13.8/13_8.sce
new file mode 100755
index 000000000..f5391597a
--- /dev/null
+++ b/1967/CH13/EX13.8/13_8.sce
@@ -0,0 +1,21 @@
+clc
+//initialisation of variables
+clear
+s1= 44.5 //cal deg^-1 mole^-1
+s2= 49 //cal deg^-1 mole^-1
+s3= 51.06 //cal deg^-1 mole^-1
+s4= 16.75 //cal deg^-1 mole^-1
+h1= -17.9 //kcal mole^-1
+h2= 0 //kcal mole^-1
+h3= -94 //kcal mole^-1
+h4= -68.3 //kcal mole^-1
+T= 25 //C
+n= 2
+//CALCULATIONS
+dS= s3+2*s4-s1-n*s2
+dH= h3+n*h4-h1-n*h2
+dF= -0.001*(273.2+T)*dS+dH
+//RESULTS
+printf ('Entropy Change = %.1f E.U',dS)
+printf ('\n Enthalpy Change = %.1f E.U',dH)
+printf ('\n Standard free energy = %.1f kcal',dF)
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)