This article presents a novel computational approach to solving models with both uninsurable idiosyncratic and aggregate risk that uses projection methods, simulation and perturbation. The approach is shown to be both as efficient and as accurate as existing methods on a model based on Krusell and Smith (1998), for which prior solutions exist. The approach has the advantage of extending straightforwardly, and with reasonable computational cost, to models with a greater range of diversity between agents, which is demonstrated by solving both a model with heterogeneity in discount-rates and a lifecycle model with incomplete markets
