Is there a general formula for solving 4th degree equations (quartic)?

There is a general formula for solving quadratic equations, namely the Quadratic Formula:

$$x = \frac{ -b \pm \sqrt{ b^2 - 4ac } }{ 2a }$$

For third degree equations of the form $$ax^3+bx^2+cx+d=0$$, there is a set of three equations: one for each root.

Is there a general formula for solving equations of the following form?

$$ax^4 + bx^3 + cx^2 + dx + e = 0$$

How about for higher degrees? If not, why not?

• Could you change the title to "Is there a general formula for solving 4th degree polynomial equations" or "Is there a general formula for solving quartic equations?" – user126 Aug 5 '10 at 23:40
• Put a*x^4+b*x^3+c*x^2+d*x +e = 0 in wolfram alpha (wolframalpha.com) and then sit back and watch the fireworks! – John Alexiou Aug 16 '11 at 2:00
• The general quartic can be reduced to just solving a quadratic. See this simple answer. – Tito Piezas III Apr 4 '15 at 12:50
• – Watson Dec 24 '16 at 17:26

There is, in fact, a general formula for solving quartic (4th degree polynomial) equations. As the cubic formula is significantly more complex than the quadratic formula, the quartic formula is significantly more complex than the cubic formula. Wikipedia's article on quartic functions has a lengthy process by which to get the solutions, but does not give an explicit formula.

Beware that in the cubic and quartic formulas, depending on how the formula is expressed, the correctness of the answers likely depends on a particular choice of definition of principal roots for nonreal complex numbers and there are two different ways to define such a principal root.

There cannot be explicit algebraic formulas for the general solutions to higher-degree polynomials, but proving this requires mathematics beyond precalculus (it is typically proved with Galois Theory now, though it was originally proved with other methods). This fact is known as the Abel-Ruffini theorem.

Also of note, Wolfram sells a poster that discusses the solvability of polynomial equations, focusing particularly on techniques to solve a quintic (5th degree polynomial) equation. This poster gives explicit formulas for the solutions to quadratic, cubic, and quartic equations.

edit: I believe that the formula given below gives the correct solutions for x to $ax^4+bx^3+cx^2+d+e=0$ for all complex a, b, c, d, and e, under the assumption that $w=\sqrt{z}$ is the complex number such that $w^2=z$ and $\arg(w)\in(-\frac{\pi}{2},\frac{\pi}{2}]$ and $w=\sqrt[3]{z}$ is the complex number such that $w^3=z$ and $\arg(w)\in(-\frac{\pi}{3},\frac{\pi}{3}]$ (these are typically how computer algebra systems and calculators define the principal roots). Some intermediate parameters $p_k$ are defined to keep the formula simple and to help in keeping the choices of roots consistent.

Let: \begin{align*} p_1&=2c^3-9bcd+27ad^2+27b^2e-72ace \\\\ p_2&=p_1+\sqrt{-4(c^2-3bd+12ae)^3+p_1^2} \\\\ p_3&=\frac{c^2-3bd+12ae}{3a\sqrt[3]{\frac{p_2}{2}}}+\frac{\sqrt[3]{\frac{p_2}{2}}}{3a} \end{align*} $\quad\quad\quad\quad$

\begin{align*} p_4&=\sqrt{\frac{b^2}{4a^2}-\frac{2c}{3a}+p_3} \\\\ p_5&=\frac{b^2}{2a^2}-\frac{4c}{3a}-p_3 \\\\ p_6&=\frac{-\frac{b^3}{a^3}+\frac{4bc}{a^2}-\frac{8d}{a}}{4p_4} \end{align*}

Then: \begin{align} x&=-\frac{b}{4a}-\frac{p_4}{2}-\frac{\sqrt{p_5-p_6}}{2} \\\\ \mathrm{or\ }x&=-\frac{b}{4a}-\frac{p_4}{2}+\frac{\sqrt{p_5-p_6}}{2} \\\\ \mathrm{or\ }x&=-\frac{b}{4a}+\frac{p_4}{2}-\frac{\sqrt{p_5+p_6}}{2} \\\\ \mathrm{or\ }x&=-\frac{b}{4a}+\frac{p_4}{2}+\frac{\sqrt{p_5+p_6}}{2} \end{align}

(These came from having Mathematica explicitly solve the quartic, then seeing what common bits could be pulled from the horrifically-messy formula into parameters to make it readable/useable.)

• The computational effort is simplified greatly if you "depress" the quartic (i.e. remove the cubic term through a change of variables) by appealing to Vieta's formulae (the mean of the roots is -b/(4a)) first. If you will notice, all the roots of the original equation have a -b/(4a) term added in to compensate for this preliminary translation. – J. M. isn't a mathematician Aug 6 '10 at 0:43
• @J. Mangaldan: Absolutely. My goal (in the edit) was to create a fully-general formula that could be applied straight away; it is not at all illustrative of how to get such a formula. (Your observation is true for the nth-degree polynomial equation formula for n=2, 3, and 4: each formula has a -b/(na) term common to every solution corresponding to the depression.) – Isaac Aug 6 '10 at 0:48
• The formula given above is not correct. One can try some examples: x^4-1,x^4+x^2 and x^4+5*x^2+4. – user14620 Aug 15 '11 at 15:02
• Wikipedia's article has an explicit formula now. – Lee Sleek May 4 '13 at 18:34
• This is not a general solution. (x+2)^4 = x^4 + 8x^3 + 24x^2 + 32x + 16 gives p1 = 0 and p2 = 0, so p3 has zero in the denominator. – Olathe Aug 25 '16 at 10:53

