Modern mathematics can obtain the area using the methods of integral calculus or its more sophisticated offspring, real analysis. However the area of a disk was studied by the Ancient Greeks. Eudoxus of Cnidus in the fifth century B.C. had found that the area of a disk is proportional to its radius squared. Archimedes used the tools of Euclidean geometry to show that the area inside a circle is equal to that of a right triangle whose base has the length of the circle's circumference and whose height equals the circle's radius in his book Measurement of a Circle. The circumference is 2πr, and the area of a triangle is half the base times the height, yielding the area πr^{2} for the disk. Prior to Archimedes, Hippocrates of Chios was the first to show that the area of a disk is proportional to the square of its diameter, as part of his quadrature of the lune of Hippocrates, but did not identify the constant of proportionality.
A variety of arguments have been advanced historically to establish the equation
A
=
π
r
2
of varying degrees of mathematical rigor. The most famous of these is Archimedes' method of exhaustion, one of the earliest uses of the mathematical concept of a limit, as well as the origin of Archimedes' axiom which remains part of the standard analytical treatment of the real number system. The original proof of Archimedes is not quite rigorous by modern standards, because it assumes that we can compare the length of arc of a circle to the length of a secant and a tangent line, and similar statements about the area, as geometrically evident.
The area of a regular polygon is half its perimeter times the apothem. As the number of sides of the regular polygon increases, the polygon tends to a circle, and the apothem tends to the radius. This suggests that the area of a disk is half the circumference of its bounding circle times the radius.
Following Archimedes & c. 260 BCE, compare the area enclosed by a circle to a right triangle whose base has the length of the circle's circumference and whose height equals the circle's radius. If the area of the circle is not equal to that of the triangle, then it must be either greater or less. We eliminate each of these by contradiction, leaving equality as the only possibility. We use regular polygons in the same way.
Suppose that the area C enclosed by the circle is greater than the area T = ^{1}⁄_{2}cr of the triangle. Let E denote the excess amount. Inscribe a square in the circle, so that its four corners lie on the circle. Between the square and the circle are four segments. If the total area of those gaps, G_{4}, is greater than E, split each arc in half. This makes the inscribed square into an inscribed octagon, and produces eight segments with a smaller total gap, G_{8}. Continue splitting until the total gap area, G_{n}, is less than E. Now the area of the inscribed polygon, P_{n} = C − G_{n}, must be greater than that of the triangle.
E
=
C
−
T
>
G
n
P
n
=
C
−
G
n
>
C
−
E
P
n
>
T
But this forces a contradiction, as follows. Draw a perpendicular from the center to the midpoint of a side of the polygon; its length, h, is less than the circle radius. Also, let each side of the polygon have length s; then the sum of the sides, ns, is less than the circle circumference. The polygon area consists of n equal triangles with height h and base s, thus equals ^{1}⁄_{2}nhs. But since h < r and ns < c, the polygon area must be less than the triangle area, ^{1}⁄_{2}cr, a contradiction. Therefore, our supposition that C might be greater than T must be wrong.
Suppose that the area enclosed by the circle is less than the area T of the triangle. Let D denote the deficit amount. Circumscribe a square, so that the midpoint of each edge lies on the circle. If the total area gap between the square and the circle, G_{4}, is greater than D, slice off the corners with circle tangents to make a circumscribed octagon, and continue slicing until the gap area is less than D. The area of the polygon, P_{n}, must be less than T.
D
=
T
−
C
>
G
n
P
n
=
C
+
G
n
<
C
+
D
P
n
<
T
This, too, forces a contradiction. For, a perpendicular to the midpoint of each polygon side is a radius, of length r. And since the total side length is greater than the circumference, the polygon consists of n identical triangles with total area greater than T. Again we have a contradiction, so our supposition that C might be less than T must be wrong as well.
Therefore, it must be the case that the area enclosed by the circle is precisely the same as the area of the triangle. This concludes the proof.
Following Satō Moshun (Smith & Mikami 1914, pp. 130–132) and Leonardo da Vinci (Beckmann 1976, p. 19), we can use inscribed regular polygons in a different way. Suppose we inscribe a hexagon. Cut the hexagon into six triangles by splitting it from the center. Two opposite triangles both touch two common diameters; slide them along one so the radial edges are adjacent. They now form a parallelogram, with the hexagon sides making two opposite edges, one of which is the base, s. Two radial edges form slanted sides, and the height is h (as in the Archimedes proof). In fact, we can assemble all the triangles into one big parallelogram by putting successive pairs next to each other. The same is true if we increase to eight sides and so on. For a polygon with 2n sides, the parallelogram will have a base of length ns, and a height h. As the number of sides increases, the length of the parallelogram base approaches half the circle circumference, and its height approaches the circle radius. In the limit, the parallelogram becomes a rectangle with width πr and height r.
There are various equivalent definitions of the constant π. The conventional definition in precalculus geometry is the ratio of the circumference of a circle to its diameter:
π
=
C
D
.
However, because the circumference of a circle is not a primitive analytical concept, this definition is not suitable in modern rigorous treatments. A standard modern definition is that π is equal to twice the least positive root of the cosine function or, equivalently, the halfperiod of the sine (or cosine) function. The cosine function can be defined either as a power series, or as the solution of a certain differential equation. This avoids any reference to circles in the definition of π, so that statements about the relation of π to the circumference and area of circles are actually theorems, rather than definitions, that follow from the analytical definitions of concepts like "area" and "circumference".
The analytical definitions are seen to be equivalent, if it is agreed that the circumference of the circle is measured as a rectifiable curve by means of the integral
C
=
2
∫
−
R
R
R
d
x
R
2
−
x
2
=
2
R
∫
−
1
1
d
x
1
−
x
2
.
The integral appearing on the right is an abelian integral whose value is a halfperiod of the sine function, equal to π. Thus
C
=
2
π
R
=
π
D
is seen to be true as a theorem.
Several of the arguments that follow use only concepts from elementary calculus to reproduce the formula
A
=
π
r
2
, but in many cases to regard these as actual proofs, they rely implicitly on the fact that one can develop trigonometric functions and the fundamental constant π in a way that is totally independent of their relation to geometry. We have indicated where appropriate how each of these proofs can be made totally independent of all trigonometry, but in some cases that requires more sophisticated mathematical ideas than those afforded by elementary calculus.
Using calculus, we can sum the area incrementally, partitioning the disk into thin concentric rings like the layers of an onion. This is the method of shell integration in two dimensions. For an infinitesimally thin ring of the "onion" of radius t, the accumulated area is 2πt dt, the circumferential length of the ring times its infinitesimal width (one can approximate this ring by a rectangle with width=2πt and height=dt). This gives an elementary integral for a disk of radius r.
A
r
e
a
(
r
)
=
∫
0
r
2
π
t
d
t
=
2
π
[
t
2
2
]
0
r
=
π
r
2
.
It is rigorously justified by the multivariate substitution rule in polar coordinates. Namely, the area is given by a double integral of the constant function 1 over the disk itself. If D denotes the disk, then the double integral can be computed in polar coordinates as follows:
A
r
e
a
(
r
)
=
∬
D
1
d
(
x
,
y
)
=
∬
D
t
d
t
d
θ
=
∫
0
r
∫
0
2
π
t
d
θ
d
t
=
∫
0
r
[
t
θ
]
0
2
π
d
t
=
∫
0
r
2
π
t
d
t
which is the same result as obtained above.
An equivalent rigorous justification, without relying on the special coordinates of trigonometry, uses the coarea formula. Define a function
ρ
:
R
2
→
R
by
ρ
(
x
,
y
)
=
x
2
+
y
2
. Note ρ is a Lipschitz function whose gradient is a unit vector

