## A New Visual Approach to Pendulum Period Determination

R. Sánchez-Martínez , 1 ∗ E. Heredia-Muñoz 2

1 Instituto de Física, Universidad Nacional Autónoma de México. Coyoacán, Cd. de México. México. 2 Facultad de Matemáticas, Universidad Veracruzana. Xalapa, Ver. México.

August 2, 2024

## Abstract

The period of oscillation of a simple pendulum ( T = 2 π √ l/g ) is a familiar formula to the average firstyear physics student. However, deriving this expression from first principles involves solving a non-linear differential equation using the small angle approximation, which may appear obscure to a student in the early stages of learning calculus. Therefore, we propose an alternative approach to the derivation of this formula employing geometry, algebra, and physical intuition. Our method follows the ideas of the fathers of calculus, replacing the circular path of the pendulum with a successive collection of infinitesimal inclined planes and summing the travel times of each plane in the limit as the number of planes becomes very large. Notably, the evaluation of this limit relies solely on geometry, making it accessible to any student, even those not yet familiar with calculus techniques.

1

## Introduction and Motivation

The simple pendulum is one of the most extensively studied problems in physics[1]. It consists of a point mass attached to a fixed string, enabling it to swing freely in a back-and-forth motion under the influence of gravity. This simple device captivated the leading figures of the 17th century such as Galileo, Huygens and Newton, and played a crucial role in several key areas: establishing the laws of collision and conservation of energy and momentum, determining the value of the acceleration due to gravity, and, perhaps most importantly, providing evidence for Newtonian synthesis of terrestrial and celestial mechanics[2].

It is well known that, for small initial angular displacements, the oscillating motion of a pendulum can be accurately described by a simple harmonic motion with period the amount of mass attached? Why do longer pendulums swing more slowly than shorter ones? And-most importantly for this work-what is the origin of the 2 π factor in the formula? At first glance, one might think it is directly related to the circular motion of the pendulum. However, it is evident that the mass does not complete a full circle during an oscillation-which would account for a factor of 2 πl or something similar. Instead, it moves along small segments of arcs in a non-uniform manner. Hence, if there is indeed a relationship to a circle, it is not the circle traced by the pendulum's trajectory. In fact, the appearance of this factor is due to the periodic relationship of the trigonometric functions used to describe its motion. Nonetheless, we would like to emphasize that even without considering any periodicity in the description of the system-as we will demonstrate in this work-a factor of π emerges!

$$T = 2 \pi \sqrt { \frac { l } { g } }. \quad \quad \ \ ( 1. 1 ) \quad \ \.$$

This is undoubtedly one of the first formulas a physics student learns in their introductory year. When it is first presented in the classroom, some questions may arise: Why doesn't the period depend on

∗ rodrigo\_smtz@comunidad.unam.mx

Deriving equation (1.1) from first principles involves writting down the equations of motion of the mass. This leads to a second-order non-linear differential equation in terms of the angle, with no analytical solution for arbitrary initial conditions. Only for sufficiently small initial angular displacements, this differ-

ential equation can be approximated to that of a simple harmonic oscillator, which has a well-known solution. As a homogeneous linear equation with constant coefficients, we can solve it using the so-called auxiliary equation method[3], resulting in the familiar periodic solutions in terms of sine and cosine, from which we can extract eq. (1.1). While this procedure would serve as an instructive exercise in a basic course on differential equations, for a student who is just learning the basics of differential calculus, this method may seem abstract and lacking in much physical interest.

The question now arises: is there an alternative procedure to derive expression (1.1) without the use of calculus? We claim that there is.

In this work, we propose a novel approach that relies in the physics and geometry of the problem. The idea is simple: following in the footsteps of Newton and Leibniz, we seek to apply the concept of partitioning a complex problem into infinitesimally small, manageable portions to the study of the pendulum's trajectory. Our approach involves replacing the circular path of the pendulum with a series of infinitesimal inclined planes due to the simplicity of the kinematic equations governing the motion on each plane. Then, by summing the travel times of these planes and taking the limit as their number approaches infinity, we can determine the descent time of the pendulum and, consequently, its period of oscillation. 1

