We present an algorithm and software routines for computing nth order Taylor series approximate solutions to dynamic, discrete-time rational expectations models around a nonstochastic steady state.