From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 249/CH16/EX16.1/16_01.sce | 24 ++++++++++++++++++++++++
 1 file changed, 24 insertions(+)
 create mode 100755 249/CH16/EX16.1/16_01.sce

(limited to '249/CH16/EX16.1/16_01.sce')

diff --git a/249/CH16/EX16.1/16_01.sce b/249/CH16/EX16.1/16_01.sce
new file mode 100755
index 000000000..5d218f839
--- /dev/null
+++ b/249/CH16/EX16.1/16_01.sce
@@ -0,0 +1,24 @@
+clear
+clc
+Co=1;k=1;t=1;//given
+//Scheme A
+//For mixed flow reactor
+//t=(Co-C1)/KC1^2
+C1=(-1+sqrt(1-4*t*(-Co)))/2*t;
+//For the plug flow reactor
+//t=1/k(1/C2-1/C1)
+C2=C1/(1+k*t*C1);
+printf("\n Conversion for flow scheme A is %f",C2)
+//Scheme B
+//For plug flow
+C3=Co/(1+k*t*Co);
+//For mixed flow reactor
+C4=(-1+sqrt(1-4*t*(-C3)))/2*t;
+printf("\n Conversion for flow scheme B is %f",C4)
+//Scheme C,D,E
+//Using exit age distribution fn for 2 equal plug-mixed flow reactor system,using fig 12.1
+t_bar=2;
+in=1000;
+C5=integrate('(2/t_bar)*(exp(1-2*t/t_bar))/(1+Co*k*t)','t',t_bar/2,in);
+printf("\n Conversion for flow scheme C,D,E is %f",C5)
+
-- 
cgit