Sedrakyan's area formula:
A Unified Distance-Based Area Formula for Convex Polygons

Hayk Sedrakyan, Nairi Sedrakyan, Christophe Pinarbasi

Download the full paper (PDF)

Abstract

This work fills a 2000-year-old gap and is the structural “missing link”: Heron's formula dates back to ∼ ⁣60\sim\!60 AD and Bretschneider's formula to 1842, yet no such unified formula existed for polygons with more than 4 sides. This work shows that these two formulas—which looked completely different for centuries—are particular cases of a single general formula. The remarkable thing is not that it uses only side lengths and certain distances, but rather that one single unified formula computes the area of any convex nn-gon for any number of sides n≥3n \geq 3 (which is the key novelty), previously thought impossible for n≥5n \geq 5. This novel formula is overcoming the “impossibility” of rigidity for n≥5n \geq 5. Applications include distance geometry, knot theory, drone networks, robotics—where coordinates are unavailable, only pairwise distances are known, and area calculations are mission-critical. The geometrically constructed point B2B_2 can be viewed as a representative point associated with the figure’s area distribution—analogous in spirit to the center of mass, which serves as a representative point for the volume of a three-dimensional body.

Introduction

In computer graphics, geographic information systems, and structural engineering, calculating the area of irregular polygons usually requires breaking them down into a mesh of triangles or using Cartesian coordinates. In many real-world scenarios (flexible structures, robotics, distance-only sensors) coordinates are genuinely unavailable, making this formula not just a niche alternative but the only clean option. This newly discovered formula allows for the area calculation of a flexible or moving convex structure using only distance constraints (sides and specific internal diagonal lengths), bypassing the need to calculate coordinate positions or interior angles. The formula is expressed via a beautifully alternating sum of squared side lengths and a geometrically constructed point B2B_2, which elegantly bypasses the structural flexibility of the polygon by introducing a single reference point relative to the diagonal A1A3A_1A_3. The distance A2B2A_2B_2 accounts for the degrees of freedom of the polygon without requiring trigonometry, coordinates, or triangulation.

