gen_impulse_reg function

Generate impulse regressors.