Dear prof. Pfeifer,
May I ask something in your code implementation for Gali (2015) chapter 3? In there, you have the following definition of natural output:
[name='Definition natural output, eq. (20)']
y_nat=psi_n_ya*a;
Is seems it should be y_nat=psi_n_ya*a + psi_y;
. In the book, we have y_t^n = \psi_{ya} a_t + \psi_y, \;\;\psi_y>0. Or \psi_y is very small to matter?