Fifth-order perturbation solution to DSGE models

This paper derives a fifth-order perturbation solution to DSGE models. The paper develops a new notation that reduces the notational complexity of high-order solutions and yields a faster code. The new notation consists of new matrix forms of high-order multivariate chain rules and a new representation of the model as a function of one vector variable. The algorithm that implements the new notation is between 3 and 55 times faster than Dynare++, depending on model size and solution order.