Quartic function

From testwiki
Jump to navigation Jump to search

Template:Short description Template:Distinguish Template:About Template:Redirect Template:Use dmy dates

Graph of a polynomial of degree 4, with 3 critical points and four real roots (crossings of the x axis) (and thus no complex roots). If one or the other of the local minima were above the x axis, or if the local maximum were below it, or if there were no local maximum and one minimum below the x axis, there would only be two real roots (and two complex roots). If all three local extrema were above the x axis, or if there were no local maximum and one minimum above the x axis, there would be no real root (and four complex roots). The same reasoning applies in reverse to polynomial with a negative quartic coefficient.

In algebra, a quartic function is a function of the form

f(x)=ax4+bx3+cx2+dx+e,Template:Ref

where a is nonzero, which is defined by a polynomial of degree four, called a quartic polynomial.

A quartic equation, or equation of the fourth degree, is an equation that equates a quartic polynomial to zero, of the form

ax4+bx3+cx2+dx+e=0,

where Template:Nowrap.[1] The derivative of a quartic function is a cubic function.

Sometimes the term biquadratic is used instead of quartic, but, usually, biquadratic function refers to a quadratic function of a square (or, equivalently, to the function defined by a quartic polynomial without terms of odd degree), having the form

f(x)=ax4+cx2+e.

Since a quartic function is defined by a polynomial of even degree, it has the same infinite limit when the argument goes to positive or negative infinity. If a is positive, then the function increases to positive infinity at both ends; and thus the function has a global minimum. Likewise, if a is negative, it decreases to negative infinity and has a global maximum. In both cases it may or may not have another local maximum and another local minimum.

The degree four (quartic case) is the highest degree such that every polynomial equation can be solved by radicals, according to the Abel–Ruffini theorem.

History

Lodovico Ferrari is credited with the discovery of the solution to the quartic in 1540, but since this solution, like all algebraic solutions of the quartic, requires the solution of a cubic to be found, it could not be published immediately.[2] The solution of the quartic was published together with that of the cubic by Ferrari's mentor Gerolamo Cardano in the book Ars Magna.[3]

The proof that four is the highest degree of a general polynomial for which such solutions can be found was first given in the Abel–Ruffini theorem in 1824, proving that all attempts at solving the higher order polynomials would be futile. The notes left by Évariste Galois prior to dying in a duel in 1832 later led to an elegant complete theory of the roots of polynomials, of which this theorem was one result.[4]

Applications

Each coordinate of the intersection points of two conic sections is a solution of a quartic equation. The same is true for the intersection of a line and a torus. It follows that quartic equations often arise in computational geometry and all related fields such as computer graphics, computer-aided design, computer-aided manufacturing and optics. Here are examples of other geometric problems whose solution involves solving a quartic equation.

In computer-aided manufacturing, the torus is a shape that is commonly associated with the endmill cutter. To calculate its location relative to a triangulated surface, the position of a horizontal torus on the Template:Math-axis must be found where it is tangent to a fixed line, and this requires the solution of a general quartic equation to be calculated.[5]

A quartic equation arises also in the process of solving the crossed ladders problem, in which the lengths of two crossed ladders, each based against one wall and leaning against another, are given along with the height at which they cross, and the distance between the walls is to be found.[6]

In optics, Alhazen's problem is "Given a light source and a spherical mirror, find the point on the mirror where the light will be reflected to the eye of an observer." This leads to a quartic equation.[7][8][9]

Finding the distance of closest approach of two ellipses involves solving a quartic equation.

The eigenvalues of a 4×4 matrix are the roots of a quartic polynomial which is the characteristic polynomial of the matrix.

The characteristic equation of a fourth-order linear difference equation or differential equation is a quartic equation. An example arises in the Timoshenko-Rayleigh theory of beam bending.[10]

Intersections between spheres, cylinders, or other quadrics can be found using quartic equations.

Inflection points and golden ratio

Letting Template:Mvar and Template:Mvar be the distinct inflection points of the graph of a quartic function, and letting Template:Mvar be the intersection of the inflection secant line Template:Mvar and the quartic, nearer to Template:Mvar than to Template:Mvar, then Template:Mvar divides Template:Mvar into the golden section:[11]

FGGH=1+52=φ(the golden ratio).

Moreover, the area of the region between the secant line and the quartic below the secant line equals the area of the region between the secant line and the quartic above the secant line. One of those regions is disjointed into sub-regions of equal area.

Solution

Nature of the roots

Given the general quartic equation

ax4+bx3+cx2+dx+e=0

