We show how the quantum trajectory method for describing the dissipative evolution of conditioned states can be implemented in a numerically efficient higher-order unraveling that goes beyond the Euler method usually employed in such simulations. We demonstrate its applicability through an analysis of spontaneous emission and find that improvements by two orders of magnitude are possible for the specific example addressed.