Binary fluid mixtures display a wide array of fluid phase behavior ranging from simple vapor-liquid equilibrium diagrams to more complicated topology like azeotropy and heteroazeotropy. As a first step towards understanding the phenomena of equilibrium of binary molecular systems, the properties of binary monoatomic fluids have been studied in this work. An additional simplification made in this work is that of a symmetrical binary system where all similar molecules interact via the hard-core Yukawa (HCY) potential, u(r), while the dissimilar molecules interact via a potential, δu(r) where δ is a scalar parameter. The value of δ here is 0.75, which leads to the dissimilar molecules showing a tendency to dislike one another and the resulting phase diagram is such that the mixing-demixing line intersects the vapor-liquid equilibrium curve away from the vapor-liquid critical point. The effect of the variable potential range parameter of the HCY potential on the topology of the phase diagrams is investigated using grand canonical transition matrix Monte Carlo simulations and the conditions of temperature and pressure at which the system exhibits azeotropy and heteroazeotropy are ascertained. We find that the densities of the mixtures, as predicted by our simulation at equimolar concentrations, are in close agreement with the self-consistent Ornstein-Zernike approximation results of Schöll-Paschinger and co-workers.