offers some annotations on the usual modern derivation.

We hope that this work is of great utility and interest to both students and teachers interested in solving mechanics problems using geometric techniques. In particular, we encourage teachers to incorporate this derivation into their basic physics classes, inviting them to reproduce our diagrams, since we believe that choosing the appropriate visual elements is an important part of solving a problem and generating meaningful learning.

## 2 Our New Derivation

## 2.1 Procedure

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-e36ca03f-d606-493b-bc37-68cde674097e.jpg)

B

Figure 2.1: Geometry of a simple pendulum.

The structure of the work is as follows: First, we describe the partition of the circular path of the pendulum into small inclined planes and calculate the time it takes for a particle to traverse each of them using basic kinematic formulas. By summing these individual time intervals, we arrive to a general expression for the approximation to the total descent time in terms of the number of inclined planes used. The solution involves taking the limit as this number tends to infinity, which could be done, in principle, with calculus techniques. However, in the spirit of maintaining a calculus-free solution, we have devised a visual method to evaluate this limit using a little analytical geometry. This evaluation is provided in the second part of this section. In Section 3, we present our final remarks, followed by two appendices. The first provides a historical note on the derivation of equation (1.1) in the 17th century, while the second

Let's examine the motion of a pendulum, as depicted in Fig. 2.1, with length l and mass m . It starts from rest at point A and falls to point B , covering an arc length &gt; AB = ∆ l ϕ . At any point P of the trajectory, it is sufficient to know the angle ϕ = ∠ BOP to determine its position completely.

As depicted in Fig. 2.2, our strategy is to divide the trajectory &gt; AB into N inclined planes, each of equal length δx and characterized by an inclination angle θ n . With this approach, we aim to break down the overall motion into manageable segments for analysis. Given that the mass experiences a constant acceleration while traversing each plane, it becomes straightforward to formulate an expression for the time interval δ n t required for the body to traverse the n -th plane. Consequently, an approximation to the total falling time from A to B is obtained by summing these intervals, and the precise value emerges as we consider an increasingly large number of inclined planes.

1 For those already familiar with calculus, this may sound contradictory, as it seems we've just described solving a problem using standard integral calculus. However, although the idea is practically the same, we will demonstrate that the expression for the total descent time of the pendulum cannot be written as a Riemann sum and, therefore, cannot be evaluated as a definite integral. This makes our proposed method fundamentally different.

= 1

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-2abd1de4-8aba-4da8-af1a-36c96f74ac2d.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-43e01e9e-d662-4e7d-80c6-fc93c6afed9d.jpg)

= 2

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-f1243a65-0cc3-450d-8f62-aff3787082f3.jpg)

= 3

Figure 2.2: Approximation of the circular path followed by the pendulum by a succesion of inclined planes.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-ec6e09cf-13ce-4e23-ad63-1f35cd8ca5b4.jpg)

Figure 2.3: Geometry of the n -th inclined plane.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-feec19fd-6887-464b-8f88-1d9b7930513e.jpg)

The procedure is as follows. Let { x , x 0 1 , . . . , x N } be a partition of the arc &gt; AB , where x 0 = A and x N = B for a positive integer N . With the help of the diagram shown in Figure 2.3, let's analyze the motion of the mass while descending through the n -th inclined plane, represented by the line x n -1 x n for n = 1 2 , , . . . , N . Each plane has a different inclination angle θ n but subtends the same central angle δϕ = ∆ ϕ/N from O and has the same length δx = 2 sin( l δϕ/ 2) . Also, note that the points x n are located at an angular distance ϕ n = ( N -n δϕ ) from the vertical, so that ϕ 0 = ∆ ϕ and ϕ N = 0 .

