4 files changed, 56 insertions, 0 deletions
diff --git a/249/CH14/EX14.1/14_01.sce b/249/CH14/EX14.1/14_01.sce
new file mode 100755
index 000000000..1c1cc15f2
--- /dev/null
+++ b/249/CH14/EX14.1/14_01.sce
@@ -0,0 +1,8 @@
+clear
+clc
+//Original and new length(m)
+L1=32;L2=50;
+sigma1=8;
+// For small deviaqtion from plug flow,sigma_sqr is directly proportional to L
+sigma2=sigma1*sqrt(L2/L1);
+printf("\n No of bottles of rose expected is %f",sigma2)
+\ No newline at end of file
diff --git a/249/CH14/EX14.2/14_02.sce b/249/CH14/EX14.2/14_02.sce
new file mode 100755
index 000000000..3f3271202
--- /dev/null
+++ b/249/CH14/EX14.2/14_02.sce
@@ -0,0 +1,8 @@
+clear
+clc
+sigma1=14;sigma2=10.5;
+L1=119;
+//spread of curve is directly proportional to sqrt of distance from origin
+L=sigma1^2*L1/(sigma1^2-sigma2^2);
+printf("\n The dumping of toxic phenol must have occured within %f",L)
+printf("miles upstream of cincinnati")
+\ No newline at end of file
diff --git a/249/CH14/EX14.3/14_03.sce b/249/CH14/EX14.3/14_03.sce
new file mode 100755
index 000000000..d08581cb6
--- /dev/null
+++ b/249/CH14/EX14.3/14_03.sce
@@ -0,0 +1,20 @@
+clear
+clc
+vo=1;
+t1=1/6;
+t2=1;
+t3=11/6;
+w=1/10;
+//Ratio of areas of the first 2 peaks
+A2_by_A1=0.5;
+R=A2_by_A1/(1-A2_by_A1);
+//From the location of 1st peak
+V1=(R+1)*vo*t1;
+//From the time between peaks
+V2=(R*vo)*((t2-t1)-(t1));
+//From fig 14.3
+N=1+(2*(t1/w))^2;
+printf("\n The reflux ratio is %f",R)
+printf("\n The volume of 1st tank is %f", V1)
+printf("\n The volume of 2nd tank is %f",V2)
+printf("\n The number of tanks are %f ",N)
+\ No newline at end of file
diff --git a/249/CH14/EX14.4/14_04.sce b/249/CH14/EX14.4/14_04.sce
new file mode 100755
index 000000000..2f9ddf7db
--- /dev/null
+++ b/249/CH14/EX14.4/14_04.sce
@@ -0,0 +1,20 @@
+clear
+clc
+//from fig E14.4a
+t2=280;t1=220;
+sigma1_sqr=100;sigma2_sqr=1000;
+dt=t2-t1;
+dsigma_sqr=sigma2_sqr-sigma1_sqr;
+N=dt^2/dsigma_sqr;
+for t=1:200
+//For N tank in series
+E(t)=((t^(N-1))*(N^N)*exp(-t*N/dt))/((factorial(N-1))*(dt^N));
+end
+for i=1:200
+    t(i)=i;
+end
+plot(t,E)
+xlabel('time(sec)')
+ylabel('E(sec^-1)')
+title('Shape of E curve')
+