Before revision:  Finally, we extend our modified Next Reaction Method to systems with delays and compare its efficiency with those of the recent algorithms designed for systems with delays.
After revision:  Finally, we extend our modified Next Reaction Method to systems with time dependent propensities as well as to systems with delays.