% Problem 4.7

A=[-5 0 1 2 
   0 -4 2 -1
   1 0 -3 0
   0 -1 1 -5]

Q=eye(4)

P=lyap(A',Q)

[V,d]=eig(P)