diff --git a/matlab/dr_block.m b/matlab/dr_block.m
index c729d6b00..2bca7bbe3 100644
--- a/matlab/dr_block.m
+++ b/matlab/dr_block.m
@@ -51,7 +51,7 @@ function [dr,info,M_,options_,oo_] = dr_block(dr,task,M_,options_,oo_)
 % along with Dynare.  If not, see <http://www.gnu.org/licenses/>.
 
 info = 0;
-verbose = 0;
+verbose = options_.verbosity;
 if options_.order > 1
     error('2nd and 3rd order approximation not implemented with block option')
 end
diff --git a/matlab/dynare_estimation_1.m b/matlab/dynare_estimation_1.m
index eac92a3e7..158558472 100644
--- a/matlab/dynare_estimation_1.m
+++ b/matlab/dynare_estimation_1.m
@@ -287,7 +287,6 @@ if ~isequal(options_.mode_compute,0) && ~options_.mh_posterior_mode_estimation
         H0 = 1e-4*eye(nx);
         crit = 1e-7;
         nit = 1000;
-        verbose = 2;
         numgrad = options_.gradient_method;
         epsilon = options_.gradient_epsilon;
         % Change some options.
diff --git a/matlab/ep_residuals.m b/matlab/ep_residuals.m
index 66593ae50..56b554de8 100644
--- a/matlab/ep_residuals.m
+++ b/matlab/ep_residuals.m
@@ -34,7 +34,7 @@ function r = ep_residuals(x, y, ix, iy, steadystate, dr, maximum_lag, endo_nbr)
 % You should have received a copy of the GNU General Public License
 % along with Dynare.  If not, see <http://www.gnu.org/licenses/>.
 
-global oo_
+global oo_ options_
 
 persistent k1 k2 weight
 
@@ -45,7 +45,7 @@ if isempty(k1)
     weight = 0.0;
 end
 
-verbose = 0;
+verbose = options_.ep.verbosity;
 
 % Copy the shocks in exo_simul.
 oo_.exo_simul(maximum_lag+1,ix) = exp(transpose(x));
diff --git a/matlab/global_initialization.m b/matlab/global_initialization.m
index 0cdd0f123..c14443ec8 100644
--- a/matlab/global_initialization.m
+++ b/matlab/global_initialization.m
@@ -626,7 +626,6 @@ options_.endogenous_prior_restrictions.moment={};
 % OSR Optimal Simple Rules
 options_.osr.tolf=1e-7;
 options_.osr.maxit=1000;
-options_.osr.verbose=2;
 
 % use GPU
 options_.gpu = 0;
diff --git a/matlab/maximize_prior_density.m b/matlab/maximize_prior_density.m
index 09cb3d774..01e9ffc8d 100644
--- a/matlab/maximize_prior_density.m
+++ b/matlab/maximize_prior_density.m
@@ -36,7 +36,6 @@ number_of_estimated_parameters = length(iparams);
 H0 = 1e-4*eye(number_of_estimated_parameters);
 crit = 1e-7;
 nit = 1000;
-verbose = 2;
 gradient_method = 2;
 gradient_epsilon = 1e-6;
 
diff --git a/matlab/osr1.m b/matlab/osr1.m
index a8d393148..6f2b08daf 100644
--- a/matlab/osr1.m
+++ b/matlab/osr1.m
@@ -76,7 +76,6 @@ inv_order_var = oo_.dr.inv_order_var;
 H0 = 1e-4*eye(np);
 crit=options_.osr.tolf;
 nit=options_.osr.maxit;
-verbose=options_.osr.verbose;
 
 %% do initial checks
 [loss,vx,info,exit_flag]=osr_obj(t0,i_params,inv_order_var(i_var),weights(i_var,i_var));
diff --git a/matlab/particle/gaussian_filter.m b/matlab/particle/gaussian_filter.m
index d82fbc7cb..da3840c74 100644
--- a/matlab/particle/gaussian_filter.m
+++ b/matlab/particle/gaussian_filter.m
@@ -55,8 +55,6 @@ persistent init_flag mf0 mf1
 persistent nodes2 weights2 weights_c2 number_of_particles
 persistent sample_size number_of_state_variables number_of_observed_variables
 
-verbose = 1;
-
 % Set default
 if isempty(start)
     start = 1;