% Solve model, generate model IRFs
[zdatalinear, zdatapiecewise, zdatass, oobase_, Mbase_] = ...
    solve_one_constraint('chap_3', 'gali_chapter_3_zlb', ...
                         'inot <= -(1/betta - 1)', 'inot > -(1/betta - 1)', ...
                         shockssequence, irfshock, 15, 20);



                          
% unpack the IRFs                          
for i=1:Mbase_.endo_nbr
  eval([deblank(Mbase_.endo_names(i,:)),'_uncdifference=zdatalinear(:,i);']);
  eval([deblank(Mbase_.endo_names(i,:)),'_difference=zdatapiecewise(:,i);']);
  eval([deblank(Mbase_.endo_names(i,:)),'_ss=zdatass(i);']);
end


nparams = size(Mbase_.param_names,1);

for i = 1:nparams
  eval([Mbase_.param_names(i,:),'= Mbase_.params(i);']);
end
