initial commit / add all books

1 files changed, 41 insertions, 0 deletions

diff --git a/1922/CH9/EX9.10.b/9_10b.sce b/1922/CH9/EX9.10.b/9_10b.sce
new file mode 100755
index 000000000..d27c6bd4b
--- /dev/null
+++ b/1922/CH9/EX9.10.b/9_10b.sce
@@ -0,0 +1,41 @@
+clc

+clear

+//Initialization of variables

+kp=1.09

+kp2=0.154

+feed=[ 1 5 0 0 0 ]

+//calculations

+x=poly(0,"x")

+vec=roots(kp/4^4 /4 *(1-x)*(5-2*x)^2 *(6+2*x)^2 -x^5)

+x=vec(5)

+pro=feed- [x 2*x -x -4*x 0]

+y=poly(0,"y")

+vec2=roots(kp2*(0.273-y)*(0.727-y)*(7.454+2*y)^2 - 4*y^2 *(2.908+2*y)^2 *4)

+y=vec2(4)

+pro2=pro-[ y 0 y -2*y -2*y]

+z=poly(0,"z")

+vec3= roots(kp*(0.189-z)*(3.546-2*z)^2 *(7.622+2*z)^2 -(0.643+z)*(3.076+4*z)^4 *4)

+z=vec3(5)

+pro3=pro2 - [z 2*z -z -4*z 0]

+w=poly(0,"w")

+vec4=roots(kp2*(0.229-w)*(0.603-w)*(7.542+2*w) - (2.916+2*w)^2 *(0.168+2*w)^2 *4)

+w=vec4(4)

+w=0.01

+pro4=pro3-[w 0 w -2*w -2*w]

+//results

+disp("feed = ")

+format('v',6);feed

+disp(feed)

+disp("After reactor 1,")

+format('v',6);pro

+disp(pro)

+disp("After reactor 2,")

+format('v',6);pro2

+disp(pro2)

+disp("After reactor 3,")

+format('v',6);pro3

+disp(pro3)

+disp("After reactor 4")

+format('v',6);pro4

+disp(pro4)

+disp("The answers are a bit different due to rounding off error in textbook")