Abstract
We discuss the numerical solution of the Lyapunov equation AHX+XA=−C, C=CH and propose a variant of the Bartels-Stewart algorithm that allows the Cholesky factor of X to be found, without first finding X, when A is stable and C is non-negative definite.