Cubic equation

From Wikipedia, the free encyclopedia

This article discusses cubic equations in one variable. For a discussion of cubic equations in two variables, see elliptic curve.
Graph of a cubic polynomial:y = x3/4 + 3x2/4 − 3x/2 − 2 = (1/4)(x + 4)(x + 1)(x − 2)
Graph of a cubic polynomial:
y = x3/4 + 3x2/4 − 3x/2 − 2
= (1/4)(x + 4)(x + 1)(x − 2)

In mathematics, a cubic equation is a polynomial equation in which the highest occurring power of the unknown is the third power. The standard form of a cubic equation is

ax^3+bx^2+cx+d=0. \,

Usually, the coefficients a,..., d are real numbers. However, most of the theory is also valid if they belong to a field of characteristic other than two or three. We will always assume that a is non-zero (otherwise it is a quadratic equation).

Solving a cubic equation amounts to finding the roots of a cubic function.

Contents

[edit] History

Cubic equations were discovered by the ancient civilizations more or less independently. The Jaina mathematicians in ancient India are sometimes credited with the discovery sometime between 400 BCE and 200 CE. In the Greek world, it was known by the fifth century BCE with problems such as doubling the cube.

According to some books on the history of mathematics, Menaechmus and Archimedes invented the method of solving cubic equations by intersecting conic sections, though historians such as Reviel Netz disputed whether the Greeks were thinking about cubic equations or just problems that can lead to cubic equations. The Persian mathematician Omar Khayyám (10481131) re-discovered the method of constructing solutions of cubic equations by intersecting a parabola with a circle. He showed how this geometric solution could be used to get a numerical answer by consulting trigonometric tables.

In the early 16th century, the Italian mathematician Scipione del Ferro (1465-1526) found a method for solving a class of cubic equations, namely those of the form x3 + mx = n. In fact, all cubic equations can be reduced to this form if we allow m and n to be negative, but negative numbers were not known at that time. Del Ferro kept his achievement secret until just before his death, when he told his student Antonio Fiore about it.

In 1530, Niccolò Tartaglia (1500-1557) received two problems in cubic equations from Zuanne da Coi and announced that he could solve them. He was soon challenged by Fiore, which led to a famous contest between the two. Each contestant had to put up a certain amount of money and to propose a number of problems for his rival to solve. Whoever solved more problems within 30 days would get all the money.

Tartaglia received questions in the form x3 + mx = n, for which he had worked out a general method. Fiore received questions in the form x3 + mx2 = n, which proved to be too difficult for him to solve, and Tartaglia won the contest.

Later, Tartaglia was persuaded by Gerolamo Cardano (1501-1576) to reveal his secret for solving cubic equations. Tartaglia did so only on the condition that Cardano would never reveal it. A few years later, Cardano learned about Ferro's prior work and broke the promise by publishing Tartaglia's method in his book Ars Magna (1545) with credit given to Tartaglia. This led to another competition between Tartaglia and Cardano, for which the latter did not show up but was represented by his student Lodovico Ferrari (1522-1565). Ferrari did better than Tartaglia in the competition, and Tartaglia lost both his prestige and income.

Cardano noticed that Tartaglia's method sometimes required him to extract the square root of a negative number. He even included a calculation with these complex numbers in Ars Magna, but he did not really understand it. Rafael Bombelli studied this issue in detail and is therefore often considered as the discoverer of complex numbers.

[edit] The nature of the roots

Every cubic equation with real coefficients has at least one solution x among the real numbers; this is a consequence of the intermediate value theorem. We can distinguish several possible cases using the discriminant,

\Delta = 4b^3d - b^2c^2 + 4ac^3 - 18abcd + 27a^2d^2. \,