An expression for the inclination angles of the planes can be found using elementary geometry. Let ¯ x n be the midpoint between x n -1 and x n and consider the auxiliar points ξ n and ξ ′ n , both lying on the vertical from O , as shown in Fig. 2.3. Since the triangles △ x ξ n n ξ ′ n and △ Ox ξ ¯ n ′ n are similar, we conclude that the angles ∠ ξ n x ξ n ′ n and ∠ ξ ′ n Ox ¯ n are equal to each other. Therefore,

$$\theta _ { n } = \phi _ { n } + \frac { 1 } { 2 } \delta \phi = \left ( N - \frac { 2 n - 1 } { 2 } \right ) \delta \phi. \quad \ ( 2. 1 )$$

It is well-known that the acceleration experienced by the mass during its descent along each plane is given by

$$a _ { n } = g \sin \theta _ { n }. \text{ \quad \ \ } ( 2. 2 )$$

This allows us to calculate the change in speed after each plane using basic kinematic formulas. Let v n be the instantaneous speed acquired after moving past the n -th plane, with v 0 = 0 since the motion starts from rest. We can calculate the value of the speed after traversing n planes as follows

$$\text{ after traversing $n$ planes as follows} \\ v _ { 0 } ^ { 2 } = 0 \\ v _ { 1 } ^ { 2 } = v _ { 0 } ^ { 2 } + 2 a _ { 1 } \delta x = 2 a _ { 1 } \delta x \\ v _ { 2 } ^ { 2 } = v _ { 1 } ^ { 2 } + 2 a _ { 2 } \delta x = 2 ( a _ { 1 } + a _ { 2 } ) \delta x \\ v _ { 3 } ^ { 2 } = v _ { 2 } ^ { 2 } + 2 a _ { 3 } \delta x = 2 ( a _ { 1 } + a _ { 2 } + a _ { 3 } ) \delta x \\ \vdots \\ v _ { n } ^ { 2 } = v _ { n - 1 } ^ { 2 } + 2 a _ { n } \delta x = 2 ( a _ { 1 } + \cdots + a _ { n } ) \delta x. \quad ( 2. 3 ) \\ \text{ Thus, the value of $n$... demands denends on our abil-}$$

Thus, the value of v n depends depends on our ability to evaluate the sum ∑ n k a k = g ∑ n k sin θ k , which is not an easy task for arbitrary angles θ k . However, an analytical closed form of the sum can be obtained if the inclination angles are small enough to allow the approximation sin θ k ≈ θ k . We will show that this is precisely our case. Up to this point, we have described the problem in the most general terms without specifying any initial conditions for the movement. Yet, it is important to remember that the expression of the period (1.1) is valid only under the small-angle approximation, 0 &lt; ∆ ϕ ≪ 1 . Now, turning back to the evaluation of the sum, since the inclination angles change from plane to plane, it suffices to demonstrate that the largest angle θ n becomes negligible as ∆ ϕ approaches zero, allowing us to use θ n itself instead of its

sine. According to eq. (2.1), the maximum inclination angle will be θ 1 = (1 -1 2 N )∆ ϕ , which is always less or equal than ∆ ϕ . Therefore, for sufficiently small initial angular displacements, the small-angle approximation for every θ k is also valid and (2.2) becomes

$$a _ { n } \approx g \theta _ { n } = g \left ( N - \frac { 2 n - 1 } { 2 } \right ) \delta \phi. \quad \text{(2.4)}$$

Now, we are able to evaluate the sum and obtain

$$\sum _ { k = 1 } ^ { n } a _ { k } = g \sum _ { k = 1 } ^ { n } \theta _ { k } = g \left ( N n - \frac { n ^ { 2 } } { 2 } \right ) \delta \phi. \quad \text{(2.5)}$$

Additionally, in this limit, the length of each plane simplifies to δx = lδϕ . Substituting this result and eq. (2.5) into (2.3) gives

From the definition of a motion under constant acceleration, we can calculate the time required for the mass to traverse the n -th plane using

