%shock.m
% We use Option 2 in relaxsetting, line 39
% 
% if t==-1        % At time -1 the parameters should have the values for the old steady state
%     alpha=aalpha_0;
% else
%     alpha=aalpha_f;
% end