The following cases need to be considered.

  • If Δ < 0, then the equation has three distinct real roots.
  • If Δ > 0, then the equation has one real root and a pair of complex conjugate roots.
  • If Δ = 0, then (at least) two roots coincide. It may be that the equation has a double real root and another distinct single real root; alternatively, all three roots coincide yielding a triple real root. A possible way to decide between these subcases is to compute the resultant of the cubic and its second derivative: a triple root exists if and only if this resultant vanishes.


See also: multiplicity of a root of a polynomial

[edit] Cardano's method

The solutions can be found with the following method due to Scipione del Ferro and Tartaglia, published by Gerolamo Cardano in 1545.

We first divide the standard equation by the leading x3 coefficient to arrive at an equation of the form

x^3 + ax^2 + bx +c = 0. \quad (1)

The substitution x = t - {a\over 3}eliminates the quadratic term; in fact, we get the equation

t^3 + pt + q = 0, \quad\mbox{where } p = b - \frac{a^2}3 \quad\mbox{and}\quad q = c + \frac{2a^3-9ab}{27}. \qquad (2)

This is called the depressed cubic.

Suppose that we can find numbers u and v such that

u^3-v^3 = q \quad\mbox{and}\quad uv = \frac{p}{3}. \quad (3)

A solution to our equation is then given by

t = v - u, \,

as can be checked by directly substituting this value for t in (2), as a consequence of the third order binomial identity

(v-u)^3+3uv(v-u)+(u^3-v^3)=0  \ .

The system (3) can be solved by solving the second equation for v, which gives

v = \frac{p}{3u}.

Substituting this into the first equation in (3) yields

u^3 - \frac{p^3}{27u^3} = q.

Moving the q to the other side and multiplying by 27u3 yields

27u^6 - 27qu^3 - p^3 = 0\,.

This can be seen as a quadratic equation for u3. If we solve this equation, we find that

u^{3}={q\over 2}\pm \sqrt{{q^{2}\over 4}+{p^{3}\over 27}}
u=\sqrt[3]{{q\over 2}\pm \sqrt{{q^{2}\over 4}+{p^{3}\over 27}}}. \quad (4)

Since t = vu, t = x + a/3, and v = p/3u, we find

x=\frac{p}{3u}-u-{a\over 3}.

