A Renaissance Drama: The Race to Solve the Cubic

A Renaissance Drama: The Race to Solve the Cubic


Latest

A Renaissance Drama: The Race to Solve the Cubic

Author: Christoph Würsch, Institute for Computational Engineering, ICE, OST

Cardano’s Method for Cubic Equations

Girolamo Cardano (1501–1576) published the general solution of the cubic in Ars Magna (1545). The depressed-cubic case was known earlier to Scipione del Ferro (c. 1465–1526) and was independently rediscovered by Niccolò Tartaglia (1500–1557). Cardano, after learning of del Ferro’s prior discovery, included the method in his treatise and credited del Ferro. The modern streamlined derivation below follows the spirit of Cardano’s presentation.

The story of the cubic formula isn’t one of quiet contemplation in a university library; it’s a tale of secrets, ambition, betrayal, and a dramatic intellectual duel fought in the public squares of 16th-century Italy.

Girolamo Cardano (1501-1576)

Our story begins in Bologna with Scipione del Ferro, a mathematics professor with a secret. In an era where academic reputations were made and broken in public challenges, a unique discovery was a weapon to be guarded. Del Ferro had found the holy grail of his time: a method for solving the “depressed” cubic equation, those of the form x3+mx=nx^3 + mx = n. For years, he kept this powerful knowledge to himself, a secret key that could unlock problems no one else could solve. Only on his deathbed, around 1526, did he pass it on, whispering the formula to his student, the mediocre but ambitious Antonio Fior. Fior now possessed a treasure he hadn’t earned, a single, powerful trick he believed would make him famous.

The Challenger from the Gutter

Meanwhile, in Venice, a brilliant, self-taught mathematician named Niccolò Fontana was making a name for himself. Scarred on the face by a French soldier’s sword as a child, the wound left him with a permanent stammer, earning him the nickname “Tartaglia” (the Stammerer), which he defiantly adopted. Tartaglia was a mathematical prodigy who rose from poverty to become one of Italy’s most respected problem-solvers.

Nicolò Tartaglia (1500-1557)

In 1535, news of Tartaglia’s prowess reached Fior, who, seeing a chance for glory, made a fatal mistake: he challenged the Stammerer to a public mathematics contest. The rules were simple: each man would set 30 problems for the other. The winner would take home the prize money and, more importantly, the glory. Fior, confident in his secret weapon, crafted all his problems based on the x3+mx=nx^3 + mx = n form, certain that Tartaglia would be unable to solve them.

Tartaglia, receiving the list of problems, realized they were all of the same, unsolved type. Panic set in. With the clock ticking, he threw himself at the problem with ferocious intensity. In a stunning display of genius, on the night of February 13, 1535, just days before the deadline, Tartaglia had a breakthrough. He rediscovered the solution for himself! He swiftly solved all 30 of Fior’s problems in a matter of hours. When the day of reckoning came, Tartaglia was triumphant. Fior, on the other hand, couldn’t solve a single one of the varied problems Tartaglia had set for him. His one-trick pony had failed, and he faded into historical obscurity.

The Oath and the Betrayal

Tartaglia’s stunning victory made him a celebrity. News of his secret method reached Milan and the ears of the brilliant, eccentric, and insatiably curious polymath Gerolamo Cardano. A renowned physician, astrologer, and compulsive gambler, Cardano was also a mathematician of the first rank and was writing a comprehensive book on algebra, Ars Magna (“The Great Art”). The cubic solution was the one missing piece.

Cardano begged Tartaglia to share the secret. Tartaglia, wary and protective of his hard-won knowledge, repeatedly refused. But Cardano was persistent. He flattered Tartaglia, invited him to Milan, and promised him patronage and introductions. Finally, he cornered Tartaglia and swore a sacred oath:

“I swear to you, by God’s holy Gospels, and as a true man of honour, not only never to publish your discoveries, if you teach them to me, but I promise you, and I pledge my faith as a true Christian, to note them down in code, so that after my death no one will be able to understand them.”