Assuming the four roots $$x_1$$, $$x_2$$, $$x_3$$, $$x_4$$ to be real, they are respectively$$x_1=-\frac{b}{4 a}-\frac{1}{2} \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}-\frac{1}{2} \sqrt{-\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{2 a^2}-\frac{4 c}{3 a}-\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}-\frac{-\frac{b^3}{a^3}+\frac{4 c b}{a^2}-\frac{8 d}{a}}{4 \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}}}$$

$$x_2=-\frac{b}{4 a}-\frac{1}{2} \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}+\frac{1}{2} \sqrt{-\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{2 a^2}-\frac{4 c}{3 a}-\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}-\frac{-\frac{b^3}{a^3}+\frac{4 c b}{a^2}-\frac{8 d}{a}}{4 \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}}}$$

$$x_3=-\frac{b}{4 a}+\frac{1}{2} \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}-\frac{1}{2} \sqrt{-\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{2 a^2}-\frac{4 c}{3 a}-\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}+\frac{-\frac{b^3}{a^3}+\frac{4 c b}{a^2}-\frac{8 d}{a}}{4 \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}}}$$

$$x_4=-\frac{b}{4 a}+\frac{1}{2} \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}+\frac{1}{2} \sqrt{-\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{2 a^2}-\frac{4 c}{3 a}-\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}+\frac{-\frac{b^3}{a^3}+\frac{4 c b}{a^2}-\frac{8 d}{a}}{4 \sqrt{\frac{\sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}{3 \sqrt[3]{2} a}+\frac{b^2}{4 a^2}-\frac{2 c}{3 a}+\frac{\sqrt[3]{2} \left(c^2-3 b d+12 a e\right)}{3 a \sqrt[3]{2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e+\sqrt{\left(2 c^3-9 b d c-72 a e c+27 a d^2+27 b^2 e\right)^2-4 \left(c^2-3 b d+12 a e\right)^3}}}}}}$$

• I think purpose of maths is to make less effort to arrive the same result. this is the opposite case – Fennekin Feb 9 '17 at 17:31
• I've condensed the post by using $\gamma$. (The page was taking a long time to load.) – Tito Piezas III Jul 2 '17 at 5:02

Edited in response to Quonux's comments.

Yes. As an answer I will use a shorter version of this Portuguese post of mine, where I deduce all the formulae. Suppose you have the general quartic equation (I changed the notation of the coefficients to Greek letters, for my convenience):

$$\alpha x^{4}+\beta x^{3}+\gamma x^{2}+\delta x+\varepsilon =0.\tag{1}$$

If you make the substitution $x=y-\frac{\beta }{4\alpha }$, you get a reduced equation of the form

$$y^{4}+Ay^{2}+By+C=0\tag{2},$$

with

$$A=\frac{\gamma }{\alpha }-\frac{3\beta ^{2}}{8\alpha ^{2}},$$

$$B=\frac{\delta }{\alpha }-\frac{\beta \gamma }{2\alpha ^{2}}+\frac{\beta^3 }{ 8\alpha^3 },$$

$$C=\frac{\varepsilon }{\alpha }-\frac{\beta \delta }{4\alpha ^{2}}+\frac{ \beta ^{2}\gamma}{16\alpha ^{3}}-\frac{3\beta ^{4}}{256\alpha ^{4}}.$$

After adding and subtracting $2sy^{2}+s^{2}$ to the LHS of $(2)$ and rearranging terms, we obtain the equation

$$\underset{\left( y^{2}+s\right) ^{2}}{\underbrace{y^{4}+2sy^{2}+s^{2}}}-\left[ \left( 2s-A\right) y^{2}-By+s^{2}-C\right] =0. \tag{2a}$$

Then we factor the quadratic polynomial $$\left(2s-A\right) y^{2}-By+s^{2}-C=\left(2s-A\right)(y-y_+)(y-y_-)$$ and make $y_+=y_-$, which will impose a constraint on $s$ (equation $(4)$). We will get:

$$\left( y^{2}+s+\sqrt{2s-A}y-\frac{B}{2\sqrt{2s-A}}\right) \left( y^{2}+s- \sqrt{2s-A}y+\frac{B}{2\sqrt{2s-A}}\right) =0,$$ $$\tag{3}$$

where $s$ satisfies the resolvent cubic equation

$$8s^{3}-4As^{2}-8Cs+\left( 4AC-B^{2}\right) =0.\tag{4}$$

The four solutions of $(2)$ are the solutions of $(3)$:

$$y_{1}=-\frac{1}{2}\sqrt{2s-A}+\frac{1}{2}\sqrt{-2s-A+\frac{2B}{\sqrt{2s-A}}}, \tag{5}$$

$$y_{2}=-\frac{1}{2}\sqrt{2s-A}-\frac{1}{2}\sqrt{-2s-A+\frac{2B}{\sqrt{2s-A}}} ,\tag{6}$$

$$y_{3}=\frac{1}{2}\sqrt{2s-A}+\frac{1}{2}\sqrt{-2s-A-\frac{2B}{\sqrt{2s-A}}} ,\tag{7}$$

$$y_{4}=\frac{1}{2}\sqrt{2s-A}-\frac{1}{2}\sqrt{-2s-A-\frac{2B}{\sqrt{2s-A}}} .\tag{8}$$

Thus, the original equation $(1)$ has the solutions $$x_{k}=y_{k}-\frac{\beta }{4\alpha }.\qquad k=1,2,3,4\tag{9}$$

Example: $x^{4}+2x^{3}+3x^{2}-2x-1=0$

$$y^{4}+\frac{3}{2}y^{2}-4y+\frac{9}{16}=0.$$