$$\delta _ { n } t = \frac { v _ { n } - v _ { n - 1 } } { a _ { n } }. \quad \quad \quad ( 2. 7 )$$

In conclusion, the time taken for the actual descent from point A to B is determined by

$$t _ { A B } = \lim _ { N \to \infty } \sum _ { n = 1 } ^ { N } \delta _ { n } t. \text{ \quad \ \ } \text{ (2.8)}$$

$$v _ { n } = \sqrt { \left ( N ^ { 2 } - ( N - n ) ^ { 2 } \right ) g l } \, \delta \phi. \quad \ \ ( 2. 6 ) \quad \ \ d u \,.$$

Substituting (2.4) and (2.6) into this expression yields our final result for the total descent time of the pendulum, displayed in eq. (2.9) below.

$$t _ { A B } = 2 \lim _ { N \to \infty } \underbrace { \left ( \sum _ { n = 1 } ^ { N } \frac { \sqrt { N ^ { 2 } - ( N - n ) ^ { 2 } } - \sqrt { N ^ { 2 } - ( N - ( n - 1 ) ) ^ { 2 } } } { 2 N - ( 2 n - 1 ) } \right ) } _ { F ( N ) } \sqrt { \frac { l } { g } }.$$

Hence, the calculation of the period of the pendulum depends on evaluating of the limit of the sum enclosed in parentheses in (2.9). For our result to be consistent with the one well-known formula (1.1), the following limit must be true naturally in the proof, serving as a reminder that the resulting equations that describe the motion of a physical system are unique and independent of the formalism applied to derive them.

$$\lim _ { N \to \infty } F ( N ) = \frac { \pi } { 4 }. \quad \quad ( 2. 1 0 ) \quad \ \ g _ { i } ^ { \text{ ti} } \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The evaluation of complex limits is a broad and fascinating subject within calculus courses, often involving specialized techniques like L'Hôpital's rule. In our commitment to making the solution accessible to firstyear students, we will now proceed with an illustrative geometrical proof. 2

## 2.2 Proof of (2.10)

Before starting the proof, let's recall that we arrived to equation (2.9) under the only assumption that the pendulum's fall time can be determined by dividing its trajectory into increasingly smaller inclined planes. At no point have we introduced the notion that this is a periodic motion -described by a sine or consine function-, and therefore, there is no apparent reason for the factor π to appear. Yet, it emerges

While the expression of F N ( ) in (2.9) may initially appear complex, the limit (2.10) can be elegantly evaluated through basic analytical geometry and the properties of angles in the circunference. With this in mind, we proceed to the proof.

Let O (0 0) , be the origin of coordinates in the plane, N a positive integer, and Q 0 ( N, 0) and

$$\text{ve} \quad P _ { n } \left ( - ( N - n + 1 ), \sqrt { N ^ { 2 } - ( N - n + 1 ) ^ { 2 } } \right ) \ ( 2. 1 1 \text{a} )$$

$$Q _ { n } \left ( N - n, \sqrt { N ^ { 2 } - ( N - n ) ^ { 2 } } \right ) \quad \quad \quad ( 2. 1 1 b )$$

for n = 1 2 , , . . . , N , points on the circunference of radius N and centered at O . Then, the points P n and Q n -1 are symmetric with respect to the vertical axis, such that P Q n n -1 is an horizontal segment. The general representation of such points in the circunference is shown in Fig. (2.4), where the specific cases for N = 1 2 , and 3 are depicted for illustration.

From the definition of the points in (2.11a) and (2.11b), the slope s n of each segment P Q n n will be

2 Now, it is evident that the limit (2.10) cannot be evaluated using integral calculus, as the sum F N ( ) does not take the form of a Riemann sum. This is primarly because each term of in sum depends on both n and N .

Figure 2.4: Visual aid for the proof of (2.10).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-c68c7ce5-4ea6-4c82-a9f5-c983ff1bed35.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-c2b88002-f75d-49e5-9ed7-3448fd4fbd6e.jpg)

given by

