This paper develops a method to compute second-order perturbations of discretetime heterogeneous agent models. It addresses the three main tasks to make secondorder approximations tractable: state reduction, generating sufficient smoothness, and fast computation of the quadratic terms in the...