Importance of Conic Section “size” in the Derivation of Propositions X-X VI in Newton’s Principia Book I Gopi Krishna Vijaya Reciprocal System Research Society, 1680 E Atkin Ave, Salt Lake City, UT-84106 Abstract: The “direct problem” of the inverse-square law is examined, by retracing the derivations of Propositions X-XVI in Newton’s Principia. It is found that transfers of constants of proportionality are inconsistent between the propositions, leading to errors in the derived relations for multiple orbits. Proposition X has to be corrected for the relations of time period to the force law by including the area constants for conic sections. Propositions XI-XIII assume the expression 2h?/L to be constant, that in turn implicitly assume a modified version of Kepler’s third law, while Propositions XIV-XVI apply this to multiple orbits leading to further “size” inconsistencies with the original form of Kepler’s third law. It has been shown how this relevance of “size” has been missed in the Newtonian literature, leading to a necessity of revising the existence of an exact inverse-square law. Introduction Newton’s force laws including the inverse-square law have received a thoroughgoing treatment not only in the 18" and 19" century, but also recently. In the past thirty years, the issue of the rigor of the Newtonian proof for the inverse argument of the inverse-square law — that the presence of an inverse-square law implies the existence of orbits which are conic sections — has received renewed attention. In addition to the arguments raised by Weinstock [1] [2] regarding the validity of this proof in the Principia, a series of papers — such as [3] [4] [5] [6] and [7] — have clarified the details of Newton’s proof and justified the rigor of the derivation of conics from the inverse-square law. They have also shown that the required proof not only was possible for Newton, but was most likely accomplished by him using the integral tables available to him at the time [3]. Although there has been extensive commentary on the functional dependences of the laws, such as the proportionality between the force F and the radial dependency 1/r’, the author is not aware of any recent detailed work on the corresponding constants of proportionality for these laws. This has been understandable since the original Principia was written in the language of geometric proportionality, and the later derivations using modern methods of calculus such as [8] utilize constants without a full examination of their interrelationships between propositions (described further in Sections 1-5). Therefore, the analysis in this paper proceeds from a different foundation: that of the use of the constants of proportionality or — simplistically called “sizes” in this work. Prior to an analysis of the propositions related directly to the inverse-square law, viz. Propositions XI — XVI (Book I), Proposition X will be examined in section 1. Following this, the same concepts are applied for the remaining propositions sequentially. 1. Proposition X, Book I In this proposition, the purpose is, according to Newton [9]: Prop. X. Prob. V. Let a body revolve in an ellipse; it is required to find the law of the centripetal force tending toward the center of the ellipse. The relevant diagram is included from the Principia below. Fig. 1: Proposition X, with C as center of forces, and particle at P It must be noted that the word “force” as used by Newton in these sections of the Principia refers to the force of acceleration, and therefore have dimensions of acceleration. The result of the proposition is to deduce that: 2 2 Therefore (by prop. 6, corol. 5), the centripetal force is as —— inversely, that is (because 2BC2 x CA? is given), as = inversely, that is, as the distance PC directly. For convenience of notation, the following variables are used: a = CA; b = BC; r = PC; F = centripetal force Therefore, the proposition states that: ee ee (1) 2BC2xCA2 = 2a?b? At this point it is important to note that the last proportionality (w.r.t ‘r’) only holds when the same ellipse is involved, i.e. where a and b can be considered constants. Where this no longer holds, e.g. between different ellipses, the proportionality of the force may differ. However, Corollary 2, Proposition X states that: COROLLARY 2. And the periodic times of the revolutions made in all ellipses universally around the same center will be equal. For in similar ellipses those times are equal (by prop. 4, corols. 3 and 8), while in ellipses having a common major axis they are to one another as the total areas of the ellipses directly and the particles of the areas described in the same time inversely; that is, as the minor axes directly and the velocities of bodies in their principal vertices inversely; that is, as those minor axes directly and the ordinates to the same point of the common axis inversely; and therefore (because of the equality of the direct and inverse ratios) in the ratio of equality. Here, Proposition IV is invoked to justify the equality of times for ellipses. However, in Proposition IV Corollary 3, it is explained only for circles that: t, =t,°oFar (2) where t,and t, are the time periods of two circular orbits. In Corollary 8, this relationship is extended to “any similar figures whatever,” with the caveat that the relationships are not directly, but only ultimately true (see [10], pg. 173). Newton highlights this disclaimer by mentioning that one must apply “uniform descriptions of areas for uniform motion,” which is a natural result of taking the limit of a curve as being a circle. Any elliptical orbit can be approximated by a circular orbit at the limit, but the converse does not hold. Since it does not follow that a curve where the limit can be approximated by a circle (with the restrictions applied by the scholium to Lemma 11) can have the same relationships between parameters as that of the circle itself, equation (2) does not automatically hold for the ellipse, and will have to be derived from Proposition X itself. Starting once more with equation (1): r F « Fa2p2 (1a) Here, the denominator is related to the square of the area, since the area of an ellipse is mab. By substituting the areal velocity (conventionally taken to be h/2) and the time period (7) conjointly for the area, one gets: a rT = rT a rT 2a7b2 2(h/2n)*T2 h2T2 (3) Where equal time periods are compared (T,=T,), as in Proposition X Corollary 2, the relationship between two ellipses becomes: F. h3 eo Bh te Fp 1 «hy # - (unless hy = hz) (4) 2 Since for the areal velocities, there is no prior reason to assume that they are equal, Newton’s statement that “... For in similar ellipses those times are equal (by prop. 4, corols. 3 and 8)...” is no longer valid. From equation (4), when time periods are equal, force is proportional to the distance if and only if the areal velocities of two ellipses are equal, and not otherwise. Hence conversely, the times will be equal only under the same condition. It is therefore incorrect to represent the proportionality of equation (la) as that between the force and distance alone. Corollary 3 of Proposition IV is not applicable directly to ellipses and Corollary 8, applying only at the limit, does not support the proportionality either. As to the latter part of Corollary 2, where the major axes are equal, the following figure (Fig. 2) will be used for reference. Here b, and bz are the minor axes of the two ellipses with the major axis a being common to both of them. The velocities at the principal vertices are v, and v, in the larger and the smaller orbit respectively, here taken in opposite directions for ease of representation (note that the ratios depend on the magnitudes, not the directions.) U2 Fig. 2: Two ellipses, with common major axis “a” and minor axes b; > bo. According to the definition of areal velocity, the relation between time periods is given by: Ty _ Ay mn dAz/dat ee Tab, Fe av2 —_ biV2 T2 _ A2 dA,/dat a Ttabz avy b2vz (5) Here A; is the area of the i-th ellipse, and dA;/dt is its rate of change. From this point onwards, Corollary 2 claims that the ratio in equation (4) is one of equality, by comparing ordinates at the limit. However, there is a straightforward way of proceeding, described as follows. Since the radius of curvature R of an ellipse at the major axis vertex is given by b*/a, for the two ellipses, the following relation holds: Ba = bi/a = (2)" (6) R, -b3/a bz From Proposition IV (F « v?/R), and combining equations (6) and (1a) at the principal vertices on major axis: Fy vi/Ra ey Re (vex is et ya PAE 2h a= POLL ' F, =v /Ro V2 Ry byVv2 2 Fy a a (72) = / > SS Tb i) apt! Gap? by ( ) Therefore from equation (7a), (7b) and (5), it is straightforward to obtain: Ty byv2n _ by Fo = = — ay 8 T2 b2vz b2 Fy ( ) Hence, the velocities at the principal vertices are equal but the time periods are not, contrary to what is deduced in Corollary 2 of Proposition X. Time periods vary as the minor axes. Since the derivation just described does not depend on any other assumptions with regard to the force law and all the constants are accounted for, it does not suffer from the shortcomings of assuming a dependence of the force solely on the distance. Following up with the constants throughout the derivation is hence seen to be of utmost importance. It is quite surprising that in common references, such as [8] (pg. 90, (b) (i)) and [10] (pg. 218-224), this dependence of force on the constants of the ellipse has been completely neglected. In fact, the question is posed to the general student audience in [10] (pg. 216) that: This corollary claims that, no matter how the sides of the ellipse are squished in or stretched out, the periodic times will remain the same. Does this seem surprising to you? Would you expect a body to take as long to go around a little bitty orbit as around a huge one? Are you curious how Newton would prove such a thing? Can we see how the force law of I.10 would have that result? Although the question is posed, the answer is not pursued systematically. Since the physical intuition is evidently surprised by the lack of effect of such a “size” difference, that relation has been followed up, and corrected for the neglect of this “size” (which is — in this particular case — a*b?.) 2. Proposition XI, XII, XIII (Book I) In these propositions, the force law for an object moving in a conic section is described. Since the expression from Proposition VI corollary 1 is used in these proofs, it is important to clarify the constants in it. The relevant diagram from Principia is reproduced below for clarity. ww Fig. 3: Movement of an object at P on a curve, Proposition VI, Corollary 1, Book I The expression for the force is ( [8], pg. 77): 2h2QR SP2OF? Fa —*_orF= limg..p SP2QT2 ) Where h/2 is the areal velocity, as usual. In the case of elliptic motion (Proposition XI), hyperbolic motion (Proposition XII), and parabolic motion (Proposition XIII) the following expression is obtained: 5 BNF LS p2 (10) Here, L is the principal latus rectum of the conic section in question, which is true at the limit of equation (9). The last line of each of the propositions ends with: 1 Fea =p (11) In the transition from equation (10) to (11), it is assumed that 2h?/L is a constant of proportion. While this is true as long as the treatment is restricted to one conic section, there is no necessity that it remains a constant when comparing different orbits and different bodies. The constant h representing areal velocity is related to the radius vector and tangential velocity at a point, and is independent of the latus rectum. If the two arbitrary constants, h and L, are to have a constant relation between them, that relationship has to be proved from the dynamics of the particle itself when applied to different orbits. An additional consequence of treating h?/L as a constant can be demonstrated by the example of its use in an ellipse with major axis a, minor axis b and time period T: 2h he 2(2mab/T)* 9 a? L a(1-e7) 2b7/a T (12) Assuming 2h?/L to be a constant is exactly the same as assuming the validity of Kepler’s third law, for a mean distance of a except that it stands for the major axis here (for more on this concept, see Section 4). The force F, as a function F (h, L, SP) can be reduced to F (SP) only under this condition. A physical or mathematical interpretation for the orbit requires a demonstration that 2h?/L is necessarily a constant, prior to moving ahead with the next proposition that utilizes F as a function of SP alone. Hence the important fact to notice is that 2h?/L is NOT proved to be a constant when comparing different orbits. 3. Proposition XIV, Book I This proposition states that ( [9], pg. 467): Prop. XIV, Thm. 6: If several bodies revolve about a common center and the centripetal force is inversely as the square of the distance of places from the center, I say that the principal latera recta of the orbits are as the squares of the areas which the bodies describe in the same time by radii drawn to the center. If several bodies are given revolving about a common center, then both the areal velocities h and their corresponding latera recta L, are both assumed to be given. It is on this basis that Newton proceeds to establish a relationship regarding latera recta. However, notwithstanding the lack of proof for the constancy of 2h*/L for several orbits, in this Proposition it is assumed that the centripetal force is varying only with respect to SP”, as a hypothesis: But the minimally small line QR is in a given time as the generating centripetal force, that is (by hypothesis), inversely as SP”. This means that, by assuming 2h?/L to be a constant, one form of Kepler’s Third Law is set up as a hypothesis in Proposition XIV. If the bodies are moving in ellipses for example, for each individual body, Proposition XI holds. But that does not mean that 2h?/L is the same for each one. There is no prior mathematical or physical requirement that ensures that the hypothesis is true, except for the actual empirical data organized by Kepler. Further in this Proposition, it is shown that: h«vL (13) However, this is just a restatement of 2h?/L being a constant, something that holds true for one orbit alone, generalized for “several bodies revolving around a center.” Hence, the actual accomplishment of this Proposition is to introduce a form of Kepler’s third law into the development of the argument, without however making it explicit. Kepler’s law is implicitly assumed — not proved or justified. The corollary to Proposition XIV adds nothing new to the proposition, but merely restates equation (13) as: AxvL:T (14) This is done by relating areal velocity to A/T, where A is the area of the ellipse and T its periodic time. 4. Proposition XV-XVI The previous section gives a different view of Proposition XV, which is where Kepler’s Law is written out for ellipses in the Principia. In Proposition XV, the relationships of Proposition XIV are developed further, to show that: Axa:b«avalL (15) This is from the definition of L. Combining equation (15) with (14) results in: ava «Tora? « T? (16) However, this is simply a restatement of equation (13) in the exact same way equation (12) was obtained through an assumption of Kepler’s third law. Therefore, one simply recovers from the equations what had initially been introduced into them, and the sequence from equation (10) to (16) is equivalent to a tautology: Kepler’s third law implies Kepler’s third law. Further proof of this deduction of an inverse-square law from Kepler’s laws, rather than the other way around, is even found in Newton’s letters ( [11], pg. 143): ... from Kepler’s rule of the periodical times of the planets being in a sesquialternate proportion of their distances from the centres of their orbs, I deduced that the forces which keep the planets in their orbs must be reciprocally as the squares of the distances from the centres about which they revolve... If indeed Kepler’s laws are the source of the Propositions XIV and XV, it is interesting to note that the Principia does not mention Kepler even once in Book I. If the above passage is indeed accurate, the implicit application of Kepler’s third law must be more closely examined. Proposition XV says ( [9], pg. 468): Prop. XV Thm. 7: Under the same suppositions as in prop. 14, I say that the squares of the periodic times in ellipses are as the cubes of the major axes. However, Kepler’s third law is originally stated thus [12]: But it is absolutely certain and exact that the proportion between the periodic times of any two planets is precisely the sesquialternate proportion of their mean distances, that is, of their actual spheres, though with this in mind, that the arithmetic mean between the two diameters of the elliptical orbit is a little less than the longer diameter. There is a clear discrepancy between “... of the major axes” of Newton and “of their [arithmetic] mean distances” derived by Kepler. Thus, from equation (12) onwards, the form of Kepler’s Law being implicitly assumed by Newton is different from the actual empirical law. Even if one grants that, as described by him above in his letter and as has been shown in the course of the derivations, Newton assumes Kepler’s Law and derives the inverse-square law from it, it does not suffice to explain why the mean distance was replaced by the major axis, since that is not the correct form of the third law. In other words, extending equation (12), it can be said that: +b 2h? _ 2(2mab/T)? | a® 4 (3 " Sense mee ee wes sue ss L 2b7/a ce Te ee Hence, Propositions XIV and XV are both compromised due to this approximation, as it turns out that the constant of proportionality between orbits (2h*/L) is no longer a constant — only an approximation that holds where a ~ b. Hence the inverse-square law results from a distortion of Kepler’s law that excludes the importance of the minor axis of an ellipse, which is one of the measures of its “size”. In effect, a relationship that holds true for circles is extrapolated to ellipses. This discrepancy has immediate consequences in the corollary. The corollary to Proposition XV states that: Therefore the periodic times in ellipses are the same as in circles whose diameters are equal to the major axes of the ellipses. Hence, one has the surprising result that there is no dependence of the periodic times on the minor axis of the ellipse, making a circle and an ellipse equivalent as long as the major axis matches up with the diameter. Just as in Proposition X, the implication of this corollary has also been noticed in Densmore’s analysis ( [10], pg. 266-267): It was intriguing, in fact downright strange, to see in Proposition I.10 Corollary 2 that you could squash the sides of the ellipse all you want and yet the periodic time would remain the same. The different ellipses could have different major axes as well as different minor axes... But notice that here in Proposition 15 Corollary we have something else even more astonishing. What we see is that this apparent indifference to stretching doesn’t just hold for ellipses with center of forces at the center. According to this proposition and corollary, it also holds when the center of forces is at a focus and we have an inverse- square law. Once it is understood that the dependence on major axis has been introduced by fiat, through an altered form of Kepler’s third law that applies for circles only, it becomes clear that the derivation is not just “intriguing”, “strange” or “astonishing”, but erroneous, and the strangeness to physical intuition follows from this fact. Proposition XVI proceeds to establish a relationship between the velocity and the perpendicular distance between the velocity vector and the radius vector from the center of forces. Once more, it is a different rearrangement of equation (13): havVL > SY:-|[vBlxVvVL 3 SY «= (13a) Here SY is the perpendicular let fall upon the velocity vector v (see Fig. 3). Corollaries to Proposition XVI serve either the purpose of algebraic rearrangements (Corollaries 1, 2 and 5), comparison between generic conic sections (Corollaries 6 and 7) and comparison between circles and generic conic sections (Corollaries 3, 4, 8, 9). Of these, the comparison with the circle makes use of the equivalence of major axis with radius that was already examined in the corollary of Proposition XV, and hence the same error is propagated. Discussion Propositions XIV, XV and XVI hence reckon with a version of Kepler’s Law that is modified from the original empirical law, and it is quite possible to miss the mathematical difference between equations (12) and (12a) when it is compared with real astronomical data because in reality, eccentricities are small and the condition a ~ b holds. It is hence something relatively easy to miss. The actual equation, if the correct dependence is taken into account is: Lr? r2 T2 r2 -T2 a+b r2 \at+b 2h? 42a — am? (44)8 ( 2a ) Ke ( 2a ) = =— a ee | = Zz (17) Here, K, is Kepler’s constant, derived empirically, a and b are the major and minor axes of an ellipse respectively (restricting analysis to only ellipses for this paper). Since 2a/(a + b)? is not a constant across different ellipses, the inverse-square law is an approximation at best, for ellipses of small eccentricity. The other alternative is to suggest that each orbit has its own particular inverse-square law — an approach that compromises the generality of a ‘law’. It is worth examining how the successors of Newtonian theory treated these constants. Euler, in his treatise Mechanica where he derives analytically several propositions from the Principia, says [13]: 718. Let the given curve again be an ellipse, but with the centre of force C placed in the other focus. The transverse axis of this is put equal to A and the latus rectum equal to L, and from the nature of the ellipse we find that 4p? = “ Hence of differentiating, we 2 272,,2 have 8pdp = ae Truly since 16p* = on then a = & and consequently: _ 4ch? ee Therefore the centre of force is inversely proportional to the square of the distance of the body from the centre of force C. Here P is the force. The book carries no further treatment about whether or not 4ch*/L is a constant across different conic sections. Laplace [14] describes the same process by setting up a differential equation for identifying the general force law for a general orbit, and showing that for conic sections the force law takes the form (a full quote is included here since this is the most common approach in modern notation): [380] ese es {oo a-(1-e?) r? [380’] therefore, the orbits of the planets and comets being conic sections, the force @ will be inversely proportional to the square of the distance of the centre of the planet from the centre of the sun. We also perceive, that if the force @ be inversely proportional to the square of the distance, or expressed by h/r?, h being a constant coefficient, the preceding equation of the conic sections, will satisfy the differential equation (4) [377] between r and v, which gives the, expression of the force, when we change ¢ into h/r?. We shall then have c2 [380”] h= EES which forms an equation of condition between the two arbitrary constant quantities a and e of the equation of conic sections... [3 80"] Hence it follows, that if the described curve be a conic section, the force will be in the inverse ratio of the square of the distance; and conversely, if the force be in the inverse ratio of the square of the distance, the described curve will be a conic section. These equations are identical to equation (12) with c being the areal velocity. Laplace therefore does not hesitate to assume his / to be a constant, even though the relationship is between two arbitrary constants, and hence not a given across orbits. This constancy is not tested. Further on, he describes: 3. The intensity of the force @, relative to each planet and comet, depends on the 2 coefficient aoe the laws of Kepler furnish the means of determining it... (1-e?) [382'] With respect to the planets, the law of Kepler, according to which the squares of the times of their revolutions, are as the cubes of the transverse axes of their ellipses, gives T? = k?- a3, k being the same for all the planets; therefore we shall have __ 2mVa(1—-e?) [383] C Hence Laplace, following Newton, takes Kepler’s Law with respect to the “cubes of the transverse axes” only, and not the mean distance as originally described by Kepler. The law has been modified in this way even until recently (see [8], ch. 4), and is necessary for the derivation of the inverse-square law and conic sections from one another, whether through geometry or calculus. It must also be noted that Laplace later uses mean distances (in [385 vai without thereby modifying the inverse-square relation — by simply modifying Kepler’s law back to its original form. It is necessary to clarify the preceding discussion in the light of the recent controversy over the derivation of the converse of the inverse-square law. Weinstock’s position ( [1] [2]) was that he agreed with the derivation from conic section orbit (CSOF) to inverse-square force (isf) law, and argues that Newton did not demonstrate the converse. Using his notation, he claims that (isf) =>(CSOF) was not convincingly demonstrated in the Principia. However, subsequent research (particularly [3]) has shown that even if the proof was not clearly demonstrated in the Principia, Newton had the tools to derive it correctly from other propositions. At the time of Newton, Kepler’s laws for elliptic orbits were well-known, but the inverse-square law was yet to be demonstrated. Therefore, any proof of the inverse-square law had to satisfy both, in this order of priority: (a) (CSOF) >(isf) (b) (isf) >(CSOF) While Weinstock and others, from the time of Bernoulli, have restricted their discussion to (b), the present discussion focuses exclusively on the validity of (a). Unless it is convincingly demonstrated that formulation of an inverse-square force for a series of orbits is justified, the converse, i.e. part (b) becomes unnecessary. Conclusion The importance of following through with the constants of proportionality, which mostly determine the “sizes” of the conics in one form or another, has been demonstrated by systematically examining the derivations in Propositions X, XI, XU, XIII, XIV, XV and XVI of the Principia, Book I. The results of this investigation are as follows: (i) Proposition X: An error has been introduced in the force dependence equation as applied to multiple orbits by neglecting the constant relating to the area. Reference to Proposition IV is shown to be invalid. By correcting the equations, the time periods of ellipses on the same major axis but different minor axes are seen to be proportional to the minor axes, and not equal as claimed. (ii) Propositions XI, XII, XIII: It is clarified that these propositions show only a conditional dependence on the inverse-square of the distance, since they assume that the relationship of two arbitrary constants (2h*/L) is also a constant — a fact that has to be proved. Assumption of this constancy is shown to be equivalent to assuming Kepler’s third law, in one specific form with major axes only. (iii) Proposition XIV, XV, XVI: Proposition XIV assumes the constancy of 2h?/L for multiple orbits, without a proof or a demonstration of its truth, and Proposition XV rediscovers Kepler’s law mentioned in (11) after having implicitly assumed it. It is shown that the form of the Kepler’s law should depend on both axes, not only on the major axis, making both the implicit assumption and its subsequent rediscovery in Proposition XV invalid. Proposition XVI rearranges the relations of the previous two, without adding anything new. It is thus clear from (1), (ii) and (iii) that the primary purpose of these propositions viz. showing that conic section orbits imply inverse-square law — (CSOF) =(isf) — is not satisfied, and is inconsistent with Kepler’s third law. It is also seen that the same error has been repeated unacknowledged in all of the major works on the subject. The development in this paper shows that this series of propositions have to be rewritten more accurately to clarify the real principles of celestial movements. Funding Statement The author would like to acknowledge the funding received by the Rudolf Steiner Charitable Trust Advised Fund of RSF Social Finance for the support of this work. References [1] R. Weinstock, "Dismantling a centuries-old myth: Newton’s Principia and inverse-square," American journal of Physics, vol. 50, pp. 610-617, 1982. [2] R. Weinstock, "Inverse-Square Orbits in Newton’s Principia and Twentieth-Century Commentary," Archive for History of the Exact Sciences, vol. 55, pp. 137-162, 2000. [3] J.B. Brackenridge, "Newton’s Easy Quadratures "Omitted for the Sake of Brevity"," Archive for History of Exact Sciences, vol. 57, pp. 313-336, 2003. [4] [5] [6] [7] [8] [9] H. Erlichson, "The Visualization of Quadratures in the Mystery of Corollary 3 to Proposition 41 of Newton's Principia," Historia Mathematica, vol. 21, pp. 148-161, 1994. M. Nauenberg, "The early application of the calculus to the inverse-square force problem," Archive for History of Exact Sciences, vol. 64, pp. 269-300, 2010. N. Guicciardini, "Johann Bernoulli, John Keill and the Inverse Problem of Central Forces," Annals of Science, vol. 52, pp. 537-575, 1995. N. Guicciardini, "Conceptualism and contextualism in the recent historiography of Newton's Principia," Historia Mathematica, vol. 30, pp. 407-431, 2003. S. Chandrasekhar, Newton's Principia for the Common Reader, New York: Oxford University Press, 1995. I. Newton, The Principia, Mathematical Principles of Natural Philosophy, 3rd ed. (1726), translated by I. Bernard Cohen and Anne Whitman, Berkeley: University of California Press, 1999. [10] D. Densmore, Newton's Principia, The Central Argument, Santa Fe, New Mexico: Green Lion Press, 2003. [11] R. S. Westfall, Never at Rest, a Biography of Isaac Newton, New York City: Cambridge University Press, 1980. [12] J. Kepler, "(Vol 209) Book 5, Chapter 3," in Harmony of the World, Translated by Aiton, Duncan & Field, Philadelphia, American Philosophical Society, 1997, p. 411. [13] L. Euler, "Chapter 5c, Section 718," in Mechanica, Vol. I, St. Petersburg, Academy of Sciences, 1736. [14] P. S. Laplace, "Book 2, Section 3," in Mécanique céleste Vol. I (translated by N. Bowditch), Boston, Hillard, Gray, Little and Wilkins, 1829, pp. 246-247.