diff --git a/matlab/dynare_solve.m b/matlab/dynare_solve.m
index d048a9936ddf8c278a7c4ff08e1c946429e4c921..401946da0c065ab2549ccc4175f5fb1e6b39677e 100644
--- a/matlab/dynare_solve.m
+++ b/matlab/dynare_solve.m
@@ -57,8 +57,8 @@ if options_.solve_algo == 0
         % Under Octave, use a wrapper, since fsolve() does not have a 4th arg
         func2 = str2func(func);
         func = @(x) func2(x, varargin{:});
-        % The Octave version of fsolve do not converge when it starts from the solution
-        [fvec,fjac] = feval(func,x,varargin{:});
+        % The Octave version of fsolve does not converge when it starts from the solution
+        fvec = feval(func,x);
         if max(abs(fvec)) >= options_.solve_tolf
             [x,fval,exitval,output] = fsolve(func,x,options);
         else