Dear Musketeers,

I am trying to use simult_ function to simulate model that is solved using order=2 (and order=3). I start the simulation from the steady-state (oo_.dr.ys), input shocks that has value zero.

For order=1, everything works as expected: the deviation from oo_.dr.ys is zero.

For order=2, and order=3, the deviation is not anymore zero (in the attached example of magnitude 10^-3).

What is wrong? The attached code builds on jermann98.mod file using the instructions I found the forum.

Best regards,

Antti

non0irf.mod (3.49 KB)