{
  "id": "c0fe91d091f86357bc8955e92eea16f4cf90068d",
  "text": "VIII. On the Theory of the Moon. By John William Lubbock, Esq. V.P. and Treas. R.S.\n\nReceived and Read March 13, 1834.\n\nWhen I commenced the investigations relating to the theory of the moon which I have had the honour to communicate to the Society, I proposed to show how, by a different but more direct method, the numerical results given by M. Damoiseau might be obtained. The approximations were in fact carried much further by M. Damoiseau than had been done before, and the details which accompany M. Damoiseau's work evince at once the immense labour of the undertaking, and inspire confidence in the accuracy of the results offered. But the state of the question is now changed by the appearance of M. Plana's admirable work, entitled \"Théorie du Mouvement de la Lune,\" in which, although M. Plana employs the same differential equations as those used by M. Damoiseau, and obtains in the same manner finally the expressions for the coordinates of the moon, in terms of the mean longitude by the reversion of series, yet M. Plana's expressions have a very different analytical character and importance, from the circumstance that the author develops all the quantities introduced by integration, according to powers of the quantity called $m$, which expresses the ratio of the sun's mean motion to that of the moon. In this form of the expression the coefficients of the different powers of $m$, of the eccentricity, &c., are determinate, as are, for example, the numerical coefficients in the expression for the sine in terms of the arc, and other similar series. An inestimable advantage results from this procedure, which more than compensates for the great increase of labour it occasions, by diminishing the danger of neglecting any terms of the same order as those taken into account, and by affording the means of verifying many terms long before final and complete results shall have been obtained independently by myself or any other person. By treating the differential equations in which the time is the independent variable, as I have proposed, similar results to those of M. Plana may be obtained directly; but the calculations which are required in either method are so prodigiously irksome and laborious, that until identical expressions have actually been obtained independently, to the extent of every sensible term, the theory of the moon cannot, I think, be considered complete. It might, indeed, be supposed that already, through the labours of mathematicians, from Clairaut to the present time, the numerical values of the coefficients of the different inequalities were ascertained with sufficient accuracy for practical purposes, and that any further researches connected with the subject would be more likely to gratify curiosity than to lead to any useful result.\nAstronomical observations are now made with so great precision, that the numerical values of the coefficients are wanted to at least the tenth of a second of space: very few, however, of the coefficients of MM. Damoiseau and Planck agree so nearly, and some differ much more, as may be seen in the following comparison of the numerical values of the coefficients of some of the arguments in the expression for the true longitude of the moon in terms of her mean longitude, being indeed those which differ the most.\n\n| Argument | Damoiseau | Planck |\n|----------|-----------|--------|\n| 1 \\(2\\tau\\) | + 2370-00 | + 2370-320 |\n| 2 \\(\\xi\\) | + 22639-70 | + 22641-626 |\n| 3 \\(2\\tau - \\xi\\) | + 4589-61 | + 4585-648 |\n| 4 \\(2\\tau + \\xi\\) | + 192-22 | + 192-146 |\n| 5 \\(\\xi_1\\) | - 673-70 | - 668-644 |\n| 6 \\(2\\tau - \\xi_1\\) | + 165-56 | + 165-850 |\n| 7 \\(2\\tau + \\xi_1\\) | - 24-82 | - 23-611 |\n| 8 \\(2\\xi\\) | + 768-72 | + 769-477 |\n| 9 \\(2\\tau - 2\\xi\\) | + 211-57 | + 212-363 |\n| 10 \\(2\\tau + 2\\xi\\) | + 14-74 | + 14-119 |\n| 11 \\(\\xi + \\xi_1\\) | - 109-27 | - 111-099 |\n| 12 \\(2\\tau - \\xi - \\xi_1\\) | + 207-09 | + 209-742 |\n| 22 \\(2\\tau + 3\\xi\\) | + 1-27 | + 3-309 |\n| 24 \\(2\\tau - 2\\xi - \\xi_1\\) | + 8-99 | + 7-762 |\n| 27 \\(2\\tau - 2\\xi + \\xi_1\\) | + 2-55 | - 1-395 |\n| 64 \\(2\\tau + 2\\gamma\\) | - 5-75 | - 3-376 |\n| 65 \\(\\xi - 2\\gamma\\) | + 39-51 | + 37-191 |\n| 110 \\(\\tau - \\xi + \\xi_1\\) | + 2-05 | + 4-66 |\n| 136 \\(4\\tau - 2\\xi\\) | + 31-19 | + 34-518 |\n| 4 \\(2\\tau - 2\\xi - \\xi_1\\) | + 3-05 | + 1-197 |\n\nWhen the coefficients of the inequalities have been determined analytically, it remains to determine with corresponding precision the numerical values of the arbitrary quantities \\(m\\), \\(e\\), and \\(\\gamma\\). The quantity \\(m\\) is already accurately known, but the quantities \\(e\\) and \\(\\gamma\\) must be obtained from the coefficients of \\(\\sin \\xi\\) in the expression for the longitude, and of \\(\\sin \\eta\\) in the expression for the latitude, by the reversion of series; and it seems to me that the manner in which these arbitrary quantities are to be determined must be carefully and rigorously defined.\n\nI propose to obtain the expression for the radius vector by means of the equation,\n\n\\[\n\\frac{d^2 r^2}{dt^2} - \\frac{\\mu}{r} + \\frac{\\mu}{a} + 2 \\int d R + r \\frac{dR}{dr} = 0.\n\\]\n\nIn order to integrate this equation, I suppose\n\n\\[\n\\frac{a}{r} = 1 + r_0 + r_1 \\cos 2\\tau + e \\left(1 - \\frac{e^2}{8}\\right) \\cos \\xi + &c.\n\\]\n\nIf \\(r\\) be used to denote the terms in \\(r\\) which are found in the elliptic expression, so that\n\n\\[\n\\frac{a}{r} = 1 + e \\left(1 - \\frac{e^2}{8}\\right) \\cos \\xi + e^2 \\left(1 - \\frac{e^2}{3}\\right) \\cos 2\\xi + \\frac{9}{8} e^3 \\cos 3\\xi + \\frac{4}{3} e^4 \\cos 4\\xi + &c.,\n\\]\n\nand\n\n\\[\n\\frac{a}{r} = \\frac{a}{r} + a \\delta \\frac{1}{r}.\n\\]\n\\[\n\\frac{d^2}{dt^2} r^2 - \\frac{d^3}{dt^3} r^3 \\frac{1}{r} + \\frac{3}{2} \\frac{d^4}{dt^4} r^4 \\left( \\frac{1}{r} \\right)^2 - \\frac{2}{2} \\frac{d^5}{dt^5} r^5 \\left( \\frac{1}{r} \\right)^3 \\\\\n- \\frac{\\mu}{r} + \\frac{\\mu}{a} + 2 \\int dR + r \\frac{dR}{dr} = 0.\n\\]\n\nLet \\( r_n \\) be that part of the coefficient of the \\( n \\)th argument in the development of the quantity\n\n\\[- r^3 \\frac{1}{r} + \\frac{3}{2} r^4 \\left( \\frac{1}{r} \\right)^2 + &c.,\\]\n\nwhich corresponds to the argument of which \\( n \\) is the index, and let \\( R_n \\) be the coefficient corresponding to the argument of which \\( n \\) is the index in the development of \\( R \\), \\( R'_n \\) the corresponding coefficient in the development of that part of \\( \\delta dR \\) which is multiplied by \\( m \\), and only arises in the second approximation, with its sign changed, then the quantities \\( r_n \\) are given by equations similar to the following,\n\n\\[ r_1 \\left\\{ 1 + 3 e^2 \\left( 1 + \\frac{e^2}{8} \\right) \\right\\} \\{2 - 2m\\}^2 - 1 \\} = (2 - 2m)^2 r_1 \\\\\n- 2 \\left\\{ \\frac{2}{2 - 2m} + 1 \\right\\} m^2 R_1 + \\frac{m^3}{(2 - 2m)} R'_1 \\}\n\nPassing over terms given by M. Planck and arising from the first approximation with which I agree, I come to \\( r_{22} \\).\n\n\\[ r_{22} \\left\\{ \\{2 - 2m + 3c\\}^2 - 1 \\right\\} = (2 - 2m + 3c)^2 r_{22} \\\\\n- 2 \\left\\{ \\frac{2 + 3c}{(2 - 2m + 3c)} + 1 \\right\\} m^2 R_{22} \\]\n\n\\[ r_{22} = \\frac{3}{2} r_{10} - \\frac{1}{16} r_1 \\quad r_{10} = \\frac{7}{2} m^2 \\quad r_1 = m^2 \\quad *R_{22} = -\\frac{25}{32} \\]\n\n\\[ r_{22} = \\frac{25 \\cdot 3 \\cdot 7}{24 \\cdot 2 \\cdot 2} m^2 - \\frac{25}{24 \\cdot 16} m^2 + \\frac{2 \\cdot 2 \\cdot 25}{24 \\cdot 32} m^2 = \\frac{2125}{384} m^2 \\]\n\nM. Planck has \\(-\\frac{1175}{384} m^2\\)\n\n\\[ r_{25} \\left\\{ \\{2 - m + 2c\\}^2 - 1 \\right\\} = (2 - m + 2c)^2 r_{25} - 2 \\left\\{ \\frac{2 + 2c}{2 - m + 2c} + 1 \\right\\} m^2 R_{25} \\]\n\n\\[ r_{25} = \\frac{3}{2} r_{13} \\quad r_{13} = -\\frac{33}{32} m^2 \\quad R_{25} = \\frac{3}{8} \\]\n\n\\[ r_{25} = -\\frac{16 \\cdot 3 \\cdot 33}{15 \\cdot 2 \\cdot 32} m^2 - \\frac{2 \\cdot 2 \\cdot 3}{15 \\cdot 8} m^2 = -\\frac{7}{4} m^2. \\quad \\text{M. Planck has } \\frac{7}{2} m^2 \\]\n\n\\[ r_{45} \\left\\{ \\{2 - m - 3c\\}^2 - 1 \\right\\} = \\{2 - m - 3c\\}^2 r_{45} - 2 \\left\\{ \\frac{2 - 3c}{2 - m - 3c} + 1 \\right\\} m^2 R_{45} \\]\n\n* Wherever I have found a disagreement with the result of M. Planck, as this might arise from an error in my development of \\( R \\), I have verified the terms employed.\n\nMDCCCXXXIV.\n\\[ r_{45} = \\frac{3}{2} r_{27} - \\frac{1}{16} r_7 \\]\n\\[ r_{27} = \\frac{15}{8} m^2 \\]\n\\[ r_7 = -\\frac{m^2}{2} \\]\n\\[ R_{45} = \\frac{7}{64} \\]\n\n\\[ r_{45} = \\frac{2.2.8}{3.15} m + \\frac{m}{2.16.2} + \\frac{2.2.7}{2.64} m = \\frac{105}{64} m \\]\nM. Plana has \\( \\frac{285}{64} m \\)\n\n\\[ r_{53} \\left\\{ (c + 3m)^2 - 1 \\right\\} = (c + 3m)^2 r_{53} - 2 \\left\\{ \\frac{c}{c + 3m} + 1 \\right\\} m^2 R_{53} \\]\n\n\\[ r_{53} = \\frac{3}{2} r_{35} \\]\n\\[ r_{35} = -\\frac{53}{16} m^2 \\]\n\\[ R_{53} = \\frac{53}{32} \\]\n\n\\[ r_{53} = -\\frac{3.53}{6.2.16} m - \\frac{2.2.53}{6.32} m = -\\frac{127}{64} m \\]\nM. Plana has \\( -\\frac{53}{32} m \\)\n\n\\[ r_{56} \\left\\{ (c - 3m)^2 - 1 \\right\\} = (c - 3m)^2 r_{56} - 2 \\left\\{ \\frac{c}{c - 3m} + 1 \\right\\} m^2 R_{56} \\]\n\n\\[ r_{56} = \\frac{3}{2} r_{35} \\]\n\\[ R_{56} = \\frac{53}{32} \\]\n\n\\[ r_{56} = \\frac{3.53}{6.2.16} m + \\frac{2.2.53}{6.32} = \\frac{371}{192} m \\]\nM. Plana has \\( \\frac{53}{64} m \\)\n\nThe development of \\( R \\) which I gave*, results from the substitution of the elliptic values of the coordinates of the sun and moon in the disturbing function. The elliptic expression for the radius vector contains no term of which the argument is \\( \\xi - 2 \\eta \\), the longitude (\\( \\lambda \\)) contains the term \\( + \\frac{3}{4} e \\gamma^2 \\sin (\\xi - 2 \\eta) \\). This is changed when the disturbing function is considered.\n\n\\[ r_{65} \\left\\{ (c - 2g)^2 (1 - 3r_0) - 1 \\right\\} = (c - 2g)^2 r_{65} - 2.2.m^2 R_{65} \\]\n\n\\[ r_{65} = \\frac{3}{2} r_{62} \\]\n\\[ r_{62} = \\frac{m^2}{2} \\]\n\\[ c = 1 - \\frac{3}{4} m^2 \\]\n\\[ g = 1 + \\frac{3}{4} m^2 \\]\n\n\\[ r_0 = \\frac{m^2}{6} \\]\n\\[ R_{65} = \\frac{9}{8} + \\frac{1}{2} r_{65} \\]\n\\[ (c - 2g) 2 (1 - 3r_0) = 1 + 4m^2 \\]\n\n\\[ r_{65} = \\frac{3}{4.2.2} - \\frac{2.2}{4} \\left( \\frac{9}{8} + \\frac{1}{2} r_{65} \\right) \\]\n\n\\[ r_{65} = -\\frac{5}{8}. \\]\n\nThis term, produced by the disturbing force, although independent of \\( m \\), together with the corresponding term in \\( \\lambda \\), renders in a certain sense incomplete the coefficients of all terms in my development of \\( R \\), of which the arguments are any combinations of the quantity \\( \\xi - 2 \\eta \\).\n\n\\[ r_{73} \\left\\{ (2 - 3m - 2g)^2 - 1 \\right\\} = (2 - 3m - 2g)^2 r_{73} - 2m^2 R_{73} \\]\n\n\\[ r_{73} = 0 \\]\n\\[ R_{73} = -\\frac{21}{16} \\]\n\n\\[ r_{73} = -\\frac{2.21}{16} m^2 = -\\frac{21}{8} m^2 \\]\nM. Plana has \\( \\frac{7}{8} m - \\frac{7}{8} m^2 \\)\n\n* Philosophical Transactions, 1831, p. 263.\nThe term above, $-\\frac{5}{8} e \\gamma^2 \\cos (\\xi - 2 \\eta)$, introduces the term $+\\frac{3}{4} e \\gamma^2 \\sin (\\xi - 2 \\eta)$ in the longitude, instead of $-\\frac{e \\gamma^2}{2} \\sin (\\xi - 2 \\eta)$. The terms in $R$ produced in consequence may easily be found from the formula\n\n$$\\delta R = -a \\left( \\frac{dR}{da} \\right) r \\delta \\frac{1}{r} + \\frac{dR}{d\\tau} \\delta \\lambda$$\n\ntaking\n\n$$r \\delta \\frac{1}{r} = -\\frac{5}{8} e \\gamma^2 \\cos (\\xi - 2 \\eta), \\quad \\text{and} \\quad \\delta \\lambda = \\frac{5}{4} e \\gamma^2 \\sin (\\xi - 2 \\eta)$$\n\nand I find that $R$ contains, instead of the terms corresponding to the same arguments given in the Philosophical Transactions, 1831, p. 263.\n\n$$+\\frac{15}{32} m^2 e \\gamma^2 \\cos (2 \\tau - \\xi + 2 \\eta) - m^2 e \\gamma^2 \\cos (2 \\tau + \\xi - 2 \\eta)$$\n\n$$+\\frac{39}{32} m^2 e e_1 \\gamma^2 \\cos (\\xi + \\xi_1 - 2 \\eta) + \\frac{105}{64} m^2 e e_1 \\gamma^2 \\cos (2 \\tau - \\xi - \\xi_1 + 2 \\eta)$$\n\n$$+\\frac{33}{64} m^2 e e_1 \\gamma^2 \\cos (e \\tau + \\xi + \\xi_1 - 2 \\eta) - \\frac{15}{84} m^2 e e_1 \\gamma^2 \\cos (2 \\tau - \\xi + \\xi_1 + 2 \\eta)$$\n\n$$-\\frac{231}{64} m^2 e e_1 \\gamma^2 \\cos (2 \\tau + \\xi - \\xi_1 - 2 \\eta)$$\n\nThe coefficient of arg. 77 is easily found as follows:\n\n$$R = -\\frac{r^2}{4r_i^3} \\{1 + 3 \\cos (2 \\lambda - 2 \\lambda) - 2 s^2\\}.$$  \n\nThis term can only arise from\n\n$$-\\frac{r^2}{4r_i^3} \\{1 - 3 s^2\\}$$\n\n$$= \\frac{r^3}{2r_i^3} \\delta \\frac{1}{r} + \\frac{3r^2}{4r_i^3} s^2.$$  \n\nIn which expression it is sufficient to write for $\\delta \\frac{1}{r}$,\n\n$$-\\frac{5}{8} e r^2 \\cos (\\xi - 2 \\eta)$$\n\nand to make\n\n$$*s^2 = -\\frac{\\gamma^2}{2} \\cos 2 \\eta + \\gamma^2 e \\cos (\\xi - 2 \\eta) + \\frac{1}{4} \\gamma^2 e^2 \\cos (2 \\xi - 2 \\eta)$$\n\nwhich gives\n\n$$R_{77} = \\frac{3.5}{2.2.8} + \\frac{3}{4.4} - \\frac{3}{4} + \\frac{3}{4.2.2.2} = 0.$$  \n\n* This is not the expression for $s^2$ in the elliptic movement; the last term is altered by the disturbing force.\nWhen the elliptic values of the coordinates are substituted in the disturbing function, the term in question arises only from the expansion of the quantity\n\n\\[ + \\frac{3}{4} \\frac{r^2}{r_i^3} s^2, \\]\n\nand in the elliptic motion\n\n\\[ s^2 = -\\frac{\\gamma^2}{2} \\cos 2\\eta + \\gamma^2 e \\cos (\\xi - 2\\eta) - \\frac{3}{8} \\gamma^2 e^2 \\cos (2\\xi - 2\\eta) + \\&c. \\]\n\n\\[ R_{77} = -\\frac{3}{4} \\cdot \\frac{3}{8} - \\frac{3}{4} + \\frac{3}{4 \\cdot 2 \\cdot 2 \\cdot 2} = -\\frac{15}{16}. \\]\n\nWriting the index between brackets instead of the cosine of the corresponding argument, in order to save space.\n\n\\[ -\\int \\frac{dR}{d\\lambda} dt \\]\n\n\\[ = \\frac{3}{4} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{5}{2} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 [1] - \\frac{9}{2} \\left\\{ 1 - \\frac{13}{24} e^2 - \\frac{5}{2} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e [3] \\]\n\n\\[ + \\frac{1}{2} \\left\\{ 1 - \\frac{19}{8} e^2 - \\frac{5}{2} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e [3] + \\frac{21}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{123}{56} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e_i [6] \\]\n\n\\[ + \\frac{3}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{e_i^2}{8} - \\frac{\\gamma^2}{2} \\right\\} m^2 e_i [7] - \\frac{15}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{\\gamma^2}{2} \\right\\} m e^2 [9] \\]\n\n\\[ + \\frac{3}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{5}{2} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e^2 [10] - \\frac{63}{4} \\left\\{ 1 - \\frac{91}{128} e^2 - \\frac{123}{56} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e e_i [12] \\]\n\n\\[ - \\frac{1}{4} \\left\\{ 1 - \\frac{19}{8} e^2 - \\frac{e_i^2}{8} - \\frac{\\gamma^2}{2} \\right\\} m^2 e e_i [13] + \\frac{9}{4} \\left\\{ 1 - \\frac{13}{24} e^2 - \\frac{e_i^2}{8} - \\frac{\\gamma^2}{2} \\right\\} m^2 e e_i [15] \\]\n\n\\[ + \\frac{7}{4} \\left\\{ 1 - \\frac{19}{8} e^2 - \\frac{123}{56} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e_i^2 [16] + \\frac{51}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{115}{51} e_i^2 - \\frac{\\gamma^2}{2} \\right\\} m^2 e_i^2 [18] \\]\n\n\\[ + \\frac{5}{16} m^3 - \\frac{35}{8} m e^2 e_i [24] - \\frac{3}{16} m^2 e^2 e_i [25] + \\frac{15}{8} m e^2 e_i [27] + \\frac{21}{16} m^2 e^2 e_i [28] \\]\n\n\\[ - \\frac{153}{4} m^2 e e_i^2 [30] + \\frac{17}{4} m^2 e e_i^2 [34] + \\frac{845}{64} m^2 e_i^3 [36] + \\frac{m^3 e_i^3}{64} [37] + \\frac{3}{64} m^2 e^4 [39] \\]\n\n\\[ + \\frac{9}{32} m^2 e^4 [40] + \\frac{49}{32} m^2 e^3 e_i [42] + \\frac{25}{32} m^2 e^3 e_i [43] - \\frac{7}{32} m^2 e^3 e_i [45] \\]\n\n\\[ + \\frac{35}{32} m^2 e^3 e_i [46] - \\frac{255}{32} m^2 e^2 e_i^2 [48] + \\frac{51}{16} m^2 e^2 e_i^2 [52] - \\frac{2535}{64} m^2 e e_i^3 [54] \\]\n\n\\[ - \\frac{1}{96} m^2 e e_i^3 [55] - \\frac{3}{32} m^2 e_i^3 [57] - \\frac{15}{32} m^2 e_i^3 [58] + \\frac{2453}{128} m^2 e_i^4 [60] - \\frac{741}{128} m^2 e_i^4 [70] \\]\n\nI have verified some of the terms in the expression for the reciprocal of the radius vector given by M. Planck, which depend on \\( \\frac{a^3}{a_i^4} \\), and arise from the second portion of \\( R \\); very few, however, of these can be obtained without a further development of the disturbing function, in consequence particularly of the term\n\n\\[ \\frac{5}{4} \\frac{a}{a_i} e_i \\cos (\\tau + \\xi). \\]\nwhich is independent of \\( m \\). In consequence of this term, all the terms in my development of \\( R \\), of which the arguments are any combination of the quantity \\( \\tau + \\xi \\), are incomplete.\n\nWhen the terms depending on \\( e^2 \\), and those depending on the square of the disturbing force, are neglected, the inequalities of longitude are given by the equation\n\n\\[\n\\frac{d \\lambda}{dt} = \\frac{h}{r^3} - \\frac{1}{r^3} \\int \\frac{dR}{d\\lambda} dt.\n\\]\n\nI find from this equation in \\( \\lambda_4 \\) the term\n\n\\[\n\\frac{175}{32} m e^2 = \\left\\{ \\frac{405}{32} m e^2 + \\frac{15}{4} m e^2 + \\frac{15}{8} m e^2 \\right\\} \\frac{1}{3},\n\\]\n\ninstead of \\( \\frac{195}{32} m^2 \\), according to M. Planck*.\n\n\\( \\lambda_6 \\) contains the term \\( -\\frac{1953}{128} m^2 e^2 \\) instead of \\( -\\frac{1253}{128} m^2 e^2 \\), for\n\n\\[\nR \\ldots \\ldots \\ldots + \\frac{369}{64} e^3 \\cos (2 \\tau - \\xi),\n\\]\n\n\\[\n\\frac{a}{r} \\ldots \\ldots \\ldots - \\frac{123}{16} m^2 e^3 \\cos (2 \\tau - \\xi),\n\\]\n\n\\[\n\\lambda_6 \\ldots \\ldots \\ldots \\left\\{ -\\frac{123}{8} m^2 - \\frac{369}{64} m^2 \\right\\} e^2 = -\\frac{1353}{64} m^2 e^2.\n\\]\n\nIf the numerical coefficient of the corresponding term in the quantity \\( -\\int \\frac{dR}{d\\lambda} dt \\), be called \\( R \\), then I find\n\n\\[\n\\lambda_{22} = \\left\\{ 2 r_{22} + r_{10} + r_{11} + \\frac{9}{8} r_1 + R_{22} + R_{10} + \\frac{5}{4} R_4 + \\frac{13}{8} R_1 \\right\\} \\frac{1}{(2 - 2m + 3c)} \\\\\n= \\left\\{ \\frac{2125}{192} m^2 + \\frac{7}{2} m^2 + \\frac{33}{16} m^2 + \\frac{9}{8} m^2 + \\frac{5}{16} m^2 + \\frac{3}{8} m^2 + \\frac{5}{8} m^2 + \\frac{13}{8} m^2 \\right\\} \\frac{1}{5} \\\\\n= \\frac{779}{192} m^2.\n\\]\n\nM. Planck has \\( \\frac{1093}{64} m^2 \\), and for the numerical value of the coefficient converted into sexagesimal seconds \\( 3''309 \\). M. Damoiseau has \\( 1''27 \\); I obtain \\( '77'' \\).\n\nI find\n\n\\[\n\\lambda_{25} = \\left\\{ 2 r_{25} + r_{13} + r_7 + R_{25} + R_{13} + \\frac{5}{4} R_7 \\right\\} \\frac{1}{(2 - m + 2c)} \\\\\n= \\left\\{ -\\frac{7}{2} m^2 - \\frac{33}{32} m^2 - \\frac{m^3}{2} - \\frac{3}{16} m^2 - \\frac{1}{4} m^2 + \\frac{5}{4} m^2 \\right\\} \\frac{1}{4} = -\\frac{5}{4} m^2.\n\\]\n\nM. Planck has \\( -\\frac{95}{64} m^2 \\), and for the numerical value of the coefficient converted into sexagesimal seconds \\( -'087 \\). M. Damoiseau has \\( -'19 \\); I obtain \\( -'073 \\).\n\n* The figures are indistinct in the copy before me.\nI find\n\n\\[ \\lambda_{28} = \\left\\{ \\frac{49}{2} m^2 + \\frac{231}{32} m^2 + \\frac{7}{2} m^2 + \\frac{91}{16} m^2 + \\frac{9}{4} m^2 + \\frac{5}{4} \\cdot \\frac{21}{8} m^2 \\right\\} \\frac{1}{4} = \\frac{673}{64} m^2. \\]\n\nM. Plana has \\( \\frac{665}{64} m^2 \\), and for the numerical value of the coefficient converted into sexagesimal seconds \\( \" \\cdot 607 \\). M. Damoiseau has \\( \" \\cdot 90 \\); I obtain \\( \" \\cdot 615 \\).\n\nI find\n\n\\[ \\lambda_{49} = \\left\\{ 2 r_{49} + r_{31} + r_{19} + R_{49} + R_{31} + \\frac{5}{4} R_{19} \\right\\} \\frac{1}{(2 + 2c)} = 0; \\]\n\nM. Plana has \\( -\\frac{885}{64} m \\).\n\nI find\n\n\\[ \\lambda_{59} = 2 r_{59} \\times \\frac{1}{4m} = \\frac{591}{16} m, \\quad \\text{M. Plana has } -\\frac{77}{32} m; \\]\n\nI find\n\n\\[ \\lambda_{60} = \\left\\{ 2 r_{60} + R_{60} \\right\\} \\frac{1}{(2 - 6m)} = -\\frac{55}{8} m. \\]\n\nThese discordances will appear very trifling, considering the nature of the calculations; and it is by no means impossible, after all, that M. Plana may be right, and that the mistake may be with me, notwithstanding all the pains I have used.\n\nBefore the terms in the longitude can be arrived at which depend on \\( \\gamma^2 \\), it is necessary to obtain the expression for the tangent of the latitude \\( s \\): this may be done by means of the equations\n\n\\[ \\frac{d^2 z}{dt^2} + \\frac{\\mu z}{r^3} + \\frac{m_z}{r^3} + \\frac{3 m_z r' r \\cos(\\lambda' - \\lambda)}{r^5} = 0 \\]\n\n\\[ z = \\frac{rs}{\\sqrt{1 + s^2}} \\quad \\frac{z}{r} = s - \\frac{s^3}{2} + \\frac{3}{8} s^5 - &c. \\]\n\n\\[ s = \\frac{z}{r} + \\frac{z^3}{2r^3} - &c. \\]\n\n\\[ \\tan^{-1}s = s - \\frac{s^3}{3} + \\frac{s^5}{5} - &c. \\]\n\nIt is, however, more convenient in the determination of \\( s \\), to adhere to the method of Clairaut, that is, to the method adopted by M. Plana, notwithstanding the difficulties which occur in that method, and to which I have before alluded.\n\nThe following is the differential equation employed.\n\n\\[ \\left\\{ \\frac{d^2 s}{d\\lambda'^2} + s \\right\\} \\left\\{ 1 - \\frac{2}{h^2} \\int r^2 \\frac{dR}{d\\lambda'} d\\lambda' \\right\\} \\]\n\n\\[ + \\frac{r'^2}{h^2} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - r' s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda'} \\right) \\frac{ds}{d\\lambda'} \\right\\} = 0 \\]\n\nSubstituting in this equation in the terms multiplied by \\( m \\), for \\( s \\), \\( r \\sin(g'\\lambda' - \\nu) \\), and for \\( \\frac{ds}{d\\lambda'} = r \\cos(g'\\lambda' - \\nu) \\), neglecting the square of the disturbing force and the cube of \\( s \\), I obtain\n\\[\n\\frac{d^2 s}{d \\lambda^2} + s - \\frac{3 r^4}{2 h^2 r_i^3} \\left\\{ \\sin (2 \\lambda - 2 \\lambda_i - g \\lambda + \\nu) - \\sin (g \\lambda - \\nu) \\right\\}\n- \\frac{r^5}{8 h^2 r_i^4} \\left\\{ \\frac{21}{2} \\cos (\\lambda - \\lambda_i - g \\lambda + \\nu) + \\frac{15}{2} \\cos (\\lambda - \\lambda_i + g \\lambda - \\nu)\n+ 9 \\cos (3 \\lambda - \\lambda_i - g \\lambda + \\nu) + 6 \\cos (3 \\lambda - 3 \\lambda_i + g \\lambda - \\nu) \\right\\}\n\\]\n\nThe simplest method of substituting for \\( \\lambda_i \\) in terms of \\( \\lambda \\) seems to me to be by first obtaining expressions for \\( \\cos \\lambda_i, \\sin \\lambda_i, \\cos 2 \\lambda_i, \\sin 2 \\lambda_i \\), in terms of \\( n_t \\). Having obtained these expressions, they may be reduced to terms of \\( \\lambda \\) by Lagrange's theorem*; but when the higher powers of \\( m \\) are neglected, it is sufficient to write \\( m \\lambda \\) instead of \\( n_t \\).\n\n\\[\n\\sin 2 \\lambda_i = (1 - 4 e_i^2) \\sin 2 n_t t - 2 e_i \\sin (2 n_t t - \\xi_i) + 2 e_i \\sin (2 n_t t + \\xi_i)\n+ \\frac{3}{4} e_i^2 \\sin (2 n_t t - 2 \\xi_i) + \\frac{13}{4} e_i^2 \\sin (2 n_t t + 2 \\xi_i) + \\&c.\n\\]\n\n\\[\n\\cos 2 \\lambda_i = (1 - 4 e_i^2) \\cos 2 n_t t - 2 e_i \\cos (2 n_t t - \\xi_i) + 2 e_i \\cos (2 n_t t + \\xi_i)\n+ \\frac{3}{4} e_i^2 \\cos (2 n_t t - 2 \\xi_i) + \\frac{13}{4} e_i^2 \\cos (2 n_t t + 2 \\xi_i) + \\&c.\n\\]\n\nGreat facility results in the following substitutions in consequence of the coefficients being alike in the corresponding arguments of the expressions \\( \\sin \\lambda_i, \\cos \\lambda_i; \\sin 2 \\lambda_i, \\cos 2 \\lambda_i, \\&c.; \\) so\n\n\\[\n\\sin (2 \\lambda - 2 \\lambda_i - g \\lambda) = (1 - 4 e^2) \\{ \\sin (2 \\lambda - g \\lambda) \\cos 2 m \\lambda - \\cos (2 \\lambda - g \\lambda) \\sin 2 m \\lambda \\}\n- 2 e_i \\{ \\sin (2 \\lambda - g \\lambda) \\cos (2 m \\lambda - c_i m \\lambda) - \\cos (2 \\lambda - g \\lambda) \\sin (2 m \\lambda - c_i m \\lambda) \\}\n+ 2 e_i \\{ \\sin (2 \\lambda - g \\lambda) \\cos (2 m \\lambda + c_i m \\lambda) - \\cos (2 \\lambda - g \\lambda) \\sin (2 m \\lambda + c_i m \\lambda) \\}\n\\]\n\n* By Lagrange's theorem, if\n\n\\[\nu = \\theta - f \\theta\n\\]\n\n\\[\n\\theta = u + fu + \\frac{d^2 (fu)^3}{2 du^2} + \\frac{d^3 (fu)^5}{2 \\cdot 3 du^3} + \\&c.\n\\]\n\nSo\n\n\\[\na_i^3 r_i^3 = 1 + \\frac{3}{2} e_i^2 \\left( 1 + \\frac{5}{4} e_i^2 \\right) + 3 e_i \\left( 1 + \\frac{9}{8} e_i^2 \\right) \\cos \\xi_i\n+ \\frac{9}{2} e_i^2 \\left( 1 + \\frac{7}{9} e_i^2 \\right) \\cos 2 \\xi_i + \\frac{53}{8} e_i^3 \\cos 3 \\xi_i + \\frac{77}{8} e_i^4 \\cos 4 \\xi_i\n\\]\n\n\\[\nm \\lambda = n_t + 2 m e \\left( 1 - \\frac{e^2}{8} \\right) \\sin \\xi + \\frac{5}{4} m e^3 \\left( 1 - \\frac{11}{30} e^2 \\right) \\sin 2 \\xi\n+ \\frac{13}{12} m e^3 \\sin \\xi + \\frac{103}{96} m e^4 \\sin 4 \\xi + \\&c.\n\\]\n\nHence evidently\n\n\\[\na_i^3 r_i^3 = 1 + \\frac{3}{2} e_i^2 \\left( 1 + \\frac{5}{4} e_i^2 \\right) + 3 e_i \\left( 1 + \\frac{9}{8} e_i^2 \\right) \\cos (c_i m \\lambda)\n+ \\frac{9}{2} e_i^2 \\left( 1 + \\frac{7}{9} e_i^2 \\right) \\cos (2 c_i m \\lambda) + \\frac{53}{8} e_i^3 \\cos (3 c_i m \\lambda) + \\frac{77}{9} e_i^4 \\cos (4 c_i m \\lambda)\n\\]\n\n+ terms multiplied by \\( m, c_i \\) may be considered as equal to unity.\n\\[\n\\frac{3}{4} e^2 \\left\\{ \\sin (2\\lambda - g\\lambda) \\cos (2m\\lambda - 2c,m\\lambda) - \\cos (2\\lambda - g\\lambda) \\sin (2m\\lambda - 2c,m\\lambda) \\right\\}\n\\]\n\\[\n+ \\frac{13}{4} e^2 \\left\\{ \\sin (2\\lambda - g\\lambda) \\cos (2m\\lambda + 2c,m\\lambda) - \\cos (2\\lambda - g\\lambda) \\cos (2m\\lambda + 2c,m\\lambda) \\right\\}\n\\]\n\\[\n= (1 - 4e^2) \\sin (2\\lambda - 2m\\lambda - g\\lambda)\n\\]\n\\[\n- 2e_i \\sin (2\\lambda - 2m\\lambda + c_i m\\lambda - g\\lambda) + 2e_i \\sin (2\\lambda - 2m\\lambda - c_i m\\lambda - g\\lambda)\n\\]\n\\[\n+ \\frac{3}{4} e^2 \\sin (2\\lambda - 2m\\lambda + 2c_i m\\lambda - g\\lambda) + \\frac{13}{4} e^2 \\sin (2\\lambda - 2m\\lambda - 2c_i m\\lambda - g\\lambda)\n\\]\n\n\\[\n\\frac{m_i r^4}{h^2 r_i^3} = m^2 \\left\\{ 1 + \\frac{2}{3} e^2 + 3e_i \\cos c_i m\\lambda + \\frac{9}{2} e^2 \\cos 2c_i m\\lambda \\right\\}\n\\]\n\\[\n\\{1 + 2e^2 - 4e \\cos c\\lambda + 5e^2 \\cos 2c\\lambda\\}\n\\]\n\n\\[\n\\frac{3}{2} \\frac{m_i r^4}{h^2 r_i^3} = m^2 \\left\\{ \\frac{3}{2} + \\frac{9}{4} e^2 + 3e_i \\cos c_i m\\lambda + \\frac{9}{2} e_i \\cos c_i m\\lambda + \\frac{15}{2} e^2 \\cos 2c_i m\\lambda\n\\]\n\\[\n- 9e_i \\cos (c_i m\\lambda) - 9e_i \\cos (c_i m\\lambda + c_i m\\lambda) + \\frac{27}{4} e^2 \\cos 2c_i m\\lambda \\right\\}\n\\]\n\nAll the terms which I have verified in the expression for the latitude in terms of the true longitude agree with those given by M. Planá.\n\nIf \\( \\lambda = nt + \\alpha \\), then by Taylor's theorem,\n\n\\[\ns = (s) + \\left( \\frac{ds}{d\\lambda} \\right) \\alpha + \\left( \\frac{d^2 s}{d\\lambda^2} \\right) \\frac{\\alpha^2}{2} + \\left( \\frac{d^3 s}{d\\lambda^3} \\right) \\alpha^3 + &c.\n\\]\n\n\\( (s) \\) being the quantity arising from the substitution of \\( nt \\) for \\( \\lambda \\), in the expression for \\( s \\) in terms of \\( \\lambda \\). In this manner I found the same terms as those given by M. Planá, except \\(-\\frac{15}{32} m e^2 \\gamma \\sin (2\\tau - 2\\xi + \\eta)\\) instead of \\(-\\frac{15}{64} m e^2 \\gamma \\sin (2\\tau - 2\\xi + \\eta)\\), and \\(-\\frac{3}{8} m e e_i \\gamma \\sin (2\\tau + \\xi + \\xi_i - \\eta)\\) instead of \\(\\frac{1}{2} m e e_i \\gamma \\sin (2\\tau + \\xi + \\xi_i - \\eta)\\).\n\nI next obtained \\( s^2 \\), in order to procure the terms in the longitude depending on \\( \\gamma^2 \\). The quantity \\( e \\) in my notation does not accord with that quantity in the work of M. Planá, but with \\( e \\left( 1 - \\frac{\\gamma^2}{4} \\right) \\); so that, in order to arrive at the same figures in some of the terms multiplied by \\( \\gamma^2 \\), this circumstance must be attended to.\n\nI find\n\n\\[\n\\lambda_{70} = \\left\\{ -\\frac{39}{16} m^2 - \\frac{11}{16} m^2 - \\frac{33}{32} m^2 - m^2 - \\frac{1}{4} m^2 - \\frac{3}{8} m^2 \\right\\} \\frac{1}{5} = -\\frac{37}{32} m^2\n\\]\n\ninstead of \\(-\\frac{39}{32} m^2\\), according to M. Planá; and I find\n\n\\[\n\\lambda_{75} = -\\left\\{ \\frac{35}{8} m + \\frac{7}{8} m - \\frac{35}{16} m - \\frac{7}{8} m - \\frac{7}{16} m \\right\\} = -\\frac{7}{4} m\n\\]\n\ninstead of \\(\\frac{7}{8} m\\), according to M. Planá.\nThese are the only discrepancies which I have noted in the terms multiplied by $\\gamma^2$ which I have examined. In making use of the development of $R$ before alluded to,\n\n$$1 - \\frac{\\gamma^2}{2} + \\frac{7}{16} \\gamma^4$$\n\nis to be substituted for $\\cos^4 \\frac{\\gamma}{2}$.\n\nThe relation between the constants $h$ and $a$ is to be obtained from the equations\n\n$$r^4 \\frac{d \\lambda^2}{dt^2} = h^2 - 2 \\int r^2 \\frac{dR}{d\\lambda} d\\lambda$$\n\n$$\\frac{r^2 d \\lambda^2 + dr^2}{dt^2} - \\frac{2 \\mu}{r} + \\frac{\\mu}{a} + 2 \\int dR = 0;$$\n\nwhere, for simplicity, I have neglected at present quantities depending on $\\gamma^2$.\n\n$$d \\cdot \\frac{1}{r} = - \\frac{dr}{r^2}$$\n\n$$d r^2 = r^4 \\left( \\frac{1}{r} \\right)^2 = \\left( \\frac{1}{r} \\right)^{-4} \\left( \\frac{d}{r} \\right)^2 = a^2 (1 - 4 r_0) \\left\\{ \\frac{e^2 c^2 n^2}{2} + \\frac{e^2 r_3^2 n^2}{2} \\right\\}$$\n\n$2 \\int dR$ can give no constant term *; therefore, considering the constant part only of the equation above, since the coefficients corresponding to all the arguments must be separately identical,\n\n$$h^2 \\left\\{ 1 + 2 r_0 + \\frac{e^2}{2} + \\frac{e^2 r_3^2}{2} \\right\\} = a \\mu \\left\\{ 1 + 2 r_0 - (1 - 4 r_0) \\right\\} \\frac{e^2 c^2}{2} + \\frac{e^2 r_3^2}{2}$$\n\n$r_4^2$, being already multiplied by $m^4$, is not to be taken into account in this approximation.\n\n$$c^2 \\{ 1 - 3 r_0 \\} - 1 + 2 m^2 = 0$$\n\n$$c^2 = 1 - \\frac{7}{2} m^2$$\n\n$$h^2 \\left\\{ 1 - m^2 + \\frac{e^2}{2} + \\frac{33}{128} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\} = a \\mu \\left\\{ 1 - \\frac{e^2}{2} - m^2 - \\frac{321}{128} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\}$$\n\n$$h^2 = a \\mu \\left\\{ 1 - \\frac{e^2}{2} - m^2 - \\frac{321}{128} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\}$$\n\n$$\\left\\{ 1 + m^2 - \\frac{e^2}{2} - \\frac{33}{128} m^2 e^2 + \\frac{3}{2} m^2 e_i^2 - m^2 e^2 \\right\\}$$\n\n$$= a \\mu \\left\\{ 1 - e^2 - \\frac{241}{64} m^2 e^2 \\right\\}$$\n\n$$h = a \\mu \\left\\{ 1 - \\frac{e^2}{2} - \\frac{241}{128} m^2 e^2 \\right\\}$$\n\n$$\\frac{d \\lambda}{dt} = \\frac{h}{r^3} - \\frac{1}{r^3} \\int \\frac{dR}{d\\lambda} dt.$$\n\n* This seems at variance with the equation of M. Planck, vol. i. p. 122,\n\n$$- \\frac{2a}{\\sigma} \\int d' \\Omega = - \\frac{m^2}{2} \\left\\{ 1 + \\frac{1}{2} \\gamma^2 + \\frac{3}{2} e^2 + \\frac{3}{2} e^2 \\gamma^2 \\cos (2 \\varphi - 2 \\theta) + &c. \\right\\}$$\n\nwhich equation I am unable to understand.\n\nMDCCCXXXIV.\nThe second term on the right-hand side of this equation gives no constant quantity: hence\n\n\\[ \\frac{h}{r^3} = \\sqrt{\\frac{\\mu}{a^3}} \\left\\{ 1 - \\frac{e^2}{2} - \\frac{241}{128} m^2 e^2 \\right\\} \\left\\{ 1 - m^2 + \\frac{e^2}{2} + \\frac{33}{128} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\} \\]\n\n\\[ \\lambda = \\sqrt{\\frac{\\mu}{a^3}} \\left\\{ 1 - m^2 - \\frac{9}{8} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\} t + &c. \\]\n\nIf the constant quantity which multiplies \\( t \\) be called\n\n\\[ n = \\sqrt{\\frac{\\mu}{a^3}} \\]\n\n\\[ \\sqrt{\\frac{\\mu}{a^3}} \\left\\{ 1 - m^2 - \\frac{9}{8} m^2 e^2 - \\frac{3}{2} m^2 e_i^2 \\right\\} = \\sqrt{\\frac{\\mu}{a^3}} \\]\n\n\\[ a = a \\left\\{ 1 + \\frac{2}{3} m^2 + \\frac{3}{4} m^2 e^2 + m^2 e_i^2 \\right\\} \\]\n\n\\[ \\frac{a}{r} = \\frac{a \\left\\{ 1 + \\frac{2}{3} m^2 + \\frac{3}{4} m^2 e^2 + m^2 e_i^2 \\right\\}}{r} = 1 - \\frac{m^2}{2} - \\frac{3}{4} m^2 e^2 - \\frac{3}{4} m^2 e_i^2. \\]\n\nReverting to the equation\n\n\\[ \\frac{d\\lambda}{dt} = \\frac{h}{r^2} - \\frac{1}{r^3} \\int \\frac{dR}{d\\lambda} dt. \\]\n\nThe second term of this equation gives no term multiplied by \\( \\cos\\xi \\); therefore,\n\n\\[ d\\lambda = ndt + \\frac{2h}{a^2} (1 + r_0) e \\cos\\xi dt \\]\n\n\\[ = ndt + 2n \\left( 1 - \\frac{m^2}{2} \\right) (1 + m^2) e \\cos\\xi dt \\]\n\n\\[ c = 1 - \\frac{7}{4} m^2 \\quad cn = cn \\quad c = 1 - \\frac{3}{4} m^2 \\]\n\n\\[ \\lambda = nt + 2 \\left( 1 + \\frac{3}{4} m^2 + \\frac{m^2}{2} \\right) e \\sin\\xi \\]\n\n\\[ \\lambda = nt + 2 \\left( 1 + \\frac{5}{4} m^2 \\right) e \\sin\\xi + &c. \\]\n\n\\[ \\frac{a}{r} = 1 + \\frac{m^2}{6} + \\frac{m^2 e_i^2}{4} + e \\left( 1 + \\frac{2}{3} m^2 \\right) \\cos\\xi + &c. \\]\n\nThese are the expressions for \\( \\lambda \\) and for \\( \\frac{1}{r} \\), when the quantity \\( e \\) is retained; but if the coefficient of \\( \\sin\\xi \\), in the expression for \\( \\lambda \\), be called \\( 2e \\), after the manner adopted for the planets in the first volume of the Mécanique Céleste, so that\n\n\\[ \\lambda = nt + 2e \\sin\\xi + &c., \\quad \\text{then} \\]\n\n\\[ \\frac{a}{r} = 1 + \\frac{m^2}{6} + \\frac{m^2 e_i^2}{4} + e \\left( 1 - \\frac{7}{12} m^2 \\right) \\cos\\xi + &c. \\]\nAs the preceding results do not quite agree with those of M. Planck, I shall endeavour to show how they may be obtained from the same equations which he employs.\n\n\\[\n\\frac{d^2}{d\\lambda^2} \\left( \\frac{1}{r} \\right) + \\frac{1}{r} - \\frac{1}{h^2} \\left\\{ 1 + 2 \\int r^2 \\frac{dR}{d\\lambda} d\\lambda \\right\\} - \\frac{r}{h^2} \\left\\{ r \\frac{dR}{dr} - \\frac{1}{r} \\frac{dR}{d\\lambda} \\cdot \\frac{dr}{d\\lambda} \\right\\} = 0\n\\]\n\n\\[ R = -\\frac{r^3}{4r_i^3} \\left\\{ 1 + 3 \\cos (2\\lambda - 2\\lambda_i) \\right\\}\n\\]\n\nIn order to integrate the preceding differential equation, let\n\n\\[\n\\frac{h^2}{r} = \\mu \\left\\{ 1 + r_0 + r_1 \\cos (2\\lambda - 2m\\lambda) + e \\cos (c\\lambda - \\pi) + er_3 \\cos (2\\lambda - 2m\\lambda - c\\lambda + \\pi) \\right\\}\n\\]\n\nThe letters \\( r_0, r_1, \\) &c., being now used in a somewhat different sense to heretofore, having now reference to the expression for \\( \\frac{1}{r} \\) in terms of the true longitude.\n\n\\[\n\\frac{1}{r} = \\frac{h^6}{\\mu^3} (1 + 3e^2 + &c.)\n\\]\n\nNeglecting \\( e_i^2 \\),\n\n\\[\nr^0 + \\frac{h^6}{2a_i^3 \\mu} (1 + 3e^2) = 0.\n\\]\n\nSubstituting in this equation for \\( h \\) its elliptic value which is allowable, \\( r_0 = -\\frac{m^2}{2} \\), also \\( r_3 = \\frac{15}{8} m \\),\n\n\\[\n\\left( \\frac{1}{r} \\right)^{-2} = \\frac{h^4}{\\mu^3} \\left\\{ 1 + m^2 + \\frac{225}{128} \\times 3m^2e^2 + \\frac{3e^3}{2} + 3m^2e^2 \\right\\}\n\\]\n\n\\[\n= \\frac{h^4}{\\mu^3} \\left\\{ 1 + \\frac{3e^2}{2} + m^2 + \\frac{1059}{128} m^2e^2 \\right\\}\n\\]\n\n\\[\n\\frac{dt}{d\\lambda} = \\frac{r^3}{h} = \\frac{h^3}{\\mu^3} \\left\\{ 1 + \\frac{3e^2}{2} + m^2 + \\frac{1059}{128} m^2e^2 \\right\\} = \\sqrt{\\frac{a^3}{\\mu}}\n\\]\n\n\\[\nh^2 = a \\left\\{ 1 + \\frac{3e^2}{2} + m^2 + \\frac{1059}{128} m^2e^2 \\right\\}^{-\\frac{1}{2}}\n\\]\n\n\\[\n= a \\left\\{ 1 - e^2 - \\frac{2}{3} m^2 - \\frac{739}{192} m^2e^2 \\right\\}\n\\]\n\n\\[\n= a \\left\\{ 1 - e^2 - \\frac{2}{3} m^2 - \\frac{739}{192} m^2e^2 \\right\\} \\left\\{ 1 + \\frac{2m^2}{3} + \\frac{3}{4} m^2e^2 \\right\\}\n\\]\n\n\\[\n= a \\left\\{ 1 - e^2 - \\frac{241}{64} m^2e^2 \\right\\} \\text{as before.}\n\\]\n\\[\n\\frac{h^2}{r} = \\frac{a \\left\\{ 1 - e^2 - \\frac{2}{3} m^2 - \\frac{739}{192} m^2 e^2 \\right\\}}{r}\n\\]\n\n\\[= 1 - \\frac{m^2}{2} + e \\cos (c \\lambda - \\varpi) + &c.\\]\n\n\\[\\frac{a}{r} = 1 + e^2 + \\frac{m^2}{6} + \\frac{899}{192} m^2 e^2 + e \\left(1 + \\frac{2}{3} m^2\\right) \\cos (c \\lambda - \\varpi) + &c.\\]\n\nMy letter \\(a\\) appears to correspond with \\(a (1 + p)\\), or \\(a\\) in M. Planà's notation.\n\n\\[\n\\begin{align*}\nn &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots n \\\\\nc &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots c \\\\\n\\mu &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\sigma \\\\\nR &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\Omega \\\\\n\\frac{1}{r} &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots u \\\\\n\\lambda' &\\quad \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots \\ldots v.\n\\end{align*}\n\\]\n\nPutting \\(e \\left(1 - \\frac{m^2}{2}\\right)\\) instead of \\(e\\) in the various expressions found above,\n\n\\[\\lambda = nt + 2 \\left(1 + \\frac{3}{4} m^2\\right) e \\sin \\xi,\\]\n\nwhich then so far agrees with the expression of M. Planà*, and I then find\n\n\\[\\frac{a}{r} = 1 + \\frac{m^2}{6} + \\frac{m^2 e^2}{4} + e \\left(1 + \\frac{m^2}{6}\\right) \\cos \\xi + &c.\\]\n\n\\[\\frac{a}{r} = 1 + e^2 + \\frac{m^2}{6} + \\frac{707}{192} m^2 e^2 + e \\left(1 + \\frac{m^2}{6}\\right) \\cos (c \\lambda - \\varpi).\\]\n\nM. Planà has\n\n\\[\\frac{a}{r} = 1 + \\frac{m^2}{6} - \\frac{45}{16} m^2 e^2 + \\frac{m^2 e^2}{4} + e \\left(1 + \\frac{m^2}{6}\\right) \\cos \\xi + &c.\\]\n\n\\[\\frac{a}{r} = 1 + e^2 + \\frac{m^2}{6} + \\frac{167}{192} m^2 e^2 + e \\left(1 + \\frac{m^2}{6}\\right) \\cos (c \\lambda - \\varpi) + &c.\\]\n\nwhich equations do not agree with those I have found. I am, however, well aware how difficult it is to escape error in these inquiries, and wish to be understood as not offering any of the results contained in this paper too confidently.\n\nWherever I presumed to have arrived at figures differing from those of M. Planà, I verified afresh all the steps of the process contained in previous papers, particularly the corresponding term in the development of \\(R\\). Thus I have found by means of the expressions given §, that three times the numerical coefficient of \\(e^3 \\cos (2 \\tau - \\xi)\\).\n\n---\n\n* Vol. i. p. 574. † Vol. i. p. 664. ‡ Vol. i. p. 636.\n\n§ Philosophical Transactions, 1832, p. 601.\nin the development of $R$,\n\n$$= \\frac{63}{16} + \\frac{81}{64} + \\frac{45}{16} + \\frac{27}{32} - \\frac{153}{16} + \\frac{18}{32} + \\frac{30}{8} + \\frac{30}{32} + \\frac{204}{16} = \\frac{1107}{64}. $$\n\nThe coefficient in question $= \\frac{369}{64} = \\frac{21}{8} \\times \\frac{153}{56}$.\n\nBut I have found that the following corrections are required. For\n\n$$+ \\frac{3}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - 4 e_i^2 \\right\\} \\cos^4 \\frac{1}{2} \\frac{a^2}{a_i^3} e_i \\cos (2 \\tau + \\xi) $$\n\nread\n\n$$+ \\frac{3}{8} \\left\\{ 1 - \\frac{5}{2} e^2 - \\frac{e_i^2}{8} \\right\\} \\cos^4 \\frac{1}{2} \\frac{a^2}{a_i^3} e_i \\cos (2 \\tau + \\xi). $$\n\nand for\n\n$$- \\frac{51}{64} \\frac{a^2}{a_i^3} e_i^2 \\gamma^2 \\cos (2 \\xi - 2 \\eta) - \\frac{45}{64} \\frac{a^2}{a_i^3} e_i^2 \\gamma^2 \\cos (2 \\xi + 2 \\eta) - \\frac{195}{64} \\frac{a^2}{a_i^3} e_i^2 \\gamma^2 \\cos (2 \\tau - 2 \\xi - 2 \\eta) $$\n\nread\n\n$$- \\frac{27}{16} \\frac{a^2}{a_i^3} e_i^2 \\gamma^2 \\cos (2 \\xi - 2 \\eta) - \\frac{27}{16} e_i^2 \\gamma^2 \\cos (2 \\xi + 2 \\eta) - \\frac{51}{16} \\frac{a^2}{a_i^3} e_i^2 \\gamma^2 \\cos (2 \\tau - 2 \\xi - 2 \\eta). $$",
  "source": "olmocr",
  "added": "2026-01-12",
  "created": "2026-01-12",
  "metadata": {
    "Source-File": "/home/jic823/projects/def-jic823/royalsociety/pdfs/108057.pdf",
    "olmocr-version": "0.3.4",
    "pdf-total-pages": 16,
    "total-input-tokens": 25962,
    "total-output-tokens": 16309,
    "total-fallback-pages": 0
  },
  "attributes": {
    "pdf_page_numbers": [
      [
        0,
        0,
        1
      ],
      [
        0,
        2796,
        2
      ],
      [
        2796,
        5567,
        3
      ],
      [
        5567,
        8005,
        4
      ],
      [
        8005,
        10405,
        5
      ],
      [
        10405,
        12324,
        6
      ],
      [
        12324,
        15399,
        7
      ],
      [
        15399,
        17733,
        8
      ],
      [
        17733,
        20193,
        9
      ],
      [
        20193,
        23228,
        10
      ],
      [
        23228,
        26103,
        11
      ],
      [
        26103,
        28443,
        12
      ],
      [
        28443,
        30490,
        13
      ],
      [
        30490,
        32473,
        14
      ],
      [
        32473,
        34770,
        15
      ],
      [
        34770,
        35845,
        16
      ]
    ],
    "primary_language": [
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en"
    ],
    "is_rotation_valid": [
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true
    ],
    "rotation_correction": [
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0
    ],
    "is_table": [
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false
    ],
    "is_diagram": [
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false
    ]
  },
  "jstor_metadata": {
    "identifier": "jstor-108057",
    "title": "On the Theory of the Moon [Continued]",
    "authors": "John William Lubbock",
    "year": 1834,
    "volume": "124",
    "journal": "Philosophical Transactions of the Royal Society of London",
    "page_count": 16,
    "jstor_url": "https://www.jstor.org/stable/108057"
  }
}