Revisiting Euler’s Orbital Calculations for the Comet of 1742

Show more

1. Introduction

Figure 1 shows the orbit of Halley’s Comet and its predicted location in 2024 relative to the orbits of the planets. Blue is above the plane of the ecliptic and green is below. Almost the entire Halley orbit is below the plane of the ecliptic. Further, Halley revolves around its orbit in retrograde motion (the opposite sense from planet revolution). In the preceding view, the planets revolve counter-clockwise and Halley revolves clockwise. Figure 2 shows the same thing, but from a top view.

Figure 3 shows the solar system view from the plane of the ecliptic. This is an imaginary plane used as the reference plane for the solar system. Most planets in our solar system orbit the Sun on or near this plane. The Earth’s rotational axis is inclined at almost 23.5 degrees with respect to the perpendicular of the ecliptic plane. This inclination is what produces the seasons on the Earth.

Figure 4 shows the ecliptic coordinate system with the reference direction in the vernal equinox ( $\Upsilon $ ). The coordinates are celestial (or ecliptic) longitude (λ) and celestial (or ecliptic) latitude (β).

Figure 1. Halley’s comet in 2024. Source: http://www.pas.rochester.edu/~blackman/ast104/halley.html.

Figure 2. Halley’s comet in 2024—top view. Source: http://www.pas.rochester.edu/~blackman/ast104/halley.html.

Figure 3. Solar system view from the Ecliptic plane. Source: http://www.pas.rochester.edu/~blackman/ast104/halley.html.

Figure 4. The ecliptic coordinate system.

The celestial latitude (β) of a star, etc., is its angular distance (from 0˚ to 90˚) north (counted positive) or south (counted negative) of the ecliptic; it is measured along the great circle through the body and the poles of the ecliptic. The celestial longitude (λ) of a body is its angular distance (from 0˚ to 360˚) from the vernal equinox, measured eastward along the ecliptic to the intersection of the body’s circle of longitude; it is measured in the same direction as the Sun’s apparent annual motion. Although observations are taken from the Earth’s surface the coordinates should strictly be geocentric and, tabulated as such, are universally applicable. A slight correction is therefore applied to convert surface (topocentric) observations to geocentric values.

Figure 5 shows the trajectory of Comet Y4 ATLAS. Discovered by the ATLAS (Asteroid Terrestrial-impact Last Alert System) survey on the evening of December 28th, 2019, Comet Y4 ATLAS was the last comet discovery of 2019. Though it was only at magnitude +19.6 at the time of discovery^{1}, the orbit of Y4 ATLAS is intriguing, bringing it in just 0.262 AU (39.2 million kilometers or 24.4 million miles) from the Sun in late May 2020. That’s interior to Mercury at perihelion, at 0.307 AU from the Sun. This comet also grabbed observer’s attention in January 2020 when it jumped 5 magnitudes to +12, or 100-fold in brightness.

Figure 6 shows a comet in a parabolic orbit and the associated orbital elements. Together with *T*, the time that the comet is at the perihelion, the longitude of the ascending node (Ω), the inclination of the orbit (*i*), the argument of the perihelion (ω), and the perihelion distance (*a*), make five unknowns that fully characterize the parabolic orbit of a comet. As regards to *i*, if one stands at the ascending

Figure 5. The orbit of Comet C/2019 Y4 ATLAS, shown near perihelion, with indications of the ecliptic geocentric longitude (λ) and the ecliptic geocentric latitude (β). Credit: NASA/JPL. Adapted from: https://www.universetoday.com/145036/comet-y4-atlas-in-outburst-first-good-comet-for-2020/.

Figure 6. Comet in a parabolic orbit and associated orbital elements.

node and looking towards the Sun, it is reckoned upward from the plane of the ecliptic, towards the left, in the same way as angles are ordinarily reckoned. It may be anything from 0˚ to 180˚. If it is between 0˚ and 90˚, it shows that the comet is moving around the Sun in the same direction as the Earth and other planets, or that its motion is direct; if between 90˚ and 180˚, the comet moves in the opposite direction, or retrograde.

Now let *r* be the distance of a comet from the Sun, and *v* its velocity. It can be shown the if
$v>h\sqrt{\frac{2}{r}}$ , the comet is moving in a hyperbola; if
$v=h\sqrt{\frac{2}{r}}$ , its path is a parabola; and if
$v<h\sqrt{\frac{2}{r}}$ , it travels in an ellipse. Here *h*is Gauss’ gravitational constant, which is equal to 0.01720209895.

Therefore, it appears that of an infinite number of velocities a comet might have, only one would cause it to move in a parabola, and so, the chances are infinity to one that no comet moves in that curve. It is, moreover, known from observation that of all comets whose orbits have been computed, not more than three of four appeared to be hyperbolas, and those were not certainly so, there is strong reason to believe that all comets are moving in ellipses. However, in the cases of all but a few their orbits are so nearly parabolas that their performances can be very well predicted by assuming a parabola as the true curve, and in the cases of a newly discovered comet, this always done (Truman, 1912). This is one of the reasons why we are only concerned here with Euler’s parabolic case.

2. Brief History of Comets Orbital Calculations

Orbit determination has a long history, beginning with the prehistoric discovery of the planets and subsequent attempts to predict their motions. Johannes Kepler used Tycho Brahe’s careful observations of Mars to deduce the elliptical shape of its orbit and its orientation in space, deriving his three laws of planetary motion in the process. The mathematical methods for orbit determination originated with the publication in 1687 of the first edition of Newton’s Principia, which gave a geometrical method for finding the orbit of a body following a parabolic path from three observations. Newton’s method was considered difficult to apply in practice, but it was used, with minor adaptations, by Edmund Halley to establish the orbits of various comets, including that which bears his name. The first analytical method based on three observations, that did not make use of geometrical considerations was presented by Euler in *Theoria* *motuum* *planetarum* *et* *cometarum* (Euler, 1744), and applied to the investigations of the comets of 1680 and 1744. This work was in turn generalized to elliptical and hyperbolic orbits by Lambert in 1761-1777.

