diff options
Diffstat (limited to '1394/CH3/EX3.3.1/Ex3_3_1.sce')
| -rwxr-xr-x | 1394/CH3/EX3.3.1/Ex3_3_1.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/1394/CH3/EX3.3.1/Ex3_3_1.sce b/1394/CH3/EX3.3.1/Ex3_3_1.sce new file mode 100755 index 000000000..895e4340e --- /dev/null +++ b/1394/CH3/EX3.3.1/Ex3_3_1.sce @@ -0,0 +1,23 @@ + +clc
+//initialization of variables
+// At 6 degree centigrade
+p1sat = 37 // Vapor pressure of benzene in mm Hg
+p = 760 // atmospheric pressure in mm Hg
+y1l = 0
+y10 = p1sat/p
+n1byDcbyl = log((1-y1l)/(1-y10))// because flux n1 = D*c/l * ln(1-y11/1-y10)
+n2byDcbyl = y10-y1l // Flux calculated assuming dilute solution as n1 = Dc/l*(y10-y1l)
+err1 = ((n1byDcbyl-n2byDcbyl)/n2byDcbyl)*100 // Percentage error
+printf("The error in measurement at 6 degree centigrade is %.1f percent",err1)
+// At 60 degree centigrade
+p1sat = 395 // Vapor pressure of benzene in mm Hg
+p = 760 // atmospheric pressure in mm Hg
+y1l = 0
+y10 = p1sat/p
+n1byDcbyl = log((1-y1l)/(1-y10))// because flux n1 = D*c/l * ln(1-y11/1-y10)
+n2byDcbyl = y10-y1l // Flux calculated assuming dilute solution as n1 = Dc/l*(y10-y1l)
+err1 = ((n1byDcbyl-n2byDcbyl)/n2byDcbyl)*100 // Percentage error
+printf("\n The error in measurement at 60 degree centigrade is %.1f percent",err1)
+
+
|