  I have created, for public use, a Python code allowing the simulation of
light curves with any given power spectral density and any probability density
function, following the algorithm described in Emmanoulopoulos et al. 2013. The
simulated products have exactly the same variability and statistical properties
as the observed light curves. The code and its documentation are available at:
  https://github.com/samconnolly/DELightcurveSimulation
  Note that a Mathematica code of the algorithm is given in Emmanoulopoulos et
al. [arXiv:1305.0304]
