diff --git a/matlab/PosteriorIRF.m b/matlab/PosteriorIRF.m
index fd37ff48dd62133653f33b53c550b9046663a846..8c094a24f72ecde2dc9bca2abfd882953b454aaa 100644
--- a/matlab/PosteriorIRF.m
+++ b/matlab/PosteriorIRF.m
@@ -398,7 +398,7 @@ if options_.TeX && any(strcmp('eps',cellstr(options_.graph_format)))
         figunumber = 0;
 
         for jj=1:nvar
-            if max(abs(MeanIRF(:,jj,ii))) > options_.impulse_responses.plot_threshold
+            if max(abs(MeanIRF(:,jj,ii))) >= options_.impulse_responses.plot_threshold
                 subplotnum = subplotnum+1;
                 
                 if subplotnum == 1 
diff --git a/matlab/PosteriorIRF_core2.m b/matlab/PosteriorIRF_core2.m
index 1fef8cdff4472e6795ec5b603c7214f4d2aa751e..51f4bc08ae927fc64c25ddc5517d3b04ba305171 100644
--- a/matlab/PosteriorIRF_core2.m
+++ b/matlab/PosteriorIRF_core2.m
@@ -100,7 +100,7 @@ for i=fpar:npar,
     figunumber = 0;
 
     for j=1:nvar
-        if max(abs(MeanIRF(:,j,i))) > options_.impulse_responses.plot_threshold
+        if max(abs(MeanIRF(:,j,i))) >= options_.impulse_responses.plot_threshold
             subplotnum = subplotnum+1;
             if subplotnum == 1 && options_.relative_irf
                 hh = dyn_figure(options_,'Name',['Relative response to orthogonalized shock to ' tit(i,:)]);
diff --git a/matlab/stoch_simul.m b/matlab/stoch_simul.m
index 03f3acc47cb6d2edbdb4e716232e466315b546fc..d6b6ceadfb01fada54f6ad96fa2a3794c02d1f46 100644
--- a/matlab/stoch_simul.m
+++ b/matlab/stoch_simul.m
@@ -215,7 +215,7 @@ if options_.irf
                          y(i_var(j),:)');
                 eval(['oo_.irfs.' deblank(M_.endo_names(i_var(j),:)) '_' ...
                       deblank(M_.exo_names(i,:)) ' = y(i_var(j),:);']);
-                if max(abs(y(i_var(j),:))) > options_.impulse_responses.plot_threshold
+                if max(abs(y(i_var(j),:))) >= options_.impulse_responses.plot_threshold
                     irfs  = cat(1,irfs,y(i_var(j),:));
                     if isempty(mylist)
                         mylist = deblank(var_list(j,:));