Moved by the solemnity of the oath, Tartaglia relented. He revealed his solution, possibly in the form of a cryptic poem, to keep it from being easily understood.

The Great Art and the Aftermath

Cardano was no Fior. He not only understood the method but, with the help of his brilliant young student Lodovico Ferrari, mastered and expanded it. Together, they figured out how to solve all cubic equations, including those with an x2x^2 term. During this work, Cardano stumbled upon a startling implication: the formula sometimes required taking the square root of a negative number to arrive at a real answer. These were the first whispers of complex numbers.

Years passed. Cardano, honoring his oath, kept the secret. But then, on a trip to Bologna, he and Ferrari made a shocking discovery while examining del Ferro’s old notebooks. There, plain as day, was the solution to the depressed cubic. Del Ferro had discovered it before Tartaglia.

Cardano saw his chance. His oath, he reasoned, bound him not to publish Tartaglia’s work. But this was del Ferro’s work, discovered earlier! Feeling absolved, in 1545, he published Ars Magna. In it, he presented the complete solution to the cubic and quartic (solved by Ferrari). He gave full credit to both del Ferro for the first discovery and Tartaglia for his independent rediscovery.

But for Tartaglia, it was the ultimate betrayal. His name was in the book, but the glory belonged to Cardano. Enraged, he launched a vicious public feud, accusing Cardano of theft and oath-breaking. Cardano himself refused to engage, letting his fiercely loyal student, Ferrari, fight in his place. Another public contest was arranged, this time between Tartaglia and Ferrari. The younger, more agile mind of Ferrari won the day. Tartaglia, defeated and humiliated, lost his prestige and died a few years later, impoverished and bitter.

Cardano’s name is the one forever attached to the method, a testament to a dramatic and bitter chapter in the history of science, where genius, secrecy, and ambition collided to give the world one of its great mathematical treasures.

We present a clear and self-contained derivation of the solution formulas for cubic equations in two stages: first the depressed cubic x3+px+q=0x^3+px+q=0, and then the general monic cubic x3+ax2+bx+c=0x^3+ax^2+bx+c=0 via a Tschirnhaus shift. Along the way we justify each algebraic step, derive the resolvent quadratic, discuss how to choose cube-root branches consistently, and analyze the structure of real and complex roots via the discriminant. Brief historical notes situate the method in the work of Cardano, with connections to del Ferro and Tartaglia.

The Depressed Cubic x3+px+q=0x^3+px+q=0

First, let’s tackle a simplified version of the cubic equation called the depressed cubic. This is a cubic equation that is missing the x2x^2 term. It has the general form:

x3+px+q  =  0x^3 + p\,x + q \;=\; 0

Here, pp and qq are known numbers, which can be real or complex. The entire strategy hinges on a brilliantly clever substitution.

Step 1: The substitution x=u+vx=u+v

The core idea is to replace the single unknown variable xx with the sum of two new unknown variables, uu and vv. So, we set x=u+vx = u+v. This might seem strange—why make the problem more complicated by introducing two unknowns instead of one? The reason is that this gives us extra flexibility. We can later impose a convenient condition on uu and vv that will make the equation much simpler.

Let’s substitute x=u+vx=u+v into our depressed cubic equation:

(u+v)3+p(u+v)+q=0(u+v)^3 + p(u+v) + q = 0

Expanding the (u+v)3(u+v)^3 term gives us:

(u3+3u2v+3uv2+v3)+p(u+v)+q=0(u^3+3u^2v + 3uv^2+v^3) + p(u+v) + q = 0

Now, let’s group the terms a bit differently. We can factor 3uv3uv out of the middle two terms of the expansion:

u3+v3+3uv(u+v)+p(u+v)+q=0u^3+v^3 + 3uv(u+v) + p(u+v) + q = 0

This is where the magic happens. We have two terms that contain (u+v)(u+v). Let’s factor that out:

u3+v3+(3uv+p)(u+v)+q=0u^3+v^3 + (3uv + p)(u+v) + q = 0

