I'm working on an algorithm to calculate all real roots of a quartic equation.
At the moment I have a solution using the method of Descartes and Vietas substitution. My approach is:
- Normalize quartic
- Convert to depressed quartic
- Decompose into quadratic factors $\rightarrow$ yields resolvent cubic
- Convert resolvent cubic to depressed cubic
- Solve the depressed cubic (which involves calculation of a cube root)
- Reverse further previous substitutions (which involves calculation of square roots)
At the moment I'm assuming that all variables are real. If a complex number would occur (e.g. while calculating the cube root or square roots) I'm neglecting the solution at further calculations. However, I'm not sure if there could be the possibility that one of my complex immediate results could lead to a real root of the quartic.
My question: Is this possible? Or in other words: Would it be necessary to also consider complex immediate results during the calculations?
If the answer would be yes, than my second question is: Is there (another) method that allows solving the quartic equation without the need for operating with complex numbers?