I guess I’m not being clear, a large part due to not understanding whats going on. I’m looking at Wouter de Haan’s slides: http://www.wouterdenhaan.com/numerical/slidesbayesian.pdf on slide 39 it mentions, a_{k,k}, a_{k,z} all these are functions of deep/structural parameters \Psi. Also in Solution and Estimation of DSGE models (version: http://www.hoover.org/sites/default/files/fernandez-rubio-shorfheide-hbofmacro_solveestimate_v1.pdf), on printed page number 70 and 76 in acrobat, it lists the matrices as a function of the deep parameters. I want dynare to output those matrices as shown on page 70 as a function of those (deep/structural) parameters, but preferably up to second or third order. As you probably assumed my familiarity in this area and intelligence is higher than it actually is, please keep it dumb. If I’m still asking the same thing and you would just be repeating yourself please feel free to ignore.