diff --git a/matlab/dsge_likelihood.m b/matlab/dsge_likelihood.m
index 21e17e57d3f67e8d4570124edd1ec4b0380c3e5e..fa9996f87a4a4945b318341ba7187ecb6e28b869 100644
--- a/matlab/dsge_likelihood.m
+++ b/matlab/dsge_likelihood.m
@@ -259,7 +259,7 @@ end
 [T,R,SteadyState,info,Model,DynareOptions,DynareResults] = dynare_resolve(Model,DynareOptions,DynareResults,'restrict');
 
 % Return, with endogenous penalty when possible, if dynare_resolve issues an error code (defined in resol).
-if info(1) == 1 || info(1) == 2 || info(1) == 5 || info(1) == 7 || info(1) == 8 || ...
+if info(1) == 1 || info(1) == 2 || info(1) == 5 || info(1) == 7 || info(1) == 8 || info(1) == 9 || ...
             info(1) == 22 || info(1) == 24 || info(1) == 19 || info(1) == 25 || info(1) == 10
     fval = objective_function_penalty_base+1;
     info = info(1);
diff --git a/matlab/non_linear_dsge_likelihood.m b/matlab/non_linear_dsge_likelihood.m
index 43bf42575259828bd9408c9dd660f6c62a44869a..66cc31570b75e0bd8984d9964464b8198e4b7c1f 100644
--- a/matlab/non_linear_dsge_likelihood.m
+++ b/matlab/non_linear_dsge_likelihood.m
@@ -214,7 +214,7 @@ end
 % Linearize the model around the deterministic sdteadystate and extract the matrices of the state equation (T and R).
 [T,R,SteadyState,info,Model,DynareOptions,DynareResults] = dynare_resolve(Model,DynareOptions,DynareResults,'restrict');
 
-if info(1) == 1 || info(1) == 2 || info(1) == 5 || info(1) == 25 || info(1) == 10 || info(1) == 7
+if info(1) == 1 || info(1) == 2 || info(1) == 5 || info(1) == 25 || info(1) == 10 || info(1) == 7 || info(1) == 9
     fval = objective_function_penalty_base+1;
     exit_flag = 0;
     return