function datos = ParInBin(prop)
% Parámetros de interacción binaria
% Se devuelve la matriz de datos para la propiedad pedida


% Propiedad 1
% Kij para ecuación de Redlich-Kwong-Soave
RKSKIJ=[
%      H               M               B               T               F
       1               -0.0222         0               0.39            0  % H
       -0.0222         0               0.0209          0.0978          0  % M 
       0               0.0209          0               0               0  % B
       0.39            0.0978          0               0               0  % T
       0               0               0               0               0  % F
    ];
    
% Propiedad 2    
% Parámetros para constante de Henry
% Pa 

HENRY=[
%   H-B            H-T            M-B            M-T
-49.1847        17.295799       -33.072399      84.633301
2445.6001       628.590027      789.23999       -1220.19995
11.752          0               9.5495          -12.693
-0.020799       0               -0.021042       0.031756
273.65          291.25          283.25          298.15
335.95          305.25          333.15          333.15
];

% Propiedad 3
% Parámetros disponibles para modelo UNIQUAC de coeficiente de actividad
% K
UNIQUAC=[
%  i=T,j=B               i=B,j=F
 -1.082100000            0.0            %  aij
 1.337700000             0.0            %  aji
 424.1051000             47.99790000    %  bij
 -522.8772000            -73.17830000   %  bji
 0.0                     0.0            %  cij
 0.0                     0.0            %  cji
 0.0                     0.0            %  dij
 0.0                     0.0            %  dji
 313.1500000             378.2500000    %  Tlower
 383.9000000             556.8500000    %  Tupper
];

switch prop
case 1 % Kij
   datos= RKSKIJ;
case 2 % Henry
   datos=HENRY;
case 3 % UNIQUAC
   datos=UNIQUAC;
end

   
   