1 files changed, 15 insertions, 0 deletions
diff --git a/1394/CH15/EX15.4.1/Ex15_4_1.sce b/1394/CH15/EX15.4.1/Ex15_4_1.sce
new file mode 100755
index 000000000..734245cb5
--- /dev/null
+++ b/1394/CH15/EX15.4.1/Ex15_4_1.sce
@@ -0,0 +1,15 @@
+
+clc
+//intialization of variables
+tB1 = 38 // days , breakthrough time
+tE1 = 46 // days, exhaustion time
+c = 2 // number of times flow doubled
+//Calculations
+theta1 = 2*tB1/(tB1+tE1)// in the first case
+ratio1 = 1-theta1 // ratio of unused bed length to total bed length
+ratio2 = ratio1*c
+tB2 = ((1/c)*(tB1 + 0.5*(tE1-tB1)))*ratio2//breakthrough time for second case
+tE2 = (c-ratio2)*tB2/ratio2//exhaustion time for second case
+//Results
+//answwer slightly wrong in textbook
+printf("The breakthrough time for this case is %.1f days",tB2)