Heron's formula

Example

border: 1px solid grey; float: right; margin: 3px;"

Figure of the triangle using CSS. If allnonneg is 0, then the figure will be hidden.

--scale: calc(160px / var(--calculator-maxabc,5) * var(--calculator-allnonneg,1)); --thickness: 3px; width: 240px; height: 150px; position: relative; display: inline-block;" width: calc(var(--calculator-c,5) * var(--scale)); border-top: var(--thickness) solid purple; left: 40px; top: 150px; position: absolute;" width: calc(var(--calculator-b,4) * var(--scale)); border-top: var(--thickness) solid green; transform: /* law of cosines */ rotate(calc(-1 * acos( (var(--calculator-c,5) * var(--calculator-c,5)

• var(--calculator-b,4) * var(--calculator-b,4)

• var(--calculator-a,3) * var(--calculator-a,3) ) / 2 / var(--calculator-b,4) / var(--calculator-c,5) ))); transform-origin: 0px calc(var(--thickness) / 2); /* rotate from the center of the line rather than the top */ left: 40px; top: 150px; position: absolute;" width: calc(var(--calculator-a,3) * var(--scale)); border-top: var(--thickness) solid orange; transform: rotate(calc(acos( (var(--calculator-c,5) * var(--calculator-c,5)

• var(--calculator-a,3) * var(--calculator-a,3)

• var(--calculator-b,4) * var(--calculator-b,4) ) / 2 / var(--calculator-a,3) / var(--calculator-c,5) ))); transform-origin: calc(var(--calculator-a,3) * var(--scale)) calc(var(--thickness) / 2); left: calc(40px + var(--calculator-c,5) * var(--scale) - var(--calculator-a,3) * var(--scale)); top: 150px; position: absolute;"

Let be the triangle with sides , , and . This triangle's semiperimeter is s = \tfrac12(a+b+c)= {}\tfrac12(4+13+15) = 16 therefore , , , and the area is \begin{align} A &= {\textstyle \sqrt{s(s-a)(s-b)(s-c)}} \[3mu] &= {\textstyle \sqrt{16 \cdot 12 \cdot 3 \cdot 1 \vphantom)} } \[3mu] &= 24. \end{align}

In this example, the triangle's side lengths and area are integers, making it a Heronian triangle. However, Heron's formula works equally well when the side lengths are real numbers. As long as they obey the strict triangle inequality, they define a triangle in the Euclidean plane whose area is a positive real number.

Alternate expressions

Heron's formula can also be written in terms of just the side lengths instead of using the semiperimeter, in several ways,

\begin{align} A &=\tfrac{1}{4}\sqrt{(a+b+c)(-a+b+c)(a-b+c)(a+b-c)} \[6mu] &=\tfrac{1}{4}\sqrt{2\bigl(a^2 b^2+a^2c^2+b^2c^2\bigr)-\bigl(a^4+b^4+c^4\bigr)} \[6mu] &=\tfrac{1}{4}\sqrt{\bigl(a^2+b^2+c^2\bigr)\vphantom)^2-2\bigl(a^4+b^4+c^4\bigr)} \[6mu] &=\tfrac{1}{4}\sqrt{4\bigl(a^2b^2+a^2c^2+b^2c^2\bigr)-\bigl(a^2+b^2+c^2\bigr)\vphantom)^2} \[6mu] &=\tfrac{1}{4}\sqrt{4a^2b^2-\bigl(a^2+b^2-c^2\bigr)\vphantom)^2}. \end{align}

After expansion, the expression under the square root is a quadratic polynomial of the squared side lengths , , .

