In this paper, a numerical simulation method for non-Markovian master equations is described.  This method is a generalization of the Markovian Monte Carlo wavefunction method into the non-Markovian regime.  As a simple example, the Janes-Cummings model i