diff --git a/matlab/lyapunov_solver.m b/matlab/lyapunov_solver.m
index 98d782cdf30d7ef239bd7968f9cc875b75217ba8..2bafb065ab159cefbfaa8c3a22666f0800ab69c9 100644
--- a/matlab/lyapunov_solver.m
+++ b/matlab/lyapunov_solver.m
@@ -41,7 +41,7 @@ function P=lyapunov_solver(T,R,Q,DynareOptions) % --*-- Unitary tests --*--
 % along with Dynare.  If not, see <http://www.gnu.org/licenses/>.
 
 if DynareOptions.lyapunov_fp == 1
-    P = lyapunov_symm(T,R*Q'*R',DynareOptions.lyapunov_fixed_point_tol,DynareOptions.qz_criterium,DynareOptions.lyapunov_complex_threshold, 3, DynareOptions.debug);
+    P = lyapunov_symm(T,R*Q*R',DynareOptions.lyapunov_fixed_point_tol,DynareOptions.qz_criterium,DynareOptions.lyapunov_complex_threshold, 3, DynareOptions.debug);
 elseif DynareOptions.lyapunov_db == 1
     [P, errorflag] = disclyap_fast(T,R*Q*R',DynareOptions.lyapunov_doubling_tol);
     if errorflag %use Schur-based method
@@ -183,4 +183,4 @@ end
 %$ end
 %$
 %$ T = all(t);
-%@eof:1
\ No newline at end of file
+%@eof:1