The statistics of photon emission from a single molecule under continuous wave excitation are considered. In particular, we study stochastic model systems where photon emission rates evolve in time with non-Markovian dynamics. Our calculations are based on the recently introduced generalized optical Bloch equation (GBE) formalism, but with numerical complications beyond those seen in previous Markovian stochastic models. A spectral representation is introduced to facilitate the numerical solution of the GBE equations for these more challenging cases. (C) 2004 American Institute of Physics.