Thursday, July 26, 2007

Late

I'm kind of late on my schedule now :(

Friday, July 20, 2007

Online lectures

http://www.procul.org/blog/2006/06/03/webcast-podcast-t%e1%bb%ab-berkeley-va-mit/

Tuesday, July 17, 2007

Solve cubic equation using Cardano Method

From Wikipedia

Cubic equation

From Wikipedia, the free encyclopedia

Jump to: navigation, search
This article discusses cubic equations in one variable. For a discussion of cubic equations in two variables, see elliptic curve.
Graph of a cubic function; the roots are where the curve crosses the x-axis (y = 0).
Graph of a cubic function; the roots are where the curve crosses the x-axis (y = 0).

In mathematics, a cubic equation is a polynomial equation of the third degree. The general form of a cubic equation is

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

where

a\ne 0 \,

(if a = 0, then the equation becomes a quadratic equation).

Usually, the coefficients a, b, c, 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

[hide]

[edit] History

Cubic equations were known to the ancient Indians and ancient Greeks since the 5th century BCE, and even earlier to the ancient Egyptians, who dealt with the problem of doubling the cube, and attempted to solve it using compass and straightedge constructions.[1] Hippocrates, Menaechmus and Archimedes are believed to have come close to solving this problem using intersecting conic sections,[1] though historians such as Reviel Netz dispute whether the Greeks were thinking about cubic equations or just problems that can lead to cubic equations.

In the 11th century, the poet-mathematician Omar Khayyám (10481131) made significant progress in the theory of cubic equations. In an early paper he wrote regarding cubic equations, he discovered that a cubic equation can have more than one solution, that it cannot be solved using earlier compass and straightedge constructions, and found a geometric solution which could be used to get a numerical answer by consulting trigonometric tables. In his later work, the Treatise on Demonstration of Problems of Algebra, he wrote a complete classification of cubic equations with general geometric solutions found by means of intersecting conic sections.[2][3]

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 Δ < href="http://en.wikipedia.org/wiki/Root_%28mathematics%29" title="Root (mathematics)">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 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 \tfrac{-1}{2} \pm i\tfrac{\sqrt{3}}{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 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 a, b, c, and d

The formula for finding the roots of a cubic function is fairly complicated. Therefore, it is common for some students to use the rational root test or a numerical solution instead.

If we have

f(x) = ax^3 + bx^2 + cx + d = a(x - x_1)(x - x_2)(x - x_3),\,

let

q = \frac{9abc - 27a^2d - 2b^3}{54a^3}

and

r = \sqrt{\left (\frac{3ac-b^2}{9a^2}\right )^3 + q^2}.


Now, let

s = \sqrt[3]{q + r}

and

t = \sqrt[3]{q - r}.

The solutions are

x_1 = s+t-\frac{b}{3a},
x_2=-\frac{1}{2}(s+t)-\frac{b}{3a}+\frac{\sqrt{3}}{2}(s-t)i,
x_3=-\frac{1}{2}(s+t)-\frac{b}{3a}-\frac{\sqrt{3}}{2}(s-t)i.

[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 <>t <>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 <>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] Notes

  1. ^ a b Guilbeau (1930).

    "The Egyptians considered the solution impossible, but the Greeks came nearer to a solution."

  2. ^ J. J. O'Connor and E. F. Robertson (1999), Omar Khayyam, MacTutor History of Mathematics archive.

    "Khayyam himself seems to have been the first to conceive a general theory of cubic equations."

  3. ^ Guilbeau (1930).

    "Omar Al Hay of Chorassan, about 1079 AD did most to elevate to a method the solution of the algebraic equations by intersecting conics."

[edit] References

  • W. S. Anglin; & J. Lambek (1995). "Mathematics in the Renaissance", in The heritage of Thales, Ch. 24. Springers.
  • Lucye Guilbeau (1930). "The History of the Solution of the Cubic Equation", Mathematics News Letter 5 (4), p. 8-12.
  • R.W.D. Nickalls (1993). A new approach to solving the cubic: Cardan's solution revealed, The Mathematical Gazette, 77:354–359.

[edit] External links