This is a socalledbelongs to the class of "nonsymmetric algebraic Ricatti equationequations", which can be solved following the method explained in Nonsymmetric algebraic Riccati equations and Wiener–Hopf factorization for M-matrices (2001). See also Federico Poloni's general exposition and in Federico Poloni's paper on Quadratic vector equations (2010).