diff --git a/matlab/dynare_solve.m b/matlab/dynare_solve.m
index 1c774421821f81faec5ab8bce521b30ed1c4130b..4d36ed4409900f079b0db812ed89e2d71a5e683a 100644
--- a/matlab/dynare_solve.m
+++ b/matlab/dynare_solve.m
@@ -3,6 +3,10 @@
 function [x,check] = dynare_solve(func,x,varargin)
   global gstep_ options_ debug_
   
+
+% unfinished
+  jacobian_flag = 0;   
+
   options_ = set_default_option(options_,'solve_algo',2);
   check = 0;
   func = str2func(func);
@@ -26,7 +30,7 @@ function [x,check] = dynare_solve(func,x,varargin)
 
   if options_.solve_algo == 1
     nn = size(x,1) ;
-    [x,check]=solve1(func,x,1:nn,1:nn,varargin{:});
+    [x,check]=solve1(func,x,1:nn,1:nn,jacobian_flag,varargin{:});
   elseif options_.solve_algo == 2
     nn = size(x,1) ;
     %    tolf = eps^(2/3) ;
@@ -82,4 +86,4 @@ function [x,check] = dynare_solve(func,x,varargin)
 
   % 08/28/03 MJ add a final call to solve1 for solve_algo == 1 in case
   %             initvals generates 'false' zeros in the Jacobian
-  
\ No newline at end of file
+