Consecutive and identical shocks

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;