RiccatiSolve

RiccatiSolve[{a, b}, {q, r}]

RiccatiSolve[{a, b}, {q, r, p}]