The solution of scattering problems involving low-contrast dielectric objects with three-dimensional arbitrary shapes is considered. Using the traditional forms of the surface integral equations, scattered fields cannot be calculated accurately if the contrast of the object is low. Therefore, we consider the stabilization of the formulations by extracting the nonradiating parts of the equivalent currents. We also investigate various types of stable formulations and show that accuracy can be improved systematically by eliminating the identity terms from the integral-equation kernels. Traditional and stable formulations are compared, not only for small scatterers but also for relatively large problems solved by employing the multilevel fast multipole algorithm. Stable and accurate solutions of dielectric contrasts as low as 10(-4) are demonstrated on problems involving more than 250000 unknowns.