Lyapunov Discrete Time Equation

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.