Given a symmetric positive definite matrix Q and a non-singular matrix L, find symmetric positive definite solution X such that X = Q + L (X inv) L^T.
Reference: Benner, P., Faßbender, H. On the Solution of the Rational Matrix Equation.
Benner, Faßbender (2007)