In the years that followed, other leading mathematicians and astronomers devoted themselves to the problem. In 1774, Boscovich developed a method that also made use of three observations, published under the title *De* *orbitus* *cometarum* *determinandis* (Boscovich, 1774). In 1780, Laplace (1784) gave an entirely new approach to the problem. In contrast to others previously methods, requiring three observations relatively close to each other, Laplace’s method could make use of observations separated by about 30 to 40 degrees of each other. The precision of the method improves if, indeed, it uses a more extensive set of observations, but, sometimes it is unsatisfactory when only three observations are available. Orbit determination by three observations was also employed in 1782 by the Portuguese astronomer José Monteiro da Rocha, in a method which is essentially the same as the method proposed by Olbers and published under von Zach’s sponsorship two years before, in 1797 (Figueiredo & Fernandes, 2020).

Another milestone in orbit determination was Carl Friedrich Gauss’ assistance in the “recovery” of the dwarf planet Ceres in 1801. Gauss’s method was able to use just three observations (in the form of celestial coordinates) to find the six orbital elements that completely describe an orbit.

The problem of comet’s orbit determination was a recurrent one along the years that followed, such that the theory of orbit determination has been developed to the point where today it is applied in GPS receivers as well as the tracking and cataloguing of newly observed minor planets.

3. Euler's Geometrical Calculations of the Parabolic Orbit for the Comet of 1742 (Euler, 1743)

The geometrical method, based on three observations, was used by Euler for the investigation of the comet of 1742 (Euler, 1743). As we shall see, crucial to the developments, is the determination of the comet’s distance to the Earth at the moment of the second observation, which Euler obtained by trial and error.

The mathematical formulation for the calculation of the orbital parameters is essentially in plane geometry, which, for a phenomenon that occurs in space, is rather difficult to follow from sketches of the space geometry projected in the plane of the paper. As pointed out earlier, the geometrical method is difficult to apply in practice. Moreover, Euler’s calculations are rather cumbersome by transforming the original operations into logarithmic operations, in which all minute details are presented. By considering different conic sections for the orbit, such as parabolas, ellipses and hyperboles, the calculations turned out very repetitive and tedious. This is another reason why only the parabolic orbit case will be considered here for the comet of 1742.

Only the paragraphs directly concerned with the parabolic case will be reviewed, and even so, they were abridged to what is real necessary to logically follow Euler’s calculations. The paragraphs sub-headings in what follows are those that appear in the original publication.

Euler begins by praising Newton’s theory on the motion of planets, and by commenting comets’ observations and the theory of comets. Then follows a discussion on the modeling of comets’ orbits from observations, by means of conics with the Sun at the focus, and by applying the principle of equal areas in equal times. He then asserts the objective of the investigation of developing a procedure to find the classical parameters that characterizes the orbit of a comet from three observations, with a final goal of checking his method by comparison of the calculated geocentric longitude and latitude of the Comet with a fourth observation.

In §. 6, 7, and 8 Euler develops the expression for the calculation of the true anomaly of the comet in a parabolic orbit^{2} as

$t+\frac{1}{3}{t}^{3}=\frac{n}{N}$

where
$t=\mathrm{tan}\frac{1}{2}\nu $ *,* ν is the anomaly, *n* is the number of days that it takes for the comet to go from the perihelion *P* to *M*,
$N=\frac{a\sqrt{a}}{384651.5}$ , and *a* is the perihelion-Sun distance *PS*. In Figure 7,
$\nu =P\stackrel{^}{S}M$ .

Euler then provides a table at the end of the manuscript that allows the direct determination of ν from tabulated values of the polynomial $t+\frac{1}{3}{t}^{3}$ , obtained from the ratio $\frac{n}{N}$ .

The paragraphs that follow, relate to the calculation of the parameters that fully characterize the parabolic orbit of the Comet of 1742.

• . XXIII.

Besides the observations on Mar. 3^{d}, 15^{h}, 10', let be chosen two other observations held before and after the conjunction^{3}, which are not too close and not too far apart from each other. If they are chosen very close to one another, the aberration error^{4} may be sensible: if however, they were chosen rather far apart, they bring forth the hindrance of curvature. From these considerations, the three sequence of the Comet’s position shown in Table 1 are seen most suitable for this purpose.

The time elapsed out between the first and the second observation is 3^{d}, 8^{h}, 21' = 3.34791 days, which is set = *m*. Indeed, the time elapsed out between the second and the third observation is 2^{d}, 23^{h}, 5' = 2.96180 days, which is set = *n*, and hence, the time elapsed out between the first and the last observation = *m* + *n* = 6.30971 days.

Figure 7. Geometrical elements used by Euler for the calculation of the comet’s anomaly $\nu =P\stackrel{^}{S}M$ . Adapted from Figure 1 (Euler, 1743).

Table 1. Observations made by the astronomer Delisle^{5} for the comet of 1742^{6}.

• . XXIV.

Figure 8 now represents the plane of the ecliptic, in which the Sun occupies the position *S*, and the three locations of the Earth are indicated by *A*, *B* and *C* according to these three observations, and once guided by the straight lines *AS*, *BS*, *CS*, these segments will have the following dimensions:

$AS=99482.70$

$BS=99576.55$

$CS=99660.46$

and

