Quartic equation

Consider a quartic equation expressed in the form ${\displaystyle a_{0}x^{4}+a_{1}x^{3}+a_{2}x^{2}+a_{3}x+a_{4}=0}$:

If a₄ (the constant term)  = 0, then one of the roots is x = 0, and the other roots can be found by dividing by x, and solving the resulting cubic equation,

${\displaystyle a_{0}x^{3}+a_{1}x^{2}+a_{2}x+a_{3}=0.\,}$

Call our quartic polynomial Q(x). Since 1 raised to any power is 1, ${\displaystyle Q(1)=a_{0}+a_{1}+a_{2}+a_{3}+a_{4}}$. Thus if ${\displaystyle a_{0}+a_{1}+a_{2}+a_{3}+a_{4}=0}$, Q(1) = 0 and so x = 1 is a root of Q(x). It can similarly be shown that if ${\displaystyle a_{0}+a_{2}+a_{4}=a_{1}+a_{3}}$, x = −1 is a root.

In either case the full quartic can then be divided by the factor (x − 1) or (x + 1) respectively yielding a new cubic polynomial, which can be solved to find the quartic's other roots.

If ${\displaystyle a_{1}=a_{0}k}$ , ${\displaystyle a_{2}=0}$ and ${\displaystyle a_{4}=a_{3}k}$, then x = −k is a root of the equation. The full quartic can then be factorized this way:

${\displaystyle a_{0}x^{4}+a_{0}kx^{3}+a_{3}x+a_{3}k=a_{0}x^{3}(x+k)+a_{3}(x+k)=(a_{0}x^{3}+a_{3})(x+k).\,}$

If ${\displaystyle a_{1}=a_{0}k}$ , ${\displaystyle a_{3}=a_{2}k}$ and ${\displaystyle a_{4}=0}$, x = 0 and x = −k are two known roots. Q(x) divided by x(x + k) is a quadratic polynomial.

A quartic equation where a₃ and a₁ are equal to 0 takes the form

${\displaystyle a_{0}x^{4}+a_{2}x^{2}+a_{4}=0\,\!}$

and thus is a biquadratic equation, which is easy to solve: let ${\displaystyle z=x^{2}}$, so our equation turns to

${\displaystyle a_{0}z^{2}+a_{2}z+a_{4}=0\,\!}$

which is a simple quadratic equation, whose solutions are easily found using the quadratic formula:

${\displaystyle z={{-a_{2}\pm {\sqrt {a_{2}^{2}-4a_{0}a_{4}}}} \over {2a_{0}}}\,\!}$

When we've solved it (i.e. found these two z values), we can extract x from them

${\displaystyle x_{1}=+{\sqrt {z_{+}}}\,\!}$

${\displaystyle x_{2}=-{\sqrt {z_{+}}}\,\!}$

${\displaystyle x_{3}=+{\sqrt {z_{-}}}\,\!}$

${\displaystyle x_{4}=-{\sqrt {z_{-}}}\,\!}$

If either of the z solutions were negative or complex numbers, then some of the x solutions are complex numbers.

${\displaystyle a_{0}x^{4}+a_{1}x^{3}+a_{2}x^{2}+a_{1}mx+a_{0}m^{2}=0\,}$

Steps:

1) Divide by x ².

2) Use variable change z = x + m/x.

To begin, the quartic must first be converted to a depressed quartic.

Let

