Most of us remember the quadratic formula from secondary school algebra. If you hand a student any equation of the form ax^2 + bx + c = 0, they can immediately rattle off the solution using "negative b plus or minus the square root of b squared minus 4ac, all over 2a." It is compact, deterministic, and dependable. But what happens when you step up to degree three? For the general cubic equation, ax^3 + bx^2 + c*x + d = 0, is there an equivalent closed-form master formula?
The answer is yes, though its history is full of mathematical intrigue, duels, and secrecy in 16th-century Italy involving Scipione del Ferro, Niccolò Tartaglia, and Gerolamo Cardano. Unlike the quadratic case, which relies on a simple completion of the square, solving the cubic requires a two-stage strategy: first, eliminating the quadratic term through a clever shift, and second, splitting the single unknown into two interacting variables that reveal a hidden quadratic equation underneath. Here is how the complete derivation unfolds.
Step 1: Depressing the Polynomial
We begin with the general polynomial equation ax^3 + bx^2 + c*x + d = 0, where we assume a is non-zero. Dividing the entire equation through by the leading coefficient a yields the monic cubic x^3 + (b/a)*x^2 + (c/a)*x + (d/a) = 0. The chief obstacle preventing a clean solution is the presence of both x^3 and x^2 terms. In quadratic algebra, completing the square works because the linear term can be absorbed into a binomial square. In a cubic, we can perform an analogous transformation—known historically as a Tschirnhaus transformation—to depress the polynomial by eliminating the x^2 term entirely.
To find the correct coordinate shift, let x = t - b/(3a). When we substitute this expression into x^3 + (b/a)x^2 + (c/a)x + (d/a) = 0 and expand every binomial term, the expansion of (t - b/(3a))^3 produces t^3 - (b/a)t^2 + [b^2/(3a^2)]t - b^3/(27a^3). Meanwhile, the expansion of (b/a)(t - b/(3a))^2 produces (b/a)t^2 - [2b^2/(3a^2)]t + b^3/(9a^3). Notice what happens when these two expressions are added: the - (b/a)*t^2 from the cubic term and the +(b/a)*t^2 from the quadratic term cancel each other out completely.
Gathering the remaining linear coefficients and the constant terms produces the classic "depressed cubic" equation: t^3 + pt + q = 0. In this simplified form, the parameter p is defined strictly by the original coefficients as p = (3ac - b^2) / (3a^2), while the constant term q is given by q = (2b^3 - 9abc + 27a^2d) / (27*a^3). Any method capable of solving this reduced form unlocks the general cubic, because once t is known, recovering the original root requires only computing x = t - b/(3a).
Step 2: Cardano's Auxiliary Substitution
Now we must solve t^3 + p*t + q = 0. Cardano’s ingenious insight was to replace the single unknown t with the sum of two auxiliary variables: t = u + v. At first glance, replacing one unknown with two seems counterproductive, but it introduces a degree of freedom that we can exploit to our advantage.
Substituting u + v in place of t yields (u + v)^3 + p*(u + v) + q = 0. Expanding the binomial cube using the identity (u + v)^3 = u^3 + v^3 + 3uv*(u + v), our equation reorganizes into u^3 + v^3 + (3uv + p)*(u + v) + q = 0.
Because we introduced two variables to represent one, we are mathematically entitled to impose any second condition on u and v that suits our convenience. The most elegant choice is to eliminate the middle term by setting 3uv + p = 0, which immediately implies that u*v = -p/3.
With that term vanquished, the cubic equation collapses into a remarkably simple relation: u^3 + v^3 = -q. Furthermore, if we cube both sides of our imposed constraint uv = -p/3, we find that (u^3)(v^3) = -p^3 / 27.
Step 3: The Resolvent Quadratic
Look closely at the two equations we now possess: we know the sum of two numbers, u^3 + v^3 = -q, and we know their product, (u^3)(v^3) = -p^3 / 27. According to Vieta’s formulas, two quantities with a known sum S and known product P are the roots of the quadratic equation z^2 - Sz + P = 0. Substituting our specific sum and product produces the "resolvent quadratic": z^2 + q*z - p^3 / 27 = 0.
Because this is a quadratic equation, we can invoke the familiar quadratic formula directly. The solutions for z are z = -q/2 ± sqrt((q/2)^2 + (p/3)^3).
Let us define the cubic discriminant expression Delta = (q/2)^2 + (p/3)^3. Because u and v are symmetric in our derivation, we can designate one root to be u^3 and the other to be v^3 without any loss of generality: u^3 = -q/2 + sqrt(Delta) and v^3 = -q/2 - sqrt(Delta).
Step 4: Extracting Roots and Handling Multiplicity
To determine u and v, we take the cube roots of both expressions. However, in the complex plane, every non-zero complex number possesses three distinct cube roots. If omega represents the primitive complex cube root of unity, omega = -1/2 + isqrt(3)/2, then the three algebraic values for u are u_0, u_0omega, and u_0*omega^2, where u_0 = cbrt(-q/2 + sqrt(Delta)).
We must be careful when pairing these roots with the three cube roots of v. Our earlier constraint required that uv = -p/3. If p is a real number, the product uv must also evaluate to a real number, meaning we cannot pair any arbitrary cube root of u with any arbitrary cube root of v. If we set v_0 = cbrt(-q/2 - sqrt(Delta)) such that u_0 * v_0 = -p/3, the three consistent pairings that satisfy the condition are:
t_1 = u_0 + v_0,
t_2 = u_0omega + v_0omega^2, and
t_3 = u_0omega^2 + v_0omega.
Step 5: The Complete Master Solution
Finally, we reverse our initial coordinate shift to express the roots of our original equation ax^3 + bx^2 + c*x + d = 0. Recalling that x = t - b/(3a), the three general solutions are:
First root: x_1 = cbrt(-q/2 + sqrt(Delta)) + cbrt(-q/2 - sqrt(Delta)) - b/(3a)
Second root: x_2 = omega * cbrt(-q/2 + sqrt(Delta)) + omega^2 * cbrt(-q/2 - sqrt(Delta)) - b/(3a)
Third root: x_3 = omega^2 * cbrt(-q/2 + sqrt(Delta)) + omega * cbrt(-q/2 - sqrt(Delta)) - b/(3a)
In these formulas, the auxiliary terms are fully determined by the original coefficients:
p = (3ac - b^2) / (3a^2)
q = (2b^3 - 9abc + 27a^2d) / (27a^3)
Delta = (q/2)^2 + (p/3)^3
omega = -1/2 + i*sqrt(3)/2
Interpreting the Discriminant
The behavior of the solutions depends on the sign of Delta = (q/2)^2 + (p/3)^3. When Delta is greater than zero, the square root yields a real value, producing one real root and two complex conjugate roots. When Delta equals zero, the polynomial has three real roots, at least two of which are identical.
The most fascinating scenario occurs when Delta is less than zero. Here, the formula requires taking the square root of a negative number, introducing imaginary components into both u^3 and v^3, even though the polynomial actually has three distinct, completely real roots. Known historically as casus irreducibilis ("the irreducible case"), this mathematical paradox demonstrated to Renaissance mathematicians that imaginary numbers were not mere fictions—they were an unavoidable algebraic detour required to discover bona fide real solutions
Comments
Post a Comment