$$s _ { n } & = \frac { \sqrt { N ^ { 2 } - ( N - n ) ^ { 2 } } - \sqrt { N ^ { 2 } - ( N - n + 1 ) ^ { 2 } } } { ( N - n ) - \left ( - ( N - n + 1 ) \right ) }, \quad \begin{array} { c } \text{to h} \\ \text{impl} \end{array} \\ \dots & \dots & \dots & \dots & \dots$$

But also, we know that any inscribed angle is equal to half of the central angle that spans the same arc, implying that which is the same expression as the terms in the sum of (2.9). That is,

$$F ( N ) = \sum _ { n = 1 } ^ { N } s _ { n }. \quad \quad ( 2. 1 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-dff0046f-43db-4e0d-a6bd-145f5ddb84ad.jpg)

N

Figure 2.5: Geometrical representation of eq. (2.14). The particular case for N = 2 (above) and the general case for arbitrary N (below) are shown. Note that the angle ∠ Q OQ 0 N is a right angle for all positive integer N .

On the other hand, the slopes s n can also be written in terms of the respective elevation angle α n = ∠ Q n -1 P Q n n using s n = tan α n . Now, given that inscribed angles in a circunference, embracing the same arc, are equal (see Fig. 2.5), we can consider all the inscribed angles α n sharing a common vertex at P 1 to obtain

$$\sum _ { n = 1 } ^ { N } \alpha _ { n } = \sum _ { n = 1 } ^ { N } \zeta Q _ { n - 1 } P _ { 1 } Q _ { n } = \zeta Q _ { 0 } P _ { 1 } Q _ { N }. \quad ( 2. 1 4 ) \quad \begin{pmatrix} 1. 1 \\ \zeta \\ \text{entities} \end{pmatrix}$$

$$\zeta Q _ { 0 } P _ { 1 } Q _ { N } = \frac { 1 } { 2 } \zeta Q _ { 0 } O Q _ { N } = \frac { \pi } { 4 }. \quad \quad ( 2. 1 5 )$$

It is important to note that this result is independent of the number of points in the partition of the arc, i.e., it is valid for any positive integer N .

This provides us with an expression for the sum of the angles, but we still need an expression for the sum of the slopes, ∑ N n tan α n . Once again, we cannot evaluate this sum analytically for arbitrary angles α n . However, as the angles α n approach zero with the radius N of the circumference tending to infinity (which occurs when considering a very large number of inclined planes), we can use again the small-angle approximation tan α n ≈ α n . Now, we are ready to evaluate the limit (2.10) as follows:

$$\ u a t e & \text{ the limit (2.10) as follows: } \\ \lim _ { N \to \infty } F ( N ) & = \lim _ { N \to \infty } \sum _ { n = 1 } ^ { N } \tan \alpha _ { n } \\ & = \lim _ { N \to \infty } \sum _ { n = 1 } ^ { N } \alpha _ { n } \\ & = \lim _ { N \to \infty } \frac { \pi } { 4 } \\ & = \frac { \pi } { 4 }. & ( 2. 1 6 ) \\ h \text{ this result, we complete the proof.}$$

With this result, we complete the proof.

## 3 Final Remarks

Substituting the final result of the proof, eq. (2.16), into (2.9) provides an expression for the total descent time of the pendulum. From this, we can easily calculate the period of one complete oscillation as T = 4 t AB and finally arrive to the expression in (1.1). This concludes our solution to the problem.

Some readers might argue that this method is not entirely calculus-free because we use the concept of

limits and some of their properties in (2.16). They are correct. However, we believe that in this context, introducing the concept of limits does not complicate the problem but rather represents the geometry and physics in a very intuitive way. Therefore, even students with limited knowledge of calculus can grasp what is happening at each step.

In conclusion, our geometric approach to deriving the period of a pendulum not only offers an alternative to the traditional calculus-based method but also improves our understanding of the underlying physics. We believe this approach can enhance the teaching and learning of fundamental physics concepts and encourage educators to integrate these geometric techniques into their curriculum.

