diff --git a/matlab/lmmcp/lmmcp.m b/matlab/lmmcp/lmmcp.m
index 2d6189fd0491554745cbd3db3ce280f964ed714e..aa170222c999699d70c5691d6bf08438cf1a0ea5 100644
--- a/matlab/lmmcp/lmmcp.m
+++ b/matlab/lmmcp/lmmcp.m
@@ -392,7 +392,7 @@ while (k < kmax) && (Psix > eps2)
   Fx        = Fxnew;
   Phix      = Phixnew;
   Psix      = Psixnew;
-  [~,DFx]   = feval(FUN,x,varargin{:});
+  [junk,DFx]   = feval(FUN,x,varargin{:});
   DPhix     = DPhi(x,Fx,DFx,lb,ub,lambda1,lambda2,n,Indexset);
   DPsix     = DPhix'*Phix;
   normDPsix = norm(DPsix);
diff --git a/matlab/perfect-foresight-models/private/simulation_core.m b/matlab/perfect-foresight-models/private/simulation_core.m
index 6e86671df6c1cd2cbf295c6c97e621cc6007c837..ecf5c8964dd37d0cb493e68153f6caaf255b7a8c 100644
--- a/matlab/perfect-foresight-models/private/simulation_core.m
+++ b/matlab/perfect-foresight-models/private/simulation_core.m
@@ -75,9 +75,9 @@ else
                 yT = y(:,periods+2);
                 z = y(:,2:periods+1);
                 illi = M_.lead_lag_incidence';
-                [i_cols,~,i_cols_j] = find(illi(:));
+                [i_cols,junk,i_cols_j] = find(illi(:));
                 illi = illi(:,2:3);
-                [i_cols_J1,~,i_cols_1] = find(illi(:));
+                [i_cols_J1,junk,i_cols_1] = find(illi(:));
                 i_cols_T = nonzeros(M_.lead_lag_incidence(1:2,:)');
                 [y,info] = dynare_solve(@perfect_foresight_problem,z(:),1, ...
                                  str2func([M_.fname '_dynamic']),y0,yT, ...
@@ -109,9 +109,9 @@ if nargout>1
     yy  = oo_.endo_simul(:,2:options_.periods+1);
     if ~exist('illi')
         illi = M_.lead_lag_incidence';
-        [i_cols,~,i_cols_j] = find(illi(:));
+        [i_cols,junk,i_cols_j] = find(illi(:));
         illi = illi(:,2:3);
-        [i_cols_J1,~,i_cols_1] = find(illi(:));
+        [i_cols_J1,junk,i_cols_1] = find(illi(:));
         i_cols_T = nonzeros(M_.lead_lag_incidence(1:2,:)');
     end
     if options_.block && ~options_.bytecode