The resolvent cubic is

$$8s^{3}-6s^{2}-\frac{9}{2}s-\frac{101}{8}=0.$$

Making the substitution $s=t+\frac{1}{4}$, we get

$$t^{3}-\frac{3}{4}t-\frac{7}{4}=0.$$

One solution of the cubic is

$$s_{1}=\left( -\frac{q}{2}+\frac{1}{2}\sqrt{q^{2}+\frac{4p^{3}}{27}}\right) ^{1/3}+\left( -\frac{q}{2}-\frac{1}{2}\sqrt{q^{2}+\frac{4p^{3}}{27}}\right) ^{1/3}-\frac{b}{3a},$$

where $a=8,b=-6,c=-\frac{9}{2},d=-\frac{101}{8}$ are the coefficients of the resolvent cubic and $p=-\frac{3}{4},q=-\frac{7}{4}$ are the coefficients of the reduced cubic. Numerically, we have $s_{1}\approx 1.6608$.

The four solutions are :

$$x_{1}=-\frac{1}{2}\sqrt{2s_{1}-A}+\frac{1}{2}\sqrt{-2s_{1}-A+\frac{2B}{ \sqrt{2s-A}}}-\frac{\beta }{4\alpha },$$

$$x_{2}=-\frac{1}{2}\sqrt{2s_{1}-A}-\frac{1}{2}\sqrt{-2s_{1}-A+\frac{2B}{ \sqrt{2s-A}}}-\frac{\beta }{4\alpha },$$

$$x_{3}=\frac{1}{2}\sqrt{2s_{1}-A}+\frac{1}{2}\sqrt{-2s_{1}-A-\frac{2B}{ \sqrt{2s-A}}}-\frac{\beta }{4\alpha },$$

$$x_{4}=\frac{1}{2}\sqrt{2s_{1}-A}-\frac{1}{2}\sqrt{-2s_{1}-A-\frac{2B}{ \sqrt{2s-A}}}-\frac{\beta }{4\alpha },$$

with $A=\frac{3}{2},B=-4,C=\frac{9}{16}$. Numerically we have $x_{1}\approx -1.1748+1.6393i$, $x_{2}\approx -1.1748-1.6393i$, $x_{3}\approx 0.70062$, $x_{4}\approx -0.35095$.

Another method is to expand the LHS of the quartic into two quadratic polynomials, and find the zeroes of each polynomial. However, this method sometimes fails. Example: $x^{4}-x-1=0$. If we factor $x^{4}-x-1$ as $x^{4}-x-1=\left( x^{2}+bx+c\right) \left( x^{2}+Bx+C\right)$ expand and equate coefficients we will get two equations, one of which is $-1/c-c^{2}\left( 1+c^{2}\right) ^{2}+c=0$. This is studied in Galois theory.

The general quintic is not solvable in terms of radicals, as well as equations of higher degrees.

• Great answer! Howevern, could you give an example of a quartic for which the method described in the last paragraph doesn't work? – Heinz Doofenschmirtz Jul 28 '15 at 8:19

There most certainly is, but its ugly, complicated, and not worth memorizing. People know about it and have quoted or cited it for you, but really they would never use it. If you want something actually useful for pen-and-paper solutions, you might want to understand the actual theory behind the solution. I will provide one method below.

The Quartic Formula is just the end result of this methodology, written in terms of the original coefficients. Because of that, the method is far easier to remember than the formula, which is why I find it annoying when people cite just the formula and tell you, "dont bother, use a computer instead." A pen-paper solution is not complicated, it's just time consuming.

Understanding how its done - even if you never use it - expands your brain and your understanding, allows you to implement it in programming, and allows you to recreate it whenever you may need it instead of hyper reliance on computers to always be there for you, which, in my opinion, makes for a poor mathematician.

There are three methods for solving Quartics that I both know and know of:

• Descartes' Quadratic Factorization
• Euler's Method
• Ferrari's Method

If anyone knows of more, please let me know.

Ferrari's Method is historically the first method discovered. Euler's Method looks a lot like Cardano's Method for the Cubic and was probably modeled after the same approach. But I am partial to Descartes' Quadratic Factorization technique. Its a relatively simple process to follow and is what I will be using below. If you want to see how the others work, let me know.

All of the above methods start out the same: depression (removing the $$n-1$$ degree term, in this case the cubic term) and normalization (making the lead coefficient 1, i.e. making the polynomial monic).

They all end up roughly the same place too: solving a cubic equation. So you need to be prepared for that. I recommend you brush up on that; I wont explain solving the cubic here but will only refer to it.

Depress and Normalize

Given the arbitrary quartic: $$ax^4 + bx^3 + cx^2 + dx + e =0$$ Where $$a,b,c,d,e\in\mathbb{R}$$.

We must convert this quartic into a depressed monic quartic. First we depress it by substituting $$x = z-\frac{b}{4a}$$. We arrive at the quartic: $$az^4 + Bz^2 + Cz + D =0$$ For some $$B,C,D\in\mathbb{R}$$. This quartic, now in $$z$$, is merely a horizontal shift from the original quartic in $$x$$. All points have shifted, therefore the roots have shifted, but by the same constant. Also, $$B,C,D$$ are computed from the previous $$a,b,c,d,e$$ and have no dependence on $$x,z$$. Whats important is that there is no $$z^3$$ term. Next we divide through by the lead coefficient $$a$$ to make this quartic monic (normalize). $$z^4 + pz^2 + qz + r = 0\;\;\;\;\;\;\;\;\;\;\;\;(1)$$

