diff --git a/matlab/metropolis_hastings_initialization.m b/matlab/metropolis_hastings_initialization.m
index 27e8f4c7454504dfe708f944fd88995e352aa9fe..f31111d43712e16ed6088383b290ab2271869fa9 100644
--- a/matlab/metropolis_hastings_initialization.m
+++ b/matlab/metropolis_hastings_initialization.m
@@ -113,6 +113,7 @@ if ~options_.load_mh_file && ~options_.mh_recover
     fprintf(fidlog,' \n');
     % Find initial values for the nblck chains...
     if nblck > 1% Case 1: multiple chains
+        set_dynare_seed('default');
         fprintf(fidlog,['  Initial values of the parameters:\n']);
         disp('Estimation::mcmc: Searching for initial values...')
         ix2 = zeros(nblck,npar);