The geometric approach to solving the problem is credited to physicist R. Sánchez-Martínez, while the elegant proof of (2.10) was provided by mathematician E. Heredia-Muñoz.

## A 17th Century Derivation

The modern formula for the period of a simple pendulum, as we know it today, eq. (1.1), dates back to the mid-17th century. In this brief appendix, we would like to replicate the train of thought of the physicists of that time who first deduced it. Our aim is to illustrate the original ideas using modern notation for greater clarity in the explanations while attempting to maintain the original terms whenever possible.

When discussing the systematic study of a pendulum's motion, we often associate the first scientific observations with Galileo Galilei (1564-1642). According to Ariotti[4], 'although Galileo may not be the discoverer of all the major properties of the pendular process, there is no doubt that he is the first to recognize fully its importance' and his works served as a foundation for his succesors to build new and more precise theories.

Already in 1632, in his Dialogo [5], Galileo was aware that the period of a simple pendulum was proportional to its length, but it is six yearls later, in the Discorsi [6], that he justifies the more correct version of the law-that the period depends on the square root of the length-by means of a series of experiments. However, his observations did not let him to specify a more precise expresion. On the other hand, while Galileo was correct in stating that the period of a pendulum is independent of the amount of mass attached to it, he mistakenly believed that it was also independent of the initial amplitude-a fact that holds true only for sufficiently small angles. Another closely related error was his assumption that pendulums starting from various points along the same arc of a circle require the same amount of time to reach the base of that circle. Physics had to wait around twenty years for Huygens to find the curve that truly satisfied this condition, and in passing, to state for the first time the appropriate expression for the period.

Christiaan Huygens (1629-1695) was a Dutch scholar, physicist, and mathematician. Being aware of the isochronism issue of the oscillations in a circular pendulum, he dedicated himself to find the real curve with such characteristics and apply it in constructing a better pendulum clock. His major results were reported in the work, Horologium Oscillatorium [7] (The Pendulum Clock), published in 1673. In that book, he finally proved for the first time that the path along which a pendulum exhibits isochronous motion was not a circular arc, but a cycloid. He went even beyond in showing that for a pendulum to follow a cycloidal path, its point of suspension must be attached between two inverted similar cycloids (see Figure A.1 below).

Figure A.1: Geometrical elements of a cycloidal pendulum. This diagram is a reproduction of fig. 63 from du Châtelet's Institutions [8], with adapted notation to aid our explanations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-efa1a8b1-6220-4661-87ce-3f312a67aaba.jpg)

To understand Huygens' original statement, let's consider the cycloid ABC shown in Fig. A.1 with center at o and generated by a circle with radius R , and two similar cycloidal segments with vertices at A and C , coinciding at O . He stated that a pendulum supended from O and with length l = 4 R will describe the cycloid ABC as it swings, due to the partial rolling of its thread along the cycloidal arcs &gt; OA and &gt; OC .

All these results are recorded in the second part of the Horologium : De Descensu Gravium &amp; motu eorum in Cyclode (On the Descent of Heavy Bodies and

Their Motion in a Cycloid). And it is precisely at the end of this part that he addresses the problem of the period in Proposition XXV, which reads as follows:

In a cycloid whose axis is vertical with its vertex facing downward, the descent times in which a mobile object, released from any point on it, reaches the lowest point of the vertex, are equal to each other; and they have the same ratio to the time of perpendicular fall along the entire axis of the cycloid as the semicircumference of a circle has to its diameter.

The first part of the proposition tells us that the falling time from any point P to B is equal to that of falling from any other point Q , i.e., it is independent of the amplitude. In the second part, Huygens provides a value for the descent time through the cycloidal arc &gt; PB in terms of the free 'perpendicular' fall through the vertical axis oB , stating that the ratio t &gt; PB : t oB is the same 'as the semicircumference of a circle has to its diameter', which is the 17th-century way to express simply π/ 2 . 3 To get an expression for the period of one oscillation requires only for us to know basic kinematic formulas. Using the basic free fall formula | oB | = 2 R = 1 2 gt 2 oB , we can determine the time it takes for a particle to fall from o to B , leading to the familiar expression