For some $$p,q,r\in\mathbb{R}$$. This has no effect on the root positions; all the lead coefficient did was scale up the non-zero values. There is absolutely no loss in generality with respect to the zeros. All of these constants, $$p,q,r$$, can be computed from the original coefficients, $$a,b,c,d,e$$. There is still no cubic term and now there is no lead coefficient either.

What might be interesting to note is what has happened to the polynomial. We started out with 5 arbitrary constants and have reduced it to 3, by normalizing the lead and removing the cubic term. We originally had arbitrary values of $$a,b,c,d,e\in\mathbb{R}$$, and now we have arbitrary values $$p,q,r\in\mathbb{R}$$. Although the latter three are computed from the original five, they have arbitrary values, and there is no loss in generality. This is a significant simplification of the problem. The non-existence of the cubic term will prove vital.

Everything so far has simply been the setup: writing the polynomial in reduced monic form. Recall all of the quartic methods accomplish at least this much. Next we implement Descartes Factorization method.

Descartes Factorization Method

We must assume that all of the coefficients are real, $$p,q,r\in\mathbb{R}$$. This is a required condition to make the methodology work. The reason is because now all solutions with non-zero imaginary components come in complex conjugate pairs. Big deal? It allows us to group two solutions together, even if they are purely real, into quadratic factors with real coefficients. We know that all real-valued monic quartics can be factored into: $$(z^2 + mz+n)(z^2+sz+t)=0\;\;\;\;\;\;\;\;\;\;\;\;(2)$$ Where $$m,n,s,t\in\mathbb{R}$$. Clearly the task at hand is to find the values of these constants. If we can convert (1) into (2) by determining quadratic factors that satisfies the equation, we can find the roots with the quadratic formula applied to those quadratic factors.

What we do is distribute these quadratic factors out into normal polynomial form and you get: $$z^4 + (m+s)z^3 + (t+n+ms)z^2 + (mt+ns)z + nt = 0$$ Which we compare term-by-term to the depressed monic in (1). You get this system of equations: $$m+s=0 \\ n+t+ms=p \\mt+ns=q \\ nt=r$$ Its plain to see then that $$s = -m$$ and that $$t=\frac{r}{n}$$ from the first and fourth equations. Our quadratic factors can be rewritten: $$(z^2 + mz+n)(z^2-mz+\frac{r}{n})=0\;\;\;\;\;\;\;\;\;\;\;\;(3)$$ We only have two unknowns in this factorization: $$m$$ and $$n$$. The remaining two of the four simultaneous equations (the second and the third equation) can be rewritten: $$\frac{r}{n}+n-m^2 = p \\ m(\frac{r}{n} -n)=q$$ By moving the $$m$$ to the right-hand side: $$\frac{r}{n}+n = p+m^2 \\ \frac{r}{n} -n=\frac{q}{m}$$

From here we form two new equations by adding and subtracting the previous two. By adding we get: $$2\frac{r}{n} = p+m^2 + \frac{q}{m}$$ By subtracting we get: $$2n = p+m^2 -\frac{q}{m}$$ Notice on the left hand side that both of these equations can be solved for $$n$$ and $$\frac{r}{n}$$ readily in terms of $$m$$, both of which appear in the quadratic factors of (3). These can be utilized later once we know the $$m$$ to complete the quadratic factor.

We can find the $$m$$ by taking these latest two equations and multiplying them, therefore eliminating the unknown $$n$$. Notice that $$n$$ is in the numerator in one and in the denominator in the other. Thus: $$4r = (p+m^2)^2 - (\frac{q}{m})^2 = m^4 + 2pm^2 + p^2 - \frac{q^2}{m^2}$$

Thus we are essentially down to $$m$$ as our last unknown. Everything else is known in terms of $$m$$, and $$m$$ is the only unknown in the above equation. One unknown, one remaining equation. Multiply the equation through by $$m^2$$ and rearrange: $$m^6 + 2pm^4 + (p^2-4r)m^2 - q^2 = 0$$ This is still rather ugly. It's a sextic, not a quartic, which is worse. But notice that the powers of $$m$$ are all even. We can substitute $$w=m^2$$ and arrive at: $$w^3 + 2pw^2 + (p^2 - 4r)w - q^2 =0$$

And thus we are essentially done. We are left with a cubic polynomial in $$w$$, which is solvable with its own techniques. Techniques which I only assume you already know about if you are trying to solve quartics. Just like with quartics, as you know already, there do exist cubic formulae, but I do recommend learning the methods behind them.

If you need help with cubics, I recommend Cardano's method (the original solution) or Vieta's Trigonometric Solution (my preferred). There is also Completing the Cube, a nice proof of concept but Id never use it. Feel free to ask a separate question for a cubic and I will be happy to answer.

The point is that the problem has been reduced from that of finding the roots of a quartic to that of finding the roots of a cubic. A simpler problem! That's usually how it goes. All quartic root finding methods require finding the roots of a cubic first, obvious or not. Just as finding the roots of a cubic involve solving a quadratic. Hope this works for you.

Anyway, solve the cubic in $$w$$. Then with that recall the quadratic substitution $$w=m^2$$ and solve for $$m$$. Then recall that we previously had equations for $$2n$$ and $$2\frac{r}{n}$$. And now you know all of the unknown terms in the quadratic factorization $$(z^2 + mz + n)(z^2 - mz + \frac{r}{n})=0$$.

Still not done. Each of these quadratic factors must now be solved using the quadratic formula, and you have solutions in $$z$$. This solves the depressed monic quartic we started Descartes Quadratic Factorization method with.

Finally