${\displaystyle Ax^{4}+Bx^{3}+Cx^{2}+Dx+E=0\qquad \qquad (1')}$

be the general quartic equation which it is desired to solve. Divide both sides by A,

${\displaystyle x^{4}+{B \over A}x^{3}+{C \over A}x^{2}+{D \over A}x+{E \over A}=0.}$

The first step should be to eliminate the x³ term. To do this, change variables from x to u, such that

${\displaystyle x=u-{B \over 4A}.}$

Then

${\displaystyle \left(u-{B \over 4A}\right)^{4}+{B \over A}\left(u-{B \over 4A}\right)^{3}+{C \over A}\left(u-{B \over 4A}\right)^{2}+{D \over A}\left(u-{B \over 4A}\right)+{E \over A}=0.}$

Expanding the powers of the binomials produces

${\displaystyle \left(u^{4}-{B \over A}u^{3}+{6u^{2}B^{2} \over 16A^{2}}-{4uB^{3} \over 64A^{3}}+{B^{4} \over 256A^{4}}\right)+{B \over A}\left(u^{3}-{3u^{2}B \over 4A}+{3uB^{2} \over 16A^{2}}-{B^{3} \over 64A^{3}}\right)+{C \over A}\left(u^{2}-{uB \over 2A}+{B^{2} \over 16A^{2}}\right)+{D \over A}\left(u-{B \over 4A}\right)+{E \over A}=0.}$

Collecting the same powers of u yields

${\displaystyle u^{4}+\left({-3B^{2} \over 8A^{2}}+{C \over A}\right)u^{2}+\left({B^{3} \over 8A^{3}}-{BC \over 2A^{2}}+{D \over A}\right)u+\left({-3B^{4} \over 256A^{4}}+{CB^{2} \over 16A^{3}}-{BD \over 4A^{2}}+{E \over A}\right)=0.}$

Now rename the coefficients of u. Let

${\displaystyle {\begin{aligned}\alpha &={-3B^{2} \over 8A^{2}}+{C \over A},\\\beta &={B^{3} \over 8A^{3}}-{BC \over 2A^{2}}+{D \over A},\\\gamma &={-3B^{4} \over 256A^{4}}+{CB^{2} \over 16A^{3}}-{BD \over 4A^{2}}+{E \over A}.\end{aligned}}}$

The resulting equation is

${\displaystyle u^{4}+\alpha u^{2}+\beta u+\gamma =0\qquad \qquad (1)}$

which is a depressed quartic equation.

If ${\displaystyle \beta =0\ }$ then we have a biquadratic equation, which (as explained above) is easily solved; using reverse substitution we can find our values for ${\displaystyle x}$.

If ${\displaystyle \gamma =0\ }$ then one of the roots is ${\displaystyle u=0\ }$ and the other roots can be found by dividing by ${\displaystyle u}$, and solving the resulting depressed cubic equation,

${\displaystyle u^{3}+\alpha u+\beta =0\,.}$

Using reverse substitution we can find our values for ${\displaystyle x}$.

Otherwise, the depressed quartic can be solved by means of a method discovered by Lodovico Ferrari. Once the depressed quartic has been obtained, the next step is to add the valid identity

${\displaystyle (u^{2}+\alpha )^{2}-u^{4}-2\alpha u^{2}=\alpha ^{2}\,}$

to equation (1), yielding

${\displaystyle (u^{2}+\alpha )^{2}+\beta u+\gamma =\alpha u^{2}+\alpha ^{2}.\qquad \qquad (2)}$

The effect has been to fold up the u⁴ term into a perfect square: (u² + α)². The second term, αu² did not disappear, but its sign has changed and it has been moved to the right side.

The next step is to insert a variable y into the perfect square on the left side of equation (2), and a corresponding 2y into the coefficient of u² in the right side. To accomplish these insertions, the following valid formulas will be added to equation (2),

${\displaystyle {\begin{aligned}(u^{2}+\alpha +y)^{2}-(u^{2}+\alpha )^{2}&=2y(u^{2}+\alpha )+y^{2}\ \ \\&=2yu^{2}+2y\alpha +y^{2},\end{aligned}}}$

and

${\displaystyle 0=(\alpha +2y)u^{2}-2yu^{2}-\alpha u^{2}\,}$

These two formulas, added together, produce

${\displaystyle (u^{2}+\alpha +y)^{2}-(u^{2}+\alpha )^{2}=(\alpha +2y)u^{2}-\alpha u^{2}+2y\alpha +y^{2}\qquad \qquad (y{\hbox{-insertion}})\,}$

which added to equation (2) produces

${\displaystyle (u^{2}+\alpha +y)^{2}+\beta u+\gamma =(\alpha +2y)u^{2}+(2y\alpha +y^{2}+\alpha ^{2}).\,}$

This is equivalent to

${\displaystyle (u^{2}+\alpha +y)^{2}=(\alpha +2y)u^{2}-\beta u+(y^{2}+2y\alpha +\alpha ^{2}-\gamma ).\qquad \qquad (3)\,}$

The objective now is to choose a value for y such that the right side of equation (3) becomes a perfect square. This can be done by letting the discriminant of the quadratic function become zero. To explain this, first expand a perfect square so that it equals a quadratic function:

${\displaystyle (su+t)^{2}=(s^{2})u^{2}+(2st)u+(t^{2}).\,}$

The quadratic function on the right side has three coefficients. It can be verified that squaring the second coefficient and then subtracting four times the product of the first and third coefficients yields zero:

${\displaystyle (2st)^{2}-4(s^{2})(t^{2})=0.\,}$

Therefore to make the right side of equation (3) into a perfect square, the following equation must be solved:

${\displaystyle (-\beta )^{2}-4(2y+\alpha )(y^{2}+2y\alpha +\alpha ^{2}-\gamma )=0.\,}$

Multiply the binomial with the polynomial,

${\displaystyle \beta ^{2}-4(2y^{3}+5\alpha y^{2}+(4\alpha ^{2}-2\gamma )y+(\alpha ^{3}-\alpha \gamma ))=0\,}$

Divide both sides by −4, and move the −β²/4 to the right,

${\displaystyle 2y^{3}+5\alpha y^{2}+(4\alpha ^{2}-2\gamma )y+\left(\alpha ^{3}-\alpha \gamma -{\beta ^{2} \over 4}\right)=0\qquad \qquad }$

This is a cubic equation for y. Divide both sides by 2,

${\displaystyle y^{3}+{5 \over 2}\alpha y^{2}+(2\alpha ^{2}-\gamma )y+\left({\alpha ^{3} \over 2}-{\alpha \gamma \over 2}-{\beta ^{2} \over 8}\right)=0.\qquad \qquad (4)}$

Equation (4) is a cubic equation nested within the quartic equation. It must be solved in order to solve the quartic. To solve the cubic, first transform it into a depressed cubic by means of the substitution

${\displaystyle y=v-{5 \over 6}\alpha .}$

Equation (4) becomes

${\displaystyle \left(v-{5 \over 6}\alpha \right)^{3}+{5 \over 2}\alpha \left(v-{5 \over 6}\alpha \right)^{2}+(2\alpha ^{2}-\gamma )\left(v-{5 \over 6}\alpha \right)+\left({\alpha ^{3} \over 2}-{\alpha \gamma \over 2}-{\beta ^{2} \over 8}\right)=0.}$

Expand the powers of the binomials,

${\displaystyle \left(v^{3}-{5 \over 2}\alpha v^{2}+{25 \over 12}\alpha ^{2}v-{125 \over 216}\alpha ^{3}\right)+{5 \over 2}\alpha \left(v^{2}-{5 \over 3}\alpha v+{25 \over 36}\alpha ^{2}\right)+(2\alpha ^{2}-\gamma )\left(v-{5 \over 6}\alpha \right)+\left({\alpha ^{3} \over 2}-{\alpha \gamma \over 2}-{\beta ^{2} \over 8}\right)=0.}$

Distribute, collect like powers of v, and cancel out the pair of v² terms,

${\displaystyle v^{3}+\left(-{\alpha ^{2} \over 12}-\gamma \right)v+\left(-{\alpha ^{3} \over 108}+{\alpha \gamma \over 3}-{\beta ^{2} \over 8}\right)=0.}$

This is a depressed cubic equation.

Relabel its coefficients,

${\displaystyle P=-{\alpha ^{2} \over 12}-\gamma ,}$

${\displaystyle Q=-{\alpha ^{3} \over 108}+{\alpha \gamma \over 3}-{\beta ^{2} \over 8}.}$

The depressed cubic now is

${\displaystyle v^{3}+Pv+Q=0.\qquad \qquad (5)}$

The solutions (any solution will do, so pick any of the three complex roots) of equation (5) are computed as (see Cubic equation)

${\displaystyle y=-{5 \over 6}\alpha +U+V\qquad \qquad (6)}$

where

${\displaystyle U={\sqrt[{3}]{-{Q \over 2}\pm {\sqrt {{Q^{2} \over 4}+{P^{3} \over 27}}}}}}$

and V is computed according to the two defining equations ${\displaystyle -U^{3}-V^{3}=Q}$ and ${\displaystyle -3UV=P}$, so

${\displaystyle V={\begin{cases}-{\frac {P}{3U}}&{\text{ if }}U\neq 0{\text{ and }}\\-{\sqrt[{3}]{Q}}&{\text{ if }}U=0\ .\end{cases}}}$

With the value for y given by equation (6), it is now known that the right side of equation (3) is a perfect square of the form

${\displaystyle (s^{2})u^{2}+(2st)u+(t^{2})=\left(\left({\sqrt {(s^{2})}}\right)u+{(2st) \over 2{\sqrt {(s^{2})}}}\right)^{2}}$

(This is correct for both signs of square root, as long as the same sign is taken for both square roots. A ± is redundant, as it would be absorbed by another ± a few equations further down this page.)

so that it can be folded:

${\displaystyle (\alpha +2y)u^{2}+(-\beta )u+(y^{2}+2y\alpha +\alpha ^{2}-\gamma )=\left(\left({\sqrt {(\alpha +2y)}}\right)u+{(-\beta ) \over 2{\sqrt {(\alpha +2y)}}}\right)^{2}.}$

Note: If β ≠ 0 then α + 2y ≠ 0. If β = 0 then this would be a biquadratic equation, which we solved earlier.

Therefore equation (3) becomes

${\displaystyle (u^{2}+\alpha +y)^{2}=\left(\left({\sqrt {\alpha +2y}}\right)u-{\beta \over 2{\sqrt {\alpha +2y}}}\right)^{2}.\qquad \qquad (7)}$

Equation (7) has a pair of folded perfect squares, one on each side of the equation. The two perfect squares balance each other.

If two squares are equal, then the sides of the two squares are also equal, as shown by:

${\displaystyle (u^{2}+\alpha +y)=\pm \left(\left({\sqrt {\alpha +2y}}\right)u-{\beta \over 2{\sqrt {\alpha +2y}}}\right).\qquad \qquad (7')}$

Collecting like powers of u produces

${\displaystyle u^{2}+\left(\mp _{s}{\sqrt {\alpha +2y}}\right)u+\left(\alpha +y\pm _{s}{\beta \over 2{\sqrt {\alpha +2y}}}\right)=0.\qquad \qquad (8)}$

Note: The subscript s of ${\displaystyle \pm _{s}}$ and ${\displaystyle \mp _{s}}$ is to note that they are dependent.

Equation (8) is a quadratic equation for u. Its solution is

${\displaystyle u={\pm _{s}{\sqrt {\alpha +2y}}\pm _{t}{\sqrt {(\alpha +2y)-4(\alpha +y\pm _{s}{\beta \over 2{\sqrt {\alpha +2y}}})}} \over 2}.}$

Simplifying, one gets

${\displaystyle u={\pm _{s}{\sqrt {\alpha +2y}}\pm _{t}{\sqrt {-\left(3\alpha +2y\pm _{s}{2\beta \over {\sqrt {\alpha +2y}}}\right)}} \over 2}.}$

This is the solution of the depressed quartic, therefore the solutions of the original quartic equation are

${\displaystyle x=-{B \over 4A}+{\pm _{s}{\sqrt {\alpha +2y}}\pm _{t}{\sqrt {-\left(3\alpha +2y\pm _{s}{2\beta \over {\sqrt {\alpha +2y}}}\right)}} \over 2}.\qquad \qquad (8')}$

Remember: The two ${\displaystyle \pm _{s}}$ come from the same place in equation (7'), and should both have the same sign, while the sign of ${\displaystyle \pm _{t}}$ is independent.

Given the quartic equation

${\displaystyle Ax^{4}+Bx^{3}+Cx^{2}+Dx+E=0,\,}$

its solution can be found by means of the following calculations:

${\displaystyle \alpha =-{3B^{2} \over 8A^{2}}+{C \over A},}$

${\displaystyle \beta ={B^{3} \over 8A^{3}}-{BC \over 2A^{2}}+{D \over A},}$

${\displaystyle \gamma =-{3B^{4} \over 256A^{4}}+{CB^{2} \over 16A^{3}}-{BD \over 4A^{2}}+{E \over A}.}$

If ${\displaystyle \,\beta =0,}$ then

${\displaystyle x=-{B \over 4A}\pm _{s}{\sqrt {-\alpha \pm _{t}{\sqrt {\alpha ^{2}-4\gamma }} \over 2}}\qquad {\mbox{(for }}\beta =0{\mbox{ only)}}.}$

Otherwise, continue with

${\displaystyle P=-{\alpha ^{2} \over 12}-\gamma ,}$

${\displaystyle Q=-{\alpha ^{3} \over 108}+{\alpha \gamma \over 3}-{\beta ^{2} \over 8},}$

${\displaystyle R=-{Q \over 2}\pm {\sqrt {{Q^{2} \over 4}+{P^{3} \over 27}}},}$

(either sign of the square root will do)

${\displaystyle U={\sqrt[{3}]{R}},}$

(there are 3 complex roots, any one of them will do)

${\displaystyle y=-{5 \over 6}\alpha +{\begin{cases}U=0&\to -{\sqrt[{3}]{Q}}\\U\neq 0,&\to U-{P \over 3U},\end{cases}}\quad \quad \quad }$

${\displaystyle W={\sqrt {\alpha +2y}}}$

${\displaystyle x=-{B \over 4A}+{\pm _{s}W\pm _{t}{\sqrt {-\left(3\alpha +2y\pm _{s}{2\beta \over W}\right)}} \over 2}.}$

The two ±_s must have the same sign, the ±_t is independent. To get all roots, compute x for ±_s,±_t = +,+ and for +,− and for −,+ and for −,−. This formula handles repeated roots without problem.

Ferrari was the first to discover one of these labyrinthine solutions. The equation which he solved was

${\displaystyle x^{4}+6x^{2}-60x+36=0}$

which was already in depressed form. It has a pair of solutions which can be found with the set of formulas shown above.

If the coefficients of the quartic equation are real then the nested depressed cubic equation (5) also has real coefficients, thus it has at least one real root.

Furthermore the cubic function

${\displaystyle C(v)=v^{3}+Pv+Q,}$

where P and Q are given by (5) has the properties that

${\displaystyle C\left({\alpha \over 3}\right)={-\beta ^{2} \over 8}<0\ }$ and

${\displaystyle \lim _{v\to \infty }C(v)=\infty ,}$ where α and β are given by (1).

This means that (5) has a real root greater than ${\displaystyle \alpha \over 3}$, and therefore that (4) has a real root greater than ${\displaystyle -\alpha \over 2}$.

Using this root the term ${\displaystyle {\sqrt {\alpha +2y}}}$ in (8) is always real, which ensures that the two quadratic equations (8) have real coefficients[2].

It could happen that one only obtained one solution through the seven formulae above, because not all four sign patterns are tried for four solutions, and the solution obtained is complex. It may also be the case that one is only looking for a real solution. Let x₁ denote the complex solution. If all the original coefficients A, B, C, D and E are real -- which should be the case when one desires only real solutions – then there is another complex solution x₂ which is the complex conjugate of x₁. If the other two roots are denoted as x₃ and x₄ then the quartic equation can be expressed as

${\displaystyle (x-x_{1})(x-x_{2})(x-x_{3})(x-x_{4})=0,\,}$

but this quartic equation is equivalent to the product of two quadratic equations:

${\displaystyle (x-x_{1})(x-x_{2})=0\qquad \qquad (9)}$

and

${\displaystyle (x-x_{3})(x-x_{4})=0.\qquad \qquad (10)}$

Since

${\displaystyle x_{2}=x_{1}^{\star }}$

then

${\displaystyle {\begin{aligned}(x-x_{1})(x-x_{2})&=x^{2}-(x_{1}+x_{1}^{\star })x+x_{1}x_{1}^{\star }\\&=x^{2}-2\operatorname {Re} (x_{1})x+[\operatorname {Re} (x_{1})]^{2}+[\operatorname {Im} (x_{1})]^{2}.\end{aligned}}}$

Let

${\displaystyle a=-2\operatorname {Re} (x_{1}),}$

${\displaystyle b=\left[\operatorname {Re} (x_{1})\right]^{2}+\left[\operatorname {Im} (x_{1})\right]^{2}}$

so that equation (9) becomes

${\displaystyle x^{2}+ax+b=0.\qquad \qquad (11)}$

Also let there be (unknown) variables w and v such that equation (10) becomes

${\displaystyle x^{2}+wx+v=0.\qquad \qquad (12)}$

Multiplying equations (11) and (12) produces

${\displaystyle x^{4}+(a+w)x^{3}+(b+wa+v)x^{2}+(wb+va)x+vb=0.\qquad \qquad (13)}$

Comparing equation (13) to the original quartic equation, it can be seen that

${\displaystyle a+w={B \over A},}$

${\displaystyle b+wa+v={C \over A},}$

${\displaystyle wb+va={D \over A},}$

and

${\displaystyle vb={E \over A}.}$

Therefore

${\displaystyle w={B \over A}-a={B \over A}+2\operatorname {Re} (x_{1}),}$

${\displaystyle v={E \over Ab}={E \over A\left([\operatorname {Re} (x_{1})]^{2}+[\operatorname {Im} (x_{1})]^{2}\right)}.}$

Equation (12) can be solved for x yielding

${\displaystyle x_{3}={-w+{\sqrt {w^{2}-4v}} \over 2},}$

${\displaystyle x_{4}={-w-{\sqrt {w^{2}-4v}} \over 2}.}$

One of these two solutions should be the desired real solution.

Most textbook solutions of the quartic equation require a magic substitution that is almost impossible to memorize. Here is a way to approach it that makes it easy to understand.

The job is done if we can factor the quartic equation into a product of two quadratics. Let

${\displaystyle {\begin{aligned}0=x^{4}+bx^{3}+cx^{2}+dx+e&=(x^{2}+px+q)(x^{2}+rx+s)\\&=x^{4}+(p+r)x^{3}+(q+s+pr)x^{2}+(ps+qr)x+qs\end{aligned}}}$

By equating coefficients, this results in the following set of simultaneous equations:

${\displaystyle {\begin{aligned}b&=p+r\\c&=q+s+pr\\d&=ps+qr\\e&=qs\end{aligned}}}$

This is harder to solve than it looks, but if we start again with a depressed quartic where ${\displaystyle b=0}$, which can be obtained by substituting ${\displaystyle (x-b/4)}$ for ${\displaystyle x}$, then ${\displaystyle r=-p}$, and:

${\displaystyle {\begin{aligned}c+p^{2}&=s+q\\d/p&=s-q\\e&=sq\end{aligned}}}$

It's now easy to eliminate both ${\displaystyle s}$ and ${\displaystyle q}$ by doing the following:

${\displaystyle {\begin{aligned}(c+p^{2})^{2}-(d/p)^{2}&=(s+q)^{2}-(s-q)^{2}\\&=4sq\\&=4e\end{aligned}}}$

If we set ${\displaystyle P=p^{2}}$, then this equation turns into the cubic equation:

${\displaystyle P^{3}+2cP^{2}+(c^{2}-4e)P-d^{2}=0\,}$

which is solved elsewhere. Once you have ${\displaystyle p}$, then:

${\displaystyle {\begin{aligned}r&=-p\\2s&=c+p^{2}+d/p\\2q&=c+p^{2}-d/p\end{aligned}}}$

The symmetries in this solution are easy to see. There are three roots of the cubic, corresponding to the three ways that a quartic can be factored into two quadratics, and choosing positive or negative values of ${\displaystyle p}$ for the square root of ${\displaystyle P}$ merely exchanges the two quadratics with one another.

The symmetric group S₄ on four elements has the Klein four-group as a normal subgroup. This suggests using a resolvent whose roots may be variously described as a discrete Fourier transform or a Hadamard matrix transform of the roots. Suppose r_i for i from 0 to 3 are roots of

${\displaystyle x^{4}+bx^{3}+cx^{2}+dx+e=0\qquad (1)}$

If we now set

${\displaystyle {\begin{aligned}s_{0}&={\tfrac {1}{2}}(r_{0}+r_{1}+r_{2}+r_{3}),\\s_{1}&={\tfrac {1}{2}}(r_{0}-r_{1}+r_{2}-r_{3}),\\s_{2}&={\tfrac {1}{2}}(r_{0}+r_{1}-r_{2}-r_{3}),\\s_{3}&={\tfrac {1}{2}}(r_{0}-r_{1}-r_{2}+r_{3}),\end{aligned}}}$

then since the transformation is an involution we may express the roots in terms of the four s_i in exactly the same way. Since we know the value s₀ = -b/2, we really only need the values for s₁, s₂ and s₃. These we may find by expanding the polynomial

${\displaystyle (z^{2}-s_{1}^{2})(z^{2}-s_{2}^{2})(z^{2}-s_{3}^{2})\qquad (2)}$

which if we make the simplifying assumption that b = 0, is equal to

${\displaystyle z^{6}+2cz^{4}+(c^{2}-4e)z^{2}-d^{2}\qquad (3)}$

This polynomial is of degree six, but only of degree three in z², and so the corresponding equation is solvable. By trial we can determine which three roots are the correct ones, and hence find the solutions of the quartic.

We can remove any requirement for trial by using a root of the same resolvent polynomial for factoring; if w is any root of (3), and if

${\displaystyle F_{1}=x^{2}+wx+{\frac {1}{2}}w^{2}+{\frac {1}{2}}c-{\frac {1}{2}}\cdot {\frac {c^{2}w}{d}}-{\frac {1}{2}}\cdot {\frac {w^{5}}{d}}-{\frac {cw^{3}}{d}}+2{\frac {ew}{d}}}$

${\displaystyle F_{2}=x^{2}-wx+{\frac {1}{2}}w^{2}+{\frac {1}{2}}c+{\frac {1}{2}}\cdot {\frac {w^{5}}{d}}+{\frac {cw^{3}}{d}}-2{\frac {ew}{d}}+{\frac {1}{2}}\cdot {\frac {c^{2}w}{d}}}$

then

${\displaystyle F_{1}F_{2}=x^{4}+cx^{2}+dx+e\qquad \qquad (4)}$

We therefore can solve the quartic by solving for w and then solving for the roots of the two factors using the quadratic formula.