diff --git a/matlab/dr1.m b/matlab/dr1.m
index 34d002cfbdaed27a96a2d1af29268e18df137c74..81dae7de1f8d984f4a76b26f40d51f09cf9ab229 100644
--- a/matlab/dr1.m
+++ b/matlab/dr1.m
@@ -425,12 +425,15 @@ for i=1:ykmax_-1
   k0 = k1;
   offset = offset + n1;
 end
-%C = kron(hx,hx);
-C = hx;
 D = [rhs; zeros(n-endo_nbr,size(rhs,2))];
-%x0 = sylvester3(A,B,C,D);
-%dr.ghxx = sylvester3a(x0,A,B,C,D);
-dr.ghxx = gensylv(2,A,B,C,D);
+if exist('gensylv')
+  C = hx;
+  dr.ghxx = gensylv(2,A,B,C,D);
+else
+  C = kron(hx,hx); 
+  x0 = sylvester3(A,B,C,D);
+  dr.ghxx = sylvester3a(x0,A,B,C,D);
+end
 
 %ghxu
 %rhs