% Problem 4.22

A=[0.995  0.485  0.000 -0.004  0.000 0.000
  -0.017  0.942  0.000 -0.018  0.000 0.000
   0.000 -0.003  0.995  0.485  0.000 0.000
   0.000 -0.014 -0.017  0.942  0.000 0.000
   0.000  0.000  0.000  0.000  1.000 0.500
   0.000  0.000  0.000  0.000  0.002 1.000]

Q=eye(6)

P=dlyap(A',Q)

[V,d]=eig(P)