A numerical technique is introduced that reduces exponentially the time required for Monte Carlo simulations of nonequilibrium systems. Results for the quasistationary probability distribution in two model systems are compared with the asymptotically exact theory in the limit of extremely small noise intensity. Singularities of the nonequilibrium distributions are revealed by the simulations.