diff --git a/matlab/check.m b/matlab/check.m
index 0a680ed431fa389c88fe52797152764efbe4c4b6..e6ea8aec3e6c229c0e5465fab871fdb357bcb607 100644
--- a/matlab/check.m
+++ b/matlab/check.m
@@ -49,11 +49,7 @@ global M_ options_ oo_
     oo_.exo_simul = tempex;
     
     eigenvalues_ = dr.eigval;
-    if(options_.model_mode == 1)
-        nyf = dr.nyf;  
-    else
-        nyf = nnz(dr.kstate(:,2)>M_.maximum_endo_lag+1);
-    end;
+    nyf = nnz(dr.kstate(:,2)>M_.maximum_endo_lag+1);
     [m_lambda,i]=sort(abs(eigenvalues_));
     n_explod = nnz(abs(eigenvalues_) > options_.qz_criterium);
      
diff --git a/matlab/global_initialization.m b/matlab/global_initialization.m
index 1a6be4ad360321a828ffcb93bb8d3899e0694123..63a59f89c248dae492b6fa32f89bb1211fe51b77 100644
--- a/matlab/global_initialization.m
+++ b/matlab/global_initialization.m
@@ -113,7 +113,6 @@ function global_initialization()
   options_.replic = 50;
   options_.drop = 100;
   options_.simul_algo = 0;
-  options_.model_mode = 0;
   % if mjdgges.dll (or .mexw32 or ....) doesn't exist, matlab/qz is added to the path. 
   % There exists now qz/mjdgges.m that contains the calls to the old Sims code 
   % Hence, if mjdgges.m is visible exist(...)==2, 
diff --git a/matlab/model_info.m b/matlab/model_info.m
index 7d6b986699234c53412ab19554212c8675b4146b..3f1d713dcacc5644872649a80c4354ddbe9f17e3 100644
--- a/matlab/model_info.m
+++ b/matlab/model_info.m
@@ -67,7 +67,7 @@ function model_info;
        fprintf('\n\n');
    end;
  else
-   fprintf('There is no block decomposition of the model.\nUse ''sparse'' or ''sparse_dll'' model''s option.\n');
+   fprintf('There is no block decomposition of the model.\nUse ''block'' model''s option.\n');
  end;
  
  
diff --git a/matlab/resid.m b/matlab/resid.m
index 5287b8812f0cca496a5ef8d3255472c942c386bc..d8b10b312ae3a51b348bf71b08d6c200db89d7b1 100644
--- a/matlab/resid.m
+++ b/matlab/resid.m
@@ -49,21 +49,10 @@ function resid(period)
   y =oo_.endo_simul(:);
   z = zeros(n,period);
   fh = str2func([M_.fname '_dynamic']);
-  if(options_.model_mode == 1 || options_.model_mode == 3)
-      addpath(M_.fname);
-  end;
   for it_=M_.maximum_lag+1:period+M_.maximum_lag
-    %if(options_.model_mode == 1 || options_.model_mode == 3)
-    %    z(:,it_-M_.maximum_lag) = feval(fh,oo_.endo_simul',oo_.exo_simul, M_.params, it_);
-    %else
         z(:,it_-M_.maximum_lag) = feval(fh,y(iyr0),oo_.exo_simul, M_.params, it_);
         iyr0 = iyr0 + n;
-    %end;
   end
-  if(options_.model_mode == 1 || options_.model_mode == 3)
-      rmpath(M_.fname);
-  end;
-  % disp([[1:period]' z']); 
 
   for i = 1:4
     disp(' ')
diff --git a/matlab/resol.m b/matlab/resol.m
index 944e47bfb7e0840122740b4b53a49a4a0b2948c4..6b5b601ec99a0593450aecc4f8eeb2f60bbcbeee 100644
--- a/matlab/resol.m
+++ b/matlab/resol.m
@@ -112,11 +112,7 @@ if ~isreal(dr.ys)
 end
 
 dr.fbias = zeros(M_.endo_nbr,1);
-if(options_.model_mode==1)
-    [dr,info,M_,options_,oo_] = dr1_sparse(dr,check_flag,M_,options_,oo_);
-else
-    [dr,info,M_,options_,oo_] = dr1(dr,check_flag,M_,options_,oo_);
-end
+[dr,info,M_,options_,oo_] = dr1(dr,check_flag,M_,options_,oo_);
 if info(1)
   return
 end