function Y = sumY(T);

% This function returns the vector of differences (1-sumY) at given values of 
% temperature (T), pressure (P) and liquid composition (X)

global X P NC R W V A RP Q

Y=UNIQUAC(T);
Y=1-sum(Y);