diff --git a/matlab/+pac/+estimate/nls.m b/matlab/+pac/+estimate/nls.m
index e59181246..ba52b5a2b 100644
--- a/matlab/+pac/+estimate/nls.m
+++ b/matlab/+pac/+estimate/nls.m
@@ -280,23 +280,23 @@ if nargin<5
 end
 
 if is_gauss_newton
-    [params1, SSR] = gauss_newton(resfun, params0);
+    [params1, SSR, exitflag] = gauss_newton(resfun, params0);
 elseif is_lsqnonlin
     if ismember('levenberg-marquardt', varargin)
         % Levenberg Marquardt does not handle boundary constraints.
-        [params1, SSR] = lsqnonlin(resfun, params0, [], [], optimset(varargin{:}));
+        [params1, SSR, ~, exitflag] = lsqnonlin(resfun, params0, [], [], optimset(varargin{:}));
     else
-        [params1, SSR] = lsqnonlin(resfun, params0, bounds(:,1), bounds(:,2), optimset(varargin{:}));
+        [params1, SSR, ~, exitflag] = lsqnonlin(resfun, params0, bounds(:,1), bounds(:,2), optimset(varargin{:}));
     end
 else
     % Estimate the parameters by minimizing the sum of squared residuals.
-    [params1, SSR] = dynare_minimize_objective(ssrfun, params0, ...
-        minalgo, ...
-        options_, ...
-        bounds, ...
-        parameter_names, ...
-        [], ...
-        []);
+    [params1, SSR, exitflag] = dynare_minimize_objective(ssrfun, params0, ...
+                                                      minalgo, ...
+                                                      options_, ...
+                                                      bounds, ...
+                                                      parameter_names, ...
+                                                      [], ...
+                                                      []);
 end
 
 % Revert local modifications to the options.