Harmonize plotting of IRFS

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.
(cherry picked from commit c9ebe21fafad941aa53dc3da42a23525f2d8b77e)