Find real roots for real quartic equations in C++

An Eigen example: find polynomial roots as the eigenvalues of the Companion Matrix