We report the phase coexistence properties of polarizable Stockmayer fluids of reduced permanent dipoles m(0)* = 1.0 and 2.0 and reduced polarizabilities alpha* = 0.00, 0.03, and 0.06, calculated by. series of grand canonical Monte Carlo simulations with the histogram reweighting method. In the histogram reweighting method, the distributions of density and energy calculated in Grand Canonical Monte Carlo simulations are stored in histograms and analyzed to construct the grand canonical partition function of the system. All thermodynamic properties are calculated from the grand partition function. The results are compared with Wertheim’s renormalization perturbation theory. Deviations between theory and simulation results for the coexistence envelope are near 2% for the lower dipole moment and 10% for the higher dipole moment we studied.