diff --git a/matlab/dynare_solve.m b/matlab/dynare_solve.m
index 41e72c4ae..d3bf68b29 100644
--- a/matlab/dynare_solve.m
+++ b/matlab/dynare_solve.m
@@ -266,6 +266,7 @@ elseif options.solve_algo == 3
     else
         [x,info] = csolve(func,x,[],1e-6,500,varargin{:});
     end
+    [fvec, fjac] = feval(func, x, varargin{:});
 elseif options.solve_algo == 10
     % LMMCP
     olmmcp = options.lmmcp;