Now we use the extra flexibility we gained from introducing two variables. We can impose a condition on uu and vv that simplifies this equation. Look at the term (3uv+p)(3uv + p). If we could make that whole term equal to zero, the equation would become much simpler! So, let’s do exactly that. We will require that uu and vv satisfy the following constraint:

3uv+p=0uv=p33uv + p = 0 \quad\Longleftrightarrow\quad uv = -\frac{p}{3}

By enforcing this condition, the middle part of our cubic equation simply vanishes! The equation now simplifies dramatically to:

u3+v3+q=0u3+v3=qu^3 + v^3 + q = 0 \quad\Longleftrightarrow\quad u^3+v^3 = -q

So, by substituting x=u+vx=u+v and imposing the condition uv=p/3uv = -p/3, we’ve transformed our original problem into finding two numbers, uu and vv, that satisfy two new equations:

  1. u3+v3=qu^3+v^3 = -q
  2. uv=p/3uv = -p/3

Step 2: The resolvent quadratic for u3u^3 and v3v^3

Let’s look at the two conditions we have. If we cube the second condition, we get (uv)3=(p/3)3(uv)^3 = (-p/3)^3, which is u3v3=p3/27u^3v^3 = -p^3/27.

Now, let’s make another small substitution to make things even clearer. Let s=u3s = u^3 and t=v3t = v^3. Our two conditions become:

  1. Sum: s+t=qs + t = -q
  2. Product: st=p3/27st = -p^3/27

This is a classic problem! If you know the sum and product of two numbers (ss and tt), you can find them by forming a quadratic equation. A quadratic equation with roots ss and tt can be written as (zs)(zt)=0(z-s)(z-t) = 0. Expanding this gives z2(s+t)z+st=0z^2 - (s+t)z + st = 0.

We already know (s+t)(s+t) and stst, so we can plug them right in:

z2(q)z+(p327)  =  0z^2 - (-q)\,z + \left(-\frac{p^3}{27}\right) \;=\; 0

This simplifies to what is known as the resolvent quadratic equation:

z2+qzp327  =  0z^2 + q\,z - \frac{p^3}{27} \;=\; 0

We can solve this for zz using the quadratic formula, z=B±B24AC2Az = \frac{-B \pm \sqrt{B^2-4AC}}{2A}. In our case, A=1A=1, B=qB=q, and C=p3/27C=-p^3/27. The two solutions for zz will be our values for ss and tt.

s,t=q±q24(1)(p3/27)2=q2±q24+p327s, t = \frac{-q \pm \sqrt{q^2 - 4(1)(-p^3/27)}}{2} = -\frac{q}{2} \pm \sqrt{\frac{q^2}{4} + \frac{p^3}{27}}

So, the two solutions are: s=q2+(q2)2+(p3)3,t=q2(q2)2+(p3)3s = -\frac{q}{2} + \sqrt{\left(\frac{q}{2}\right)^2 + \left(\frac{p}{3}\right)^3}, \qquad t = -\frac{q}{2} - \sqrt{\left(\frac{q}{2}\right)^2 + \left(\frac{p}{3}\right)^3}

Let’s define the expression inside the square root as DD, which is closely related to the discriminant: D  :=  (q2)2+(p3)3D \;:=\; \left(\frac{q}{2}\right)^2 + \left(\frac{p}{3}\right)^3 With this definition, our solutions for ss and tt become much cleaner: s=q/2+Ds = -q/2 + \sqrt{D} and t=q/2Dt = -q/2 - \sqrt{D}.

Step 3: Taking matched cube roots and solving for xx

We’ve found ss and tt. Remember that we defined s=u3s=u^3 and t=v3t=v^3. To find uu and vv, we just need to take the cube roots:

u=s3=q2+D3u = \sqrt[3]{s} = \sqrt[3]{-\frac{q}{2}+\sqrt{D}} v=t3=q2D3v = \sqrt[3]{t} = \sqrt[3]{-\frac{q}{2}-\sqrt{D}}

Since our original substitution was x=u+vx=u+v, the solution for xx appears to be:

