diff --git a/doc/manual/source/the-model-file.rst b/doc/manual/source/the-model-file.rst
index 5bd4bb2e6e9912fc04b7c3f1d882ee1aa9ad3fd8..6c1a320046bd87a62cedc6beceb0ad020d7073bd 100644
--- a/doc/manual/source/the-model-file.rst
+++ b/doc/manual/source/the-model-file.rst
@@ -6102,7 +6102,7 @@ block decomposition of the model (see :opt:`block`).
                 routine (available under MATLAB if the Optimization Toolbox is
                 installed; available under Octave if the `optim
                 <https://octave.sourceforge.io/optim/>`__ package from
-                Octave-Forge is installed).
+                Octave-Forge is installed). Only supported for ``method_of_moments``. 
 
            ``101``
 
diff --git a/matlab/optimization/dynare_minimize_objective.m b/matlab/optimization/dynare_minimize_objective.m
index 06ecb9be9c59e8e3ba1eff62f1fa399c834d6995..62f7019ae7017264386a189a85031aca008b3654 100644
--- a/matlab/optimization/dynare_minimize_objective.m
+++ b/matlab/optimization/dynare_minimize_objective.m
@@ -542,6 +542,9 @@ switch minimizer_algorithm
     opt_par_values = opt_par_values(:);
   case 13
     % Matlab's lsqnonlin (Optimization toolbox needed).
+    if ~isfield(options_,'mom')
+        error('Option mode_compute=13 is available only for method_of_moments estimation.')
+    end
     if isoctave && ~user_has_octave_forge_package('optim')
         error('Option mode_compute=13 requires the optim package')
     elseif ~isoctave && ~user_has_matlab_license('optimization_toolbox')