Dont forget about the original quartic we had at the very beginning, prior to the depression and normalization. We had introduced a horizontal shift of $$x = z-\frac{b}{4a}$$. Doing this last bit will solve the original quartic in terms of $$x$$, which is the solution you want.

You are going to arrive at a set of solutions when done. Be sure to check your answers. You may have redundant or superfluous solutions. Some redundant solutions may be written in very different algebraic ways, but will represent the same numeric value.

If you express the final answer of $$x$$ in terms of the original $$a,b,c,d,e$$ you will just have the same "quartic formulae" that other people are citing you. The expression will of course be slightly different depending on which of the quartic methods you employ.

Concerns

If youre concerned about the assumption that coefficients $$p,q,r$$ are real, dont be. All it means is that $$a,b,c,d,e$$ are real, which is usually a good assumption. In fact, we can generalize. The values $$p,q,r$$ can be made complex, implying only that the original quartic has complex $$a,b,c,d,e$$. It also means you'll have to solve a cubic with complex coefficients. This is doable and the math still works fine.

We can reduce the problem of solving the general quartic to merely solving a quadratic. Given,

$$x^4+ax^3+bx^2+cx+d=0$$

Then the four roots are,

$$x_{1,2} = -\frac{a}{4}+\frac{\color{red}\pm\sqrt{u}}{2}\color{blue}+\frac{1}{4}\sqrt{3a^2-8b-4u+\frac{-a^3+4ab-8c}{\color{red}\pm\sqrt{u}}}\tag1$$

$$x_{3,4} = -\frac{a}{4}+\frac{\color{red}\pm\sqrt{u}}{2}\color{blue}-\frac{1}{4}\sqrt{3a^2-8b-4u+\frac{-a^3+4ab-8c}{\color{red}\pm\sqrt{u}}}\tag2$$

where,

$$u = \frac{a^2}{4}-\frac{2b}{3} +\frac{1}{3}\left(v_1^{1/3}\zeta_3+\frac{b^2 - 3 a c + 12 d}{v_1^{1/3}\zeta_3}\right)\tag{3a}$$

or alternatively,

$$u = \frac{a^2}{4}-\frac{2b}{3} +\frac{1}{3}\left(\color{blue}{v_1^{1/3}+v_2^{1/3}}\right)\tag{3b}$$

with $v_1$ any non-zero root of the quadratic,

$$v^2 + (-2 b^3 + 9 a b c - 27 c^2 - 27 a^2 d + 72 b d)v + (b^2 - 3 a c + 12 d)^3 = 0$$

and a chosen cube root of unity $\zeta_3^3 = 1$ such that $u$ is also non-zero. (Normally, use $\zeta_3=1$, but not when $a^3-4ab+8c = 0$.)

P.S. This is essentially the method used by Mathematica, though much simplified for aesthetics.

• I noticed that by defining $x=(v_1+v_2)/2$ and $y=(v_1-v_2)/2$ one can rewrite $v_1^{1/3}+v_2^{1/3}=2(b^2-3ac+12d)^{1/2}\cos\left(\frac{2}{3}\arctan\left(((x^2+y^2)^{1/2}-x)/y\right)\right)$ - which is true at least if the square root in $v_i$ turns out imaginary. Now I wonder if there is any trigonometric identity that can resolve the $\cos(2/3\arctan(...))$ into something rational involving roots? – Kagaratsch Sep 11 '15 at 1:06
• @Status: No, it does not imply $x^4=0$. I've clarified the post. If $v_1, v_2$ are zero, then the expression in blue in $(3b)$ is also zero, but one can still use this $u$ to find the roots you mentioned. (The case when $u, v_1, v_2$ are all zero is degenerate and the quartic has 4 identical roots.) To explain the two forms, the advantage of $(3a)$ is there is only one complex cube root $v_1^{1/3}$, while the advantage of $(3b)$ is there is no variable denominator (hence no division by zero). – Tito Piezas III Sep 21 '17 at 12:00
• Solving the quartic requires the resolution of a cubic (itself requiring the resolution of a quadratic and possibly trigonometric functions), then of two quadratics. This is provably unavoidable. – Yves Daoust Oct 26 '17 at 18:42

Yes, there is a quartic formula.

There is no such solution by radicals for higher degrees. This is a result of Galois theory, and follows from the fact that the symmetric group $S_5$ is not solvable. It is called Abel's theorem. In fact, there are specific fifth-degree polynomials whose roots cannot be obtained by using radicals from $\mathbb{Q}$.

What has not been mentioned thus far is that one can in fact use any number of "auxiliary cubics" in the solution of the quartic equation. Don Herbison-Evans, in this page (Wayback link, as the original page has gone kaput), which was adapted from his technical report, mention five such possible auxiliary cubics.

Given the quartic equation

$$x^4 + ax^3 + bx^2 + cx + d = 0$$

the five possible auxiliary cubics are referred to in the document as

Christianson-Brown:

$$y^3 +\frac{4a^2b - 4b^2 - 4ac + 16d - \frac34a^4}{a^3 - 4ab + 8c}y^2 + \left(\frac3{16}a^2 - \frac{b}{2}\right)y - \frac{1}{64}(a^3 + 4a b - 8c) = 0$$

Descartes-Euler-Cardano:

$$y^3 + \left(2b - \frac34 a^2\right)y^2 + \left(\frac3{16}a^4 - a^2b + ac + b^2 - 4d\right)y - \frac{1}{64}(a^3 + 4a b - 8c)^2 = 0$$

Ferrari-Lagrange

$$y^3 + by^2 + (ac - 4d)y + a^2d + c^2 - 4bd = 0$$

