shock_matrix = zeros(options_.irf,M_.exo_nbr); %create shock matrix with number of time periods in columns
shock_matrix = zeros(options_.irf,M_.exo_nbr);
shock_matrix(1,strmatch('e0',M_.exo_names,'exact')) = 1;
shock_matrix(2,strmatch('e0',M_.exo_names,'exact')) = 1;
shock_matrix(3,strmatch('e0',M_.exo_names,'exact')) = 1;