Dear all, I need to solve the following Lyapunov Matrix equation: C=ACA' + B, with A and B given square symmetric matrices. Does anyone knows of a package that can solve the lyapunov matrix equation in R? Or even a C/Fortran implementation? I did not find one on netlib. Thank you.