x := ( 0.0001 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.

7 0.75 0.8 0.85 0.9 0.95 0.99999 )

T
x0 := x x1 := 1 - x0

long := length ( x0 ) length(v), cuenta el numero de elementos del vector v

0 
dx :=   0 
1  dy :=   PmmHg := 760 Pbar := 1.01325
1 

HEXANO (0)  507.5   30.1 

Tc :=   Pc :=  
 562.2   48.9 
BENCENO (1)

 -7.46765  Vpb :=  1.44211   -3.28222  Vpd :=  -2.50941 

Vpa :=     Vpc :=    
 -6.98273   1.33213   -2.62863   -3.33399 

 

( b ) :=  1 -
b 

 Tc 


  1  ( Vpa ( b) + Vpb ( b)
1.5 3
+ Vpc ( b) + Vpd ( b)
6
)

Pvp( b ) := Pc e
1- ( b) 


 0.0513   0.108   77.95 

 0.1216   0.216   75.95 
     
 0.2189   0.334   73.85 
 0.3108   0.422 
     72.45 
 0.3811   0.483   71.55 
xexp :=  0.503  yexp :=  0.579   
texp :=  70.45 
 0.6031   0.657 
     69.75 
 0.7171   0.746   
 0.7982   0.814   69.15 
 68.75 
 0.8982   0.902   68.65 
     
 0.9402   0.942 
 68.55 

NSG := 3 NC := 2 i := 1 .. 2 z := 10 m := 1 .. NSG k := 1 .. NSG

R := 0.9011 Q := 0.848 v := 2 v := 0
1 1 1, 1 1, 2

R := 0.674 Q := 0.54 v := 4 v := 0
2 2 2, 1 2, 2

R := 0.5313 Q := 0.4 v := 0 v := 6
3 3 3, 1 3, 2
a := 0 a := 0 a := 61.13
1, 1 1, 2 1, 3

a := 0 a := 0
2, 1 2, 2 a := 61.13
2, 3
a := -11.12 a := -11.12
3, 1 3, 2 a := 0
3, 3
PARTE COMBINATORIAL

LnC( b , c) := Tb
x  x0
1 c
x 1-x
2 1
for i  1 .. NC
NSG
r 
i (vk , iRk)
k =1
NSG
q 
i (vk , iQk)
k =1
for i  1 .. NC
r x
i i

i NC
(rj xj)
j =1
q x
i i

i NC
(qj xj)
j =1

( i i) ( i )
z
l   r -q - r -1
i 2
for i  1 .. NC

 i   NC
Ln C  ln   q  ln (xj lj)
z i i
+ +l - 
i  xi  2 i  i  i xi
    j =1
re1  Ln C
0, 0 1
re1  Ln C
0, 1 2
re1
LnR ( b , c) := x  x0
1 c
x 1 -x
2 1
Tb
for m  1 .. NSG
for i  1 .. NC
v
m, i
X 
m, i NSG
v
n, i
n =1
for m  1 .. NSG
for i  1 .. NC
X Q
m, i m

m, i NSG
(Xn , i Qn)
n =1
for m  1 .. NSG
for n  1 .. NSG
- am , n
T
e ( T)  e
 e ( T)
m, n
for k  1 .. NSG
for i  1 .. NC
  NSG  NSG   
 Q  1 - ln 
m, i k, m
Ln ´
k, i k   ( 
m, i
-
m , k)

 NSG 
m = 1 


m=1 
n = 1
( n, i

n , m)


  
for m  1 .. NSG
NC
(xj  vm , j )
j =1

m NSG NC
(xj  vn , j )
n =1 j =1
for m  1 .. NSG
Q
m m

m NSG
( n
Q
n )
n =1
for k  1 .. NSG
  NSG  NSG   
 Q  1 - ln
k, m
( m , k)
-  
m
Ln 
k k   m  NSG 
m = 1  m=1



n = 1
( n

n , m)


  
for i  1 .. 2
NSG
Ln R 
i vk , i ( Ln k
- Ln ´
k, i 
)
k =1
re2  Ln R
0, 0 1
re2  Ln R
0, 1 2
re2
 1 ( b , c) := re3  LnC ( b , c) + LnR( b , c)
0, 0 0, 0
re3
e

2 ( b , c) := re4  LnC ( b , c) + LnR( b , c)

0, 1 0, 1
re4
e

n := 373.15

Tteorica( n ) := for c  0 .. long - 1

x0  1 ( b , c)  Pvp( b )
c 0
y0 ( b , c) 
Pbar
x1  2 ( b , c)  Pvp( b )
c 1
y1 ( b , c) 
Pbar
bn
a  root ( y0 ( b , c) + y1 ( b , c) - 1 , b )
c
a

Tt := Tteorica( n ) - 273.15

TC := Tteorica( n )
0 
0 
dx :=  
dy :=  
1  1 

y0 := for c  0 .. long - 1 y1 := for c  0 .. long - 1

t  TC t  TC
c c

(
1 TC , c  x0 ) (
2 TC , c  x1)
( ) ( )
c c c c
b  Pvp TC 0 b  Pvp TC 1
c c Pbar c c Pbar
b b

0 := for c  0 .. long - 1 1 := for c  0 .. long - 1

b  TC b  TC
c c
B  1 ( b , c) B  2 ( b , c)
c c
B B

comprobacion := y0 + y1
 0 0 0 0
0
0 80.164 0 2.52·10-4 0 1.786 0 1
0 1·10 -4
1 78.018 1 0.109 1 1.651 1 1.002
1 0.05
2 76.366 2 0.194 2 1.539 2 1.008
2 0.1
3 75.06 3 0.263 3 1.446 3 1.017
3 0.15
4 74 4 0.322 4 1.368 4 1.029
4 0.2
5 73.124 5 0.373 5 1.303 5 1.044
5 0.25
6 72.385 6 0.419 6 1.247 6 1.062
6 0.3
7 71.755 7 0.461 7 1.201 7 1.082
7 0.35
8 71.212 8 0.502 8 1.162 8 1.104
8 0.4
Tt = 9 70.741 y0 = 9 0.54 0= 9 1.129 1= 9 1.128
10 70.331 x0 = 9 0.45
10 0.578 10 1.101 10 1.153
10 0.5
11 69.975 11 0.616 11 1.078 11 1.181
11 0.55
12 69.667 12 0.654 12 1.059 12 1.21
12 0.6
13 69.406 13 0.692 13 1.043 13 1.241
13 0.65
14 69.187 14 0.731 14 1.03 14 1.273
14 0.7
15 69.011 15 0.771 15 1.02 15 1.307
15 0.75
16 68.876 16 0.813 16 1.013 16 1.342
16 0.8
17 68.782 17 0.856 17 1.007 17 1.379
17 0.85
18 68.731 18 0.902 18 1.003 18 1.416
18 0.9
19 68.724 19 0.95 19 1.001 19 1.455
19 0.95
20 68.762 20 1 20 1 20 1.495
20 1

0 0 0
0 77.95 0 0.051 0 0.108
1 75.95 1 0.122 1 0.216
2 73.85 2 0.219 2 0.334
3 72.45 3 0.311 3 0.422

texp = 4 71.55 xexp = 4 0.381 yexp = 4 0.483

5 70.45 5 0.503 5 0.579
6 69.75 6 0.603 6 0.657
7 69.15 7 0.717 7 0.746
8 68.75 8 0.798 8 0.814
9 68.65 9 0.898 9 0.902
10 68.55 10 0.94 10 0.942