Question: Solving Lyapunov equations

When executing "P, k := LyapunovSolve(A, C)" command, I got the following warning: 

"Warning, Matrices have common or very close eigenvalues; perturbed values were used to solve the equation"

Although Maple can give the solution, P and k don't satisfy the equation: 

A.P+P.A^%T-C = 0        (1)

Are there any other ways to solve the Lyapunov equation (1)?

Thank you for your help.

LyapunovSolve.mw

Please Wait...