Theorem (Sedrakyan's area formula).

The area of a convex nn-gon A1A2…AnA_1A_2\ldots A_n is:

144A1A32⋅A2B22−(a12−a22+⋯+(−1)n−1an2)2\frac{1}{4}\sqrt{4A_1A_3^{2}\cdot A_2B_2^{2}-\Big(a_1^{2}-a_2^{2}+\cdots+(-1)^{n-1}a_n^{2}\Big)^{2}}

where AiAi+1=aiA_iA_{i+1}=a_i (indices mod nn) and point B2B_2 is described in the Theorem provided in the main result.

Remark. This formula has been presented at international conferences as Sedrakyan's area formula. The formula was initially realized by Hayk Sedrakyan, experimentally verified with Christophe Pinarbasi, and formally proved by Nairi Sedrakyan and Hayk Sedrakyan. The author thanks Christophe Pinarbasi for stimulating discussions. Having recognized the result but lacking a formal proof at the time, the authors have, since 2022, shared this topic with several people; however, none of them was able to prove the result. Nevertheless, the authors would like to thank them for their efforts.

Particular cases (Heron's and Bretschneider's formulas)

Theorem (n=3n = 3, Heron's formula). The area ATA_T of a triangle A1A2A3A_1A_2A_3 with side lengths aa, bb, cc is (in the section main results and proofs we show that for a triangle point B2B_2 coincides with vertex A1A_1):

AT=144⋅A1A32⋅A2B22−(a12−a22+a32)2=144a2c2−(a2−b2+c2)2A_T=\dfrac{1}{4}\sqrt{4\cdot A_1A_3^{2}\cdot A_2B_2^{2}-(a_1^{2}-a_2^{2}+a_3^{2})^{2}}=\frac{1}{4}\sqrt{4a^{2}c^{2}-(a^{2}-b^{2}+c^{2})^{2}}

In the next section we show that this formula of ATA_T is actually Heron's formula.

Theorem (n=4n = 4, Bretschneider's formula). The area AQA_Q of a quadrilateral A1A2A3A4A_1A_2A_3A_4 with side lengths aa, bb, cc, dd and diagonals ee, ff is (in the section main results and proofs we show that for a quadrilateral point B2B_2 coincides with A4A_4):

AQ=144⋅A1A32⋅A2B22−(a12−a22+a32−a42)2=144e2f2−(a2−b2+c2−d2)2A_Q=\dfrac{1}{4}\sqrt{4\cdot A_1A_3^{2}\cdot A_2B_2^{2}-(a_1^{2}-a_2^{2}+a_3^{2}-a_4^{2})^{2}}=\dfrac{1}{4}\sqrt{4e^{2}f^{2}-(a^{2}-b^{2}+c^{2}-d^{2})^{2}}

In the next section we show that this formula of AQA_Q is actually Bretschneider's formula.

Remark. In applied sciences, we often encounter situations where we can measure distances, but we have absolutely no way to map out an (x,y)(x, y) coordinate grid. For example:

Distance-Only Sensor Networks: Imagine a swarm of drones or localized robotic nodes. They can ping each other to measure their exact mutual distances (the side lengths and select diagonals of an nn-gon), but they do not know their orientation in global space. They have no coordinate origin.

Flexible Structural Mechanics: Imagine an engineering framework or a architectural mesh under stress. The beams (sides) are shifting and flexing. In these scenarios, you cannot use the coordinate-based formula because you have no coordinates. This area formula provides the only direct mathematical pathway to compute or track the changing area purely out of a distance matrix.

Proof of equivalence for the cases n=3n = 3 and n=4n = 4

The case n=3n = 3 (Heron's formula). Let us prove that the formula of ATA_T is actually Heron's formula. By difference of squares formula, we have

AT=144a2c2−(a2−b2+c2)2=14(2ac−(a2−b2+c2))(2ac+(a2−b2+c2))A_T=\frac{1}{4}\sqrt{4a^{2}c^{2}-(a^{2}-b^{2}+c^{2})^{2}}=\frac{1}{4}\sqrt{\Big(2ac-(a^{2}-b^{2}+c^{2})\Big)\Big(2ac+(a^{2}-b^{2}+c^{2})\Big)}
=14(b2−(a2−2ac+c2))((a2+2ac+c2)−b2)=14(b2−(a−c)2)((a+c)2−b2)=\frac{1}{4}\sqrt{\Big(b^{2}-(a^{2}-2ac+c^{2})\Big)\Big((a^{2}+2ac+c^{2})-b^{2}\Big)}=\dfrac{1}{4}\sqrt{\Big(b^{2}-(a-c)^{2}\Big)\Big((a+c)^{2}-b^{2}\Big)}
=14(b−(a−c))(b+(a−c))(a+c−b)(a+c+b).=\frac{1}{4}\sqrt{(b-(a-c))(b+(a-c))(a+c-b)(a+c+b)}.

Thus, it follows that

AT=144a2c2−(a2−b2+c2)2=(b−a+c)(b+a−c)(a+c−b)(a+b+c)4.A_T=\frac{1}{4}\sqrt{4a^{2}c^{2}-(a^{2}-b^{2}+c^{2})^{2}}=\frac{\sqrt{(b-a+c)(b+a-c)(a+c-b)(a+b+c)}}{4}.

The case n=4n = 4 (Bretschneider's formula). Let us prove that the formula for AQA_Q is actually Bretschneider's formula. To find the area AQA_Q of any quadrilateral Bretschneider's formula states:

AQ=(s−a)(s−b)(s−c)(s−d)−abcdcos⁡2 ⁣α+γ2,s=a+b+c+d2,A_Q=\sqrt{(s-a)(s-b)(s-c)(s-d)-abcd\cos^2\!\tfrac{\alpha+\gamma}{2}},\quad s=\tfrac{a+b+c+d}{2},

where α,γ\alpha, \gamma are opposite angles. Let θ\theta be the angle between diagonals ee and ff. Let us prove that

144e2f2−(a2−b2+c2−d2)2=12efsin⁡θ=(s−a)(s−b)(s−c)(s−d)−abcdcos⁡2 ⁣α+γ2.\frac{1}{4}\sqrt{4e^{2}f^{2}-(a^{2}-b^{2}+c^{2}-d^{2})^{2}}=\dfrac{1}{2}ef\sin\theta=\sqrt{(s-a)(s-b)(s-c)(s-d)-abcd\cos^2\!\tfrac{\alpha+\gamma}{2}}.

Step 1: Our formula equals 12efsin⁡θ\frac{1}{2}ef\sin\theta. Place the quadrilateral in the plane with A1A_1 at the origin and let e⃗=A1A3→\vec{e}=\overrightarrow{A_1A_3}, u⃗=A1A2→\vec{u}=\overrightarrow{A_1A_2}, v⃗=A1A4→\vec{v}=\overrightarrow{A_1A_4}. We get

a2−b2+c2−d2=∣u⃗∣2−∣u⃗−e⃗∣2+∣v⃗−e⃗∣2−∣v⃗∣2=2(u⃗−v⃗)⋅e⃗=2efcos⁡θ.a^2-b^2+c^2-d^2=|\vec{u}|^2-|\vec{u}-\vec{e}|^2+|\vec{v}-\vec{e}|^2-|\vec{v}|^2=2(\vec{u}-\vec{v})\cdot\vec{e}=2ef\cos\theta.

Thus

144e2f2−(a2−b2+c2−d2)2=144e2f2−4e2f2cos⁡2θ=14⋅2efsin⁡θ=12efsin⁡θ.\frac{1}{4}\sqrt{4e^2f^2-(a^2-b^2+c^2-d^2)^2}=\frac{1}{4}\sqrt{4e^2f^2-4e^2f^2\cos^2\theta}=\frac{1}{4}\cdot 2ef\sin\theta=\frac{1}{2}ef\sin\theta.

Step 2: Bretschneider's formula also equals 12efsin⁡θ\frac{1}{2}ef\sin\theta. Splitting the quadrilateral along diagonal ee into two triangles gives:

AQ=12absin⁡α+12cdsin⁡γ,A_Q=\tfrac{1}{2}ab\sin\alpha+\tfrac{1}{2}cd\sin\gamma,

and by the law of cosines e2=a2+b2−2abcos⁡α=c2+d2−2cdcos⁡γe^2=a^2+b^2-2ab\cos\alpha=c^2+d^2-2cd\cos\gamma. Using the identity 2sin⁡αsin⁡γ=cos⁡(α−γ)−cos⁡(α+γ)2\sin\alpha\sin\gamma=\cos(\alpha-\gamma)-\cos(\alpha+\gamma) and expanding (s−a)(s−b)(s−c)(s−d)(s-a)(s-b)(s-c)(s-d), one verifies by direct algebra that:

16(s−a)(s−b)(s−c)(s−d)−16abcdcos⁡2 ⁣α+γ2=4e2f2sin⁡2θ,16(s-a)(s-b)(s-c)(s-d)-16abcd\cos^2\!\tfrac{\alpha+\gamma}{2}=4e^2f^2\sin^2\theta,

and hence Bretschneider's formula gives AQ=12efsin⁡θA_Q=\frac{1}{2}ef\sin\theta.

Conclusion. Both formulas equal 12efsin⁡θ\frac{1}{2}ef\sin\theta, so they are equivalent.

The main result

Theorem (Sedrakyan's area formula).

The area of a convex nn-gon A1A2…AnA_1A_2\ldots A_n is:

144A1A32⋅A2B22−(a12−a22+⋯+(−1)n−1an2)2\frac{1}{4}\sqrt{4A_1A_3^{2}\cdot A_2B_2^{2}-\Big(a_1^{2}-a_2^{2}+\cdots+(-1)^{n-1}a_n^{2}\Big)^{2}}

where AiAi+1=aiA_iA_{i+1}=a_i (indices mod nn), and B2=ℓ2∩ℓ2′B_2=\ell_2\cap\ell_2', lines ℓ2\ell_2 and ℓ2′\ell_2' are defined:

  • ℓ2\ell_2: the line perpendicular to A1A3A_1A_3 passing through the point (x0,0)(x_0, 0), where
    x0=a2−a32+a42−⋯+(−1)nan22ax_0=\frac{a^2-a_3^2+a_4^2-\cdots+(-1)^n a_n^2}{2a}
    and a=A1A3a=A_1A_3, we take A1=(0,0)A_1=(0,0), A3=(a,0)A_3=(a,0).
  • ℓ2′\ell_2': the line parallel to A1A3A_1A_3 passing through the point (0,y0)(0, y_0), where
    y0=2((A1A2A3)−(A1A2…An))ay_0=\frac{2\bigl((A_1A_2A_3)-(A_1A_2\ldots A_n)\bigr)}{a}
    lying on the side of A1A3A_1A_3 opposite to A2A_2, and a=A1A3a=A_1A_3.

Important (constructing the point B2B_2). Looking at the form of y0y_0 in the previous equation, it initially appears that the area of (A1A2…An)(A_1A_2\dots A_n) is required to construct line ℓ2′\ell_2' and, consequently, to locate point B2B_2. This seems to introduce a circular argument, implying that we must already know the area of (A1A2…An)(A_1A_2\dots A_n) to find it. However, in Remark 1 of the next section, we demonstrate that the value of this area is not actually needed to construct ℓ2′\ell_2' and determine the location of B2B_2.

The meaning of the formula

At first glance, the introduction of the point B2B_2 may seem somewhat unusual. However, consider the analogous concept of the center of mass of a three-dimensional body. The center of mass is a single point that represents the distribution of mass throughout the body, and it is well known that this point may lie outside the body itself. In two dimensions, instead of volume, we consider area. The concept of volume in three dimensions corresponds to the concept of area in two dimensions. In a similar way, the point B2B_2 can be viewed as a representative point associated with the area distribution of the figure.

The proof of the main result

Let us denote the area of a convex nn-gon A1A2…AnA_1A_2\ldots A_n by (A1A2…An)(A_1A_2\ldots A_n). Let A1A2=a1A_1A_2=a_1, A2A3=a2A_2A_3=a_2, ..., An−1An=an−1A_{n-1}A_n=a_{n-1}, AnA1=anA_nA_1=a_n.

Problem 1. Find the locus of all points MM in the plane such that

MA32−MA12=a32−a42+⋯+(−1)n−1an2.MA_3^{2}-MA_1^{2}=a_3^{2}-a_4^{2}+\cdots+(-1)^{n-1}a_n^{2}.

Solution of Problem 1. Choose a coordinate system as shown in the figure.

Coordinate system for the locus in Problem 1

We have

MA32−MA12=a32−a42+⋯+(−1)n−1an2.MA_3^{2}-MA_1^{2}=a_3^{2}-a_4^{2}+\cdots+(-1)^{n-1}a_n^{2}.

Thus

(x−a)2+y2−x2−y2=a32−a42+⋯+(−1)n−1an2.(x-a)^{2}+y^{2}-x^{2}-y^{2}=a_3^{2}-a_4^{2}+\cdots+(-1)^{n-1}a_n^{2}.
x=a2−a32+⋯+(−1)nan22a.x=\dfrac{a^{2}-a_3^{2}+\cdots+(-1)^{n}a_n^{2}}{2a}.

Let us prove that the locus of all points MM in the plane is line l2l_2 which passes through the point

(a2−a32+⋯+(−1)nan22a, 0)\Big(\dfrac{a^{2}-a_3^{2}+\cdots+(-1)^{n}a_n^{2}}{2a},\,0\Big)

and is perpendicular to line A1A3A_1A_3.

If n=3n = 3, then a=a3a=a_3 and line l2l_2 passes through the point A1(0,0)A_1(0, 0). If n=4n = 4, then line l2l_2 is passing through point A4A_4, as

A4A32−A4A12=a32−a42.A_4A_3^{2}-A_4A_1^{2}=a_3^{2}-a_4^{2}.

Problem 2. Find the locus of all points MM in the plane such that (MA1A2A3)=(A1A2…An)(MA_1A_2A_3)=(A_1A_2\ldots A_n) and that points MM and A2A_2 are not on the same side of line A1A3A_1A_3.

Solution of Problem 2. Choose a coordinate system as shown in the solution to problem 1. Let M(x,y)M(x,y), where y≤0y\leq 0, we have

(−y)⋅a2=(A1A2…An)−(A1A2A3).\frac{(-y)\cdot a}{2}=(A_1A_2\ldots A_n)-(A_1A_2A_3).

Thus

y=2((A1A2…An)−(A1A2A3))(−a)=2((A1A2A3)−(A1A2…An))a.y=\frac{2\Big((A_1A_2\ldots A_n)-(A_1A_2A_3)\Big)}{(-a)}=\frac{2\Big((A_1A_2A_3)-(A_1A_2\ldots A_n)\Big)}{a}.

Let us prove that the locus of all points MM is line l2′l_2^{\prime} passing through point

(0, 2((A1A2A3)−(A1A2…An))a)\Big(0,\,\frac{2\Big((A_1A_2A_3)-(A_1A_2\ldots A_n)\Big)}{a}\Big)

and is perpendicular to line l2l_2.

If n=3n = 3, then line l2′l_2^{\prime} coincides with line A1A3A_1A_3, as

2((A1A2A3)−(A1A2…An))a=0.\frac{2\Big((A_1A_2A_3)-(A_1A_2\ldots A_n)\Big)}{a}=0.

If n=4n = 4, then line l2′l_2^{\prime} passes through point A4A_4, as

2((A1A2A3)−(A1A2A3A4))a=−2(A1A4A3)a.\frac{2\Big((A_1A_2A_3)-(A_1A_2A_3A_4)\Big)}{a}=-\frac{2(A_1A_4A_3)}{a}.

Let us denote by B2B_2 the intersection point of lines l2l_2 and l2′l_2^{\prime}, in a similar way denote by BiB_i the intersection point of lines lil_i and li′l_i^{\prime}, where i=1,…,ni = 1,\ldots,n.

If n≥5n\geq 5, then let us start with the following remark.

Remark 1. Let a convex nn-gon be given. Construct point MnM_n on line C2C3C_2C_3 such that

(MnC1C2)=(C1C2…Cn).(M_nC_1C_2)=(C_1C_2\ldots C_n).

Let us construct point MnM_n be mathematical induction.

The base case. If n=4n = 4, then from point C4C_4 let us draw a line parallel to line C1C3C_1C_3. Denote by M4M_4 the intersection point of that line and ray C2C3C_2C_3 (see figure).

Base case construction of point M4

Indeed, we have

(M4C1C2)=(C1C2C3)+(M4C1C3)=(C1C2C3)+(C4C1C3)=(C1C2C3C4).(M_4C_1C_2)=(C_1C_2C_3)+(M_4C_1C_3)=(C_1C_2C_3)+(C_4C_1C_3)=(C_1C_2C_3C_4).

The inductive hypothesis. Assume that for n=kn = k the point MkM_k is already constructed, then let us show how point Mk+1M_{k+1} is constructed. Consider the convex polygon C1C2C4…Ck+1C_1C_2C_4\ldots C_{k+1}. According to the inductive hypothesis on ray C2C4C_2C_4 we can construct point MkM_k such that

(MkC1C2)=(C1C2C4…Ck+1).(M_kC_1C_2)=(C_1C_2C_4\ldots C_{k+1}).

Let us draw a line from point MkM_k parallel to line C1C4C_1C_4 (see figure). Denote by NkN_k the intersection point of that line and line C3C4C_3C_4.

Inductive step construction of point N_k

We have

(NkC1C4)=(MkC1C4).(N_kC_1C_4)=(M_kC_1C_4).

Thus

(NkC1C2C4)=(MkC1C4)+(C1C2C4)=(MkC1C2)=(C1C2C4…Ck+1).(N_kC_1C_2C_4)=(M_kC_1C_4)+(C_1C_2C_4)=(M_kC_1C_2)=(C_1C_2C_4\ldots C_{k+1}).

Thus, we get that

(NkC1C2C3)=(C1C2C3…Ck+1).(N_kC_1C_2C_3)=(C_1C_2C_3\ldots C_{k+1}).
(1)

It is left for quadrilateral C1C2C3NkC_1C_2C_3N_k on ray C2C3C_2C_3 to construct point Mk+1M_{k+1}, such that

(Mk+1C1C2)=(C1C2C3Nk)(M_{k+1}C_1C_2)=(C_1C_2C_3N_k)
(2)

From (1) and (2), we get

(Mk+1C1C2)=(C1C2C3…Ck+1).(M_{k+1}C_1C_2)=(C_1C_2C_3\ldots C_{k+1}).

This ends the proof for n=k+1n = k+1. Thus, this ends the proof for any n≥4n\geq 4.

To construct line l2′l_2^{\prime} it is sufficient to construct point Mn−1M_{n-1} for the polygon A1A3A4…A4A_1A_3A_4\ldots A_4 and to draw a line passing through point Mn−1M_{n-1} and parallel to line A1A3A_1A_3.

Let us prove the following theorem.

Theorem. The area of a convex nn-gon A1A2…AnA_1A_2\ldots A_n is:

144A1A32⋅A2B22−(a12−a22+⋯+(−1)n−1an2)2.\frac{1}{4}\sqrt{4A_1A_3^{2}\cdot A_2B_2^{2}-\Big(a_1^{2}-a_2^{2}+\cdots+(-1)^{n-1}a_n^{2}\Big)^{2}}.

To prove this theorem, we use the following lemma.

Lemma (the Trapezoid Law [Sedrakyan et al.], see [1], chapter 5). Let sides ABAB and CDCD of quadrilateral ABCDABCD be parallel to each other, then

2⋅AB⋅CD=AC2+BD2−AD2−BC2.2\cdot AB\cdot CD=AC^{2}+BD^{2}-AD^{2}-BC^{2}.

Let us continue the proof of the theorem, consider the the following figure.

Figure for the proof of the main theorem

We have

A2′A32−A2′A12=a32−⋯+(−1)n−1an2.A_2^{\prime}A_3^{2}-A_2^{\prime}A_1^{2}=a_3^{2}-\cdots+(-1)^{n-1}a_n^{2}.

Moreover

(A1A2…An)=(B2A1A2A3)=(B2A1A3)+(A1A2A3)=(A_1A_2\ldots A_n)=(B_2A_1A_2A_3)=(B_2A_1A_3)+(A_1A_2A_3)=
=(B2A1A3)+(A1A2′A3)=(B2A1A2′A3)=12A1A3⋅A2′B2.=(B_2A_1A_3)+(A_1A_2^{\prime}A_3)=(B_2A_1A_2^{\prime}A_3)=\dfrac{1}{2}A_1A_3\cdot A_2^{\prime}B_2.

On the other hand, according to the Pythagorean theorem, we have

12A1A3⋅A2′B2=144⋅A1A32⋅A2′B22=144A1A32⋅(A1A32⋅(A2B22−A2′A22)).\frac{1}{2}A_1A_3\cdot A_2^{\prime}B_2=\frac{1}{4}\sqrt{4\cdot A_1A_3^{2}\cdot A_2^{\prime}B_2^{2}}=\dfrac{1}{4}\sqrt{4A_1A_3^{2}\cdot(A_1A_3^{2}\cdot(A_2B_2^{2}-A_2^{\prime}A_2^{2}))}.
(1)

According to the lemma, we have

4A1A32⋅A2′A22=(A2′A32−A2′A12−A2A32)2.4A_1A_3^{2}\cdot A_2^{\prime}A_2^{2}=(A_2^{\prime}A_3^{2}-A_2^{\prime}A_1^{2}-A_2A_3^{2})^{2}.
(2)

Point A2′A_2^{\prime} is on line l2l_2, thus according to Problem 1, we have

A2′A32−A2′A12=a32−⋯+(−1)n−1an2.A_2^{\prime}A_3^{2}-A_2^{\prime}A_1^{2}=a_3^{2}-\cdots+(-1)^{n-1}a_n^{2}.
(3)

From (1), (2), (3), we get that

12A1A3⋅A2′B2=4A1A32⋅A2B22−(a12−a22+a32−⋯+(−1)n−1an2)2.\frac{1}{2}A_1A_3\cdot A_2^{\prime}B_2=\sqrt{4A_1A_3^{2}\cdot A_2B_2^{2}-\Big(a_1^{2}-a_2^{2}+a_3^{2}-\cdots+(-1)^{n-1}a_n^{2}\Big)^{2}}.

Thus

(A1A2…An)=144A1A32⋅A2B22−(a12−a22+a32−⋯+(−1)n−1an2)2.(A_1A_2\ldots A_n)=\frac{1}{4}\sqrt{4A_1A_3^{2}\cdot A_2B_2^{2}-\Big(a_1^{2}-a_2^{2}+a_3^{2}-\cdots+(-1)^{n-1}a_n^{2}\Big)^{2}}.
(4)

This ends the proof of the theorem.

Conclusion. Triangles (n=3n=3) are inherently rigid because fixing their side lengths completely locks their interior angles and area into place. Quadrilaterals, pentagons, and higher-order polygons (n≥4n\geq 4) are flexible, meaning their sides can deform and flex like a mechanical linkage, constantly changing their angles and overall area. This fundamental fluidness is why a unified, side-length-based area formula was long considered impossible for shapes with five or more sides. Conversely, in geometry, you cannot find the area of an irregular nn-gon purely from its sides because the shape can bend. By capturing that exact architectural fluidness inside a single point B2B_2 and an alternating sum of squared sides, this area formula accomplishes something structurally revolutionary that did not have an existing. It proves that the chaotic flexibility of any arbitrary polygon can be structurally bound and unified under a single algebraic law, opening the door to solving distance-only constraints that were previously considered mathematically impossible.

Example showing how the formula works

Example. Let us find the area of a regular hexagon A1A2…A6A_1A_2\ldots A_6, where A1A2=1A_1A_2=1.

Solution. Line l2l_2 is determined by all points MM satisfying the following equation

MA32−MA12=a32−a42+a52−a62=0.MA_3^{2}-MA_1^{2}=a_3^{2}-a_4^{2}+a_5^{2}-a_6^{2}=0.

So, it will be the perpendicular bisector of line segment A1A3A_1A_3 (see figure).

Now, let us construct line l2′l_2^{\prime} and find the intersection point B2B_2 of lines l2l_2 and l2′l_2^{\prime}. Using the theorem we have that l2′l_2^{\prime} is parrallel to A1A3A_1A_3 and is passing through the point (0,y0)(0, y_0), where

y0=2((A1A2A3)−(A1A2A3A4A5A6))A1A3=−2(A1A3A4A5A6)A1A3.y_0=\dfrac{2\Big((A_1A_2A_3)-(A_1A_2A_3A_4A_5A_6)\Big)}{A_1A_3}=\dfrac{-2(A_1A_3A_4A_5A_6)}{A_1A_3}.

Notice that y0<0y_0 < 0, so −y0>0-y_0 > 0. We get

(A1A3A4A5A6)=(−y0)⋅A1A32=B2H⋅A1A32=(B2A1A3).(A_1A_3A_4A_5A_6)=\dfrac{(-y_0)\cdot A_1A_3}{2}=\dfrac{B_2H\cdot A_1A_3}{2}=(B_2A_1A_3).
Regular hexagon with constructed point B2

Therefore, point B2B_2 is a point on line l2l_2, such that

(B2A1A3)=(A1A3A4A5A6).(B_2A_1A_3)=(A_1A_3A_4A_5A_6).

When A1A5∥A6B2A_1A_5\parallel A_6B_2 and A3A5∥A4B2A_3A_5\parallel A_4B_2, then (A1MA6)=(B2MA5)(A_1MA_6)=(B_2MA_5) and (A3NA4)=(B2NA5)(A_3NA_4)=(B_2NA_5). Hence

(B2A1A3)=(A1A3A4A5A6),(B_2A_1A_3)=(A_1A_3A_4A_5A_6),

as A1A3NA5MA_1A_3NA_5M is in common for triangle B2A1A3B_2A_1A_3 and pentagon A1A3A4A5A6A_1A_3A_4A_5A_6. We obtain

A2B2=A2A5+A5B2=2+A1A6=3,  A1A3=3.A_2B_2=A_2A_5+A_5B_2=2+A_1A_6=3,\ \ A_1A_3=\sqrt{3}.

Applying the area formula (note that a1=a2=a3=a4=a5=a6a_1=a_2=a_3=a_4=a_5=a_6), we get

(A1A2A3A4A5A6)=144⋅3⋅9−(a12−a22+a32−a42+a52−a62)2=144⋅3⋅9−02=332.(A_1A_2A_3A_4A_5A_6)=\frac{1}{4}\sqrt{4\cdot 3\cdot 9-\Big(a_1^2-a_2^2+a_3^2-a_4^2+a_5^2-a_6^2\Big)^{2}}=\frac{1}{4}\sqrt{4\cdot 3\cdot 9-0^{2}}=\frac{3\sqrt{3}}{2}.

Verification of the answer. A regular hexagon with side length 1 can be divided into six equilateral triangles of side length 1. The area of one equilateral triangle with side length 1 is

A△=34.A_{\triangle}=\frac{\sqrt{3}}{4}.

Therefore, the area of the hexagon is

(A1A2A3A4A5A6)=6⋅34=332.(A_1A_2A_3A_4A_5A_6)=6\cdot\frac{\sqrt{3}}{4}=\frac{3\sqrt{3}}{2}.

Applications

Application 1. A triangle has side lengths 3\sqrt{3}, 5\sqrt{5}, and 7\sqrt{7}. Find its area.

Application 2. Consider a trapezoid with bases of lengths 5 and 9. If the sum of the squares of its diagonals is 196, find the sum of the squares of its legs.

Application 3. The area of triangle with sides aa, bb, cc is equal to the area of triangle with sides aa, bb, dd, where d≠cd\neq c, and that a2+b2=16a^{2}+b^{2}=16, find c2+d2c^{2}+d^{2}.

Application 4. Each of two non-congruent triangles has two sides of lengths 5 and 7, and the triangles have equal circumradii. Find the product of the lengths of their third sides.

Application 5. Let a1a_1, a2a_2, a3a_3 be the side lengths of a triangle, prove that

1).   2a12a22+2a22a32+2a12a32>a14+a24+a34.1).\ \ \ 2a_1^{2}a_2^{2}+2a_2^{2}a_3^{2}+2a_1^{2}a_3^{2}>a_1^{4}+a_2^{4}+a_3^{4}.
2).   a12a2+a22a1+a22a3+a32a2+a12a3+a32a1>a13+a23+a33.2).\ \ \ a_1^{2}a_2+a_2^{2}a_1+a_2^{2}a_3+a_3^{2}a_2+a_1^{2}a_3+a_3^{2}a_1>a_1^{3}+a_2^{3}+a_3^{3}.

Application 6. Let ABCDABCD be a convex quadrilateral, such that AB=aAB=a, BC=bBC=b, CD=cCD=c, AD=dAD=d, AC=eAC=e, BD=fBD=f. Prove that

∣a2−b2+c2−d2∣<2ef.\vert a^{2}-b^{2}+c^{2}-d^{2}\vert<2ef.

Application 7. Let A1A2…A2kA_1A_2\ldots A_{2k} be a convex polygon such that A1A3=A2A4=⋯=A2k−2A2k=A2k−1A1=A2kA2A_1A_3=A_2A_4=\cdots=A_{2k-2}A_{2k}=A_{2k-1}A_1=A_{2k}A_2. Prove that

A1B1=A2B2=⋯=A2kB2k.A_1B_1=A_2B_2=\cdots=A_{2k}B_{2k}.

Solutions of applications

Application 1. A triangle has side lengths 3\sqrt{3}, 5\sqrt{5}, and 7\sqrt{7}. Find its area.

Solution. According to this area formula, we have

Area=144⋅(3)2⋅(7)2−((3)2−(5)2+(7)2)2=144⋅3⋅7−(3−5+7)2=594.\mathrm{Area}=\dfrac{1}{4}\sqrt{4\cdot(\sqrt{3})^{2}\cdot(\sqrt{7})^{2}-\Big((\sqrt{3})^{2}-(\sqrt{5})^{2}+(\sqrt{7})^{2}\Big)^{2}}=\dfrac{1}{4}\sqrt{4\cdot 3\cdot 7-(3-5+7)^{2}}=\dfrac{\sqrt{59}}{4}.

Application 2. Consider a trapezoid with bases of lengths 5 and 9. If the sum of the squares of its diagonals is 196, find the sum of the squares of its legs.

Solution. Let the diagonals be ee and ff (see the figure). Given that e2+f2=196e^{2}+f^{2}=196. It is a known result that if e2+f2=(a+b)2e^{2}+f^{2}=(a+b)^{2}, where aa and bb are the bases then this condition is equivalent to the condition that the diagonals of the trapezoid must be perpendicular to each other. In this case the bases are 5 and 9, given that e2+f2=196=(5+9)2e^{2}+f^{2}=196=(5+9)^{2}. This means that the diagonals of this trapezoid must be perpendicular to each other. Then, we have

Area of trapezoid=e⋅f2.\text{Area of trapezoid}=\frac{e\cdot f}{2}.

On the other hand (by the considered area formula), we have

Area of trapezoid=144e2f2−(a2−52+c2−92)2.\text{Area of trapezoid}=\frac{1}{4}\sqrt{4e^{2}f^{2}-(a^{2}-5^{2}+c^{2}-9^{2})^{2}}.

Thus, we get

ef2=144e2f2−(a2−52+c2−92)2.\frac{ef}{2}=\frac{1}{4}\sqrt{4e^{2}f^{2}-(a^{2}-5^{2}+c^{2}-9^{2})^{2}}.

Squaring both side of this equation, we get

e2f24=116(4e2f2−(a2−52+c2−92)2).\frac{e^{2}f^{2}}{4}=\frac{1}{16}\big(4e^{2}f^{2}-(a^{2}-5^{2}+c^{2}-9^{2})^{2}\big).

This means that

(a2−52+c2−92)2=0.(a^{2}-5^{2}+c^{2}-9^{2})^{2}=0.

So, we get

a2+c2=52+92=106.a^{2}+c^{2}=5^{2}+9^{2}=106.

Application 3. The area of triangle with sides aa, bb, cc is equal to the area of triangle with sides aa, bb, dd, where d≠cd\neq c, and that a2+b2=16a^{2}+b^{2}=16, find c2+d2c^{2}+d^{2}.

Solution. As the areas of both triangles are equal, then we have

144a2b2−(a2+b2−c2)2=144a2b2−(a2+b2−d2)2.\frac{1}{4}\sqrt{4a^2b^2-(a^2+b^2-c^2)^2}=\frac{1}{4}\sqrt{4a^2b^2-(a^2+b^2-d^2)^2}.

Squaring both sides and canceling the expression 4a2b24a^{2}b^{2} on both side, we get

(a2+b2−c2)2=(a2+b2−d2)2.(a^2+b^2-c^2)^2=(a^2+b^2-d^2)^2.

Therefore,

(a2+b2−c2)2−(a2+b2−d2)2=0.(a^2+b^2-c^2)^2-(a^2+b^2-d^2)^2=0.

Note that this can be factored as follows:

(a2+b2−c2−(a2+b2−d2))⋅(a2+b2−c2+(a2+b2−d2))=0.\Big(a^2+b^2-c^2-(a^2+b^2-d^2)\Big)\cdot\Big(a^2+b^2-c^2+(a^2+b^2-d^2)\Big)=0.

Simplifying, we get

(c2−d2)⋅(2(a2+b2)−(c2+d2))=0.(c^2-d^2)\cdot\big(2(a^2+b^2)-(c^2+d^2)\big)=0.

Given that a2+b2=16a^{2}+b^{2}=16, hence

(c2−d2)⋅(2⋅16−(c2+d2))=0.(c^2-d^2)\cdot(2\cdot 16-(c^{2}+d^{2}))=0.

Given also that c≠dc\neq d, thus c2−d2≠0c^{2}-d^{2}\neq 0. Therefore,

32−(c2+d2)=0.32-(c^{2}+d^{2})=0.

We obtain

c2+d2=32.c^{2}+d^{2}=32.

Application 4. Each of two non-congruent triangles has two sides of lengths 5 and 7, and the triangles have equal circumradii. Find the product of the lengths of their third sides.

Solution of Application 3. Denote their areas by A1A_1 and A2A_2, denote the length of their third sides by aa and bb, then by this area formula, we have

A1=144⋅52⋅72−(52−a2+72)2,A_1=\dfrac{1}{4}\sqrt{4\cdot 5^{2}\cdot 7^{2}-(5^{2}-a^{2}+7^{2})^{2}},
A2=144⋅52⋅72−(52−b2+72)2.A_2=\dfrac{1}{4}\sqrt{4\cdot 5^{2}\cdot 7^{2}-(5^{2}-b^{2}+7^{2})^{2}}.

Denote by R1R_1 and R2R_2 their circumradii, then by the law of sines, we have

A1=5⋅7a4R1,A_1=\dfrac{5\cdot 7a}{4R_1},
A2=5⋅7b4R2.A_2=\dfrac{5\cdot 7b}{4R_2}.

As R1=R2R_1=R_2, we get

5⋅7a4A1=5⋅7b4A2.\dfrac{5\cdot 7a}{4A_1}=\dfrac{5\cdot 7b}{4A_2}.

Thus, it follows that

a4⋅52⋅72−(52−a2+72)2=b4⋅52⋅72−(52−b2+72)2.\dfrac{a}{\sqrt{4\cdot 5^{2}\cdot 7^{2}-(5^{2}-a^{2}+7^{2})^{2}}}=\dfrac{b}{\sqrt{4\cdot 5^{2}\cdot 7^{2}-(5^{2}-b^{2}+7^{2})^{2}}}.

Squaring both sides, we get

a2(4⋅52⋅72−(52−b2+72)2)=b2(4⋅52⋅72−(52−a2+72)2).a^{2}\Big(4\cdot 5^{2}\cdot 7^{2}-(5^{2}-b^{2}+7^{2})^{2}\Big)=b^{2}\Big(4\cdot 5^{2}\cdot 7^{2}-(5^{2}-a^{2}+7^{2})^{2}\Big).

Simplifying and using difference of squares formula, we get

(a2−b2)⋅4⋅52⋅72=(a(52−b2+72)−b(52−a2+72))(a(52−b2+72)+b(52−a2+72)).(a^{2}-b^{2})\cdot 4\cdot 5^{2}\cdot 7^{2}=\Big(a(5^{2}-b^{2}+7^{2})-b(5^{2}-a^{2}+7^{2})\Big)\Big(a(5^{2}-b^{2}+7^{2})+b(5^{2}-a^{2}+7^{2})\Big).

We obtain

(a2−b2)⋅4⋅52⋅72=(a−b)(52+ab+72)(a+b)(52−ab+72).(a^{2}-b^{2})\cdot 4\cdot 5^{2}\cdot 7^{2}=(a-b)(5^{2}+ab+7^{2})(a+b)(5^{2}-ab+7^{2}).

As triangles are different, then a≠ba\neq b. Canceling a2−b2a^{2}-b^{2} and (a−b)(a+b)(a-b)(a+b), we get

4⋅52⋅72=(52+72)2−(ab)2.4\cdot 5^{2}\cdot 7^{2}=(5^{2}+7^{2})^{2}-(ab)^{2}.

Hence

(ab)2=(52+72)2−4⋅52⋅72=(52+72−2⋅5⋅7)(52+72+2⋅5⋅7)=(25+49−70)(25+49+70)=4⋅144.(ab)^{2}=(5^{2}+7^{2})^{2}-4\cdot 5^{2}\cdot 7^{2}=(5^{2}+7^{2}-2\cdot 5\cdot 7)(5^{2}+7^{2}+2\cdot 5\cdot 7)=(25+49-70)(25+49+70)=4\cdot 144.

Therefore

ab=2⋅12=24.ab=2\cdot 12=24.

Application 5. Let a1a_1, a2a_2, a3a_3 be the side lengths of a triangle, prove that

1).   2a12a22+2a22a32+2a12a32>a14+a24+a34.1).\ \ \ 2a_1^{2}a_2^{2}+2a_2^{2}a_3^{2}+2a_1^{2}a_3^{2}>a_1^{4}+a_2^{4}+a_3^{4}.
2).   a12a2+a22a1+a22a3+a32a2+a12a3+a32a1>a13+a23+a33.2).\ \ \ a_1^{2}a_2+a_2^{2}a_1+a_2^{2}a_3+a_3^{2}a_2+a_1^{2}a_3+a_3^{2}a_1>a_1^{3}+a_2^{3}+a_3^{3}.

