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?