During my research, I was analyzing a linear system with characteristic polynomial of six order, $p(x)=x^6+3ax^5+3a^2x^4+a^3x^3 - a^3b_1b_3b_3$. I used Mathematica and it found all roots as the following:
Image may be NSFW.
Clik here to view.
I was quite impressed since it is not obvious for me how to factor this equation. I am curious how Mathematica find these solutions? What is the algorithm?
Thank you for reading!
