diff --git a/matlab/metropolis.m b/matlab/metropolis.m
index 9497b3ac2b49c206bd21ee8ecca4b433f34ebbe5..18da87a73c95824aaf744e761ff69b5d20d2382b 100644
--- a/matlab/metropolis.m
+++ b/matlab/metropolis.m
@@ -2542,8 +2542,7 @@ function metropolis(xparam1,vv,gend,data,rawdata,mh_bounds)
 	tit(lgx_orig_ord_,:) = lgx_;
 	for i = 1:exo_nbr
 	  if SS(i,i) > 1e-13
-	    y=irf(dr_,cs(lgx_orig_ord_,i),nirfs,options_.drop, ...
-		  options_.replic, options_.order);
+	    y=irf(dr_,cs(lgx_orig_ord_,i),nirfs,options_.drop,options_.replic, options_.order);
 	    if options_.relative_irf
 	      y = 100*y/cs(i,i); 
 	    end
@@ -2587,7 +2586,7 @@ function metropolis(xparam1,vv,gend,data,rawdata,mh_bounds)
 	  end
 	  % Get rotation
 	  if dsge_prior_weight > 0
-	    Atheta(dr_.order_var,:) = dr_.ghu(:,lgx_orig_ord_)*sqrt(Sigma_e_(lgx_orig_ord_,lgx_orig_ord_));%%%
+	    Atheta(dr_.order_var,lgx_orig_ord_) = dr_.ghu*sqrt(Sigma_e_);
 	    A0 = Atheta(bayestopt_.mfys,:);
 	    [OMEGAstar,SIGMAtr] = qr2(A0');
 	  end
@@ -2709,7 +2708,7 @@ function metropolis(xparam1,vv,gend,data,rawdata,mh_bounds)
 	  end
 	  % Get rotation
 	  if dsge_prior_weight > 0
-	    Atheta(dr_.order_var,:) = dr_.ghu(:,lgx_orig_ord_)*sqrt(Sigma_e_(lgx_orig_ord_,lgx_orig_ord_));%%%
+	    Atheta(dr_.order_var,lgx_orig_ord_) = dr_.ghu*sqrt(Sigma_e_);%%%
 	    A0 = Atheta(bayestopt_.mfys,:);
 	    [OMEGAstar,SIGMAtr] = qr2(A0');
 	  end
diff --git a/matlab/simult_.m b/matlab/simult_.m
index 5e0cc603ebf73859e75120b2dedb6b32dcc31365..4434c5146e92aa1e6b131bdbe07e79f86fb957e8 100644
--- a/matlab/simult_.m
+++ b/matlab/simult_.m
@@ -1,6 +1,5 @@
 % Copyright (C) 2001 Michel Juillard
 %
-
 function y_ = simult_(y0,dr,ex_,iorder)
 global endo_nbr ykmin_ xkmin_ it_ options_ iy_
   iter = size(ex_,1)-xkmin_;