$$\frac { t } { t _ { \overline { P B } } } = \frac { \pi } { 2 } \quad \Rightarrow \quad T \coloneqq 4 t _ { \widehat { P B } } = 2 \pi \sqrt { \frac { 4 R } { g } }. \quad ( \text{A.} 1 ) \quad \text{B}$$

as Descartes, Newton, and Leibniz, accompanied by beautiful, clear, and useful illustrations. Our Figure A.1 is indeed a reproduction of Figure 63 of the Institutions with notation adapted to our needs. Using our introduced notation, du Châtelet's paragraph § 474 reads

It follows from the equality of the small arc &gt; PB of the circle A BC ′ ′ and this portion of the cycloid ABC that the time during which a body makes an oscillation in a very small arc of a circle is to the time of the perpendicular fall over half the length of the pendulum, as the circumference of the circle is to its diameter, because the time of an oscillation in a cycloid follows such a proportion.

That is, when a circular pendulum describes very small oscillations, their period approximates to Huygens' formula, eq. (A.1), because the circular arc can be approximated by a cycloidal arc under these conditions. We find fascinating that the circular arc, typically considered a simple geometric form employed to study more complex curves like the cycloid, behaves differently in this context. Here, the motion on a 'more complex' cycloidal arc was used to approximate the motion on a simpler circular arc.

## B Modern Derivation

This is the exact expression for the period of a cycloidal, isochrone pendulum. But what about the simple circular pendulum? It seems to be clear to the scientists at the time, that this expression was also valid for very small oscillations, which would be approximated by a cycloid. In the fifth part of the Horologium , Huygens himself presents a theorem (Theo. IX) stating that the time of oscillation for a conical pendulum exhibits the same behavior as that of a corresponding cycloidal one, as long as it describes 'the smallest revolutions' ( circuitus minimos ). We can find little information about this approximation in future texts. In our view, one of the best explanations of the simple pendulum's period is provided by Émilie du Châtelet[8]. In her textbook Institutions de Physique (1740), we find excellent explanations of the most recent ideas of great scientists such

For completeness, we will outline the modern derivation of the pendulum's period to highlight concepts that might be obscure to students in their early stages of learning calculus or that lack clear physical meaning.

As we said in the introduction of this work, the evolution of the angle with respect to time, ϕ t ( ) , completely describes the motion of this system. The differential equation governing the dynamics of the system is

$$\ddot { \phi } + \frac { g } { l } \sin \phi = 0. \text{ \quad \ \ } ( B. 1 )$$

It can be determined using various methods: analyzing the forces acting on the mass and applying Newton's laws, formulating the Lagrangian of the system and solving the Euler-Lagrange equations, or employing the Hamiltonian function and reducing the Hamilton equations. Indeed, deriving the equations of mo-

3 The first mathematician to use the symbol π to denote the ratio of a circle's circumference to its diameter was William Jones in 1706. He introduced this notation in his work Synopsis Palmariorum Matheseos , although it was later popularized by Leonhard Euler by using it systematically in his influential works, such as Mechanica (1736) and Introductio in analysin infinitorum (1748).

tion for a simple pendulum is typically one of the initial examples used to illustrate the practical application of a newly introduced classical mechanics formalism (see for example [9] or [10]). For this illustration, we will employ the Newtonian approach.

Figure B.1: Free-body diagram of the forces acting on the pendulum mass.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00201v1-figure-d3820f86-b434-474e-8d30-df481f85fb47.jpg)

As shown in the free-body diagram, Fig. B.1, the only forces acting on the mass are its weight mg and the tension F l generated by the cord that constrains the pendulum to move in a circular path. Applying Newton's second law, ∑ F = m a , in the tangential direction, ˆ e ϕ , yields