$\text{angle}ASB=\text{3}\u02da,\text{19}\text{'},\text{53}\text{'}\text{'}$

$\text{angle}BSC=\text{2}\u02da,\text{56}\text{'},\text{22}\text{'}\text{'}$

$\text{angle}ASC=\text{6}\u02da,\text{16}\text{'},\text{15}\text{'}\text{'}$ .^{11}

Now that the segment *BS* shows the distance of the Comet to the Sun in the intermediate observation, may be drawn the lines *A*α, and *C*γ which designate the Comet’s distances to the Earth in the first and last observation, whose positions are to be determined such that

$\text{angle}\text{\hspace{0.17em}}SA\alpha =\text{42}\u02da,\text{1}\text{'}$ , $\text{angle}\text{\hspace{0.17em}}AmB=\text{45}\u02da,\text{20}\text{'},\text{53}\text{'}\text{'}$

$\text{angle}\text{\hspace{0.17em}}SC\gamma =\text{45}\u02da,2\text{4}\text{'}$ , $\text{angle}\text{\hspace{0.17em}}CnB=\text{48}\u02da,\text{20}\text{'},\text{22}\text{'}\text{'}$

seconds were disregarded in the location of the Comet, because these observations are not exact beyond minutes. Then, consequently, the segments *A*α & *C*γ

Figure 8. Triangles involved in the calculations considered in §. XXIV. Adapted from Figure 3 (Euler, 1743).

intercepts the segment *BS* in the points *m* & *n*, which should be firstly defined from the following formulas, which result from the resolution of the triangles^{12}

$Sm=\frac{AS\cdot \mathrm{sin}SA\propto}{\mathrm{sin}AmB}=93603.40;\text{\hspace{0.17em}}\text{\hspace{0.17em}}Am=\frac{AS\cdot \mathrm{sin}ASB}{\mathrm{sin}AmB}=8126.40$

$Sn=\frac{CS\cdot \mathrm{sin}SC\gamma}{\mathrm{sin}CnB}=94982.25;\text{\hspace{0.17em}}\text{\hspace{0.17em}}Cn=\frac{CS\cdot \mathrm{sin}BSC}{\mathrm{sin}CnB}=6840.66$

$Bm=BS-Sm=99576.55-93603.40=5973.15$

$Bn=BS-Sn=99576.55-94982.25=4594.30$

• . XXV.

If we now set *apc* as the portion of the Comet’s orbit described during the interval of the assumed observations, then *a* will be its position in the first observation, from which the perpendicular *a*α is dropped onto the ecliptic plane, defining the segment *A*α on this plane; and the same with *c*, which is the position of the Comet in the third observation, then, defining the segment *C*γ on the ecliptic; and also *p* is the location in the intermediate observation, with the perpendicular *p*π falling onto the point π on the line *BS*. On account that the areas described around the Sun are proportional to the times, then the ratio of area *aSp* to the area *pSc* is equal to the ratio of *m* to *n*: and if this ratio holds for the Comet’s orbit projected onto the plane of the ecliptic απγ, then, also the ratio of the area α*S*π to the area γ*S*π is as *m* to *n*. Moreover, it is also clear that the ratio of the sectors α*S*π and γ*S*π come near to approach the ratio of the triangles α*S*β & β*S*γ, such that the error is certainly imperceptible: since the ratio between these areas is almost the same, then because the angles α*S*β and β*S*γ, as the following will show, are rather small to such an extent for the curvature not to upset in any way this convenience. Then, since the ratio of the triangles α*S*β and β*S*γ will be as *m* to *n*, also
$\alpha \beta :\beta \gamma =m:n$ , an thus, the ratio of the segments of the cords αγ will be known by means of factors of *SB*: whence, if the point β in the line *BS* becomes known, through which the cord αγ transits, at the same time, the position of the cord or angle *B*βγ can be defined. So, let us set the angle
$B\beta \gamma =\phi $ , then, in the triangle α*m*β, the angle
$\alpha m\beta =AmB$ &
$m\alpha \beta =\phi -AmB$ ; in the triangle *βnγ*, indeed the angle
$\beta n\gamma =CnB$ , &
$\beta n\gamma ={180}^{\circ}-\phi -BnC$ . Then, it follows that^{13}

$\alpha \beta =\frac{\beta m\cdot \mathrm{sin}AmB}{\mathrm{sin}\left(\phi -AmB\right)}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\beta \gamma =\frac{\beta n\cdot \mathrm{sin}BnC}{\mathrm{sin}\left(\phi +BnC\right)}$

since we have that $\alpha \beta :\beta \gamma =m:n$ , and because $\mathrm{sin}\left(\phi -AmB\right)=\mathrm{sin}\phi \mathrm{cos}AmB-\mathrm{cos}\phi \mathrm{sin}AmB$ & $\mathrm{sin}\left(\phi +BnC\right)=\mathrm{sin}\phi \mathrm{cos}BnC+\mathrm{cos}\phi \mathrm{sin}BnC$ , by dividing both denominators by $\mathrm{cos}\phi $ , gives

$\frac{n\cdot \beta m}{\mathrm{tan}\phi \cdot \mathrm{cot}AmB-1}=\frac{m\cdot \beta n}{\mathrm{tan}\phi \cdot \mathrm{cot}BnC+1}$

or

$\mathrm{tan}\phi =\frac{\frac{m}{n}\cdot \frac{\beta n}{\beta m}+1}{\frac{m}{n}\cdot \frac{\beta n}{\beta m}\mathrm{cot}AmB-\mathrm{cot}BnC}$

thus, from the knowledge of the point β, the position of the cord αγ becomes known, and further, from the given latitudes, the points *a* and *c* in the Comet’s orbit.

• . XXVI.

Therefore, this matter has been reduced to the question of defining the point β on the line *B*β, such that the solution of the problem lies just on one unknown. Yet, it is not fit, as is usually done in other circumstances, to employ some character to this unknown, and to determine it by an equation, because of the too complicated expressions to which we would have arrive. Then, it will be employed some tentative values for *B*β, indeed ruling out the false ones, such that the most problematic undefined letters are avoided in the trigonometric calculation (Figure 9).

Let then be assumed some arbitrary values for the interval *B*β, and for whatever convenience the angle
$\phi =B\beta \gamma $ will be defined, and two of the Comet’s position *a *and *c* will be held, as well as their distances to the Sun *AS*, *CS*, and they will be able to assign both the angle *aSc* and the cord *ac*. Then, from the developments of the preceding sections, having been established that the Comet’s orbit is a parabola or an ellipse, no matter how oblong, let be theoretically

Figure 9. Triangles involved in the calculations considered in §. XXVI. Adapted from Figure 3 (Euler, 1743).

found the time that the Comet should had taken to travel the distance *ac*, and how this time compares with the observed time of 6.30971 days, and this will be the evidence that the true value for *B*β has been assumed, if however, it disagrees, the approximate true value is found from the comparison of two additional hypothesis for *B*β. Let then be set again two values for *B*β itself, not too much different from the one that was found, such that neither of them are not too much in disagreement with the reality, for the approximate true value of *B*β to be found.

• . XXVII.

Therefore, two values 8000 and 9000 will be attributed to *B*β itself, since I already know from other rather crude attempts that the true value lies between these limits; and I will first determine the angle
$\phi =B\beta \gamma $ .

From §. XXV. We know that

$\phi ={\mathrm{tan}}^{-1}\left[\frac{\frac{m}{n}\cdot \frac{\beta n}{\beta m}+1}{\frac{m}{n}\cdot \frac{\beta n}{\beta m}\mathrm{cot}AmB-\mathrm{cot}BnC}\right]$

where

$\frac{m}{n}=\frac{3.34791}{2.96180}=1.130$

for

$B\beta =8000$ $\{\begin{array}{l}\beta m=B\beta -Bm=8000-5973.15=2026.85\\ \beta n=B\beta -Bn=8000-4594.30=3405.70\end{array}\Rightarrow \phi =71\u02da12\text{'}22\text{'}\text{'}$

for

$B\beta =9000$ $\{\begin{array}{l}\beta n=B\beta -Bn=9000-4594.30=3026.85\\ \beta m=B\beta -Bm=9000-5973.15=4405.70\end{array}\Rightarrow \phi =74\u02da27\text{'}29\text{'}\text{'}$

Finally

for

$B\beta =8000$ $\{\begin{array}{l}\alpha m\beta =AmB=45\u02da20\text{'}53\text{'}\text{'}\\ A\propto \beta =m\propto \beta =\phi -AmB=25\u02da51\text{'}29\text{'}\text{'}\\ \beta n\gamma =CnB=48\u02da20\text{'}22\text{'}\text{'}\\ C\gamma \beta ={180}^{\circ}-\phi -\beta n\gamma =60\u02da27\text{'}16\text{'}\text{'}\end{array}$

for

$B\beta =9000$ $\{\begin{array}{l}\alpha m\beta =AmB={45}^{\circ}20\text{'}53\text{'}\text{'}\\ A\propto \beta =m\propto \beta =\phi -AmB={29}^{\circ}6\text{'}36\text{'}\text{'}\\ \beta n\gamma =CnB={48}^{\circ}20\text{'}22\text{'}\text{'}\\ C\gamma \beta ={180}^{\circ}-\phi -\beta n\gamma ={58}^{\circ}12\text{'}9\text{'}\text{'}\end{array}$

• . XXVIII.

Now, in the triangles α*m*β, β*n*γ all angles and the sides β*m* and β*n* are given, and the remaining sides will be defined by the following formulas

$\alpha \beta =\frac{\beta m\cdot \mathrm{sin}\alpha m\beta}{\mathrm{sin}A\alpha \beta}=\{\begin{array}{l}3305.98\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 4426.15\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

$\alpha m=\frac{\beta n\cdot \mathrm{sin}\phi}{\mathrm{sin}A\alpha \beta}=\{\begin{array}{l}4399.44\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 5994.33\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

$\beta \gamma =\frac{\beta n\cdot \mathrm{sin}\beta n\gamma}{\mathrm{sin}C\gamma \beta}=\{\begin{array}{l}2924.70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 3827.71\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

$\gamma n=\frac{\beta n\cdot \mathrm{sin}\phi}{\mathrm{sin}C\gamma \beta}=\{\begin{array}{l}3706.04\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 4994.15\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

Having been found α*m* and γ*n*, then

$\begin{array}{l}A\alpha =Am+\alpha m=\{\begin{array}{l}12525.84\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 14120.73\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}\\ C\gamma =Cn+n\gamma =\{\begin{array}{l}10546.70\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 11834.81\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}\end{array}$

• . XXIX.

Now, from the latitudes that has been observed, the segments *a*α and *c*γ will be obtained from
$a\alpha =A\propto \cdot \mathrm{tan}lat.$ and
$c\gamma =C\gamma \cdot \mathrm{tan}lat.$ (Figure 10). Indeed, the latitude in the first observation = 70˚, 30', and the latitude in the third observation = 77˚, 37'. Then,

$a\alpha =A\alpha \cdot \mathrm{tan}70\u02da,30\text{'}=\{\begin{array}{l}35371.88\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 39875.71\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

and

$c\gamma =C\gamma \cdot \mathrm{tan}77\u02da,37\text{'}=\{\begin{array}{l}48035.68\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 53902.47\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

Additionally, having been found *a*α and *c*γ, may be guided the segment *aK* parallel to αγ, then

Figure 10. Triangle involved in the calculations considered in §. XXIX. Adapted from Figure 3 (Euler, 1743).

$aK=\alpha \beta +\beta \gamma =\{\begin{array}{l}6230.68\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 8298.86\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

and

$cK=c\gamma -a\alpha =\{\begin{array}{l}12663.80\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 14026.76\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

then, since
$\mathrm{tan}cak=\frac{cK}{ak}$ , the angle *cak* is found; and out of this angle, we have additionally that the cord

$ac=\frac{cK}{\mathrm{sin}caK}=\{\begin{array}{l}14113.58\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 16297.89\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

• . XXX.

May be guided the segments *S*α and *S*γ, and in the disclosed cord αγ may be guided the normal *SQ*: then, in the right triangle β*SQ* the side
$S\beta =SB-B\beta $ , which together with the given angle will result in
$SQ=S\beta \mathrm{sin}\phi $ ,
$\beta Q=S\beta \mathrm{cos}\phi $ (Figure 11). Hence,
$Q\alpha =Q\beta -\alpha \beta $ and
$Q\gamma =Q\beta +\beta \gamma $ : and also
$\frac{QS}{Q\alpha}=\mathrm{tan}S\alpha Q$ and
$\frac{QS}{Q\gamma}=\mathrm{tan}S\gamma Q$ , and once these angles are known, then,
$S\alpha =\frac{QS}{\mathrm{sin}S\alpha Q}$ and
$S\gamma =\frac{QS}{\mathrm{sin}S\gamma Q}$ .

• . XXXI.

Once *S*α and *S*γ are known, and since the perpendiculars *a*α and *c*γ are given, then
$\mathrm{tan}aS\alpha =\frac{a\alpha}{S\alpha}$ and
$\mathrm{tan}cS\gamma =\frac{c\gamma}{S\gamma}$ , producing the heliocentric longitudes

*aS*α and *cS*γ, and then, the Comet’s distances to the Sun are found out from the following formulas:

Figure 11. Triangles involved in the calculations considered in §. XXX. Adapted from Figure 3 (Euler, 1743).

$Sa=\frac{a\alpha}{\mathrm{sin}aS\alpha}=\{\begin{array}{l}97228.33\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 97974.17\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

$Sc=\frac{c\gamma}{\mathrm{sin}cS\gamma}=\{\begin{array}{l}104282.54\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =8000\\ 106360.48\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}B\beta =9000\end{array}$

• . XXXII.

Having been now found the two Comet’s distances to the Sun *Sa* and *Sc*, along with the cord *ac*, it will be possible to define the time required by the Comet to complete this distance under the hypothesis that the Comet’s orbit is either parabolic or elliptic: indeed, first, if the Comet’s orbit is established as a parabola, then set
$Sa=f,Sc=g$ and
$ac=k$ , the time through *ac* will be

$=\frac{{\left(\frac{f+g+k}{2}\right)}^{3/2}-{\left(\frac{f+g-k}{2}\right)}^{3/2}}{3m}$ days^{14}, knowing that
$m=384651.5$ , which is a constant that depends only on the Earth-Sun distance (Figure 12).

This expression gives 5.82218 days for
$B\beta =8000$ , and 6.70979 days for
$B\beta =9000$ . The actual time that the Comet took to travel to the arc *ac* was found to be equal to 6.30971 days. In §. XXXIII, Euler found by interpolation that the value for *B*β should be around 8514. Later, in §. XXXVIII, Euler refined the value for *B*β to 8471.61.

• . XLII.

In this paragraph, Euler calculates the Comet’s distances to the Earth in the 1^{st} and 3^{rd} observations from (Figure 13)

Figure 12. Geometrical elements used by Euler for the calculation of the time that it takes for the comet to travel the arc *ac*. Adapted from Figure 2 (Euler, 1743).

Figure 13. Triangles involved in the calculations considered in §. XLII. Adapted from Figure 3 (Euler, 1743).

$Aa=\frac{a\alpha}{\mathrm{sin}lat.}=39777.5$

and

$Cc=\frac{c\gamma}{\mathrm{sin}lat.}=52134.6$

and from these, a new value for $ac=15245,34$ .

• . XLIII.

When the cord *ac* is extended out, it meets the prolonged segment γα in *N*, and the point *N* will be simultaneously in the orbit of the Comet and in the ecliptic; and if the segment *SN* is drawn, then this segment will belong to the nodal line^{15}. Since the angle
$aN\alpha =caK$ , then
$\frac{a\alpha}{N\alpha}=\mathrm{tan}caK$ , and therefore,

$\alpha N=\frac{a\alpha}{\mathrm{tan}caK}$ . Moreover, once the normal *SQ* is drawn from *S* to αγ, the segments *S*α and *S*γ will be defined, and the angles *BS*α and *BS*γ are given, the heliocentric longitude of the Comet will be obtained in the 1^{st} and 3^{rd} observation. The results are as follows:

• $\alpha N=\frac{a\alpha}{\mathrm{tan}caK}$ .

• Since $\beta Q=S\beta \cdot \mathrm{cos}\phi $ , then, $\alpha Q=\beta Q-\alpha \beta $ and $\gamma Q=\beta Q+\beta \gamma $ .

• Since $SQ=S\beta \cdot \mathrm{sin}\phi $ , then, $S\stackrel{^}{\alpha}Q={\mathrm{tan}}^{-1}\frac{SQ}{\alpha Q}=75\u02da20\text{'}48\text{'}\text{'}$ , $S\stackrel{^}{\gamma}Q={\mathrm{tan}}^{-1}\frac{SQ}{\gamma Q}=70\u02da59\text{'}48\text{'}\text{'}$ .

• $\alpha \stackrel{^}{S}\beta =S\stackrel{^}{\alpha}Q-Q\stackrel{^}{\beta}S=S\stackrel{^}{\alpha}Q-\phi =2\u02da,22\text{'},12\text{''}$ .

• Heliocentric longitude of the Comet in the 1^{st} observation = longitude of *B* +
$\alpha \stackrel{^}{S}\beta $ = 5^{S}, 24˚, 24', 26'' + 2˚, 22', 12'' = 5^{S}, 26˚, 46', 38''.

• $\beta \stackrel{^}{S}\gamma =Q\stackrel{^}{\beta}S-\beta \stackrel{^}{S}\gamma =\phi -S\stackrel{^}{\gamma}Q=1\u02da,58\text{'},48\text{''}$ .

• Heliocentric longitude of the Comet in the 3^{rd} observation = longitude of *B *-
$\beta \stackrel{^}{S}\gamma $ = 5^{S}, 24˚, 24', 26'' - 1˚, 58', 48'' = 5^{S}, 22˚, 25', 38''.

• Heliocentric latitude of the Comet in the 1^{st} observation =
${\mathrm{tan}}^{-1}\frac{a\alpha}{S\alpha}$ = 22˚, 36', 29''.

• Sun-Comet distance in the 1^{st} observation
$Sa=\frac{a\alpha}{\mathrm{sin}aS\alpha}=97537.59$ .

• Heliocentric latitude of the Comet in the 3^{rd} observation =
${\mathrm{tan}}^{-1}\frac{c\gamma}{S\gamma}$ = 28˚, 55', 44''.

• Sun-Comet distance in the 3^{rd} observation
$Sc=\frac{c\gamma}{\mathrm{sin}cS\gamma}=105270.13$ .

• . XLIV.

Having been found these lines, then
$QN=\alpha Q-\alpha N$ . Hence,
$\mathrm{tan}QSN=\frac{QN}{QS}$ ; therefore,
$\beta \stackrel{^}{S}N=\beta \stackrel{^}{S}Q-Q\stackrel{^}{S}N={90}^{\circ}-\text{\phi}-Q\stackrel{^}{S}N={15}^{\circ},18\text{'},41\text{''}$ , thus the position of the nodal line becomes known, and also the heliocentric longitude of the ascending node *N* as

$N=\beta \stackrel{^}{S}N+\text{longitude}\text{\hspace{0.17em}}\text{of}\text{\hspace{0.17em}}B=15\u02da,18\text{'},41\text{''}+{5}^{\text{S}},24\u02da,24\text{'},26\text{''}={6}^{\text{S}},9\u02da,43\text{'},7\text{''}$ . Then, from *α*, the perpendicular α*R* to *SN* is drawn, giving
$\alpha R=S\alpha \mathrm{sin}\alpha SR$ , knowing that
$\alpha \stackrel{^}{S}R=Q\stackrel{^}{\alpha}S-Q\stackrel{^}{S}N$ and
$SR=S\alpha \mathrm{cos}\alpha SR$ ; and once *aR* is drawn, which is normal to, then it will be found out that
$\mathrm{tan}aR\alpha =\frac{a\alpha}{\alpha R}$ ; the angle
$a\stackrel{^}{R}\alpha ={\mathrm{tan}}^{-1}\frac{a\alpha}{\alpha R}=61\u02da,43\text{'},44\text{''}$ will be the inclination of the Comet’s orbit in relation to the ecliptic. Hence, indeed,
$aR=\frac{a\alpha}{\mathrm{sin}aR\alpha}$ and
$\mathrm{sin}aSN=\frac{aR}{Sa}$ , and the angle
$a\stackrel{^}{S}N={\mathrm{sin}}^{-1}\frac{aR}{Sa}=25\u02da,52\text{'},49\text{''}$ will be the heliocentric distance of the Comet to the node in the first observation. Finally, then, the angle
$a\stackrel{^}{S}c$ will be found from
$\mathrm{cos}aSc=\frac{S{a}^{2}+S{c}^{2}-a{c}^{2}}{2\cdot Sa\cdot Sc}$ , and hence, the Comet’s heliocentric distance to the node in the third observation
$\begin{array}{l}=a\stackrel{^}{S}N+a\stackrel{^}{S}c=a\stackrel{^}{S}N+2\cdot {\mathrm{cos}}^{-1}\sqrt{\left(\frac{Sa+Sc+ac}{2\cdot Sa}\right)\cdot \left(\frac{Sa+Sc-ac}{2\cdot Sc}\right)}\\ =25\u02da,52\text{'},49\text{''}+7\u02da,26\text{'},3\text{''}=33\u02da,18\text{'},52\text{''}\end{array}$ .

• . XLV.

Having been now found both Comet-Sun distances $Sa=f$ ; $Sc=g$ , together with the cord $ac=k$ or the angle $a\stackrel{^}{S}c=q$ , it will be possible to define the parabolic and the elliptical orbit of the Comet. So, be considered the angle $P\stackrel{^}{S}a=\nu $ , and if it is set $B=g\mathrm{cos}q-f$ and $C=\frac{g-f}{1-\alpha}$ , then,

$\mathrm{sin}\frac{1}{2}\nu =\frac{\sqrt{\left(k+B\right)\left(k+C\right)}-\sqrt{\left(k-B\right)\left(k-C\right)}}{2k}$ , where for the parabola $\alpha =0$ and for the ellipse $\alpha =\frac{1}{50}$ . Thus, having been found the angle ν, then

$PS=a=\frac{f}{2-\alpha}+\frac{f\left(1-\alpha \right)}{2-\alpha}\mathrm{cos}\nu $ .^{16} For the parabolic orbit, the numerical value found by Euler for the perihelion distance is
$PS=a=75210.10$ , and for the anomaly
$a\stackrel{^}{S}p=\nu =57\u02da,10\text{'},5\text{''}$ . Then, by subtracting, from this angle, the angle
$a\stackrel{^}{S}N=25\u02da,52\text{'},49\text{''}$ , the distance of the perihelion to the node was found as
$\omega =31\u02da,17\text{'},16\text{''}$ .

• . XLVI.

Having been found the orbit of the Comet, it remains to be defined the time in which the Comet transits by the perihelion, which will be found out, if we determine the time that it will take for the Comet to move from the perihelion *P* to the position *a*. Be set
$\mathrm{tan}a\stackrel{^}{S}P=\mathrm{tan}\frac{1}{2}\nu =t$ , then the time that it takes for the Comet to arrive at *a *from *P*, expressed in days
$=\frac{a\sqrt{a}}{384651.5}\left(t+\frac{1}{3}{t}^{3}\right)$ , and by putting
$\frac{a\sqrt{a}}{384651.5}=N$ , and
$T=+\frac{1}{3}{t}^{3}$ , gives the time for going through
$Pa=NT$ . These give
$T=0.598773$ ,
$N=53.623$ ,
$NT=32.10771={32}^{\text{d}},{2}^{\text{h}},35\text{'}$ . Yet, the Comet is seen at *a* in Feb. 28^{d}, 6^{h}, 49', therefore, the Comet had been at the perihelion on Jan. 27^{d}, 4^{h}, 14'.

• . XLVII.

In this paragraph, Euler gives the calculated parabolic orbital parameters to be applied to the calculations that follow in the next paragraphs.

• Distance from the perihelion to the Sun $a=75210.10$ .

• Semi latus rectum $b=150420.20$ .

• Date in January that the Comet transits by the perihelion: 27^{d}, 4^{h}, 14'.

• Distance of the perihelion to the node = 31˚, 17', 16''.

• Heliocentric longitude of the ascending node = 6^{S}, 9˚, 43', 7''.

• Inclination of the Comet’s orbit = 61˚, 43', 44''.

• . XLVIII.

Let us calculate for the parabolic orbit case, and for the time of the last observation on March 17^{d}, 7^{h}, 50', for which the observed longitude of the Comet was 2^{S}, 18˚, 52', with a northern latitude = 63˚, 13’. Indeed, for this time the Earth’s position^{17} was drawn = 6^{S}, 7˚, 57', 30'' and the Earth-Sun distance = 99970.05. For these conditions, and following the procedures of the previous paragraphs Euler obtained the following: true anomaly ν = 74˚, 56', 29''; Comet’s distance to the node = 43˚, 39’, 13’’.

• . XLIX.

Be considered in the spherical triangle ☊*c*γ shown in Figure 14, the ecliptic side ☊γ, and the orbit of the Comet ☊*c*, then the angle *c*☊γ will be the inclination in relation to the ecliptic.

Then, by assuming that *c*☊ is the distance from the Comet to the node, hence, *c*, indeed, will be the position of the Comet in its orbit; whence if the normal *c*γ to ☊*γ* is drawn, then the arc *c*γ will be the heliocentric longitude of the Comet, and if ☊*γ* is subtracted from the longitude of the node, the result will give the heliocentric longitude of the Comet. Indeed, from the nature of [right] spherical triangles, we have that sin*c*γ = sin*c*☊∙sin*c*☊γ, and tan☊γ = cos*c*☊γ∙tan*c*☊.^{18} From these relations, the following numerical results were found by Euler: Heliocentric Latitude of the Comet *c*γ = 43˚, 39', 13''; Heliocentric Longitude of the Comet 5^{S}, 15˚, 24', 3''.

• . L.

Then, in Figure 15 let the Sun be set at *S* and *T* be the position of the Earth, then, outside the plane of the ecliptic, emerges the Comet at *c*, whence in the ecliptic, falls the perpendicular *c*γ, and *S*γ can be drawn, being the angle
$c\stackrel{^}{S}\text{\gamma}$ the heliocentric latitude of the Comet. Then,
$c\gamma =Sc\mathrm{sin}cS\gamma $ , and
$S\gamma =Sc\mathrm{cos}cS\gamma $ . Now, from γ, be drawn the perpendicular γ*p* to *ST*, then, from the given angle
$\gamma \stackrel{^}{S}T$ , which is certainly the difference between the heliocentric longitudes of the Earth and the Comet
$\gamma p=S\gamma \mathrm{sin}\gamma ST$ and
$Sp=S\gamma \mathrm{cos}\gamma ST$ , with
$Tp=ST-Sp$ . From these, it is found out that
$\mathrm{tan}ST\gamma =\frac{p\gamma}{pT}$ , which is the angle that when added to the longitude of the Sun, gives the geocentric longitude of the Comet (Figure 16).

Besides, we have that $T\gamma =\frac{p\gamma}{\mathrm{sin}ST\gamma}$ ; and $\mathrm{tan}cT\gamma =\frac{c\gamma}{T\gamma}$ ; will be the geocentric latitude of the Comet, and, finally, $\frac{c\gamma}{\mathrm{sin}cT\gamma}=Tc$ is the distance from the Comet to the Earth.

Figure 14. Spherical triangle ☊*c*γ for the calculation of the heliocentric longitude of the Comet in the last observation. Adapted from Figure 4 (Euler, 1743).

Figure 15. Auxiliary figure provided by Euler to determine the geocentric latitude and longitude of the comet.^{19} Adapted from Figure 5 (Euler, 1743).

Figure 16. Geocentric longitude of the comet.

From these expressions, Euler calculates
$S\stackrel{^}{T}\gamma ={2}^{\text{S}},11\u02da,8\text{'},21\text{''}$ , which when added to the longitude of the Sun = 7˚, 57', 30''
$\left[{\lambda}_{Sun}={6}^{\text{S}},7\u02da,57\text{'},30\text{''}+{180}^{\circ}=7\u02da,57\text{'},30\text{''}\left(\mathrm{mod}360\u02da\right)\right]$ , gives 2^{S}, 19˚, 5', 51'' for the geocentric longitude of the Comet. For the geocentric latitude of the Comet, Euler obtains 62˚, 6', 10''.

In §. LI, Euler compares the observed geocentric longitude of the Comet of 2^{S}, 18˚, 52' with the calculated value of 2^{S}, 19˚, 5', 51'', saying that they do not too much disagree, whereas the observed geocentric latitude of 63˚, 13' much disagrees with the calculated value of 62˚, 6', 10'', and concludes that this disagreement can be attributed to error in the observation.

4. Conclusion

Euler’s calculations for investigating the parabolic orbit for the Comet of 1742 were revisited. Besides the true anomaly, Euler succeeded in finding the classical orbital parameters from three observations of the Comet taken from the Earth: longitude of the ascending node, the inclination of the orbit, the argument of the perihelion, and the perihelion distance, which fully characterized the orbit of that Comet. By the end of the manuscript, Euler uses his method to compute the geocentric longitude and latitude of the Comet in the last observation, and compared them with the measurements. The geocentric longitude agreed quite well, however, the geocentric latitude showed some deviation, which Euler attributed to an error in the observation.

NOTES

^{1}The visual brightness of comets is measured by the visual magnitude scale. It can sometimes be referred to as the “apparent” magnitude scale as it is always based on an object’s brightness as viewed from the Earth. This is the same scale that is used for all astronomical objects from asteroids and meteors to planets, moons, galaxies and spacecraft. There are a few peculiar characteristics of the visual magnitude scale. The first thing that may seem a little odd for a scale is that the lower the number, the brighter the object is. Compared to many other measuring scales, this seems backward. The second characteristic is that each increase in integer represents an increase in brightness of 2.5 times. For example, an object of magnitude 5 is two and a half times the brightness of an object of magnitude 6. Source: http://www.cometwatch.co.uk/comet-info/the-magnitude-scale/.

^{2}The same developments appear in another publication by Euler namely E840—*De motu cometarum in orbitis parabolicis, solem in foco habentibus*. A translation of this manuscript can be found at: https://scholarlycommons.pacific.edu/euler-works/840/.

^{3}In astronomy, a conjunction occurs when any two astronomical objects (such as comets, asteroids, moons, planets, and stars) appear to be close together in the sky, as observed from Earth. If two objects have the same ecliptic longitude, they are considered to be in conjunction with one another. The conjunction under consideration here is that of the Comet and the Sun as seen from the Earth.

^{4}A small periodical change of position in the stars and other heavenly bodies, due to the combined effect of the motion of light and the motion of the observer.

^{5}Joseph Nicholas Delisle (1688-1768), French astronomer who graduated and later taught mathematics and astronomy at the College de France. From 1710 to 1715, he set up an observatory on the Palais de Luxembourg. In 1725 he went to St. Petersburg, where he stayed for 21 years. He was a member of the Imperial Academy of Sciences as well as director of the St. Petersburg Observatory.

^{6}These observations were personally communicated to Euler by Delisle who took them in the St. Petersburg Observatory. Coincidentally, 1742 is the year of the death of the British astronomer Edmund Halley.

^{7}1 Signe (1^{s}) => 30˚. From antiquity through the 18th century, ecliptic longitude was commonly measured using twelve zodiacal signs, each of 30˚ longitude, a practice that continues in modern astrology. The signs approximately corresponded to the constellations crossed by the ecliptic.

^{8}Solar longitude (the ecliptic longitude of the Sun) in effect describes the position of the earth in its orbit, being zero at the moment of the vernal equinox. Solar longitude calculator https://kikimoreau.shinyapps.io/SollongCalc/.

^{9}Roughly, the distance from Earth to the Sun is equal to about 150 million kilometers (93 million miles), which is called Astronomical Unit (AU).

^{10}This is the conjunction time, which was not actually observed, but was estimated by Euler instead.

^{11}These results come from the angular differences in the solar longitude for the different observations shown in Table 1.

^{12}These expressions come from the application of The Law of Sines to the triangles indicated in Figure 8.

^{13}These expressions come from the application of The Law of Sines to the triangles indicated in Figure 9, which are drawn on the ecliptic plane.

^{14}This expression was developed twice by Euler, in §. XIII, and in §. XIV. It was also independently developed by I. H. Lambert in *Insigniores**Orbitae Cometarvm Proprietates* (Most Remarkable Properties of Comets’ Orbits), 1761, and became later known as The Euler-Lambert Equation for Parabolic Motion.

^{15}This happens because the nodal line is the interception of both planes, the ecliptic plane and the plane of the orbit, and then, since the cord *ac* belongs to the plane of the orbit, it necessarily intercepts the plane of the ecliptic in the nodal line.

^{16}This expression was developed in §. XV.

^{17}This is actually the longitude of the Earth as seen from the Sun.

^{18}These relations come from the application of Napier’s rules for this right spherical triangle.

^{19}The ecliptical longitude λ (lambda) is the position along the ecliptic, relative to the vernal equinox (so relative to the stars). If you look at the Sun from the Earth, then you’re looking in exactly the opposite direction than if you look at the Earth from the Sun, so those directions are 180˚ apart. Hence, the ecliptic longitude of the Sun, as seen from the Earth, is equal to λ_{Sun} = λ_{Earth} + 180˚. The value of λ_{Sun} determines when the (astronomical) seasons begin: when λ_{Sun} = 0, then spring begins in the northern hemisphere, and autumn in the southern hemisphere. Each next multiple of 90˚ brings the start of the next season.

References

[1] Boscovich, R. J. (1774). De orbitus cometarum determinandis ope trium observationem parum a se invicem remotarum. Mémoires de Mathématique et de Physique, présentes à Académie Royale des Sciences par divers Savants, 6, 198-215.

[2] Euler, L. (1743). Determinatio orbitae cometae qui mense Martio huius anni 1742 potissimum fuit observatus. Miscellanea Berolinensia, 7, 1-90.

[3] Euler, L. (1744). Theoria motuum planetarum et cometarum. Berlin: Ambrose Haude.

[4] Figueiredo, F. B., & Fernandes, J. M. (2020). José Monteiro da Rocha (1734-1819) and His 1782 Work on the Determination of Comet Orbits. Journal for the History of Astronomy, 51, 461-481.

https://doi.org/10.1177/0021828620947844

[5] Laplace, P. S. (1784). Mémoire sur la Détermination des Orbites des Comètes. Mémoires de l'Académie Royale des Sciences de Paris, 10, 93-146.

https://gallica.bnf.fr/ark:/12148/bpt6k775981/f1.item

[6] Truman, O. H. (1912). Computation of a Comet’s Orbit. MS (Master of Science) thesis, Iowa: State University of Iowa.

https://doi.org/10.17077/etd.c9h10h63