The reaction field factors appear in the generalization of the Kirkwood model of solvation for a molecule imbedded in a cavity surrounded by a dielectric continuum. A general algorithm is proposed to compute these factors for a distributed multipole analysis of the charge distribution of the solute.