Problem: Check the stability (in Lyapunov sense) for the state coefficient matrix
We solve the Lyapunov equation using lyap() function in MATLAB and LyapunovSolve[] function in Mathematica, and then check that the solution is positive definite (i.e. all its eigenvalues are positive).