x  =  q2+D3  +  q2D3x \;=\; \sqrt[3]{-\frac{q}{2}+\sqrt{D}} \;+\; \sqrt[3]{-\frac{q}{2}-\sqrt{D}}

This is Cardano’s famous formula for the depressed cubic! However, there’s a very important subtlety here. Any non-zero complex number has three distinct cube roots. So when we write s3\sqrt[3]{s}, there are three possible values for uu. Similarly, there are three possible values for vv. This gives 3×3=93 \times 3 = 9 possible combinations for the pair (u,v)(u, v). Which one is correct?

This is where our constraint from Step 1, uv=p/3uv = -p/3, becomes critical again. We must choose the cube root for uu and the cube root for vv such that their product is exactly p/3-p/3. This requirement “matches” the roots and reduces the nine possible pairs to just three valid ones, which correspond to the three roots of the cubic equation.

Branch matching with cubic roots. To understand this matching, we use the complex roots of unity. The three cube roots of 1 are 11, ω=e2πi/3=12+i32\omega = e^{2\pi i/3} = -\frac{1}{2}+i\frac{\sqrt{3}}{2}, and ω2=e4πi/3=12i32\omega^2 = e^{4\pi i/3} = -\frac{1}{2}-i\frac{\sqrt{3}}{2}. Notice that 1+ω+ω2=01+\omega+\omega^2=0.

If u0u_0 is one specific cube root of ss, then the other two are ωu0\omega u_0 and ω2u0\omega^2 u_0. Similarly, if v0v_0 is one cube root of tt, the others are ωv0\omega v_0 and ω2v0\omega^2 v_0.

Let’s say we pick a pair (u0,v0)(u_0, v_0) such that their product u0v0=p/3u_0 v_0 = -p/3. What about the other combinations?

  • (ωu0)×(ωv0)=ω2(u0v0)(\omega u_0) \times (\omega v_0) = \omega^2 (u_0 v_0) --- This product is wrong.
  • (ωu0)×(ω2v0)=ω3(u0v0)=1×(u0v0)(\omega u_0) \times (\omega^2 v_0) = \omega^3 (u_0 v_0) = 1 \times (u_0 v_0) --- This product is right!
  • (ω2u0)×(ωv0)=ω3(u0v0)=1×(u0v0)(\omega^2 u_0) \times (\omega v_0) = \omega^3 (u_0 v_0) = 1 \times (u_0 v_0) --- This is also right.

To keep the product uvuv constant, if we multiply u0u_0 by a root of unity (say, ωk\omega^k), we must multiply v0v_0 by its inverse (ωk\omega^{-k}). This gives us the three valid pairs that satisfy the constraint: (uk,vk):=(ωku0, ωkv0),k=0,1,2(u_k,v_k) := \big(\omega^k u_0,\ \omega^{-k} v_0\big), \qquad k=0,1,2 The three sums xk=uk+vkx_k=u_k+v_k give the three roots of the cubic equation.

Step 4: The discriminant and the nature of roots

The value D=(q/2)2+(p/3)3D = (q/2)^2 + (p/3)^3 tells us about the nature of the roots. It’s directly related to the polynomial discriminant, Δ\Delta, by the formula:

Δ  =  4p327q2  =  108D\Delta \;=\; -4p^3 - 27q^2 \;=\; -108\,D

The sign of Δ\Delta (or equivalently, the sign of DD) separates the solutions into three cases for equations with real coefficients pp and qq:

  • D>0D > 0 (Δ<0\Delta < 0): In this case, D\sqrt{D} is a real number. The two cube roots in the formula are of real numbers. This gives one real root and two complex conjugate roots. Cardano’s formula works perfectly here and directly gives you the real root.

  • D=0D = 0 (Δ=0\Delta = 0): Here, D=0\sqrt{D}=0, so u=v=q/23u=v=\sqrt[3]{-q/2}. The equation has multiple real roots. Specifically, if pp and qq are not both zero, there is a double root and a simple root. If p=q=0p=q=0, there is a triple root at x=0x=0.

  • D<0D < 0 (Δ>0\Delta > 0): This is the most interesting and historically challenging case, known as the casus irreducibilis (the irreducible case). Here, D\sqrt{D} is an imaginary number. This means that to find the roots, you must take the cube root of complex numbers. The paradox is that in this case, the equation has three distinct real roots. The formula forces you to go through the complex numbers to find real solutions, which baffled mathematicians for centuries.