Neumark

$$y^3 - 2by^2 + (ac + b^2 - 4d)y + a^2d - abc + c^2 = 0$$

Yacoub-Fraidenraich-Brown

$$(a^3 - 4ab + 8c)y^3 + (a^2b - 4b^2 + 2ac + 16d)y^2 + (a^2c - 4bc + 8ad)y + a^2d - c^2 = 0$$

See the page for how to obtain the quadratics that will yield the solutions to the original quartic equation from a root of the auxiliary cubic.

Let me also mention this old ACM algorithm in Algol. Netlib has a C implementation of that algorithm.

• Interesting content in that article. Thanks. Im wondering if you are aware that some of these are not really unique. If you depress the quartic and make $a=0$ then the YFB and CB equations are identical, and the DEC and N are horizontal reflections of one another. After a simple depression there are only three unique equations to concern yourself with. – CogitoErgoCogitoSum Jul 1 '17 at 2:39

Regarding the inability to solve the quintic, this is sort-of true and sort-of false. No, there is no general solution in terms of $+$, $-$, $\times$ and $\div$, along with $\sqrt[n]{}$. However, if you allow special theta values (a new operation, not among the standard ones!) then yes, you can actually write down the solutions of arbitrary polynomials this way. Also, you can do construct lengths equal to the solutions by intersecting lower degree curves (for a quintic, you can do so with a trident and a circle.)

• "However, if you allow special theta values (a new operation, not among the standard ones!) then yes, you can actually write down the solutions of arbitrary polynomials this way." Do you know of a reference which expands on this point? – Zach Conn Dec 8 '10 at 20:35
• @Zach: Umemura shows how to use (multidimensional) theta functions for representing roots of polynomials. – J. M. isn't a mathematician Aug 16 '11 at 4:47
• On the topic of solving the quintic I recommend the paper by Peter Doyle and Curt McMullen, "Solving the quintic by iteration", Acta Math. 163 (1989), no. 3-4, 151–180. – Lee Mosher Jul 1 '17 at 23:31
• @Lee, as the paper's title itself notes, the method presented is an iterative method, and would not furnish a symbolic solution. – J. M. isn't a mathematician Nov 14 '17 at 4:16
• In the spirit of this answer, I'm only pointing out (and hence it was just a comment) another method in the topic of "solving" polynomial equations. – Lee Mosher Nov 14 '17 at 13:17

To your question, "if not, why not", I could say this (which I think goes a little further that the previous answers). A polynomial $p(x) \in \mathbb{Q}[x]$ of degree $n$ has a Galois group $G = G(p)$ attached to it. This is a subgroup of the symmetric group $S_{n},$ and this identification comes about because $G$ is a group of permutations of the roots of $p(x)$. The equation $p(x)$ is solvable by radicals if and only if $G$ is a solvable group, which is a key theorem of Galois theory. When $n \leq 4$, the symmetric group $S_n$ is itself a solvable group, so all its subgroups are solvable, and the group $G$ must be solvable. As remarked in earlier answers, when $n \geq 5,$ the group $S_n$ is never solvable. This does not mean that the polynomial $p(x)$ is never solvable by radicals, but (depending what its Galois group is), we can not be sure that it is (and there are explicit examples when it is not for each $n \geq 5$). I do make the point though that the solvability or otherwise of $p(x)$ by radicals really depends on the factorization of $p(x)$ as a product of irreducible polynomials, rather than just the degree of $p(x)$. If $p(x)$ is a product of irreducible polynomials which each have degree at most 4, then its Galois group is solvable, so $p(x)$ is solvable by radicals.

Contrary to common opinion, the resolution of the quartic with real coeficients is not so difficult.

I take for granted that depletion and reduction of the equation by a linear change of variable

$$\alpha x^4+\beta x^3+\gamma x^2+\delta x+\epsilon=0\to x^4+px^2+qx+r=0$$ are known.

Then we factor this depleted quadrinomial as

$$(x^2+ax+b)(x^2-ax+c)=x^4+(-a^2+b+c)x^2+a(c-b)x+bc$$ and identify the coefficients,

$$-a^2+b+c=p,\\a(c-b)=q,\\bc=r.$$

We can express $$c\pm b$$ in terms of $$a$$, solve for $$b,c$$ then form the product $$4a^2bc$$:

$$c-b=\frac qa,\\c+b=p+a^2,$$ hence $$4a^2bc=\left(a(p+a^2)-q\right)\left(a(p+a^2)+q\right)=a^2(p^2+a^2)^2-q^2=4a^2r.$$

This is a cubic equation in $$a^2$$, and fortunately it always has a positive solution.

Then after resolution, from $$a$$ we obtain $$b$$ and $$c$$, and the roots of the two quadratic trinomials.

• @ Yves: Thanks for the updated notation. A linear change represents here the substitution $x=y-b/4$? Does it lead to the next substitution $z=a^2$, where $f(z)=4r$ is the cubic for $z$? I am having some problems with the quartic factorization in my question. Maybe, you will find this question interesting :-). Thanks. – justik Aug 13 '18 at 17:08

Yes, there is such a formula as the "Quartic Formula." However, it is too complicated for practical use, so if someone needs to find a root to one of these equations, one may use the Rational Root Theorem, Factor Theorem. If the roots are approximated, then one would use Newton's Method (studied in calculus.) Newton's method is how commercial and scientific computer programs approximate roots of all equations (not just polynomials).

I'm making a algorithm for my app. So I've tried many methods, including the general and direct formula described in Wikipedia (Quartic function).