with real coefficients and Template:Math the nature of its roots is mainly determined by the sign of its discriminant

Δ=256a3e3192a2bde2128a2c2e2+144a2cd2e27a2d4+144ab2ce26ab2d2e80abc2de+18abcd3+16ac4e4ac3d227b4e2+18b3cde4b3d34b2c3e+b2c2d2

This may be refined by considering the signs of four other polynomials:

P=8ac3b2

such that Template:Math is the second degree coefficient of the associated depressed quartic (see below);

R=b3+8da24abc,

such that Template:Math is the first degree coefficient of the associated depressed quartic;

Δ0=c23bd+12ae,

which is 0 if the quartic has a triple root; and

D=64a3e16a2c2+16ab2c16a2bd3b4

which is 0 if the quartic has two double roots.

The possible cases for the nature of the roots are as follows:[12]

There are some cases that do not seem to be covered, but in fact they cannot occur. For example, Template:Math, Template:Mvar = 0 and Template:Mvar ≤ 0 is not one of the cases. In fact, if Template:Math and Template:Mvar = 0 then Template:Mvar > 0, since 16a2Δ0=3D+P2; so this combination is not possible.

General formula for roots

Solution of x4+ax3+bx2+cx+d=0 written out in full. This formula is too unwieldy for general use; hence other methods, or simpler formulas for special cases, are generally used.

The four roots Template:Math, Template:Math, Template:Math, and Template:Math for the general quartic equation

ax4+bx3+cx2+dx+e=0

with Template:Mvar ≠ 0 are given in the following formula, which is deduced from the one in the section on Ferrari's method by back changing the variables (see Template:Slink) and using the formulas for the quadratic and cubic equations.

x1,2 =b4aS±124S22p+qSx3,4 =b4a+S±124S22pqS

where Template:Mvar and Template:Mvar are the coefficients of the second and of the first degree respectively in the associated depressed quartic

p=8ac3b28a2q=b34abc+8a2d8a3

and where

S=1223 p+13a(Q+Δ0Q)Q=Δ1+Δ124Δ0323

(if Template:Math or Template:Math, see Template:Slink, below)

with

Δ0=c23bd+12aeΔ1=2c39bcd+27b2e+27ad272ace

and

Δ124Δ03=27Δ , where Δ is the aforementioned discriminant. For the cube root expression for Q, any of the three cube roots in the complex plane can be used, although if one of them is real that is the natural and simplest one to choose. The mathematical expressions of these last four terms are very similar to those of their cubic counterparts.

Special cases of the formula

  • If Δ>0, the value of Q is a non-real complex number. In this case, either all roots are non-real or they are all real. In the latter case, the value of S is also real, despite being expressed in terms of Q; this is casus irreducibilis of the cubic function extended to the present context of the quartic. One may prefer to express it in a purely real way, by using trigonometric functions, as follows:
S=1223 p+23aΔ0cosφ3
where
φ=arccos(Δ12Δ03).
  • If Δ0 and Δ0=0, the sign of Δ124Δ03=Δ12 has to be chosen to have Q0, that is one should define Δ12 as Δ1, maintaining the sign of Δ1.
  • If S=0, then one must change the choice of the cube root in Q in order to have S0. This is always possible except if the quartic may be factored into (x+b4a)4. The result is then correct, but misleading because it hides the fact that no cube root is needed in this case. In fact this caseTemplate:Clarify may occur only if the numerator of q is zero, in which case the associated depressed quartic is biquadratic; it may thus be solved by the method described below.
  • If Δ=0 and Δ0=0, and thus also Δ1=0, at least three roots are equal to each other, and the roots are rational functions of the coefficients. The triple root x0 is a common root of the quartic and its second derivative 2(6ax2+3bx+c); it is thus also the unique root of the remainder of the Euclidean division of the quartic by its second derivative, which is a linear polynomial. The simple root x1 can be deduced from x1+3x0=b/a.
  • If Δ=0 and Δ00, the above expression for the roots is correct but misleading, hiding the fact that the polynomial is reducible and no cube root is needed to represent the roots.

Simpler cases

Reducible quartics

Consider the general quartic

Q(x)=a4x4+a3x3+a2x2+a1x+a0.

It is reducible if Template:Math, where Template:Math and Template:Math are non-constant polynomials with rational coefficients (or more generally with coefficients in the same field as the coefficients of Template:Math). Such a factorization will take one of two forms:

Q(x)=(xx1)(b3x3+b2x2+b1x+b0)

or

Q(x)=(c2x2+c1x+c0)(d2x2+d1x+d0).