∇
ρ

=
1
(almost everywhere). Let D be the disc
ρ
<
1
in
R
2
. We will show that
L
2
(
D
)
=
π
, where
L
2
is the twodimensional Lebesgue measure in
R
2
. We shall assume that the onedimensional Hausdorff measure of the circle
ρ
=
r
is
2
π
r
, the circumference of the circle of radius r. (This can be taken as the definition of circumference.) Then, by the coarea formula,
L
2
(
D
)
=
∬
D

∇
ρ

d
L
2
=
∫
R
H
1
(
ρ
−
1
(
r
)
∩
D
)
d
r
=
∫
0
1
H
1
(
ρ
−
1
(
r
)
)
d
r
=
∫
0
1
2
π
r
d
r
=
π
.
Similar to the onion proof outlined above, we could exploit calculus in a different way in order to arrive at the formula for the area of a disk. Consider unwrapping the concentric circles to straight strips. This will form a right angled triangle with r as its height and 2πr (being the outer slice of onion) as its base.
Finding the area of this triangle will give the area of the disk
Area
=
1
2
⋅
base
⋅
height
=
1
2
⋅
2
π
r
⋅
r
=
π
r
2
The opposite and adjacent angles for this triangle are respectively in degrees 9.0430611..., 80.956939... and in radians 0.1578311... A233527, 1.4129651... A233528.
Explicitly, we imagine dividing up a circle into triangles, each with a height equal to the circle's radius and a base that is infinitesimally small. The area of each of these triangles is equal to
1
/
2
⋅
r
⋅
d
u
. By summing up (integrating) all of the areas of these triangles, we arrive at the formula for the circle's area:
A
r
e
a
(
r
)
=
∫
0
2
π
r
1
2
r
d
u
=
[
1
2
r
u
]
0
2
π
r
=
π
r
2
.
It too can be justified by a double integral of the constant function 1 over the disk by reversing the order of integration and using a change of variables in the above iterated integral:
A
r
e
a
(
r
)
=
∬
D
1
d
(
x
,
y
)
=
∬
D
t
d
t
d
θ
=
∫
0
2
π
∫
0
r
t
d
t
d
θ
=
∫
0
2
π
1
2
r
2
d
θ
Making the substitution
u
=
r
θ
,
d
u
=
r
d
θ
converts the integral to
∫
0
2
π
r
1
2
r
2
r
d
u
=
∫
0
2
π
r
1
2
r
d
u
which is the same as the above result.
The triangle proof can be reformulated as an application of Green's theorem in fluxdivergence form (i.e. a twodimensional version of the divergence theorem), in a way that avoids all mention of trigonometry and the constant π. Consider the vector field
r
=
x
i
+
y
j
in the plane. So the divergence of r is equal to two, and hence the area of a disc D is equal to
A
=
1
2
∬
D
div
r
d
A
.
By Green's theorem, this is the same as the outward flux of r across the circle bounding D:
A
=
1
2
∮
∂
D
r
⋅
n
d
s
where n is the unit normal and ds is the arc length measure. For a circle of radius R centered at the origin, we have

