10 files changed, 112 insertions, 0 deletions
diff --git a/1652/CH13/EX13.1/13_1.sce b/1652/CH13/EX13.1/13_1.sce
new file mode 100755
index 000000000..d588e21ab
--- /dev/null
+++ b/1652/CH13/EX13.1/13_1.sce
@@ -0,0 +1,8 @@
+clc
+//Initialization of variables
+n1=10 //mol
+n2=12 //mol
+//calculations
+dn=n1-n2
+//results
+printf("dHp = dEv- %d *RT",dn)
diff --git a/1652/CH13/EX13.10/13_10.sce b/1652/CH13/EX13.10/13_10.sce
new file mode 100755
index 000000000..7e4df4d01
--- /dev/null
+++ b/1652/CH13/EX13.10/13_10.sce
@@ -0,0 +1,10 @@
+clc
+//Initialization of variables
+x=0.5
+P=0.468 //atm
+//calculations
+P1=x*P
+P2=x*P
+Kp=P1*P2
+//results
+printf("Final pressure after quilibrium = %.4f atm^2",Kp)
diff --git a/1652/CH13/EX13.2/13_2.sce b/1652/CH13/EX13.2/13_2.sce
new file mode 100755
index 000000000..5e2e41710
--- /dev/null
+++ b/1652/CH13/EX13.2/13_2.sce
@@ -0,0 +1,8 @@
+clc
+//Initialization of variables
+Ht1=-22063 //cal
+T=298.15 //K
+//calculations
+H=Ht1 +0.5293*T + 0.3398*10^-3 *T^2 - 2.039*10^-7 *T^3
+//results
+printf("Enthalpy = %d cal",H)
diff --git a/1652/CH13/EX13.3/13_3.sce b/1652/CH13/EX13.3/13_3.sce
new file mode 100755
index 000000000..a31a5211b
--- /dev/null
+++ b/1652/CH13/EX13.3/13_3.sce
@@ -0,0 +1,7 @@
+clc
+//Initialization of variables
+Cp=0.797 //cal/deg/mol
+//calculations
+S=Cp/3
+//results
+printf("Entropy = %.3f eu/mol",S)
diff --git a/1652/CH13/EX13.4/13_4.sce b/1652/CH13/EX13.4/13_4.sce
new file mode 100755
index 000000000..1478c4661
--- /dev/null
+++ b/1652/CH13/EX13.4/13_4.sce
@@ -0,0 +1,10 @@
+clc
+//Initialization of variables
+T1=77.32 //K
+P=1 //atm
+T2=126 //K
+Pc=33.5 //atm
+//calculations
+dS=27/32 *1.987*P/Pc *(T2/T1)^3
+//results
+printf("Change in entropy = %.2f eu/mol",dS)
diff --git a/1652/CH13/EX13.5/13_5.sce b/1652/CH13/EX13.5/13_5.sce
new file mode 100755
index 000000000..e36ef4fc0
--- /dev/null
+++ b/1652/CH13/EX13.5/13_5.sce
@@ -0,0 +1,19 @@
+clc
+//Initialization of variables
+S1=57.47
+S2=50.34
+S3=49
+H1=8.09
+H2=21.06
+H3=0
+F1=12.39
+F2=20.72
+F3=0
+//calculations
+dS=S1-S2-0.5*S3
+dH=H1-H2-0.5*H3
+dF=F1-F2-0.5*F3
+//results
+printf("Change in entropy = %.2f eu",dS)
+printf("\n Change in enthalpy = %.2f kcal",dH)
+printf("\n Change in free energy = %.2f kcal",dF)
diff --git a/1652/CH13/EX13.6/13_6.sce b/1652/CH13/EX13.6/13_6.sce
new file mode 100755
index 000000000..4bb5de92e
--- /dev/null
+++ b/1652/CH13/EX13.6/13_6.sce
@@ -0,0 +1,13 @@
+clc
+//Initialization of variables
+P1=0.01
+P2=0.1
+P3=0.01
+dF0=-54640 //cal
+T=298.15 //K
+R=1.987 //cal/deg
+//calculations
+Qp=P1/(P2*P3^0.5)
+dF=dF0+R*T*log(Qp)
+//results
+printf("change in free energy = %d cal",dF)
diff --git a/1652/CH13/EX13.7/13_7.sce b/1652/CH13/EX13.7/13_7.sce
new file mode 100755
index 000000000..fb8d80aac
--- /dev/null
+++ b/1652/CH13/EX13.7/13_7.sce
@@ -0,0 +1,11 @@
+clc
+//Initialization of variables
+disp("From table 13.4 ")
+logKfwater=40.04724
+logKfH2=0
+logKfO2=0
+//calculations
+logK=logKfwater-logKfH2-0.5*logKfO2
+K=10^logK
+//results
+printf("Equilibrium constant = %.4e",K)
diff --git a/1652/CH13/EX13.8/13_8.sce b/1652/CH13/EX13.8/13_8.sce
new file mode 100755
index 000000000..5cfd955d9
--- /dev/null
+++ b/1652/CH13/EX13.8/13_8.sce
@@ -0,0 +1,10 @@
+clc
+//Initialization of variables
+Kp=1.1*10^40 //atm^-0.5
+dn=-0.5
+R=0.08206 //lt atm/deg mol
+T=298.15 //K
+//calculations
+Kc=Kp*(R*T)^(-dn)
+//results
+printf("Kc = %.1e (mol/lt)^-0.5",Kc)
diff --git a/1652/CH13/EX13.9/13_9.sce b/1652/CH13/EX13.9/13_9.sce
new file mode 100755
index 000000000..607bf83fe
--- /dev/null
+++ b/1652/CH13/EX13.9/13_9.sce
@@ -0,0 +1,16 @@
+clc
+//Initialization of variables
+Kp=0.141 //atm
+P=1 //atm
+nu=2
+R=0.08206 //lt atm/deg mol
+T=298.15 //K
+M=92.02 //g/mol
+//calculations
+a=poly(0,"a");
+p=Kp*a^2 +4*a^2*P -Kp
+z=roots(p)
+alpha=z(1)
+wbyV=P*M/(R*T*(1+(nu-1)*alpha))
+//results
+printf("Density of the equilibrium mixture = %.2f g/lt",wbyV)