Development of a theoretical treatment of multiple coherent light scattering in an ultracold atomic gas is reported. Specific application is made to coherent backscattering of a weak-radiation field from realistically modeled samples of ultracold atomic 85Rb. Comprehensive Monte Carlo simulations of the spatial, spectral, and polarization dependence of the backscattering line shape are made and compared with available experimental results