Note that there are six possibilities in computing u with (4), since there are two solutions to the square root (\pm), and three complex solutions to the cubic root (the principal root and the principal root multiplied by {-1\over 2} \pm i{\sqrt{3}\over 2}. However, the sign of the square root (plus or minus) does not affect the final resulting x, although care must be taken in two special cases to avoid divisions by zero. First, if p = 0, then one should choose the positive square root so that u does not equal zero, i.e., u = +\sqrt[3]{q}. Second, if p = q = 0, then we have the triple real root x = −a/3.

In summary, for the cubic equation

x^3 + ax^2 + bx +c = 0,\

the solutions for x are given by

x=\frac{p}{3u}-u-{a\over 3}

where

p = b - \frac{a^2}3
q = c + \frac{2a^3-9ab}{27}
u=\sqrt[3]{{q\over 2}\pm \sqrt{{q^{2}\over 4}+{p^{3}\over 27}}}.

[edit] Lagrange resolvents

The symmetric group S3 of order three has the cyclic group of order three as a normal subgroup, which suggests making use of the discrete Fourier transform of the roots, an idea due to Lagrange. Suppose that r0, r1 and r2 are the roots of equation (1), and define \zeta = (-1+i\sqrt{3})/2, so that ζ is a primitive third root of unity. We now set

s_0 = r_0 + r_1 + r_2,\,
s_1 = r_0 + \zeta r_1 + \zeta^2 r_2,\,
s_2 = r_0 + \zeta^2 r_1 + \zeta r_2.\,

The roots may then be recovered from the three si by inverting the above linear transformation, giving

r_0 = (s_0 + s_1 + s_2)/3,\,
r_1 = (s_0 + \zeta^2 s_1 + \zeta s_2)/3,\,
r_2 = (s_0 + \zeta s_1 + \zeta^2 s_2)/3.\,

We already know the value s0 = −a, so we only need to seek values for the other two. However, if we take the cubes, a cyclic permutation leaves the cubes invariant, and a transposition of two roots exchanges s13 and s23, hence the polynomial

(z-s_1^3)(z-s_2^3) \qquad (5)

is invariant under permutations of the roots, and so has coefficients expressible in terms of (1). Using calculations involving symmetric functions or alternatively field extensions, we can calculate (5) to be

{z}^{2}+ \left( -9\,ba+2\,{a}^{3}+27\,c \right) z+ \left( {a}^{2}-3\,b\right)^{3}.

The roots of this quadratic equation are

\frac92\,ab-{a}^{3}- \frac{27}{2}\,c \pm \frac32\,\sqrt{3\Delta},

where Δ is the discriminant defined above. Taking cube roots give us s1 and s2, from which we can recover the roots ri of (1).

[edit] Factorization

If r is any root of (1), then we may factor using r to obtain

(xr)(x2 + (a + r)x + b + ar + r2) = x3 + ax2 + bx + c.

Hence if we know one root we can find the other two by solving a quadratic equation, giving

\frac12 \left(-a-r \pm \sqrt{-3r^2-2ar+a^2-4b}\right)

for the other two roots.

[edit] Solution in terms of Chebyshev radicals

If we have a cubic equation which is already in depressed form, we may write it as \,x^3 - 3px - q = 0. Substituting x = \sqrt{p} z we obtain z^3 - 3z - p^{-\frac{3}{2}}q = 0 or equivalently

z^3 - 3z = p^{-\frac{3}{2}}q \ .

From this we obtain solutions to our original equation in terms of the Chebyshev cube root C_{1\over3} as

r_0 = \sqrt{p}\,C_{1\over3}(p^{-\frac{3}{2}}q),\,
r_1 = -\sqrt{p}\,C_{1\over3}(-p^{-\frac{3}{2}}q),\,
r_2 = -r_0 - r_1 \ .

If now we start from a general equation

x^3 + ax^2 + bx +c = 0 \qquad (1)

and reduce it to the depressed form under the substitution x = ta/3, we have \, p = (a^2-3b)/9 and \, q = -(2a^3-9ab+27c)/27, leading to

t_{a;b;c} = p^{-\frac{3}{2}}q = -\frac{2a^3-9ab+27c}{(a^2-3b)^{3/2}}.

This gives us the solutions to (1) as

r_0 = \sqrt{p}\,C_{1\over3}(t_{a;b;c})-{a\over 3} ,\,
r_1 = -\sqrt{p}\,C_{1\over3}(-t_{a;b;c})-{a\over 3},\,
r_2 = -r_0 - r_1 - a \ .

[edit] The case of a cubic equation with real coefficients

Suppose the coefficients of (1) are real. If s is the quantity q/r from the section on real roots, then s = t2; hence 0 < s < 4 is equivalent to −2 < t < 2, and in this case we have a polynomial with three distinct real roots, expressed in terms of a real function of a real variable, quite unlike the situation when using cube roots. If s > 4 then either t > 2 and C_{1\over3}(t) is the sole real root, or t < −2 and -C_{1\over3}(-t) is the sole real root. If s < 0 then the reduction to Chebyshev polynomial form has given a t which is a pure imaginary number; in this case iC_{1\over3}(-it)-iC_{1\over3}(it) is the sole real root. We are now evaluating a real root by means of a function of a purely imaginary argument; however we can avoid this by using the function

S_{1\over3}(t) = iC_{1\over3}(-it)-iC_{1\over3}(it) = 2 \operatorname{sinh}\left(\operatorname{arcsinh}\left({t\over2}\right)/3\right),\,

which is a real function of a real variable with no singularities along the real axis. If a polynomial can be reduced to the form x3 + 3xt with real t, this is a convenient way to solve for its roots.

[edit] See also

[edit] External links

[edit] References