This paper develops a fast method of computing arbitrary order perturbation approximations to bond prices in DSGE models. The procedure is implemented to third order where it can shorten the approximation process by more than 100 times. In a consumption-based endowment model with habits, it is...