diff --git a/matlab/dyn_first_order_solver.m b/matlab/dyn_first_order_solver.m
index 54f5566351340550f32ecbdc5cbf5384f121eac2..cb5e407a833c14cf36940598f853feb3b2d819eb 100644
--- a/matlab/dyn_first_order_solver.m
+++ b/matlab/dyn_first_order_solver.m
@@ -274,7 +274,8 @@ else
     opts.UT = true;
     opts.TRANSA = true;
     hx1 = linsolve(tt(indx_stable_root, indx_stable_root),Z11,opts)';
-    opts.TRANSA = false; % needed by Octave 4.0.0
+    opts.UT = false;      % needed by Octave 4.0.0
+    opts.TRANSA = false;  % needed by Octave 4.0.0
     hx2 = linsolve(Z11,ss(indx_stable_root, indx_stable_root)',opts)';
     hx =  hx1*hx2;
     ghx = [hx(k1,:); gx(k2(nboth+1:end),:)];