Proof of Application 5. We have

4a32a12−(a12−a22+a32)2>0,4a_3^{2}a_1^{2}-(a_1^{2}-a_2^{2}+a_3^{2})^{2}>0,

and

(a2+a3−a1)(a2+a1−a3)(a1+a3−a2)>0.(a_2+a_3-a_1)(a_2+a_1-a_3)(a_1+a_3-a_2)>0.

This ends the proof of both inequalities of Problem 3.

If n=4n = 4, then from (4) and from the fact that point B2B_2 coincides with point A4A_4 we get

(A1A2A3A4)=144A1A32⋅A2A42−(a12−a22+a32−a42)2.(A_1A_2A_3A_4)=\dfrac{1}{4}\sqrt{4A_1A_3^{2}\cdot A_2A_4^{2}-\Big(a_1^{2}-a_2^{2}+a_3^{2}-a_4^{2}\Big)^{2}}.
(5)

Problem 5. Let ABCDABCD be a convex quadrilateral, such that AB=aAB=a, BC=bBC=b, CD=cCD=c, AD=dAD=d, AC=eAC=e, BD=fBD=f. Prove that

∣a2−b2+c2−d2∣<2ef.\vert a^{2}-b^{2}+c^{2}-d^{2}\vert<2ef.

Proof of Application 5. The proof of Problem 4 follows from (5).

