plotG(vx,vy,char,size,color)

9886fc08-03cd-4650-a1ac-f8f10b3ff521 26 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

xL 0 0.0917 0.1842 0.2779 0.3732 0.4703 0.5699 0.6721 0.7776 0.8867 1111 mat :

xV 0 0.1543 0.2943 0.4188 0.5315 0.6327 0.7236 0.8049 0.8776 0.9424 1111 mat :

data xL xV augment :

4 appVersion 0.98.6179.21440 ≡

Benzene purity in distillate
Benzene purity in bottom
Benzene fraction in feed
'q' value
Reflux ratio xD 97 % * : xB 2 % * : xF 40 % * : q 1.5 : R 3.5 : 51 sys

x qline x q q 1 - / * xF q 1 - / - : x top_op x R 1 R + / * xD 1 R + / + : 21 sys

System data given by Valentino [20201119]

x bottom_op x R xF xB - (* q xD xB - (* + q xD xB - (* R xF xB - (* + xD - xF + / * xB xF xD - (* q xD xB - (* R xF xB - (* + xD - xF + / + :

Thiele(X,Y,0) "K" .... Thiele(X,Y,1) "Table"

X Y Δ Thiele ε 1018 -^: n X rows : i 1 n range M i 1 el Y i el : for j 2 n range i j n range dM M i j 1 - el M j 1 - j 1 - el - : dM dM abs ε < ε M i j el X i el X j 1 - el - dM / : if : 21 line for for i 1 n range c i 1 el X i el : c i 2 el M i i el : 21 line for Δ 0 ≡ c M if 61 line :

data List collects index list select 111 range : ct 1 : i select a ct el data i 1 el : b ct el data i 2 el : ct ct 1 + : 31 line for a b augment 51 line :

select the support points XY data :

length Thiele j XY rows : 11

Thiele on length of support points 'j' K X x Cfr N j 1 + el K j 1 + el : c j el N j el : a j el K j el : jj j 1 - (0 range c jj el a jj 1 + el : a jj el K jj el x X jj el - c jj el / + : 21 line for a 1 el 51 line :

Cfr(K,X,x) Expansion

symbolic expansion of Thiele continued fraction algo style K X x Cfr K 1 el x X 1 el - K 2 el x X 2 el - K 3 el x X 3 el - K 4 el x X 4 el - K 5 el x X 5 el - K 6 el x X 6 el - K 7 el x X 7 el - K 8 el x X 8 el - K 9 el x X 9 el - K 10 el x X 10 el - K 11 el / + / + / + / + / + / + / + / + / + / +

recast the Cf on the system  [K,X,x] x Cf K X x Cfr :

data set vx data 1 col : vy data 2 col : 21 sys

support points X XY 1 col : Y XY 2 col : 21 sys

Thiele K's coeff's
@ support points X X Y 0 Thiele 00 0.0917 0.5943 0.1842 2.9277 0.2779 0.3913 - 0.3732 0.4692 - 0.4703 1.1534 - 0.5699 4.3393 - 0.6721 0.0613 0.7776 2.5357 0.8867 0.8674 1 0.2305 112 mat

Thiele K's continued fraction K X Y 0 Thiele 2 col : 0 0.5943 2.9277 0.3913 - 0.4692 - 1.1534 - 4.3393 - 0.0613 2.5357 0.8674 0.2305 111 mat

Export fitted data set Cont x 01 0.01 range : i 1 x rows range v i el x i el Cf eval : for x v augment 31 line :

jnct xD q 1 - (* xF R 1 + (* + q R + / : 0.457

x operating_line x bottom_op x jnct < x top_op jnct x < cases :

x y axes x y / - :

x xLine 1 x * 0 + :

u 0.45 :

User Cf x Cf u - x 01 solve : op x operating_line u - x 01 solve : 21 sys :

pairwise sol's @ 'u' User 0.3031 0.361 21 sys

plot x Cf data x operating_line axes jnct 0.75 jnct 0.5 22 mat Cf u op u 22 mat x xLine 71 sys :

plot

x Th 10 x ≤ (x 1 ≤ (& cases :

x Op 10 x ≤ (x 1 ≤ (& cases :

Export the solved  binary McCabe data set McCabe i 1 u rows range sol i el Cf x Cf u i el - x 01 solve : op x operating_line u i el - x 01 solve : 21 mat : for sol u augment 21 line :

u xB 0.1543 0.2943 0.4188 0.5315 0.6327 0.7236 0.8049 0.8776 0.9424 xD 111 mat :

McCabe McCabe :

data UnestRow v data 1 el : j 2 data rows range v v data j el stack : for v 31 line :

i 1 u rows range duplicate i el u i el u i el 21 mat : for

Tray observer j 5 :

Unest_solutions McCabe 1 col UnestRow : m Unest_solutions : n duplicate UnestRow : n n : ZigZag m n augment : Data data . 14 black plot : 61 sys

dCf McCabe j 1 el 1 el : dop McCabe j 1 el 2 el : 21 sys

product_1 dCf 0 + 40 red 15 mat : product_2 dop 0 + 40 red 15 mat : 21 sys

TRAY xL j el xV j el + 100 black 15 mat :

TRAY 0.3732 0.5315 + 100 black 15 mat

plot ZigZag x Cf x Th * x operating_line x Op * Data product_1 product_2 TRAY 71 sys :

plot

McCabe j 1 el McCabe j 2 el 12 mat 0.3732 0.4257 21 mat 0.5315 12 mat

dCf dop 12 mat 0.3732 0.4257 12 mat

Smoe kind of step construct

plotG(vx,vy,char,size,color)

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

xL 0 0.0917 0.1842 0.2779 0.3732 0.4703 0.5699 0.6721 0.7776 0.8867 1111 mat :

xV 0 0.1543 0.2943 0.4188 0.5315 0.6327 0.7236 0.8049 0.8776 0.9424 1111 mat :

data xL xV augment :

Data data . 12 black plot :

m 0 0.0917 0.0917 0.1842 0.1842 0.2779 0.2779 0.3732 0.3732 0.4703 0.4703 0.5699 0.5699 0.6721 0.6721 0.7776 0.7776 0.8867 0.8867 11211 mat :

n 00 0.1543 0.1543 0.2943 0.2943 0.4188 0.4188 0.5315 0.5315 0.6327 0.6327 0.7236 0.7236 0.8049 0.8049 0.8776 0.8776 0.9424 0.9424 1211 mat :

step m n augment :

x Cf x Th * step x operating_line x Op * Data 4 1 sys

0.0113 0.02 21 mat 0.02 0.0917 0.1265 21 mat 0.1543 0.1842 0.2375 21 mat 0.2943 0.2779 0.3363 21 mat 0.4188 0.3732 0.4257 21 mat 0.5315 0.4703 0.5363 21 mat 0.6327 0.5699 0.6532 21 mat 0.7236 0.6721 0.7577 21 mat 0.8049 0.7776 0.8512 21 mat 0.8776 0.8867 0.9345 21 mat 0.9424 0.9387 0.97 21 mat 0.97 112 mat