function y=bubal(alpha,x)
% Equilibrium vapor over ideal, constant relative binary mixture
% function y=bubal(alpha,x)
% Argument  gives
%  alpha    realtive volatility
%    x      liquid mol fraction
% Returns y the vapor mol fraction
y=alpha*x/((alpha-1)*x+1);