The same relation can be expressed using the Cayley–Menger determinant,{{cite journal

-16A^2 = \begin{vmatrix} 0 & a^2 & b^2 & 1 \ a^2 & 0 & c^2 & 1 \ b^2 & c^2 & 0 & 1 \ 1 & 1 & 1 & 0 \end{vmatrix}.

History

The formula is credited to Heron (or Hero) of Alexandria ( 60 AD),{{cite journal

A formula equivalent to Heron's was discovered by Chinese mathematician Qin Jiushao:

A = \frac1{2}\sqrt{a^2 c^2 - \left(\frac{a^2 + c^2 - b^2}{2}\right)^2},

published in Mathematical Treatise in Nine Sections in 1247.

Proofs

There are many ways to prove Heron's formula, for example using trigonometry as below, or the incenter and one excircle of the triangle, or as a special case of De Gua's theorem (for the particular case of acute triangles), or as a special case of Brahmagupta's formula (for the case of a degenerate cyclic quadrilateral).

Trigonometric proof using the law of cosines

A modern proof, which uses algebra and is quite different from the one provided by Heron, follows.{{cite book | url-access=registration Let be the sides of the triangle and the angles opposite those sides. Applying the law of cosines we get

\cos \gamma = \frac{a^2+b^2-c^2}{2ab}

From this proof, we get the algebraic statement that

\sin \gamma = \sqrt{1-\cos^2 \gamma} = \frac{\sqrt{4a^2 b^2 -\bigl(a^2 +b^2 -c^2\bigr)\vphantom)^2 }}{2ab}.

The altitude of the triangle on base has length , and it follows

\begin{align} A &= \tfrac12 (\mbox{base}) (\mbox{altitude}) \[6mu] &= \tfrac12 ab\sin \gamma \[6mu] &= \frac{ab}{4ab}\sqrt{4a^2 b^2 -\bigl(a^2 +b^2 -c^2\bigr)\vphantom)^2} \[6mu] &= \tfrac14\sqrt{-a^4 - b^4 - c^4 + 2 a^2 b^2 + 2 a^2 c^2 + 2 b^2 c^2} \[6mu] &= \tfrac14\sqrt{(a + b + c)(- a + b + c)(a - b + c)(a + b - c)} \[6mu] &= \sqrt{ \left(\frac{a + b + c}{2}\right) \left(\frac{- a + b + c}{2}\right) \left(\frac{a - b + c}{2}\right) \left(\frac{a + b - c}{2}\right)} \[6mu] &= \sqrt{s(s-a)(s-b)(s-c)}. \end{align}

Algebraic proof using the Pythagorean theorem

The following proof is very similar to one given by Raifaizen.{{Cite journal By the Pythagorean theorem we have b^2 = h^2 + d^2 and a^2 = h^2 + (c - d)^2 according to the figure at the right. Subtracting these yields a^2 - b^2 = c^2 - 2cd. This equation allows us to express in terms of the sides of the triangle: d = \frac{-a^2 + b^2 + c^2}{2c}. For the height of the triangle we have that h^2 = b^2 - d^2. By replacing with the formula given above and applying the difference of squares identity we get \begin{align} h^2 &= b^2-\left(\frac{-a^2 + b^2 + c^2}{2c}\right)^2 \ &= \frac{(2bc - a^2 + b^2 + c^2)(2bc + a^2 - b^2 - c^2)}{4c^2} \ &= \frac{\big((b + c)^2 - a^2\big)\big(a^2 - (b - c)^2\big)}{4c^2} \ &= \frac{(b + c - a)(b + c + a)(a + b - c)(a - b + c)}{4c^2} \ &= \frac{2(s - a) \cdot 2s \cdot 2(s - c) \cdot 2(s - b)}{4c^2} \ &= \frac{4s(s - a)(s - b)(s -c )}{c^2}. \end{align}

We now apply this result to the formula that calculates the area of a triangle from its height: \begin{align} A &= \frac{ch}{2} \ &= \sqrt{\frac{c^2}{4} \cdot \frac{4s(s - a)(s - b)(s - c)}{c^2}} \ &= \sqrt{s(s - a)(s - b)(s - c)}. \end{align}

Trigonometric proof using the law of cotangents

If is the radius of the incircle of the triangle, then the triangle can be broken into three triangles of equal altitude and bases and Their combined area is A = \tfrac12ar + \tfrac12br + \tfrac12cr = rs, where s = \tfrac12(a + b + c) is the semiperimeter.

The triangle can alternately be broken into six triangles (in congruent pairs) of altitude and bases and of combined area (see law of cotangents) \begin{align} A &= r(s-a) + r(s-b) + r(s-c) \[2mu] &= r^2\left(\frac{s - a}{r} + \frac{s - b}{r} + \frac{s - c}{r}\right) \[2mu] &= r^2\left(\cot{\frac{\alpha}{2}} + \cot{\frac{\beta}{2}} + \cot{\frac{\gamma}{2}}\right) \[3mu] &= r^2\left(\cot{\frac{\alpha}{2}} \cot{\frac{\beta}{2}} \cot{\frac{\gamma}{2}}\right)\[3mu] &= r^2\left(\frac{s - a}{r} \cdot \frac{s - b}{r} \cdot \frac{s - c}{r}\right) \[3mu] &= \frac{(s-a)(s-b)(s-c)}{r}. \end{align}

The middle step above is \cot{\tfrac{\alpha}{2}} + \cot{\tfrac{\beta}{2}} + \cot{\tfrac{\gamma}{2}} = {}\cot{\tfrac{\beta}{2}}\cot{\tfrac{\gamma}{2}},}} the triple cotangent identity, which applies because the sum of half-angles is \tfrac\alpha2 + \tfrac\beta2 + \tfrac\gamma2 = \tfrac\pi2.

Combining the two, we get A^2 = s(s - a)(s - b)(s - c), from which the result follows.

Numerical stability

Heron's formula as given above is numerically unstable for triangles with a very small angle, causing round-off error when computing with limited precision such as when using floating-point arithmetic. Such triangles have one or two sides whose length is very close to the semiperimeter, leading to catastrophic cancellation. A stable alternative involves arranging the lengths of the sides so that a \ge b \ge c and computing A = \tfrac14 \sqrt{\big(a + (b + c)\big) \big(c - (a - b)\big) \big(c + (a - b)\big) \big(a + (b - c)\big)}. The extra parentheses indicate the order of operations required to achieve numerical stability in the evaluation.

Similar triangle-area formulae

Three other formulae for the area of a general triangle have a similar structure as Heron's formula, expressed in terms of different variables.

First, if and are the medians from sides and respectively, and their semi-sum is , then A = \tfrac43 \sqrt{\sigma (\sigma - m_a)(\sigma - m_b)(\sigma - m_c)}.

Next, if , , and are the altitudes from sides and respectively, and semi-sum of their reciprocals is , then A^{-1} = 4 \sqrt{H\bigl(H-h_a^{-1}\bigr)\bigl(H-h_b^{-1}\bigr)\bigl(H-h_c^{-1}\bigr)}.

Finally, if and are the three angle measures of the triangle, and the semi-sum of their sines is , then \begin{align} A &= D^{2} \sqrt{S(S-\sin \alpha)(S-\sin \beta)(S-\sin \gamma)} \[5mu] &= \tfrac12 D^{2} \sin \alpha,\sin \beta,\sin \gamma, \end{align}

where is the diameter of the circumcircle, D = a/{\sin \alpha} = b/{\sin \beta} = c/{\sin \gamma}. This last formula coincides with the standard Heron formula when the circumcircle has unit diameter.

Generalizations

Brahmagupta's formula gives the area of a cyclic quadrilateral whose sides have lengths as

K=\sqrt{(s-a)(s-b)(s-c)(s-d)}

where s = \tfrac12(a + b + c + d) is the semiperimeter.

Heron's formula is also a special case of the formula for the area of a trapezoid or trapezium based only on its sides. Heron's formula is obtained by setting the smaller parallel side to zero.

Expressing Heron's formula with a Cayley–Menger determinant in terms of the squares of the distances between the three given vertices, A = \frac{1}{4} \sqrt{- \begin{vmatrix} 0 & a^2 & b^2 & 1 \ a^2 & 0 & c^2 & 1 \ b^2 & c^2 & 0 & 1 \ 1 & 1 & 1 & 0 \end{vmatrix} } illustrates its similarity to Tartaglia's formula for the volume of a three-simplex.

Another generalization of Heron's formula to pentagons and hexagons inscribed in a circle was discovered by David P. Robbins.

Degenerate and imaginary triangles

If one of three given lengths is equal to the sum of the other two, the three sides determine a degenerate triangle, a line segment with zero area. In this case, the semiperimeter will equal the longest side, causing Heron's formula to equal zero.

If one of three given lengths is greater than the sum of the other two, then they violate the triangle inequality and do not describe the sides of a Euclidean triangle. In this case, Heron's formula gives an imaginary result. For example if and , then . This can be interpreted using a triangle in the complex coordinate plane , where "area" can be a complex-valued quantity, or as a triangle lying in a pseudo-Euclidean plane with one space-like dimension and one time-like dimension.

Volume of a tetrahedron

If are lengths of edges of the tetrahedron (first three form a triangle; opposite to and so on), then \text{volume} = \frac {\sqrt {,( - a + b + c + d),(a - b + c + d),(a + b - c + d),(a + b + c - d)}} {192,u,v,w}

where \begin{align} a &= \sqrt {\vphantom{y}xYZ}, \qquad b = \sqrt {XyZ}, & c &= \sqrt {\vphantom{y}XYz}, \qquad d = \sqrt {\vphantom{X}xyz}, \[4mu] X &= (-U + v + w),(U + v + w), & x &= (U - v + w),(U + v - w), \ Y &= (-V + w + u),(V + w + u), & y &= (V - w + u),(V + w - u), \ Z &= (-W + u + v),(W + u + v), & z &= (W - u + v),(W + u - v). \end{align}

Spherical and hyperbolic geometry

L'Huilier's formula relates the area of a triangle in spherical geometry to its side lengths. For a spherical triangle with side lengths and , semiperimeter , and area ,{{cite book | editor-last1 = Gamkrelidze | editor-first1 = R. V. | editor-last2 = Vinberg | editor-first2 = E. B. | publication-date = 1993 \tan^2 \frac S 4 = \tan \frac s2 \tan\frac{s-a}2 \tan\frac{s-b}2 \tan\frac{s-c}2

For a triangle in hyperbolic geometry the analogous formula is \tan^2 \frac S 4 = \tanh \frac s2 \tanh\frac{s-a}2 \tanh\frac{s-b}2 \tanh\frac{s-c}2.

References

1. The formula used here is the [[#Numerical stability. numerically stable formula]] (relabeled for {{tmath. a \leq b \leq c), not simply {{tmath. \textstyle ~!!\sqrt{s(s-a)(s-b)(s-c)}! . For example, with {{tmath. 0.205 but the naive implementation produces {{tmath. 0.000 instead. {{pb The area is reported as "Not a triangle" when the side lengths fail the [[triangle inequality]]. When the area is equal to zero, the three side lengths specify a [[degenerate triangle]] with three colinear points.
2. "Heron's Formula".
3. 秦, 九韶. (1773). "數學九章 (四庫全書本)".
4. (15 December 1997). "Personal email communication between mathematicians John Conway and Peter Doyle".
5. Lévy-Leblond, Jean-Marc. (2020-09-14). "A Symmetric 3D Proof of Heron's Formula". The Mathematical Intelligencer.
6. (1983). "Proceedings of the American Statistical Association, Statistical Computing Section".
7. Bényi, Árpád. (July 2003). "A Heron-type formula for the triangle". Mathematical Gazette.
8. Mitchell, Douglas W.. (November 2005). "A Heron-type formula for the reciprocal area of a triangle". Mathematical Gazette.
9. Mitchell, Douglas W.. (2009). "A Heron-type area formula in terms of sines". Mathematical Gazette.
10. (2009). "Disentangling a triangle". American Mathematical Monthly.
11. Robbins, D. P.. (1994). "Areas of Polygons Inscribed in a Circle". Discrete & Computational Geometry.
12. Schwartz, Mark. (2007). "Review of Conics". The American Mathematical Monthly.
13. Kahan, William. (3 April 2012). "What has the Volume of a Tetrahedron to do with Computer Programming Languages?".