1 files changed, 37 insertions, 0 deletions
diff --git a/249/CH18/EX18.5/18_05.sce b/249/CH18/EX18.5/18_05.sce
new file mode 100755
index 000000000..deab73fbb
--- /dev/null
+++ b/249/CH18/EX18.5/18_05.sce
@@ -0,0 +1,37 @@
+clear
+clc
+CAo=0.1;
+eA=3;
+rA=[3.4;5.4;7.6;9.1];
+CA=[0.039;0.0575;0.075;0.092];
+XA=zeros(4,1);
+inv_rA=zeros(4,1);
+for i=1:4
+XA(i)=(1-CA(i)/CAo)/(1+eA*CA(i)/CAo);
+inv_rA(i)=1/rA(i);
+end
+//W=FAo*integral(dXA/-rA) from 0 to 0.35
+//Using Trapezoidal rule to find area,XA must be in increasing order
+//Sorting XA and accordingly inv_rA
+for i=1:4
+    small=XA(i);
+    for j=i:4
+        next=XA(j);
+        if small>next
+            temp=XA(i);
+            XA(i)=XA(j);
+            XA(j)=temp;
+            temp1=inv_rA(i);
+            inv_rA(i)=inv_rA(j);
+            inv_rA(j)=temp1;
+        end
+    end
+end
+plot(XA,inv_rA)
+xlabel('XA');ylabel('-1/rA');
+//extending points to include XA=0.35
+XA(5)=0.35;inv_rA(5)=0.34;
+Area=inttrap(XA,inv_rA);
+W=Area*2000;
+printf("Amount of catalyst needed(kg) is %f",W)
+disp('The answer slightly differs from those given in book as trapezoidal rule is used for calculating area')