To solve or simulate the Holstein equation, the Monte Carlo method is commonly employed. An absorption coefficient is calculated for an experiment with a certain opacity, atomic species, Doppler-broadened lineshape, etc., and then a test is made to see whether the photon escapes after n flights through the atomic vapor (see Figure 1 in the reference).