With this notebook you can calculate the following quantities for a 
BINARY (or one-component) mixture of STICKY Hard Spheres:

--> yHsBmcsl[i,j,nco,eeta] =  the contact value of the cavity function for HS given by the BMCSL formula.
-->yShsPy[i,j,nco,con,eeta,taus] = the contact value of the cavity function for SHS given by the PY approximation.
--> yHsPy[i,j,nco,eeta] = the contact value of the cavity function for HS given by the PY approximation.
--> yProCon[i,j,nco,con,eeta,taus] = the contact value of the cavity function for SHS given by equation (3.1) of the paper  A. Santos, S. B. Yuste and M. López de Haro, (J. Chem. Phys. vol 109, num. 16, 22 October 1998.
--> 1/ chim1HsBmcsl[nco,eeta] = the isothermal susceptibility given by the BMCSL formula for the HS mixture.
--> chiShsPy[nco,eeta,taus] = the isothermal susceptibility given by the PY approximation for the SHS mixture.
--> 1/ chim1HsPy[nco,eeta] = the isothermal susceptibility given by the PY approximation for the HS mixture.
 --> chiPro[nco,eeta,taus] = the isothermal susceptibility of the SHS mixture  given by the approximate formula (3.2) of the paper  A. Santos, S. B. Yuste and M. López de Haro, (J. Chem. Phys. vol 109, num. 16, 22 October 1998).

The mixture is defined by:
-> nco=number of species
-> eeta=volume packing fraction
-> con=list with the molar fractions of each species,e.g, {x_1,x_2}, 
where  x_1 is the molar fraction of the species with the largest diameter
-> ratio=list with the diameter ratios, e.g., {1,ratio_2} 
 where ratio_2=sigma_2/sigma_1=sigma_2
-> taus=list with the tau parameters characterizing the adhesiveness of the 
intermolecular potential , e.g., {{tau_11,tau_12},{tau_21,tau_22}}.

How to use this notebook:
(1) Introduce the  nco, con  and ratio  (the values of eta and taus can be 
given later)  of your hard-sphere mixture in the corresponding cell (see 
below) , and 
(2) Evaluate the next outest cell.

This cell must be evaluated for each mixture  with different values of nco, 
con and/or ratio

Here you have to introduce the values of nco, eeta, con , ratio and  taus 
corresponding to your sticky-hard-sphere mixture

For example, these are the values for our Examples  1 and 2  (binary mixtures)

nco=2;
 con={1/2,1/2};
ratio={1,1};

For the  Example 3  (binary mixture) we used this values:

nco=2;
 con={2/3,1/3};
ratio={1,1/2};

For the  Example 4  (one-component mixture) we used this values:

nco=2;   
con={1/2,1/2};
ratio={1,1};

The core of the program...

EXAMPLES

Example 1. A binary mixture with
 nco=2;   
con={1/2,1/2};
ratio={1,1};
eeta=0.4;
taus={{1,10},{10,1}};

eeta=0.4;
taus={{1,10},{10,1}};

{chim1HsBmcsl[2,eeta],1/chiShsPy[nco,eeta,taus],chim1HsPy[nco,eeta],
  1/chiPro[nco,eeta,taus]}

{Table[yHsBmcsl[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yShsPy[i,j,nco,con,eeta,taus],{i,nco},{j,nco}],
	Table[yHsPy[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yProCon[i,j,nco,con,eeta,taus],{i,nco},{j,nco}]}

Example 2 A binary mixture with
 nco=2;   
con={1/2,1/2};
ratio={1,1};
eeta=0.4;
taus={{2,10},{10,2}};

eeta=0.4;
taus={{2,10},{10,2}};

{chim1HsBmcsl[2,eeta],1/chiShsPy[nco,eeta,taus],chim1HsPy[nco,eeta],
  1/chiPro[nco,eeta,taus]}

{Table[yHsBmcsl[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yShsPy[i,j,nco,con,eeta,taus],{i,nco},{j,nco}],
	Table[yHsPy[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yProCon[i,j,nco,con,eeta,taus],{i,nco},{j,nco}]}

Example 3. A binary mixture with
 nco=2;
 con={2/3,1/3};
ratio={1,1/2};
eeta=0.4;
taus={{1,10},{10,1}};

eeta=0.4;
taus={{1,10},{10,1}};

{chim1HsBmcsl[2,eeta],1/chiShsPy[nco,eeta,taus],chim1HsPy[nco,eeta],
  1/chiPro[nco,eeta,taus]}

{Table[yHsBmcsl[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yShsPy[i,j,nco,con,eeta,taus],{i,nco},{j,nco}],
	Table[yHsPy[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yProCon[i,j,nco,con,eeta,taus],{i,nco},{j,nco}]}//TableForm

Example 4. A one-component mixture with eeta=0.49 and tau=1. We evaluate this 
example as a binary mixture with
nco=2;
con={1/2,1/2};
ratio={1,1}; 
eeta=0.49;
taus={{1,1},{1,1}};

eeta=0.49;
taus={{1,1},{1,1}};

{chim1HsBmcsl[nco,eeta],1/chiShsPy[nco,eeta,taus],chim1HsPy[nco,eeta],
  1/chiPro[nco,eeta,taus]}

{Table[yHsBmcsl[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yShsPy[i,j,nco,con,eeta,taus],{i,nco},{j,nco}],
	Table[yHsPy[i,j,nco,eeta],{i,nco},{j,nco}],
	Table[yProCon[i,j,nco,con,eeta,taus],{i,nco},{j,nco}]}