% Problem 4.7

A=[0.012 0.047 0.097 0.072 -0.019
   0.014 0.056 0.115 0.085 -0.023
   0.066 0.252 0.581 0.431 -0.116
   0.028 0.105 0.240 0.178 -0.049
   0.001 0.003 0.003 0.002 -0.001]

Q=eye(5)

P=dlyap(A',Q)

[V,d]=eig(P)