addpath C:\dynare\4.6.0\matlab

clear all;
clc;

dynare Themodel_EK.mod 


initial_condition_states = repmat(oo_.dr.ys,1,M_.maximum_lag);
shock_matrix = zeros(options_.irf,M_.exo_nbr); %create shock matrix with number of time periods in columns
% set shocks for pure news 
shock_matrix(1,strmatch('eTauhat',M_.exo_names,'exact')) = 0.03; %set tax shock (use any shock size you want)
shock_matrix(1+4,strmatch('eIhat',M_.exo_names,'exact')) = 0.25; %
y2 = simult_(M_,options_,initial_condition_states,oo_.dr,shock_matrix,2);
IRFmodel025c = y2(:,M_.maximum_lag+1:end)-repmat(oo_.dr.ys,1,options_.irf); %deviation from steady states

save IRFmodel025c;


