Solving Jermann Quadrini 2012 in the correct way, with minimal modification to model