I had problems with the general formula. Not all equations can be solved. The same happens with the first and third answer here in this question, so I've followed the Ferrari solution in Wikipedia and finally it works smoothly:

Notice that I use cubic root below. The result here is always the principal value. It can be a complex number. In general, a cubic root has 3 different roots. The square root has two results, here is always the positive one.

The starting point is

$$a_1x^4 + a_2x^3+a_3x^2 + a_4x + a_5$$

Below, I explaind the method with exeptions, including the 3rd degree auxiliar equation, step by step:

a) Divide the equation for $$a_1$$ getting the new coefficients $$b$$, $$c$$ , $$d$$ and $$e$$

$$\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad x^4 + bx^3 + cx^2 + dx + e \qquad\qquad\qquad$$ (1)

b) Obtain the below quartic depressed (quartic without $$x^3$$ term, so $$b$$ will be $$0$$). As the sum of roots is represented by $$-b$$, for eliminating the term with $$b$$, the new equation should have sum of roots equal to $$0$$. It's enough to replace $$x$$ by $$x'-b/4$$, where $$b/4$$ is the fourth part of symmetric of the sum of the roots. So the new equation will have its roots increased by $$b/4$$. After will be necessary reduces all roots by $$-b/4$$. Using $$x$$ and not $$x'$$ for simplicity, it results in:

$$x^4 + px^2 + qx + r$$

... if one calculates $$p$$, $$q$$ and $$r$$

$$\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad p = \dfrac {8c - 3b^2}{8}\qquad\qquad\qquad\qquad$$(2)

$$\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad q = \dfrac {b^3 - 4bc + 8d}{8}\qquad\qquad\qquad$$(3)
$$r = \dfrac {-3b^4 + 256e - 64bd + 16b^2c}{256}$$

c) If $$q=0$$ is results in the following biquadratic equation, that can be easily solved:

$$x^4 + px^2 + r$$

• Do $$y = x^2$$
• Solve the 2nd degree equation $$y^2 + py + r = 0$$
• Suppose that the solution is $$y_1$$ and $$y_2$$ ($$y_1$$ can be igual to $$y_2$$)
• Finally the quartic roots will be the square root of the root minus minus the original displacement: $$x_1=\sqrt y_1 -b/4$$, $$\enspace x_2=-\sqrt y_1 -b/4$$, $$\enspace\ x_3 =\sqrt y_2 - b/4 \enspace$$ and $$\enspace x_4 = 2-\sqrt y_2 -b/4$$

d) if not biquadratic, obtain the resolvent cubic...

$$tx^3 + ux^2 + vx + w$$

... calculating $$t$$, $$u$$, $$v$$ and $$w$$

$$t = 8$$ $$u = 8p$$ $$v = 2p^2 - 8r$$ $$w = -q^2$$

e) Solve the cubic using the Cardano's formula. Before we eliminate $$x^2$$ term (with $$u$$). One should replace $$x$$ by $$x'-u/3$$. Again, $$u/3$$ is a third part of symmetric of the sum of the roots represented by $$-u$$, so the new equation will have a sum of the roots equal to $$0$$. The new roots are higher in $$u/3$$, so after will be necessary reduces all roots by $$-u/3$$. Using $$x$$ and not $$x'$$ for simplicity, it results in:

• The new depressed cubic has the form: $$x^3 + px + q= 0$$

where $$p = c- b^2/3$$ $$q = \dfrac {2b^3 - 9bc}{27} + d$$

• if $$q=0$$ the roots trivially are $$r_1=0$$, $$\enspace r_2=\sqrt {-p}\enspace$$ and $$\enspace r_3= -\sqrt {-p}\enspace$$, so the original roots will be $$x_1=-b/3$$, $$\enspace x_2=\sqrt {-p} - b/3\enspace$$ and $$\enspace x_3=-\sqrt {-p}- b/3$$

• if $$q<>0$$ use the below formula for calculate the first root, a real one

$$r_1 = \sqrt [3] {- \dfrac q 2 + \sqrt {\dfrac {q^2}{4} + \dfrac {p^3}{27} } } +\sqrt [3] {- \dfrac q 2 - \sqrt {\dfrac {q^2}{4} + \dfrac {p^3}{27} } }$$

So the first root is

$$x_1 = r_1 - b/3$$

• The other roots can be calculated easily remember that the sum of roots is 0 and the product of roots is $$-q$$ (Vieta's Formula)

$$r_1 + r_2 + r_3 = 0$$ $$r_1 \times r_2 \times r_3 = -q$$

• It result in a system of 2 equations: (call the known root $$r_1$$ of $$s$$)

$$r_2 + r_3 =-s$$ $$r_2 \times r_3 = -q/s$$

• So it obtains, in the end, a 2nd degree equation:

$$r_2 -q/(s r_2) =-s$$ $$r^2_2 -(q/s) =-s r_2$$ $$r^2_2 +s r_2 -(q/s) = 0$$

• For simmetry the solutions of this 2nd degree equation corresponds to $$r_1$$ and $$r_2$$ and it generates, finally, the final auxiliar cubic equation solution

$$x_1 = r_1 - b/3$$ $$x_2 = r_2 - b/3$$ $$x_3 = r_3 - b/3$$

f) From above, choose a nonzero cubic root $$m$$ from $$x_1$$, $$x_2$$ or $$x_3$$.

$$\qquad\qquad\qquad\qquad m$$ is a nonzero cubic root from $$\quad tx^3 + ux^2 + vx + w\qquad\qquad$$ (4)

g) So using $$b$$, $$p$$, $$q$$ and $$m$$ defined above in (1), (2), (3) and (4), the four roots are:

