In stoch_simul.m IRFs were only printed if they were not flat instead of checking whether the IRF is different from 0. Unit root responses were thus filtered out, e.g. C=0.2 for all t was not displayed.