diff --git a/matlab/dynare_solve.m b/matlab/dynare_solve.m
index 366bfa3a45235c97dba303d861f43f92c6c673b4..4dc408eb6ea8eed0472c3b6c4906f2a4057bb072 100644
--- a/matlab/dynare_solve.m
+++ b/matlab/dynare_solve.m
@@ -264,6 +264,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;