r

=
R
and
n
=
r
/
R
, so the above equality is
A
=
1
2
∮
∂
D
r
⋅
r
R
d
s
=
R
2
∮
∂
D
d
s
.
The integral of ds over the whole circle
∂
D
is just the arc length, which is its circumference, so this shows that the area A enclosed by the circle is equal to
R
/
2
times the circumference of the circle.
Note that the area of a semicircle of radius r can be computed by the integral
∫
−
r
r
r
2
−
x
2
d
x
.
By trigonometric substitution, we substitute
x
=
r
sin
θ
d
x
=
r
cos
θ
d
θ
θ
=
arcsin
(
x
r
)
=
∫
−
π
2
π
2
r
2
(
1
−
sin
2
θ
)
⋅
r
cos
θ
d
θ
=
2
r
2
∫
0
π
2
cos
2
θ
d
θ
=
π
r
2
2
.
The last step follows since, the trigonometric identity
cos
(
θ
)
=
sin
(
π
/
2
−
θ
)
implies that
cos
2
θ
and
sin
2
θ
have equal integrals over the interval
[
0
,
π
/
2
]
, using integration by substitution. But on the other hand, since
cos
2
θ
+
sin
2
θ
=
1
, the sum of the two integrals is the length of that interval, which is
π
/
2
. Consequently, the integral of
cos
2
θ
is equal to half the length of that interval, which is
π
/
4
.
Therefore, the area of a circle of radius r, which is twice the area of the semicircle, is equal to
2
⋅
π
r
2
2
=
π
r
2
.
This particular proof may appear to beg the question, if the sine and cosine functions involved in the trigonometric substitution are regarded as being defined in relation to circles. However, as noted earlier, it is possible to define sine, cosine, and π in a way that is totally independent of trigonometry, in which case the proof is valid by the change of variables formula and Fubini's theorem, assuming the basic properties of sine and cosine (which can also be proved without assuming anything about their relation to circles).
The circle is the closed curve of least perimeter that encloses the maximum area. This is known as the isoperimetric inequality, which states that if a rectifiable Jordan curve in the Euclidean plane has perimeter C and encloses an area A (by the Jordan curve theorem) then
4
π
A
≤
C
2
.
Moreover, equality holds in this inequality if and only if the curve is a circle, in which case
A
=
π
r
2
and
C
=
2
π
r
.
The calculations Archimedes used to approximate the area numerically were laborious, and he stopped with a polygon of 96 sides. A faster method uses ideas of Willebrord Snell (Cyclometricus, 1621), further developed by Christiaan Huygens (De Circuli Magnitudine Inventa, 1654), described in Gerretsen & Verdenduin (1983, pp. 243–250).
Given a circle, let u_{n} be the perimeter of an inscribed regular ngon, and let U_{n} be the perimeter of a circumscribed regular ngon. Then u_{n} and U_{n} are lower and upper bounds for the circumference of the circle that become sharper and sharper as n increases, and their average (u_{n} + U_{n})/2 is an especially good approximation to the circumference. To compute u_{n} and U_{n} for large n, Archimedes derived the following doubling formulae:
u
2
n
=
U
2
n
u
n
(geometric mean), and
U
2
n
=
2
U
n
u
n
U
n
+
u
n
(harmonic mean).
Starting from a hexagon, Archimedes doubled n four times to get a 96gon, which gave him a good approximation to the circumference of the circle.
In modern notation, we can reproduce his computation (and go further) as follows. For a unit circle, an inscribed hexagon has u_{6} = 6, and a circumscribed hexagon has U_{6} = 4√3. Doubling seven times yields
(Here (u_{n} + U_{n})/2 approximates the circumference of the unit circle, which is 2π, so (u_{n} + U_{n})/4 approximates π.)
The last entry of the table has ^{355}⁄_{113} as one of its best rational approximations; i.e., there is no better approximation among rational numbers with denominator up to 113. The number ^{355}⁄_{113} is also an excellent approximation to π, better than any other rational number with denominator less than 16604.
Snell proposed (and Huygens proved) a tighter bound than Archimedes':
n
3
sin
π
n
2
+
cos
π
n
<
π
<
n
(
2
sin
π
3
n
+
tan
π
3
n
)
.
This for n = 48 gives a better approximation (about 3.14159292) than Archimedes' method for n = 768.
Let one side of an inscribed regular ngon have length s_{n} and touch the circle at points A and B. Let A′ be the point opposite A on the circle, so that A′A is a diameter, and A′AB is an inscribed triangle on a diameter. By Thales' theorem, this is a right triangle with right angle at B. Let the length of A′B be c_{n}, which we call the complement of s_{n}; thus c_{n}^{2}+s_{n}^{2} = (2r)^{2}. Let C bisect the arc from A to B, and let C′ be the point opposite C on the circle. Thus the length of CA is s_{2n}, the length of C′A is c_{2n}, and C′CA is itself a right triangle on diameter C′C. Because C bisects the arc from A to B, C′C perpendicularly bisects the chord from A to B, say at P. Triangle C′AP is thus a right triangle, and is similar to C′CA since they share the angle at C′. Thus all three corresponding sides are in the same proportion; in particular, we have C′A : C′C = C′P : C′A and AP : C′A = CA : C′C. The center of the circle, O, bisects A′A, so we also have triangle OAP similar to A′AB, with OP half the length of A′B. In terms of side lengths, this gives us
c
2
n
2
=
(
r
+
1
2
c
n
)
2
r
c
2
n
=
s
n
s
2
n
.
In the first equation C′P is C′O+OP, length r+^{1}⁄_{2}c_{n}, and C′C is the diameter, 2r. For a unit circle we have the famous doubling equation of Ludolph van Ceulen,
c
2
n
=
2
+
c
n
.
If we now circumscribe a regular ngon, with side A″B″ parallel to AB, then OAB and OA″B″ are similar triangles, with A″B″ : AB = OC : OP. Call the circumscribed side S_{n}; then this is S_{n} : s_{n} = 1 : ^{1}⁄_{2}c_{n}. (We have again used that OP is half the length of A′B.) Thus we obtain
c
n
=
2
s
n
S
n
.
Call the inscribed perimeter u_{n} = ns_{n}, and the circumscribed perimenter U_{n} = nS_{n}. Then combining equations, we have
c
2
n
=
s
n
s
2
n
=
2
s
2
n
S
2
n
,
so that
u
2
n
2
=
u
n
U
2
n
.
This gives a geometric mean equation.
We can also deduce
2
s
2
n
S
2
n
s
n
s
2
n
=
2
+
2
s
n
S
n
,
or
2
U
2
n
=
1
u
n
+
1
U
n
.
This gives a harmonic mean equation.
When more efficient methods of finding areas are not available, we can resort to “throwing darts”. This Monte Carlo method uses the fact that if random samples are taken uniformly scattered across the surface of a square in which a disk resides, the proportion of samples that hit the disk approximates the ratio of the area of the disk to the area of the square. This should be considered a method of last resort for computing the area of a disk (or any shape), as it requires an enormous number of samples to get useful accuracy; an estimate good to 10^{−n} requires about 100^{n} random samples (Thijsse 2006, p. 273).
We have seen that by partitioning the disk into an infinite number of pieces we can reassemble the pieces into a rectangle. A remarkable fact discovered relatively recently (Laczkovich 1990) is that we can dissect the disk into a large but finite number of pieces and then reassemble the pieces into a square of equal area. This is called Tarski's circlesquaring problem. The nature of Laczkovich's proof is such that it proves the existence of such a partition (in fact, of many such partitions) but does not exhibit any particular partition.
Circles can be defined in nonEuclidean geometry, and in particular in the hyperbolic and elliptic planes.
For example, the unit sphere
S
2
(
1
)
is a model for the twodimensional elliptic plane. It carries an intrinsic metric that arises by measuring geodesic length. The geodesic circles are the parallels in a geodesic coordinate system.
More precisely, fix a point
z
∈
S
2
(
1
)
that we place at the zenith. Associated to that zenith is a geodesic polar coordinate system
(
ϕ
,
θ
)
,
0
≤
ϕ
≤
π
,
0
≤
θ
<
2
π
, where z is the point
ϕ
=
0
. In these coordinates, the geodesic distance from z to any other point
x
∈
S
2
(
1
)
having coordinates
(
ϕ
,
θ
)
is the value of
ϕ
at x. A spherical circle is the set of points a geodesic distance R from the zenith point z. Equivalently, with a fixed embedding into
R
3
, the spherical circle of radius
R
≤
π
centered at z is the set of x in
S
2
(
1
)
such that
x
⋅
z
=
cos
R
.
We can also measure the area of the spherical disk enclosed within a spherical circle, using the intrinsic surface area measure on the sphere. The area of the disk of radius R is then given by
A
=
∫
0
2
π
∫
0
R
sin
(
ϕ
)
d
ϕ
d
θ
=
2
π
(
1
−
cos
R
)
.
More generally, if a sphere
S
2
(
ρ
)
has radius of curvature
ρ
, then the area of the disk of radius R is given by
A
=
2
π
ρ
2
(
1
−
cos
(
R
/
ρ
)
)
.
Observe that, as an application of L'Hôpital's rule, this tends to the Euclidean area
π
R
2
in the flat limit
ρ
→
∞
.
The hyperbolic case is similar, with the area of a disk of intrinsic radius R in the (constant curvature
−
1
) hyperbolic plane given by
A
=
2
π
(
1
−
cosh
R
)
where cosh is the hyperbolic cosine. More generally, for the constant curvature
−
k
hyperbolic plane, the answer is
A
=
2
π
k
−
2
(
1
−
cosh
(
k
R
)
)
.
These identities are important for comparison inequalities in geometry. For example, the area enclosed by a circle of radius R in a flat space is always greater than the area of a spherical circle and smaller than a hyperbolic circle, provided all three circles have the same (intrinsic) radius. That is,
2
π
(
1
−
cos
R
)
<
π
R
2
<
2
π
(
1
−
cosh
R
)
for all
R
>
0
. Intuitively, this is because the sphere tends to curve back on itself, yielding circles of smaller area than those in the plane, whilst the hyperbolic plane, when immersed into space, develops fringes that produce additional area. It is more generally true that the area of the circle of a fixed radius R is a strictly decreasing function of the curvature.
In all cases, if
k
is the curvature (constant, positive or negative), then the isoperimetric inequality for a domain with area A and perimeter L is
L
2
≥
4
π
A
−
k
A
2
where equality is achieved precisely for the circle.
We can stretch a disk to form an ellipse. Because this stretch is a linear transformation of the plane, it has a distortion factor which will change the area but preserve ratios of areas. This observation can be used to compute the area of an arbitrary ellipse from the area of a unit circle.
Consider the unit circle circumscribed by a square of side length 2. The transformation sends the circle to an ellipse by stretching or shrinking the horizontal and vertical diameters to the major and minor axes of the ellipse. The square gets sent to a rectangle circumscribing the ellipse. The ratio of the area of the circle to the square is π/4, which means the ratio of the ellipse to the rectangle is also π/4. Suppose a and b are the lengths of the major and minor axes of the ellipse. Since the area of the rectangle is ab, the area of the ellipse is πab/4.
We can also consider analogous measurements in higher dimensions. For example, we may wish to find the volume inside a sphere. When we have a formula for the surface area, we can use the same kind of “onion” approach we used for the disk.