I try to factorize any polynomial like :
$x^4 + a.x^3 + b.x^2 + a.x + 1$ with $ a, b \in\Bbb{R}$
into :
$(x^2 + c.x + d)(x^2 + e.x + f)$ with $ c, d, e, f \in\Bbb{R}$
I also want $c(a, b)$, $d(a, b)$, $e(a, b)$, $f(a, b)$ to be continous, so I can be smooth at runtime when changing $(a, b)$.
The only way I succeeded for now is to compute all the roots and regroup them by conjugate pairs. But it's tricky because roots can be paired in many way, there are also cases where there are many solutions (when all roots are real), and I noticed roots may swap for specific values of $(a, b)$.
I would like to now if there is a simpler known method for this typical polynom.
Many thanks !