diff --git a/matlab/dynare_solve.m b/matlab/dynare_solve.m
index 9e27df3eeaf7a389a65ca024268003c002dba3a6..45463f2833340382a178eff68ac0e206fe245cc2 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