$$- m g \sin \phi = m l \ddot { \phi }, \text{ \quad \ \ } ( B. 2 )$$

where we used the relation a ϕ = lϕ ¨ for circular motion. Equation (B.2) can be easily rearranged into (B.1). This differential equation has not analytical solution for arbitrary initial conditions. However, if the amplitude of angular displacement is small enough, so the small angle approximation, sin ϕ ≈ ϕ , holds true, it can be approximated as

$$\ddot { \phi } + \frac { g } { l } \phi = 0. \text{ \quad \ \ } \text{(B.3)}$$

This new differential equation has the form of a simple harmonic oscillator . It is a second-order homogeneous linear equation, whose solution method is well known. The conventional approach involves proposing a solution in the form of ϕ t ( ) = ae bt and then substituting it into the differential equation. This substitution leads to an algebraic equation, allowing the determination of the values for a and b in relation to the system parameters and initial conditions. Although this solution method would serve as an instructive excercise in a basic course on differential equations, for a student who is just learning the basics of differential calculus, this procedure may seem abstract and lacking much physical interest.

Using ϕ (0) = ∆ ϕ and ˙ (0) ϕ = 0 , we get the periodic solution

$$\phi ( t ) = \Delta \phi \cos ( \omega t ), \text{ \quad \ \ } ( \mathbb { B }. 4 )$$

where ω := √ g/l is the natural frequency of the motion. The period T of the motion is the time taken for the pendulum to return to the exact same point. Since the motion starts at t = 0 , then we seek for a t = T such that ϕ T ( ) = ϕ (0) . Using the periodic properties of the cosine function, we can assure that this happens when ωT = 2 π . From this condition, we obtain the famous formula of the period of oscillation (1.1).

## References

- [1] G. L. Baker and J. A. Blackburn, The pendulum: a case study in physics (Oxford University Press, 2008).
- [2] M. R. Matthews, C. Gauld, and A. Stinner, 'The Pendulum: Its Place in Science, Culture and Pedagogy', Science &amp; Education 13 , 261-277 (2004).
- [3] D. G. Zill, A First Course in Differential Equations with Modeling Applications , 9th ed. (Cengage Learning, 2009).
- [4] P. E. Ariotti, 'Aspects of the conception and development of the pendulum in the 17th century', Archive for History of Exact Sciences 8 , 329-410 (1972).
- [5] G. Galilei, Dialogo sopra i due massimi sistemi del mondo (Gio. Batista Landini, 1632).
- [6] G. Galilei, Discorsi e dimostrazioni matematiche intorno a due nuove scienze attenenti alla meccanica e i movimenti locali (Elsevier, 1638).
- [7] C. Huygens, Horologium Oscillatorium Sive de Motu Pendulorum ad Horologia Aptato Demonstrationes Geometricae (F. Muguet, 1673).
- [8] Émilie du Châtelet, Institutions de Physique (Prault fils, 1740).
- [9] J. R. Taylor, Classical mechanics (University science books, 2005).
- [10] G. Goldstein, C. P. Poole, and J. L. Safko, Classical mechanics (Pearson, 2011).
Rodrigo Sánchez-Martínez Esteban Heredia-Muñoz
A new Visual approach to pendulum period determination
The period of oscillation of a simple pendulum ($T = 2\pi\sqrt{l/g}$) is a
familiar formula to the average first-year physics student. However, deriving
this expression from first principles involves solving a non-linear
differential equation using the small angle approximation, which may appear
obscure to a student in the early stages of learning calculus. Therefore, we
propose an alternative approach to the derivation of this formula employing
geometry, algebra, and physical intuition. Our method follows the ideas of the
fathers of calculus, replacing the circular path of the pendulum with a
successive collection of infinitesimal inclined planes and summing the travel
times of each plane in the limit as the number of planes becomes very large.
Notably, the evaluation of this limit relies solely on geometry, making it
accessible to any student, even those not yet familiar with calculus
techniques.