lyap — Equation de Lyapunov
[X]=lyap(A,C,flag)
X = lyap(A,C,flag) résout l'équation matricielle de
Lyapunov en temps continu ou discret
A'*X + X*A = C ( flag = 'c' ) A'*X*A - X = C ( flag = 'd' )
Une solution unique existe si A n'a pas de valeur propre
sur l'axe imaginaire (flag='c') ou si 1 n'est pas
valeur propre de A (flag='d').