1 files changed, 21 insertions, 0 deletions
diff --git a/632/CH7/EX7.7/example7_7.sce b/632/CH7/EX7.7/example7_7.sce
new file mode 100755
index 000000000..84c0bf150
--- /dev/null
+++ b/632/CH7/EX7.7/example7_7.sce
@@ -0,0 +1,21 @@
+//clc()
+//lnPas = 14.5463 - 2940.46/(T - 35.93)
+//lnPbs = 14.2724 - 2945.47 / (T - 49.15)
+//xa = (P - Pbs)/(Pas - Pbs)
+//Ya = Pas * (P - Pbs)/(P * (Pas - Pbs))
+Ya = 0.4;
+P = 65;//kPa
+//various temperature value are assumed and tried till LHS = RHS, we get
+T = 334.15;//K
+Pas = exp(14.5463 - 2940.46/(T - 35.93));
+Pbs = exp(14.2724 - 2945.47 / (T - 49.15));
+xa = (P - Pbs)/(Pas - Pbs);
+disp("K",T,"(a)The Dew point temperature at 65 kPa = ")
+disp(xa,"    Concentration of the first drop of liquid = ")
+T1 = 327;//K
+Pas1 = exp(14.5463 - 2940.46/(T1 - 35.93));
+Pbs1 = exp(14.2724 - 2945.47 / (T1 - 49.15));
+xa1 = Ya * Pbs1 / (Pas1 - Ya*(Pas1 - Pbs1));
+P1 = xa1 * Pas1 / Ya;
+disp("kPa",P1,"(b)The dew point pressure at 327 K = ") 
+disp(xa1,"    Concentration at 327K = ")