% master.m

% Set path for KWRE package
setpath

vlstrbc

% calculate RE solution

[M,PI,G,ylist,xlist,lpd,elist]=RESOLKW('linrbc','drvrbc');


% Set shock vector
mm=ROWS(M);
shock=zeros(1,mm);

sc=input('Choose 1 for A shock and 2 for G shock')

if (sc==1)
% A shock
shock(2)=1;
else
% Gshock
shock(3)=1;
end

% Compute impulse response

nir=20;

IR=IMPKW(M,PI,G,ylist,xlist,lpd,shock,nir);

% plot impulse responses;
plotrbc
Error using ==> input
Cannot call 'input' from evalc.