In either case, the roots of Template:Math are the roots of the factors, which may be computed using the formulas for the roots of a quadratic function or cubic function.

Detecting the existence of such factorizations can be done [[Resolvent cubic#Factoring quartic polynomials|using the resolvent cubic of Template:Math]]. It turns out that:

  • if we are working over Template:Math (that is, if coefficients are restricted to be real numbers) (or, more generally, over some real closed field) then there is always such a factorization;
  • if we are working over Template:Math (that is, if coefficients are restricted to be rational numbers) then there is an algorithm to determine whether or not Template:Math is reducible and, if it is, how to express it as a product of polynomials of smaller degree.

In fact, several methods of solving quartic equations (Ferrari's method, Descartes' method, and, to a lesser extent, Euler's method) are based upon finding such factorizations.

Biquadratic equation

If Template:Math then the function

Q(x)=a4x4+a2x2+a0

is called a biquadratic function; equating it to zero defines a biquadratic equation, which is easy to solve as follows

Let the auxiliary variable Template:Math. Then Template:Math becomes a quadratic Template:Math in Template:Math: Template:Math. Let Template:Math and Template:Math be the roots of Template:Math. Then the roots of the quartic Template:Math are

x1=+z+,x2=z+,x3=+z,x4=z.

Quasi-palindromic equation

The polynomial

P(x)=a0x4+a1x3+a2x2+a1mx+a0m2

is almost palindromic, as Template:Math (it is palindromic if Template:Math). The change of variables Template:Math in Template:Math produces the quadratic equation Template:Math. Since Template:Math, the quartic equation Template:Math may be solved by applying the quadratic formula twice.

Solution methods

Converting to a depressed quartic

For solving purposes, it is generally better to convert the quartic into a depressed quartic by the following simple change of variable. All formulas are simpler and some methods work only in this case. The roots of the original quartic are easily recovered from that of the depressed quartic by the reverse change of variable.

Let

a4x4+a3x3+a2x2+a1x+a0=0

be the general quartic equation we want to solve.

Dividing by Template:Math, provides the equivalent equation Template:Math, with Template:Math, Template:Math, Template:Math, and Template:Math. Substituting Template:Math for Template:Mvar gives, after regrouping the terms, the equation Template:Math, where

p=8c3b28=8a2a43a328a42q=b34bc+8d8=a334a2a3a4+8a1a428a43r=3b4+256e64bd+16b2c256=3a34+256a0a4364a1a3a42+16a2a32a4256a44.

If Template:Math is a root of this depressed quartic, then Template:Math (that is Template:Math is a root of the original quartic and every root of the original quartic can be obtained by this process.

Ferrari's solution

As explained in the preceding section, we may start with the depressed quartic equation

y4+py2+qy+r=0.

This depressed quartic can be solved by means of a method discovered by Lodovico Ferrari. The depressed equation may be rewritten (this is easily verified by expanding the square and regrouping all terms in the left-hand side) as

(y2+p2)2=qyr+p24.

Then, we introduce a variable Template:Mvar into the factor on the left-hand side by adding Template:Math to both sides. After regrouping the coefficients of the power of Template:Mvar on the right-hand side, this gives the equation Template:NumBlk which is equivalent to the original equation, whichever value is given to Template:Mvar.

As the value of Template:Mvar may be arbitrarily chosen, we will choose it in order to complete the square on the right-hand side. This implies that the discriminant in Template:Mvar of this quadratic equation is zero, that is Template:Mvar is a root of the equation

(q)24(2m)(m2+pm+p24r)=0,

which may be rewritten as

Template:NumBlk

This is the resolvent cubic of the quartic equation. The value of Template:Mvar may thus be obtained from Cardano's formula. When Template:Mvar is a root of this equation, the right-hand side of equation (Template:EquationNote) is the square

(2myq22m)2.

However, this induces a division by zero if Template:Math. This implies Template:Math, and thus that the depressed equation is bi-quadratic, and may be solved by an easier method (see above). This was not a problem at the time of Ferrari, when one solved only explicitly given equations with numeric coefficients. For a general formula that is always true, one thus needs to choose a root of the cubic equation such that Template:Math. This is always possible except for the depressed equation Template:Math.

Now, if Template:Mvar is a root of the cubic equation such that Template:Math, equation (Template:EquationNote) becomes

(y2+p2+m)2=(y2mq22m)2.

This equation is of the form Template:Math, which can be rearranged as Template:Math or Template:Math. Therefore, equation (Template:EquationNote) may be rewritten as

(y2+p2+m+2myq22m)(y2+p2+m2my+q22m)=0.

This equation is easily solved by applying to each factor the quadratic formula. Solving them we may write the four roots as

y=±12m±2(2p+2m±12qm)2,

where Template:Math and Template:Math denote either Template:Math or Template:Math. As the two occurrences of Template:Math must denote the same sign, this leaves four possibilities, one for each root.

Therefore, the solutions of the original quartic equation are

x=a34a4+±12m±2(2p+2m±12qm)2.

A comparison with the general formula above shows that Template:Math.

Descartes' solution

Descartes[14] introduced in 1637 the method of finding the roots of a quartic polynomial by factoring it into two quadratic ones. Let

x4+bx3+cx2+dx+e=(x2+sx+t)(x2+ux+v)=x4+(s+u)x3+(t+v+su)x2+(sv+tu)x+tv

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

{b=s+uc=t+v+sud=sv+tue=tv

This can be simplified by starting again with the depressed quartic Template:Math, which can be obtained by substituting Template:Math for Template:Math. Since the coefficient of Template:Math is Template:Math, we get Template:Math, and:

{p+u2=t+vq=u(tv)r=tv

One can now eliminate both Template:Mvar and Template:Mvar by doing the following:

u2(p+u2)2q2=u2(t+v)2u2(tv)2=u2[(t+v+(tv))(t+v(tv))]=u2(2t)(2v)=4u2tv=4u2r

If we set Template:Math, then solving this equation becomes finding the roots of the resolvent cubic

Template:NumBlk

which is done elsewhere. This resolvent cubic is equivalent to the resolvent cubic given above (equation (1a)), as can be seen by substituting U = 2m.

If Template:Math is a square root of a non-zero root of this resolvent (such a non-zero root exists except for the quartic Template:Math, which is trivially factored),

{s=u2t=p+u2+q/u2v=p+u2q/u

The symmetries in this solution are as follows. 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 Template:Mvar for the square root of Template:Mvar merely exchanges the two quadratics with one another.

The above solution shows that a quartic polynomial with rational coefficients and a zero coefficient on the cubic term is factorable into quadratics with rational coefficients if and only if either the resolvent cubic (Template:EquationNote) has a non-zero root which is the square of a rational, or Template:Math is the square of rational and Template:Math; this can readily be checked using the rational root test.[15]

Euler's solution

A variant of the previous method is due to Euler.[16][17] Unlike the previous methods, both of which use some root of the resolvent cubic, Euler's method uses all of them. Consider a depressed quartic Template:Math. Observe that, if

then

Therefore, Template:Math. In other words, Template:Math is one of the roots of the resolvent cubic (Template:EquationNote) and this suggests that the roots of that cubic are equal to Template:Math, Template:Math, and Template:Math. This is indeed true and it follows from Vieta's formulas. It also follows from Vieta's formulas, together with the fact that we are working with a depressed quartic, that Template:Math. (Of course, this also follows from the fact that Template:Math.) Therefore, if Template:Math, Template:Math, and Template:Math are the roots of the resolvent cubic, then the numbers Template:Math, Template:Math, Template:Math, and Template:Math are such that

{r1+r2+r3+r4=0(r1+r2)(r3+r4)=α(r1+r3)(r2+r4)=β(r1+r4)(r2+r3)=γ.

It is a consequence of the first two equations that Template:Math is a square root of Template:Math and that Template:Math is the other square root of Template:Math. For the same reason,

Therefore, the numbers Template:Math, Template:Math, Template:Math, and Template:Math are such that

{r1+r2+r3+r4=0r1+r2=αr1+r3=βr1+r4=γ;

the sign of the square roots will be dealt with below. The only solution of this system is:

{r1=α+β+γ2r2=αβγ2r3=α+βγ2r4=αβ+γ2.

Since, in general, there are two choices for each square root, it might look as if this provides Template:Math choices for the set Template:Math}, but, in fact, it provides no more than Template:Math such choices, because the consequence of replacing one of the square roots by the symmetric one is that the set Template:Math} becomes the set Template:Math}.

In order to determine the right sign of the square roots, one simply chooses some square root for each of the numbers Template:Math, Template:Math, and Template:Math and uses them to compute the numbers Template:Math, Template:Math, Template:Math, and Template:Math from the previous equalities. Then, one computes the number Template:Math. Since Template:Math, Template:Math, and Template:Math are the roots of (Template:EquationNote), it is a consequence of Vieta's formulas that their product is equal to Template:Math and therefore that Template:Math. But a straightforward computation shows that

Template:Math

If this number is Template:Math, then the choice of the square roots was a good one (again, by Vieta's formulas); otherwise, the roots of the polynomial will be Template:Math, Template:Math, Template:Math, and Template:Math, which are the numbers obtained if one of the square roots is replaced by the symmetric one (or, what amounts to the same thing, if each of the three square roots is replaced by the symmetric one).

This argument suggests another way of choosing the square roots:

Of course, this will make no sense if Template:Math or Template:Math is equal to Template:Math, but Template:Math is a root of (Template:EquationNote) only when Template:Math, that is, only when we are dealing with a biquadratic equation, in which case there is a much simpler approach.

Solving by Lagrange resolvent

The symmetric group Template:Math on four elements has the Klein four-group as a normal subgroup. This suggests using a Template:Visible anchor whose roots may be variously described as a discrete Fourier transform or a Hadamard matrix transform of the roots; see Lagrange resolvents for the general method. Denote by Template:Math, for Template:Math from Template:Math to Template:Math, the four roots of Template:Math. If we set

s0=12(x0+x1+x2+x3),s1=12(x0x1+x2x3),s2=12(x0+x1x2x3),s3=12(x0x1x2+x3),

then since the transformation is an involution we may express the roots in terms of the four Template:Math in exactly the same way. Since we know the value Template:Math, we only need the values for Template:Math, Template:Math and Template:Math. These are the roots of the polynomial

(s2s12)(s2s22)(s2s32).

Substituting the Template:Math by their values in term of the Template:Math, this polynomial may be expanded in a polynomial in Template:Math whose coefficients are symmetric polynomials in the Template:Math. By the fundamental theorem of symmetric polynomials, these coefficients may be expressed as polynomials in the coefficients of the monic quartic. If, for simplification, we suppose that the quartic is depressed, that is Template:Math, this results in the polynomial Template:NumBlk This polynomial is of degree six, but only of degree three in Template:Math, and so the corresponding equation is solvable by the method described in the article about cubic function. By substituting the roots in the expression of the Template:Math in terms of the Template:Math, we obtain expression for the roots. In fact we obtain, apparently, several expressions, depending on the numbering of the roots of the cubic polynomial and of the signs given to their square roots. All these different expressions may be deduced from one of them by simply changing the numbering of the Template:Math.

These expressions are unnecessarily complicated, involving the cubic roots of unity, which can be avoided as follows. If Template:Math is any non-zero root of (Template:EquationNote), and if we set

F1(x)=x2+sx+c2+s22d2sF2(x)=x2sx+c2+s22+d2s

then

F1(x)×F2(x)=x4+cx2+dx+e.

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

This gives exactly the same formula for the roots as the one provided by Descartes' method.

Solving with algebraic geometry

There is an alternative solution using algebraic geometry[18] In brief, one interprets the roots as the intersection of two quadratic curves, then finds the three reducible quadratic curves (pairs of lines) that pass through these points (this corresponds to the resolvent cubic, the pairs of lines being the Lagrange resolvents), and then use these linear equations to solve the quadratic.

The four roots of the depressed quartic Template:Math may also be expressed as the Template:Mvar coordinates of the intersections of the two quadratic equations Template:Math and Template:Math i.e., using the substitution Template:Math that two quadratics intersect in four points is an instance of Bézout's theorem. Explicitly, the four points are Template:Math for the four roots Template:Math of the quartic.

These four points are not collinear because they lie on the irreducible quadratic Template:Math and thus there is a 1-parameter family of quadratics (a pencil of curves) passing through these points. Writing the projectivization of the two quadratics as quadratic forms in three variables:

F1(X,Y,Z):=Y2+pYZ+qXZ+rZ2,F2(X,Y,Z):=YZX2

the pencil is given by the forms Template:Math for any point Template:Math in the projective line — in other words, where Template:Math and Template:Math are not both zero, and multiplying a quadratic form by a constant does not change its quadratic curve of zeros.

This pencil contains three reducible quadratics, each corresponding to a pair of lines, each passing through two of the four points, which can be done (42) = Template:Math different ways. Denote these Template:Math, Template:Math, and Template:Math. Given any two of these, their intersection has exactly the four points.

The reducible quadratics, in turn, may be determined by expressing the quadratic form Template:Math as a Template:Math matrix: reducible quadratics correspond to this matrix being singular, which is equivalent to its determinant being zero, and the determinant is a homogeneous degree three polynomial in Template:Math and Template:Math and corresponds to the resolvent cubic.

See also

Notes

Template:Note For the purposes of this article, e is used as a variable as opposed to its conventional use as Euler's number (except when otherwise specified).

References

Template:Reflist

Further reading

Template:Polynomials