th_autocovariances.m creates Gamma_y{nar+2} (variance decomposition)
even if there is a single shock. In that case variance decomposition
is one. It is better if Gamma_y has always the expected
size. disp_th_moments.m only displays variance decomposition if there
is more than one shock.