Solution of the Lyapunov matrix equation
Stephanie J. Barnett, Colin Storey · Electronics Letters · 1966
An estimate is given for the reduction obtained by introducing a skew-symmetric matrix in the number of Operations needed to solve the matrix equation A'L+LA = −C. Also, an iterative scheme is suggested whereby the order of the largest matrix to be inverted can be reduced to 2n−3.