diff --git a/matlab/dynare_identification.m b/matlab/dynare_identification.m
index b20574c03..ad026f9d9 100644
--- a/matlab/dynare_identification.m
+++ b/matlab/dynare_identification.m
@@ -271,7 +271,7 @@ if iload <=0,
         else
             params = prior_draw();
         end
-        [~, ideJ, ideH, ideGP, ~ , info] = ...
+        [dum1, ideJ, ideH, ideGP, dum2 , info] = ...
             identification_analysis(params,indx,indexo,options_MC,data_info, prior_exist, name_tex,0);
         if iteration==0,
             MAX_tau   = min(SampleSize,ceil(MaxNumberOfBytes/(size(ideH.siH,1)*nparam)/8));
@@ -423,7 +423,7 @@ if SampleSize > 1,
         jcrit=find(idemoments.ino);
         if length(jcrit)<SampleSize,
             if isempty(jcrit),
-                [~,jmax]=max(idemoments.cond);
+                [dum,jmax]=max(idemoments.cond);
                 fprintf('\n')
                 tittxt = 'Draw with HIGHEST condition number';
                 fprintf('\n')
@@ -436,7 +436,7 @@ if SampleSize > 1,
                 disp_identification(pdraws(jmax,:), idemodel_max, idemoments_max, name);
                 close all,
                 plot_identification(pdraws(jmax,:),idemoments_max,idehess_max,idemodel_max,idelre_max,1,tittxt,name,IdentifDirectoryName,1);
-                [~,jmin]=min(idemoments.cond);
+                [dum,jmin]=min(idemoments.cond);
                 fprintf('\n')
                 tittxt = 'Draw with SMALLEST condition number';
                 fprintf('\n')
diff --git a/matlab/plot_identification.m b/matlab/plot_identification.m
index e7b1083df..d3f370dd5 100644
--- a/matlab/plot_identification.m
+++ b/matlab/plot_identification.m
@@ -273,11 +273,11 @@ else
         eval(['print -dpdf ' IdentifDirectoryName '/' M_.fname '_ident_COND']);
         if options_.nograph, close(gcf); end
         ncut=floor(SampleSize/10*9);
-        [~,is]=sort(idelre.cond);
+        [dum,is]=sort(idelre.cond);
         [proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberLRE', 1, [], IdentifDirectoryName, 0.1);
-        [~,is]=sort(idemodel.cond);
+        [dum,is]=sort(idemodel.cond);
         [proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberModel', 1, [], IdentifDirectoryName, 0.1);
-        [~,is]=sort(idemoments.cond);
+        [dum,is]=sort(idemoments.cond);
         [proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberMoments', 1, [], IdentifDirectoryName, 0.1);
 %         [proba, dproba] = stab_map_1(idemoments.Mco', is(1:ncut), is(ncut+1:end), 'HighestCondNumberMoments_vs_Mco', 1, [], IdentifDirectoryName);
 %         for j=1:nparam,