Continuum scattering by free electrons can be significant in early type stars, while in late type stars Rayleigh scattering by hydrogen atoms or molecules may be important. Computer programs used to construct models of stellar atmospheres generally treat the scattering of the continuum radiation as isotropic and unpolarized, but this scattering has a dipole angular dependence and will produce polarization. We review an accurate method for evaluating the polarization and limb darkening of the radiation from model stellar atmospheres. We use this method to obtain results for: (i) Late type stars, based on the MARCS code models (Gustafsson et al. 2008), and (ii) Early type stars, based on the NLTE code TLUSTY (Lanz and Hubeny 2003). These results are tabulated at http://www.astro.umd.edu/~jph/Stellar_Polarization.html While the net polarization vanishes for an unresolved spherical star, this symmetry is broken by rapid rotation or by the masking of part of the star by a binary companion or during the transit of an exoplanet. We give some numerical results for these last cases.