$$x_1 = \dfrac {-b}{4} + \frac {\sqrt {2m} + \sqrt{ -2p - 2m - \dfrac{q \sqrt 2}{\sqrt m} }}{2}$$ $$x_1 = \dfrac {-b}{4} + \frac {\sqrt {2m} - \sqrt{ -2p - 2m - \dfrac{q \sqrt 2}{\sqrt m} }}{2}$$ $$x_3 = \dfrac {-b}{4} - \frac {\sqrt {2m} - \sqrt{ -2p - 2m + \dfrac{q \sqrt 2}{\sqrt m} }}{2}$$ $$x_4 = \dfrac {-b}{4} - \frac {\sqrt {2m} - \sqrt{ -2p - 2m + \dfrac{q \sqrt 2}{\sqrt m} }}{2}$$

I have used Wolfram Alpha for check some equations and its roots and it works in all cases.

If anyone is interested, there is a quartic equation calculator with a method and complete steps to solving the roots in the link here.

• Which method would be easiest to implement given that one is interested in real roots only? Does that simplify things at all? – Jim White Aug 29 '17 at 11:42

There are formulas for up to fourth-degree polynomials. From there, no. The demonstration that there can not be is Group Theory and it is something strong.

There is a general formula for solving fourth-degree equations, too.

Quartic equation:

$$t = p x + q x^2 + r x^3 + s x^4$$

Solution:

$$m = r^3 - 4 q r s + 8 p s^2$$

$$f = -(r^6 - 6 q r^4 s + 8 q^2 r^2 s^2 + 8 p r^3 s^2 - 16 p q r s^3 + 8 p^2 s^4 + 8 r^2 s^3 t)$$

$$a = 2 s (3 r^5 - 16 q r^3 s + 16 q^2 r s^2 + 20 p r^2 s^2 - 16 p q s^3 + 32 r s^3 t)$$

$$b = 4 s^2 (3 r^4 - 14 q r^2 s + 8 q^2 s^2 + 20 p r s^2 + 32 s^3 t)\\ c = 8 m s^3$$

$$X = \dfrac{-(a b + 9 c f) + \sqrt{(a b + 9 c f)^2 - 4 (b^2 - 3 a c) (a^2 + 3 b f)}}{2 (b^2 - 3 a c)}$$

$$A = a + 2 b X + 3 c X^2$$

$$F = a X + b X^2 + c X^3 - f$$

$$W = (A^3 - 27 c F^2)^{1/3}$$

$$Y = X + \dfrac{3 F}{W - A}$$

$$T = p Y + q Y^2 + r Y^3 + s Y^4 - t$$

$$P = 3 m r + 4 s^2 (q^2 - 2 p r + 4 s t) + 4 s (3 r^3 - 10 q r s + 12 p s^2) Y + 4 s^2 (3 r^2 - 8 q s) Y^2$$

$$Q = m (r + 4 s Y) (m (r + 4 s Y) - 16 s^3 T)$$

$$R = -r^2 + 2 q s - 2 r s Y - 4 s^2 Y^2$$

$$R2 = \pm \sqrt{P \pm 2 \sqrt{Q}}$$

$$Z = \dfrac{R2 - R}{4 s^2}$$

$$x = Y \pm \sqrt{Z}$$

But this solution gives eight roots, four true and four excess.

Wolfram-code for verifing:

{t, p, q, r, s} = RandomInteger[{-1000, 1000}, 5];
Print["\nEquation: ", t, " = ", p, "*x + ", q, "*x^2 + ", r, "*x^3 + ", s, "*x^4\n"];
Print["Ordinary solution:"];
Print[x /. (-t + p x + q x^2 + r x^3 + s x^4 // NSolve) // Sort, "\n"];
m = r^3 - 4 q r s + 8 p s^2;
f = -(r^6 - 6 q r^4 s + 8 q^2 r^2 s^2 + 8 p r^3 s^2 - 16 p q r s^3 + 8 p^2 s^4 + 8 r^2 s^3 t);
a = 2 s (3 r^5 - 16 q r^3 s + 16 q^2 r s^2 + 20 p r^2 s^2 - 16 p q s^3 + 32 r s^3 t);
b = 4 s^2 (3 r^4 - 14 q r^2 s + 8 q^2 s^2 + 20 p r s^2 + 32 s^3 t);
c = 8 m s^3;
X = (-(a b + 9 c f) + Sqrt[(a b + 9 c f)^2 - 4 (b^2 - 3 a c) (a^2 + 3 b f)])/(2 (b^2 - 3 a c));
A = a + 2 b X + 3 c X^2;
F = a X + b X^2 + c X^3 - f;
W = (A^3 - 27 c F^2)^(1/3);
Y = X + (3 F)/(W - A);
T = p Y + q Y^2 + r Y^3 + s Y^4 - t;
P = 3 m r + 4 s^2 (q^2 - 2 p r + 4 s t) + 4 s (3 r^3 - 10 q r s + 12 p s^2) Y + 4 s^2 (3 r^2 - 8 q s) Y^2;
Q = m (r + 4 s Y) (m (r + 4 s Y) - 16 s^3 T);
R = -r^2 + 2 q s - 2 r s Y - 4 s^2 Y^2;
R2 = {Sqrt[P + 2 Sqrt[Q]], Sqrt[P - 2 Sqrt[Q]], -Sqrt[P + 2 Sqrt[Q]], -Sqrt[P - 2 Sqrt[Q]]};
Z = (R2 - R)/(4 s^2);
Print["Not-ordinary solution:"];
Union[Y + Sqrt[Z], Y - Sqrt[Z]] // N // Sort