# Resol.m and dr1.m

**URL:** <https://forum.dynare.org/t/resol-m-and-dr1-m/164>\
**Category:** Dynare help (legacy posts)\
**Created:** [31 August 2006 16:19 UTC](https://forum.dynare.org/t/resol-m-and-dr1-m/164 "2006-08-31T16:19:54Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![krishnarao](https://forum.dynare.org/letter_avatar_proxy/v4/letter/k/ecc23a/32.png) [@krishnarao](https://forum.dynare.org/u/krishnarao)\
**Post date:** [31 August 2006 16:19 UTC](https://forum.dynare.org/t/resol-m-and-dr1-m/164/1 "2006-08-31T16:19:54Z")

</div>

My colleagues and I have a great interest in the details of the metropolis, DsgeSmoother, and DiffuseKalman routines. Though we have successfully documented many aspects of these programs, the purpose and functioning of resol.m still confuses us. At line 1656 of metropolis we see:

“% FIFTH, I update variable dr\_  
dr\_ = resol(ys\_,0);”

At first, we thought that resol.m updates dr\_ given the steady-state output of DsgeSmoother, ys. Now, we know that ys\_ for us is always empty or zero and not equal to ys. How and why do resol.m and dr1.m update dr\_.

Thanks,  
Kieran and Krishna

---

<div class="post-metadata">

**Author:** ![MichelJuillard](https://forum.dynare.org/user_avatar/forum.dynare.org/micheljuillard/32/5318_2.png) [@MichelJuillard](https://forum.dynare.org/u/MichelJuillard)\
**Post date:** [13 September 2006 19:32 UTC](https://forum.dynare.org/t/resol-m-and-dr1-m/164/2 "2006-09-13T19:32:38Z")

</div>

1. resol has two functions:

- finding the steady state corresponding to the value of the parameters
- finding the first order approximate solution of the rational expectation model  
this solution is then saved in dr\_.ghx and dr\_.ghu

Note that the fragment of code that you are quoting is in a part of the routine where we exploit previous Metropolis draws in order to compute the posterior distribution of forecasts of the model

Best

Michel