Application 7. Let A1A2…A2kA_1A_2\ldots A_{2k} be a convex polygon such that A1A3=A2A4=⋯=A2k−2A2k=A2k−1A1=A2kA2A_1A_3=A_2A_4=\cdots=A_{2k-2}A_{2k}=A_{2k-1}A_1=A_{2k}A_2. Prove that

A1B1=A2B2=⋯=A2kB2k.A_1B_1=A_2B_2=\cdots=A_{2k}B_{2k}.

Proof of Application 7. The proof of Problem 5 follows from the theorem.

Alternative proof of the lemma. Choose a coordinate system as in the solution to Problem 1.

Coordinate system for the alternative proof of the lemma

We have

AC2+BD2−AD2−BC2=(c+b)2+d2+(c−a)2+d2−c2−d2−(c+b−a)2−d2=2ab=2⋅AB⋅CD.AC^{2}+BD^{2}-AD^{2}-BC^{2}=(c+b)^{2}+d^{2}+(c-a)^{2}+d^{2}-c^{2}-d^{2}-(c+b-a)^{2}-d^{2}=2ab=2\cdot AB\cdot CD.

This ends the proof of the lemma.

References

  1. Acerbi, F., Vitrac, B., Heron of Alexandria: Metrica, Heron's formula. Serra, 2014.
  2. Bretschneider, C. A., Untersuchung der trigonometrischen Relationen des geradlinigen Viereckes. Archiv der Math. 2, 225–261, 1842.
  3. Coolidge, J. L., A Historically Interesting Formula for the Area of a Quadrilateral. Amer. Math. Monthly 46, 345–347, 1939.
  4. Heath, T. L., A History of Greek Mathematics, Vol. 2. Clarendon Press, 1921.
  5. Sedrakyan, H., Sedrakyan, N., Algebraic Inequalities. Springer International Publishing, 2018.
  6. Sedrakyan, H., Sedrakyan, N., Mathematics Competitions Must-Know Geometry: Practice and Theory Inspired by the Style of AMC and AIME Geometry (Unofficial Guide). 2023.

Download the full paper (PDF)

Explore more of the Sedrakyans' work on the home page.

© 2026 sedrakyans.org · Hayk Sedrakyan, Nairi Sedrakyan, Christophe Pinarbasi. All rights reserved.