A computational procedure is described to predict homogeneous azeotropic loci of binary mixtures with cubic equations of state and gas-liquid P-T-x critical information. The performance of the method is demonstrated on three binary systems exhibiting either negative or positive azeotropic behavior which persists up to the critical region. The method of Heidemann and Khalil was used for correlating the critical data while that of Wang and Whiting was used to predict the azeotropic lines together with a stability test applied on each estimated point. Good agreement between calculated and measured azeotropic data was obtained for the systems studied.