Step 5: Trigonometric (all-real) representation when D<0D<0

For the casus irreducibilis, where we have three real roots but the formula involves complex numbers, there is a more practical way to find the roots using trigonometry. This method avoids complex radicals. If p<0p<0 and D<0D<0, we can make the substitution:

x  =  2p3  cosθx \;=\; 2\sqrt{-\frac{p}{3}}\;\cos\theta

Plugging this into the cubic x3+px+q=0x^3+px+q=0 and using the trigonometric identity cos(3θ)=4cos3θ3cosθ\cos(3\theta)=4\cos^3\theta-3\cos\theta leads to:

2(p3)3/2cos(3θ)+q=02\left(-\frac{p}{3}\right)^{3/2} \cos(3\theta) + q = 0

We can then solve for cos(3θ)\cos(3\theta):

cos(3θ)  =  q227p3\cos(3\theta) \;=\; -\frac{q}{2}\,\sqrt{-\frac{27}{p^3}}

The value on the right-hand side will be between -1 and 1 precisely when D0D \le 0. We can find the angle 3θ3\theta by taking the arccosine. Since cosine is a periodic function, there are multiple angles that work, which is how we find all three roots. The three solutions for xx are given by:

xk  =  2p3  cos ⁣(13arccos ⁣(q227p3)+2πk3),k=0,1,2x_k \;=\; 2\sqrt{-\frac{p}{3}}\; \cos\!\left(\frac{1}{3}\arccos\!\Bigl( -\frac{q}{2}\sqrt{-\frac{27}{p^3}}\Bigr) + \frac{2\pi k}{3}\right), \qquad k=0,1,2

This formula gives all three real roots without any intermediate complex numbers.

Special cases. If p=0p=0, then x3+q=0x=q3x^3+q=0 \Rightarrow x=-\sqrt[3]{q}. If q=0q=0, then x(x2+p)=0x=0x(x^2+p)=0 \Rightarrow x=0 or x=±px=\pm\sqrt{-p}.

The General Cubic x3+ax2+bx+c=0x^3+ax^2+bx+c=0

Now, what if our cubic equation has an x2x^2 term? x3+ax2+bx+c  =  0x^3 + a\,x^2 + b\,x + c \;=\; 0 The strategy is to transform this general cubic into a depressed cubic, which we already know how to solve!

Step 1: Depress the cubic (kill the x2x^2-term)

We can eliminate the x2x^2 term by shifting the function horizontally. This is achieved with a substitution known as a Tschirnhaus transformation. We set:

x=ya3x = y - \frac{a}{3}

Why this specific value? When you substitute this into the general cubic and expand the terms, the terms involving y2y^2 will perfectly cancel out. Let’s see it: (ya/3)3+a(ya/3)2+...(y-a/3)^3 + a(y-a/3)^2 + ... The first term gives 3y2(a/3)=ay2-3y^2(a/3) = -ay^2. The second term gives +ay2+ay^2. These two terms cancel, eliminating the squared term!

Substituting x=ya/3x = y-a/3 fully and collecting terms gives us a new depressed cubic in the variable yy:

y3+py+q  =  0y^3 + p\,y + q \;=\; 0

where the new coefficients pp and qq are defined in terms of the old ones (a,b,ca,b,c):

p=ba23,q=2a327ab3+cp = b - \frac{a^2}{3}, \quad q = \frac{2a^3}{27} - \frac{ab}{3} + c

Step 2: Apply Cardano to the depressed form

We now have a depressed cubic y3+py+q=0y^3 + py + q = 0. We can solve this for yy using the exact method we developed in the first section. We calculate our value DD using the new pp and qq:

D=(q2)2+(p3)3D = \left(\frac{q}{2}\right)^2 + \left(\frac{p}{3}\right)^3 Let s=q/2+Ds = -q/2 + \sqrt{D} and t=q/2Dt = -q/2 - \sqrt{D}. Then, by taking properly matched cube roots u=s3u=\sqrt[3]{s} and v=t3v=\sqrt[3]{t} such that uv=p/3uv=-p/3, we find the solutions for yy:

y=u+v=q2+D3+q2D3y = u + v = \sqrt[3]{-\frac{q}{2}+\sqrt{D}} + \sqrt[3]{-\frac{q}{2}-\sqrt{D}}

Finally, we must remember that we solved for yy, not our original variable xx. To get the solution for xx, we just undo the substitution from Step 1: x=ya/3x = y - a/3. This gives the complete formula for the general cubic equation:

x  =  a3  +  q2+D3  +  q2D3x \;=\; -\frac{a}{3} \;+\; \sqrt[3]{-\frac{q}{2}+\sqrt{D}} \;+\; \sqrt[3]{-\frac{q}{2}-\sqrt{D}} where pp and qq are calculated from a,b,ca,b,c as shown in Step 1.

Step 3: All three roots via cubic roots of unity

Just as before, we can find all three roots. Let u0u_0 and v0v_0 be principal cube roots of ss and tt that satisfy u0v0=p/3u_0v_0=-p/3. The three solutions for yy are:

yk=ωku0+ωkv0,k=0,1,2y_k = \omega^{\,k}u_0 + \omega^{-k} v_0,\qquad k=0,1,2 The three solutions for xx are then found by undoing the shift for each yky_k:

xk  =  yka3,k=0,1,2x_k \;=\; y_k - \frac{a}{3},\qquad k=0,1,2

Discriminant and Complex-Root Analysis for the General Cubic

The horizontal shift x=ya/3x = y - a/3 does not change the nature of the roots (e.g., how many are real vs. complex). Therefore, the discriminant of the depressed cubic y3+py+q=0y^3+py+q=0 tells us everything we need to know about the roots of the original cubic x3+ax2+bx+c=0x^3+ax^2+bx+c=0.

We use the same discriminant, calculated using the pp and qq from the depressed form: Δ  =  4p327q2  =  108D\Delta \;=\; -4p^3 - 27q^2 \;=\; -108\,D

  • Δ>0\Delta > 0 (D<0D < 0): three distinct real roots.
  • Δ=0\Delta = 0 (D=0D = 0): a multiple root (either a triple root, or one simple + one double root). All roots are real.
  • Δ<0\Delta < 0 (D>0D > 0): one real root and two nonreal complex-conjugate roots.

Summary

Cardano’s method is a beautiful and powerful algorithm for solving any cubic equation. The overall strategy is:

  1. Depress the Cubic: If the equation has an x2x^2 term, use the substitution x=ya/3x = y - a/3 to eliminate it. This reduces the problem to solving a simpler “depressed” cubic.
  2. The u+vu+v Substitution: For the depressed cubic y3+py+q=0y^3+py+q=0, substitute y=u+vy = u+v. This allows you to impose the condition uv=p/3uv = -p/3, which simplifies the equation immensely.
  3. Solve the Resolvent Quadratic: The problem is now reduced to finding u3u^3 and v3v^3. Their sum is q-q and their product is p3/27-p^3/27. These two values are the roots of a simple quadratic equation.
  4. Take Matched Cube Roots: Find uu and vv by taking the cube roots. Crucially, you must choose the roots that satisfy the condition uv=p/3uv = -p/3. Using the complex roots of unity allows you to find all three valid pairs, giving the three solutions for yy.
  5. Undo the Shift: Once you have the solutions for yy, use the relation x=ya/3x = y - a/3 to find the solutions to your original equation.

The discriminant Δ=4p327q2\Delta=-4p^3-27q^2 reveals the nature of the roots (real, complex, or multiple) without having to solve the entire equation.

Further reading: G. Cardano, Ars Magna, 1545 (English trans. T. Richard Witmer, 1968); R. Bombelli, L’Algebra, 1572, for early systematic use of complex numbers in the cubic.