Computation of Equilibria in Heterogeneous Agent Models