{
  "id": "13a88f37783fda85e23b08142ee3e00a05756ef0",
  "text": "III. Researches in Physical Astronomy. By John William Lubbock, Esq., V.P. and Treas. R.S.\n\nRead December 9, 1830.\n\nIn last April I had the honour of presenting to the Society a paper containing expressions for the variations of the elliptic constants in the theory of the motions of the planets. The stability of the solar system is established by means of these expressions, if the planets move in a space absolutely devoid of any resistance*, for it results from their form that however far the approximation be carried, the eccentricity, the major axis, and the tangent of the inclination of the orbit to a fixed plane, contain only periodic inequalities, each of the three other constants, namely, the longitude of the node, the longi-\n\n* When the body moves in a medium which resists according to any power of the velocity, the contrary obtains, the major axis and eccentricity acquiring a term which varies with the time, while the longitude of the perihelion and longitude of the epoch have only periodic inequalities. This results from the equations given in the former part of this paper, Phil. Trans. Part II. 1830, page 340.\n\n\\[\nda = -2c \\left( \\frac{\\mu}{a} \\right)^{n-1} \\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n+1} \\left( \\frac{1-e \\cos v}{n} \\right) dv\n\\]\n\n\\[\nde = -2c \\left( \\frac{\\mu}{a} \\right)^{n-1} \\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n-1} \\left( \\frac{1-e^2}{2} \\right) \\cos d v\n\\]\n\n\\[\ned \\varpi = -2c \\left( \\frac{\\mu}{a} \\right)^{n-1} \\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n-1} \\left( \\frac{\\sqrt{1-e^2}}{n} \\right) \\sin v dv\n\\]\n\n\\[\nd\\varepsilon - d\\varpi = 2c \\left( \\frac{\\mu}{a} \\right)^{n-1} \\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n-1} \\left\\{ \\frac{1-e^2 \\cos v}{e} \\right\\} dv\n\\]\n\n\\[\n\\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n+1} = 1 + (n+1)e \\cos v + (n+1)^2 e^2 \\cos^2 v\n\\]\n\n\\[\n\\left\\{ \\frac{1+e \\cos v}{1-e \\cos v} \\right\\}^{n-1} = 1 + (n-1)e \\cos v + (n-1)^2 e^2 \\cos^2 v\n\\]\ntude of the perihelion, and the longitude of the epoch, contains a term which varies with the time, and hence the line of apsides and the line of nodes revolve continually in space. The stability of the system may therefore be inferred, which would not be the case if the eccentricity, the major axis, or the tangent of the inclination of the orbit to a fixed plane contained a term varying with the time, however slowly.\n\nThe problem of the precession of the equinoxes admits of a similar solution; of the six constants which determine the position of the revolving body, and the axis of instantaneous rotation at any moment, three have only periodic inequalities, while each of the other three has a term which varies with the time. From the manner in which these constants enter into the results, the equilibrium of the system may be inferred to be stable, as in the former case. Of the constants in the latter problem, the mean angular velocity of rotation\n\n\\[ da = -\\frac{2ca}{n} \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ 1 + ne \\cos v + n(n+1)e^2 \\cos^2 v + \\text{&c.} \\right\\} dv \\]\n\n\\[ de = -\\frac{2c}{n} \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ \\cos v + (n+1)e \\cos^2 v + (n+1)^2e^2 \\cos^3 v + \\text{&c.} \\right\\} (1-e^2)dv \\]\n\nneglecting the terms which are periodic\n\n\\[ da = -\\frac{2ca}{n} \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ 1 + \\frac{n(n+1)}{2}e^2 + \\text{&c.} \\right\\} dv \\]\n\n\\[ de = -\\frac{2c}{n} \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ \\frac{(n+1)}{2}e + \\text{&c.} \\right\\} (1-e^2)dv \\]\n\nThe major axis decreases perpetually, the eccentricity diminishes perpetually until it reaches zero, while the perihelion retains the same mean position, and the longitude of the epoch the same mean value. I stated inadvertently in the former part of this paper, p. 340, that the variations of the eccentricity are all periodical.\nmay be considered analogous to the mean motion of a planet, or its major axis; the geographical longitude, and the cosine of the geographical latitude of the pole of the axis of instantaneous rotation, to the longitude of the perihelion and the eccentricity; the longitude of the first point of Aries and the obliquity of the ecliptic, to the longitude of the node and the inclination of the orbit to a fixed plane; and the longitude of a given line in the body revolving, passing through its centre of gravity, to the longitude of the epoch. By the stability of the system I mean that the pole of the axis of rotation has always nearly the same geographical latitude, and that the angular velocity of rotation, and the obliquity of the ecliptic vary within small limits, and periodically. These questions are considered in the paper I now have the honour of submitting to the Society. It remains to investigate the effect which is produced by the action of a resisting medium; in this case the latitude of the pole of the axis of rotation, the obliquity of the ecliptic, and the angular velocity of rotation might vary considerably, although slowly, and the climates undergo a considerable change.\n\nThe co-efficients of the terms in the development of $R$, multiplied by the squares and products of the eccentricities, are susceptible of very great simplification, in consequence of the equations of condition which obtain between the quantities of which the general symbol is $b$. I have now given the development of $R$, as far as the terms depending upon the squares and products of the eccentricities, in its simplest form. See p. 30.\n\nI have also given methods of obtaining the inequalities of the radius vector, of longitude, and of latitude in the planetary theory. The expressions in this paper differ in form from those of Laplace, but their identity may be shown by means of equations of condition which obtain between some of the quantities involved.\n\nI have taken as a numerical example, the calculation of the co-efficients of some of the inequalities in the theory of Jupiter, disturbed by Saturn.\n\nOn the Precession of the Equinoxes.\n\nLet O be the origin of the co-ordinate axes, coinciding with some point in the interior of the mass $M$.\n\nLet $x, y, z$ be the co-ordinates of any element $dm$ parallel to three rectangular\naxes \\( O_x, O_y, O_z \\), fixed in space, \\( x_i, y_i, z_i \\), the co-ordinates of the same element parallel to three other rectangular axes \\( O_{x'}, O_{y'}, O_{z'} \\), fixed in the mass \\( M' \\) and revolving with it. Let the line \\( N_0N' \\) be the intersection of the plane \\( x_iy_i \\) with the plane \\( xy \\),\n\nLet the angle \\( NOx = \\psi, NOx_i = \\phi \\), and the inclination of the plane \\( x_iy_i \\) upon \\( xy = \\theta \\).\n\n\\[\n\\begin{align*}\nx_i &= x (\\cos \\theta \\sin \\psi \\sin \\phi + \\cos \\psi \\cos \\phi) + y (\\cos \\theta \\cos \\psi \\sin \\phi - \\sin \\psi \\cos \\phi) - z \\sin \\theta \\sin \\phi \\\\\ny_i &= x (\\cos \\theta \\sin \\psi \\cos \\phi - \\cos \\psi \\sin \\phi) + y (\\cos \\theta \\cos \\psi \\cos \\phi + \\sin \\psi \\sin \\phi) - z \\sin \\theta \\cos \\phi \\\\\nz_i &= x \\sin \\theta \\sin \\psi + y \\sin \\theta \\cos \\psi + z \\cos \\theta\n\\end{align*}\n\\]\n\nLet \\( X_i, Y_i, Z_i \\) be the accelerating forces which act upon the element \\( dm \\) in the direction of the axes \\( O_{x'}, O_{y'}, O_{z'} \\),\n\n\\[\n\\int (y_i^2 + z_i^2) \\, dm = A, \\quad \\int (x_i^2 + z_i^2) \\, dm = B, \\quad \\int (x_i^2 + y_i^2) \\, dm = C\n\\]\n\n\\[\n\\begin{align*}\np \\, dt &= \\sin \\phi \\sin \\theta \\, d\\psi - \\cos \\phi \\, d\\theta \\\\\nq \\, dt &= \\cos \\phi \\sin \\theta \\, d\\psi + \\sin \\phi \\, d\\theta \\\\\nr \\, dt &= d\\phi - \\cos \\theta \\, d\\psi\n\\end{align*}\n\\]\n\n\\[\nC \\, dr + (B - A) \\, p \\, q \\, dt = d \\, t \\int (x_i \\, Y_i - y_i \\, X_i) \\, dm\n\\]\n\n\\[\nB \\, dq + (A - C) \\, r \\, p \\, dt = d \\, t \\int (z_i \\, X_i - x_i \\, Z_i) \\, dm\n\\]\n\n\\[\nA \\, dp + (C - B) \\, q \\, r \\, dt = d \\, t \\int (y_i \\, Z_i - x_i \\, Y_i) \\, dm\n\\]\n\nIf the axis of instantaneous rotation coincides with the line \\( OI \\) at any instant\n\n\\[\n\\begin{align*}\n\\cos IOx_i &= \\frac{p}{\\sqrt{p^2 + q^2 + r^2}} \\\\\n\\cos IOy_i &= \\frac{q}{\\sqrt{p^2 + q^2 + r^2}} \\\\\n\\cos IOz_i &= \\frac{r}{\\sqrt{p^2 + q^2 + r^2}} \\\\\n\\sin IOz_i &= \\frac{\\sqrt{p^2 + q^2}}{\\sqrt{p^2 + q^2 + r^2}}\n\\end{align*}\n\\]\n\nIf \\( z_iIL \\) be a great circle cutting the plane \\( x_iy_i \\) in \\( L \\),\n\n\\[\n\\cos x_iOL = \\frac{\\cos IOx_i}{\\sin IOz_i} = \\frac{p}{\\sqrt{p^2 + q^2}}\n\\]\n\nIf the accelerating forces \\( X, Y, Z = 0 \\) and \\( B = A \\), the integrals of the preceding equations are\n\n\\[\nr = n, \\quad p = c \\cos \\frac{C - A}{A} (nt + \\nu), \\quad q = c \\sin \\frac{C - A}{A} (nt + \\gamma)\n\\]\nand neglecting $c^2$,\n\n$$\\theta = \\omega + \\frac{c}{n} \\frac{A}{C} \\sin \\left( \\frac{C}{A} n t + \\frac{C - A}{A} \\gamma \\right)$$\n\n$$\\psi = \\psi_0 - \\frac{c}{n \\sin \\omega} \\frac{A}{C} \\cos \\left( \\frac{C}{A} n t + \\frac{C - A}{A} \\gamma \\right)$$\n\n$$\\phi = \\phi_0 + n t - \\frac{c \\cos \\omega}{n \\sin \\omega} \\frac{A}{C} \\cos \\left( \\frac{C}{A} n t + \\frac{C - A}{A} \\gamma \\right)$$\n\n$n$, $c$, $\\gamma$, $\\omega$, $\\psi_0$, and $\\phi_0$ being constants. In the problem of the Precession of the Equinoxes $\\omega$ is the mean obliquity of the ecliptic, $\\psi_0$ is the longitude of the first point of Aries when $t = 0$ reckoned from some fixed line.\n\n$\\sin I O z_i = \\frac{c}{\\sqrt{n^2 + c^2}}$, hence it appears that if a body whose form is that of a figure of revolution be made to revolve, and be acted upon by no extraneous force, the axis of instantaneous rotation revolves about the axis of the figure, the latitude of the former axis remaining constant.\n\nThe angle $x_i O I_i = \\frac{C - A}{A} (n t + \\gamma)$\n\nIf the forces $X_i$, $Y_i$, $Z_i$ arise from the attractions of a distant luminary $M'$, of which the coordinates referred to the axes $O x_i$, $O y_i$, $O z_i$, are $x_i', y_i', z_i'$, the force varying inversely as the square of the distance,\n\n$$X_i = \\frac{M' (x_i' - x_i)}{\\{(x_i - x_i')^2 + (y_i - y_i')^2 + (z_i - z_i')^2\\}^{3/2}}$$\n\n$$Y_i = \\frac{M' (y_i' - y_i)}{\\{(x_i - x_i')^2 + (y_i - y_i')^2 + (z_i - z_i')^2\\}^{3/2}}$$\n\n$$Z_i = \\frac{M' (z_i' - z_i)}{\\{(x_i - x_i')^2 + (y_i - y_i')^2 + (z_i - z_i')^2\\}^{3/2}}$$\n\n$$\\int (x_i Y_i - y_i X_i) \\, dm = M' \\int \\frac{(x_i y_i' - y_i x_i')}{r'^3} \\left\\{ 1 + \\frac{3}{r^2} \\left( x_i x_i' + y_i y_i' + z_i z_i' \\right) + &c. \\right\\} \\, dm$$\n\nBy the properties of the principal axes $\\int x_i y_i \\, dm = 0$, $\\int x_i z_i \\, dm = 0$, $\\int y_i z_i \\, dm = 0$, and by the properties of the centre of gravity $\\int x_i \\, dm = 0$, $\\int y_i \\, dm = 0$, $\\int z_i \\, dm = 0$, whence\n\n$$\\int (x_i Y_i - y_i X_i) \\, dm = \\frac{3 M' x_i' y_i'}{r'^3} \\int (x_i^2 - y_i^2) \\, dm = \\frac{3 M' (B - A)}{r'^3} x_i' y_i' \\, dt$$\n\\[ Cdr + (B - A)pqdt = \\frac{3M'(B - A)}{r^5} x'_i y'_i dt \\]\n\n\\[ Bdq + (A - C)rpdt = \\frac{3M'(A - C)}{r^5} z'_i x'_i dt \\]\n\n\\[ Adp + (C - B)qr dt = \\frac{3M'(C - B)}{r^5} y'_i z'_i dt \\]\n\nSubstituting for \\( x'_i, y'_i, z'_i \\) their values from equations, p. 20, upon the supposition of \\( \\psi = 0 \\),\n\n\\[ Cdr + (B - A)pqdt = \\frac{3M'(B - A)}{2r^5} \\left\\{ \\left( y'\\cos\\theta - z'\\sin\\theta \\right)^2 - x'^2 \\right\\} \\sin 2\\varphi \\]\n\n\\[ + 2x'\\left( y'\\cos\\theta - z'\\sin\\theta \\right) \\cos^2\\varphi \\right\\} \\]\n\n\\[ Bdq + (A - C)rpdt = \\frac{3M'(A - C)}{r^5} \\left\\{ x'\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\cos\\varphi \\right. \\]\n\n\\[ + \\left( y'\\cos\\theta - z'\\sin\\theta \\right)\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\sin\\varphi \\right\\} \\]\n\n\\[ Adp + (C - B)qr dt = \\frac{3M'(C - B)}{r^5} \\left\\{ \\left( y'\\cos\\theta - z'\\sin\\theta \\right)\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\cos\\varphi \\right. \\]\n\n\\[ - x'\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\sin\\varphi \\right\\} \\]\n\nIf\n\n\\[ \\frac{3M'}{r^5} \\left\\{ \\left( y'\\cos\\theta - z'\\sin\\theta \\right)\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\right\\} = P \\]\n\n\\[ \\frac{3M'}{r^5} \\left\\{ x'\\left( y'\\sin\\theta + z'\\cos\\theta \\right) \\right\\} = P' \\]\n\n\\[ Bdq + (A - C)rpdt = (A - C)dt \\left\\{ P'\\cos\\varphi + P\\sin\\varphi \\right\\} \\]\n\n\\[ Adp + (C - B)qr dt = (C - B)dt \\left\\{ P\\cos\\varphi - P'\\sin\\varphi \\right\\} \\]\n\n\\( P \\) and \\( P' \\) may be developed according to sines and cosines of angles increasing proportionally to the time. Let \\( k\\cos(it + \\varepsilon) \\) be any term of \\( P \\), \\( k'\\sin(it + \\varepsilon) \\) the corresponding term of \\( P' \\),\n\n\\[ Bdq + (A - C)rpdt = \\frac{A - C}{2}dt \\left\\{ (k + k')\\sin(\\phi + it + \\varepsilon) + (k - k')\\sin(\\phi - it - \\varepsilon) \\right\\} \\]\n\n\\[ Adp + (C - B)qr dt = \\frac{C - B}{2}dt \\left\\{ (k + k')\\cos(\\phi + it + \\varepsilon) + (k - k')\\cos(\\phi - it - \\varepsilon) \\right\\} \\]\n\nThe equations which were given p. 21, may still be considered as afford-\ning a solution of the problem by making the constants \\( n, c, \\gamma, \\omega, \\psi_0 \\) and \\( \\phi_0 \\) vary,\n\n\\[\n\\begin{align*}\ndn &= 0 \\\\\n\\sin \\frac{C - A}{A} (nt + \\gamma) dc + c \\cos \\frac{C - A}{A} (nt + \\gamma) d\\gamma + \\frac{(A - C)}{A} nc \\cos \\frac{C - A}{A} (nt + \\gamma) dt \\\\\n&= \\frac{A - C}{2A} dt \\left\\{ (k + k') \\sin (\\phi + it + \\varepsilon) + (k - k') \\sin (\\phi - it - \\varepsilon) \\right\\} \\\\\n\\cos \\frac{C - A}{A} (nt + \\gamma) dc - c \\sin \\frac{C - A}{A} (nt + \\gamma) d\\gamma - \\frac{A - C}{A} nc \\sin \\frac{C - A}{A} (nt + \\gamma) dt \\\\\n&= -\\frac{A - C}{2A} dt \\left\\{ (k + k') \\cos (\\phi + it + \\varepsilon) + (k - k') \\cos (\\phi - it - \\varepsilon) \\right\\}\n\\end{align*}\n\\]\n\nsince \\( \\phi = \\phi_0 + nt \\) nearly\n\n\\[\ndc = \\frac{A - C}{2A} dt \\left\\{ -(k + k') \\cos \\left( \\phi_0 + nt + \\frac{C - A}{A} (nt + \\gamma) + it + \\varepsilon \\right) \\right. \\\\\n- (k - k') \\cos \\left( \\phi_0 + nt + \\frac{C - A}{A} (nt + \\gamma) - it - \\varepsilon \\right) \\right\\}\n\\]\n\n\\[\nc \\frac{C - A}{A} d\\gamma + \\frac{A - C}{A} nc dt = \\frac{A - C}{2A} dt \\left\\{ (k + k') \\sin \\left( \\phi_0 + nt + \\frac{C - A}{A} (nt + \\gamma) + it + \\varepsilon \\right) \\right. \\\\\n+ (k - k') \\sin \\left( \\phi_0 + nt + \\frac{C - A}{A} (nt + \\gamma) - it - \\varepsilon \\right) \\right\\}\n\\]\n\n\\[\nd\\omega + \\frac{dc}{n} \\frac{A}{C} \\sin \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) + \\frac{c}{n} \\frac{C - A}{C} \\cos \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) d\\gamma = 0\n\\]\n\n\\[\nd\\psi_0 - \\frac{dc}{n \\sin \\omega} \\frac{A}{C} \\cos \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) + \\frac{c \\cos \\omega}{n \\sin \\omega^2} \\frac{A}{C} \\cos \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) d\\omega \\\\\n+ \\frac{c}{n \\sin \\omega} \\frac{C - A}{A} \\sin \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) d\\gamma = 0\n\\]\n\n\\[\nd\\phi_0 - \\frac{dc}{n \\sin \\omega} \\frac{A}{C} \\cos \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) + \\frac{c}{n \\sin \\omega^2} \\frac{A}{C} \\cos \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) d\\omega \\\\\n+ \\frac{c}{n \\sin \\omega} \\frac{C - A}{A} \\sin \\left( \\frac{C}{A} nt + \\frac{C - A}{A} \\gamma \\right) d\\gamma = 0\n\\]\n\nFrom the preceding expressions it may be inferred that\n\n\\( n = \\text{constant.} \\)\n\n\\( \\frac{dc}{dt} = \\text{series of cosines without any constant quantity, unless the mean mo-} \\)\ntion of rotation is commensurate to the mean motion of revolution of the luminary $M'$.\n\n$$\\frac{d\\gamma}{dt} = \\text{series of sines without any constant quantity, except in a similar case.}$$\n\n$c$ being equal to a constant $+$ a series of sines.\n\n$$\\frac{dw}{dt} = \\text{a series of sines without any constant quantity.}$$\n\n$$\\frac{d\\psi_0}{dt} = \\text{a series of cosines} + \\text{a constant quantity.}$$\n\n$$\\frac{d\\phi_0}{dt} = \\text{a series of cosines} + \\text{a constant quantity.}$$\n\nIn the general case where $A$ is not equal to $B$, $n = \\text{constant} + \\text{series of cosines}$.\n\nThe form of the preceding expressions is not affected however, for the approximation may be carried so that except in the case of commensurability above mentioned, the mean motion of rotation being also nearly twice the mean motion of the planet $M'$ in its orbit or greater,\n\n$$n = \\text{constant} + \\text{series of cosines without any constant quantity multiplied by the time.}$$\n\n$c = \\text{constant} + \\text{series of sines or cosines without any constant quantity multiplied by the time.}$\n\n$\\omega = \\text{constant} + \\text{series of cosines without any constant quantity multiplied by the time.}$\n\n$\\gamma = \\text{constant} + \\text{series of cosines or sines} + \\text{a constant quantity multiplied by the time.}$\n\n$\\psi_0 = \\text{constant} + \\text{series of sines} + \\text{a constant quantity multiplied by the time.}$\n\n$\\phi_0 = \\text{constant} + \\text{series of sines} + \\text{a constant quantity multiplied by the time.}$\n\nThe constant quantity multiplied by the time in the value of $\\psi_0$ is the precession of the equinox.\n\nIf $z' = 0$, (which amounts to taking for the fixed plane the orbit of the planet $M'$) and $n'$ be its mean motion, then neglecting the eccentricity,\n\n$$x' = a' \\cos n't, \\quad y' = a' \\sin n't, \\quad r' = a',$$\n\n$$P = \\frac{3 M'}{2 r'^3} \\sin \\omega \\cos \\omega (1 - \\cos 2n't), \\quad P' = \\frac{3 M'}{2 r'^3} \\sin \\omega \\sin 2n't.$$\nSupposing $\\phi_0 = 0$, $\\gamma = 0$, $c_0 = 0$, $c_0$ being in fact imperceptible to observation, and neglecting $\\cos 2n't$, $\\sin 2n't$, in order to find the constant part of $\\frac{d\\psi_0}{dt}$,\n\n$$\\frac{dc_0}{dt} = -\\frac{3(A - C)}{2nAa^3} M' \\sin \\omega \\cos \\omega \\cos \\frac{C}{A} n't$$\n\n$$\\frac{cd\\gamma}{dt} = -\\frac{3(A - C)}{2nAa^3} M' \\sin \\omega \\cos \\omega \\sin \\frac{C}{A} n't$$\n\n$$\\frac{d\\psi_0}{dt} = \\frac{3(C - A)}{2nCa^3} M' \\cos \\omega = \\frac{3(C - A)}{2nC} n'^2 \\cos \\omega$$\n\nThis result agrees with that given in the Méc. Cél. vol. ii. p. 318, and with that given by M. Poisson, Mémoires de l'Academie, vol. vii. p. 247. In Laplace's notation $\\omega = h$, $m = n'$. In M. Poisson's notation $\\omega = \\theta$, $m = n'$.\n\nOn the Theory of the Motion of the Planets, continued from Part II. 1830, p. 357.\n\nFrom the general equations given in the Méc. Cél. vol. i. p. 268, the following may be inferred.\n\n$$\\frac{3}{2} \\frac{a}{a_i} b_{1,2} = \\frac{(a^2 + a_i^2)}{a_i^2} b_{1,1} - \\frac{a}{a_i} b_{1,0}$$\n\n$$\\frac{5}{2} \\frac{a}{a_i} b_{1,3} = \\frac{2(a^2 + a_i^2)}{a_i^2} b_{1,2} - \\frac{3}{2} \\frac{a}{a_i} b_{1,1}$$\n\n$$\\frac{7}{2} \\frac{a}{a_i} b_{1,4} = \\frac{3(a^2 + a_i^2)}{a_i^2} b_{1,3} - \\frac{5}{2} \\frac{a}{a_i} b_{1,2}$$\n\n$$\\frac{9}{2} \\frac{a}{a_i} b_{1,5} = \\frac{4(a^2 + a_i^2)}{a_i^2} b_{1,4} - \\frac{7}{2} \\frac{a}{a_i} b_{1,3}$$\n\n$$\\frac{a}{2a_i} b_{3,2} = \\frac{(a^2 + a_i^2)}{a_i^2} b_{3,1} - \\frac{3}{2} \\frac{a}{a_i} b_{3,0}$$\n\n$$\\frac{3}{2} \\frac{a}{a_i} b_{3,3} = \\frac{2(a^2 + a_i^2)}{a_i^2} b_{3,2} - \\frac{5}{2} \\frac{a}{a_i} b_{3,1}$$\n\n$$\\frac{5}{2} \\frac{a}{a_i} b_{3,4} = \\frac{3(a^2 + a_i^2)}{a_i^2} b_{3,3} - \\frac{7}{2} \\frac{a}{a_i} b_{3,2}$$\n\n$$\\frac{7}{2} \\frac{a}{a_i} b_{3,5} = \\frac{4(a^2 + a_i^2)}{a_i^2} b_{3,4} - \\frac{9}{2} \\frac{a}{a_i} b_{3,3}$$\n\n$$-\\frac{a}{2a_i} b_{5,2} = \\frac{(a^2 + a_i^2)}{a_i^2} b_{5,1} - \\frac{5}{2} \\frac{a}{a_i} b_{5,0}$$\n\n$$\\frac{a}{2a_i} b_{5,3} = \\frac{2(a^2 + a_i^2)}{a_i^2} b_{5,2} - \\frac{7}{2} \\frac{a}{a_i} b_{5,1}$$\n\n$$\\frac{3}{2} \\frac{a}{a_i} b_{5,4} = \\frac{3(a^2 + a_i^2)}{a_i^2} b_{5,3} - \\frac{9}{2} \\frac{a}{a_i} b_{5,2}$$\n\n$$\\frac{5}{2} \\frac{a}{a_i} b_{5,5} = \\frac{4(a^2 + a_i^2)}{a_i^2} b_{5,4} - \\frac{11}{2} \\frac{a}{a_i} b_{5,3}$$\n\n$$b_{1,0} = \\frac{a^2 + a_i^2}{a_i^2} b_{3,0} - \\frac{a}{a_i} b_{3,1}$$\n\n$$b_{1,1} = \\frac{a^2 + a_i^2}{a_i^2} b_{3,1} - \\frac{2a}{a_i} b_{3,0} - \\frac{a}{a_i} b_{3,1}$$\n\nMDCCCXXXI.\n\\[ b_{1,2} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,2} - \\frac{a}{a_i} b_{5,1} - \\frac{a}{a_i} b_{5,3} \\]\n\n\\[ b_{1,3} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,3} - \\frac{a}{a_i} b_{5,2} - \\frac{a}{a_i} b_{5,4} \\]\n\n\\[ b_{5,0} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,0} - \\frac{a}{a_i} b_{5,1} \\]\n\n\\[ b_{5,1} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,1} - 2 \\frac{a}{a_i} b_{5,0} - \\frac{a}{a_i} b_{5,2} \\]\n\n\\[ b_{5,2} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,2} - \\frac{a}{a_i} b_{5,1} - \\frac{a}{a_i} b_{5,3} \\]\n\n\\[ b_{5,3} = \\frac{a^2 + a_i^2}{a_i^2} b_{5,3} - \\frac{a}{a_i} b_{5,2} - \\frac{a}{a_i} b_{5,4} \\]\n\n\\[ b_{1,1} = \\frac{a}{a_i} \\left\\{ b_{5,0} - \\frac{1}{2} b_{5,2} \\right\\} \\]\n\n\\[ 2b_{1,2} = \\frac{a}{2a_i} \\left\\{ b_{5,1} - b_{5,3} \\right\\} \\]\n\n\\[ 3b_{1,3} = \\frac{a}{2a_i} \\left\\{ b_{5,2} - b_{5,4} \\right\\} \\]\n\n\\[ 4b_{1,4} = \\frac{a}{2a_i} \\left\\{ b_{5,3} - b_{5,5} \\right\\} \\]\n\n\\[ b_{5,1} = 3 \\frac{a}{a_i} \\left\\{ b_{5,0} - \\frac{1}{2} b_{5,2} \\right\\} \\]\n\n\\[ 2b_{5,2} = \\frac{3}{2} \\frac{a}{a_i} \\left\\{ b_{5,1} - b_{5,3} \\right\\} \\]\n\n\\[ 3b_{5,3} = \\frac{3}{2} \\frac{a}{a_i} \\left\\{ b_{5,2} - b_{5,4} \\right\\} \\]\n\n\\[ 4b_{5,4} = \\frac{3}{2} \\frac{a}{a_i} \\left\\{ b_{5,3} - b_{5,5} \\right\\} \\]\n\n\\[ 2b_{5,0} = \\frac{2 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,0} - 2 \\frac{a}{a_i} b_{5,1}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ b_{5,1} = \\frac{3 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,1} - 2 \\cdot 3 \\frac{a}{a_i} b_{5,2}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ b_{5,2} = \\frac{5 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,2} - 2 \\cdot 5 \\frac{a}{a_i} b_{5,3}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ b_{5,3} = \\frac{7 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,3} - 2 \\cdot 7 \\frac{a}{a_i} b_{5,4}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ 2b_{5,0} = \\frac{2 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,0} + \\frac{2}{3} \\frac{a}{a_i} b_{5,1}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ b_{5,1} = \\frac{5}{3} \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,1} - \\frac{2}{3} \\frac{a}{a_i} b_{5,2} \\]\n\n\\[ b_{5,2} = \\frac{7}{3} \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,2} - \\frac{2}{3} \\frac{a}{a_i} b_{5,3} \\]\n\n\\[ b_{5,3} = \\frac{9}{3} \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,3} - \\frac{2 \\cdot 5}{3} \\frac{a}{a_i} b_{5,4} \\]\n\n\\[ 2b_{5,0} = \\frac{2 \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,0} + \\frac{2 \\cdot 3}{5} \\frac{a}{a_i} b_{5,1}}{\\left( 1 - \\frac{a^2}{a_i^2} \\right)^2} \\]\n\n\\[ b_{5,1} = \\frac{7}{5} \\left( 1 + \\frac{a^2}{a_i^2} \\right) b_{5,1} + \\frac{2}{5} \\frac{a}{a_i} b_{5,2} \\]\n\\[ b_{7,2} = \\frac{9}{5} \\left(1 + \\frac{a^2}{a_i^2}\\right) b_{5,3} - \\frac{2}{5} \\frac{a}{a_i} b_{5,3}, \\quad b_{7,3} = \\frac{11}{5} \\left(1 + \\frac{a^2}{a_i^2}\\right) b_{5,3} - \\frac{2}{5} \\frac{3a}{a_i} b_{5,4} \\]\n\n\\[ r = a F(n t), \\quad \\frac{dR}{dr} = \\frac{d}{da} \\frac{R}{da} F(n t), \\quad \\frac{daFnt}{da} = Fnt, \\quad \\frac{rdR}{dr} = a \\left(\\frac{dR}{da}\\right), \\]\n\n\\[ \\left\\{a^2 - 2aa_i \\cos \\theta + a_i^2\\right\\}^{-\\frac{1}{2}} = \\frac{1}{a_i} \\left\\{b_{1,0} + b_{1,1} \\cos \\theta + b_{1,2} \\cos 2\\theta + &c.\\right\\} \\]\n\n\\[ -\\left(a - a_i \\cos \\theta\\right)^{-\\frac{1}{2}} = \\frac{1}{a_i} \\left\\{\\frac{db_{1,0}}{da} + \\frac{db_{1,1}}{da} \\cos \\theta + \\frac{db_{1,2}}{da} \\cos 2\\theta + &c.\\right\\} \\]\n\n\\[ -\\left\\{\\frac{a - a_i \\cos \\theta}{a_i^3}\\right\\} \\left\\{b_{3,0} + b_{3,1} \\cos \\theta + b_{3,2} \\cos 2\\theta + &c.\\right\\} = \\frac{1}{a_i} \\left\\{\\frac{db_{1,0}}{da} + \\frac{db_{1,1}}{da} \\cos \\theta + \\frac{db_{1,2}}{da} \\cos 2\\theta + &c.\\right\\} \\]\n\nwhence\n\n\\[ \\frac{ad.b_{1,0}}{da} = -\\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{3,0} - \\frac{1}{2} b_{3,1}\\right\\} \\]\n\n\\[ \\frac{ad.b_{1,1}}{da} = -\\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{3,1} - b_{3,0} - \\frac{1}{2} b_{3,2}\\right\\} \\]\n\n\\[ \\frac{ad.b_{1,2}}{da} = -\\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{3,2} - \\frac{1}{2} b_{3,1} - \\frac{1}{2} b_{3,3}\\right\\} \\]\n\n\\[ \\frac{ad.b_{1,3}}{da} = -\\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{3,3} - \\frac{1}{2} b_{3,2} - \\frac{1}{2} b_{3,4}\\right\\} \\]\n\nsimilarly\n\n\\[ \\frac{ad.b_{3,0}}{da} = -3 \\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{5,0} - \\frac{1}{2} b_{5,1}\\right\\} \\]\n\n\\[ \\frac{ad.b_{3,1}}{da} = -3 \\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{5,1} - b_{5,0} - \\frac{1}{2} b_{5,2}\\right\\} \\]\n\n\\[ \\frac{ad.b_{5,0}}{da} = -5 \\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{7,0} - \\frac{1}{2} b_{7,1}\\right\\} \\]\n\n\\[ \\frac{ad.b_{5,1}}{da} = -5 \\frac{a}{a_i} \\left\\{\\frac{a}{a_i} b_{7,1} - b_{7,0} - \\frac{1}{2} b_{7,2}\\right\\} \\]\n\\[\n\\frac{a^2 d^3 b_{1,0}}{da^2} = \\frac{a}{a_i} \\left[ \\frac{2a}{a_i} b_{3,0} - \\frac{1}{2} b_{3,1} \\right]\n\\]\n\n\\[\n\\frac{a^2 d^3 b_{1,1}}{da^2} = \\frac{a}{a_i} \\left[ \\frac{2a}{a_i} b_{3,1} - b_{3,2} \\right]\n\\]\n\n\\[\n\\frac{a^2 d^3 b_{1,2}}{da^2} = \\frac{a}{a_i} \\left[ \\frac{2a}{a_i} b_{3,2} + \\frac{1}{2} b_{3,1} - \\frac{3}{2} b_{3,3} \\right]\n\\]\n\n\\[\n\\frac{a^2 d^3 b_{1,3}}{da^2} = \\frac{a}{a_i} \\left[ \\frac{2a}{a_i} b_{3,3} + b_{3,2} - 2b_{3,4} \\right]\n\\]\n\n\\[\n\\frac{a^2 d^3 b_{1,4}}{da^2} = \\frac{a}{a_i} \\left[ \\frac{2a}{a_i} b_{3,4} + \\frac{3}{2} b_{3,3} - \\frac{5}{2} b_{3,5} \\right]\n\\]\n\nThe value of \\( R \\) given p. 349 of the former part of this paper is susceptible of much simplification. The first term of \\( R \\) for instance\n\n\\[\nm_i \\left\\{ -\\frac{b_{1,0}}{a_i} + \\frac{3(a^2 e^2 + a_i^2 e_i^2)}{2.2 a_i^3} b_{3,0} + \\frac{a a_i}{2 a_i^3} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{(e^2 + e_i^2)}{2} \\right) b_{3,1} \\right.\n\\]\n\n\\[\n- \\frac{3}{2.4 a_i^5} \\left( 2a^4 e^2 + 5a^2 a_i^2 (e^2 + e_i^2) + 2a_i^4 e_i^2 \\right) b_{5,0} + \\frac{3.2}{2.4 a_i^5} (a^2 e^2 + a_i^2 e_i^2) a a_i b_{5,1}\n\\]\n\n\\[\n+ \\frac{1.3.3}{2.4.2} \\frac{a^2 a_i^2}{a_i^5} (e^2 + e_i^2) b_{5,2} \\right\\}\n\\]\n\n\\[\n= m_i \\left\\{ -\\frac{b_{1,0}}{a_i} + \\frac{3(a^2 e^2 + a_i^2 e_i^2)}{2.2 a_i^3} b_{3,0} + \\frac{a a_i}{2 a_i^3} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) b_{3,1} \\right.\n\\]\n\n\\[\n- \\frac{3.2}{2.4} \\frac{(a^2 + a_i^2)}{a_i^3} \\frac{(a^2 e^2 + a_i^2 e_i^2)}{a_i^3} b_{5,0} - \\frac{3.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} (e^2 + e_i^2) b_{5,0}\n\\]\n\n\\[\n+ \\frac{3.2}{2.4} \\frac{(a^2 e^2 + a_i^2 e_i^2)}{a_i^3} \\frac{a}{a_i} b_{5,1} + \\frac{3.3}{2.4.2} \\frac{a^2 a_i^2}{a_i^5} (e^2 + e_i^2) b_{5,2} \\right\\}\n\\]\n\nand since\n\n\\[\nb_{3,0} = \\frac{(a^2 + a_i^2)}{a_i^2} b_{5,0} - \\frac{a}{a_i} b_{3,1} \\quad \\text{See p. 26}\n\\]\n\n\\[\nb_{3,1} = \\frac{3a}{a_i} \\left\\{ b_{5,0} - \\frac{1}{2} b_{5,2} \\right\\}\n\\]\n\nthis term reduces itself to\n\n\\[\nm_i \\left\\{ -\\frac{b_{1,0}}{a_i} + \\frac{a a_i}{2 a_i^3} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) b_{3,1} - \\frac{3}{2.4} \\frac{a a_i}{a_i^3} (e^2 + e_i^2) b_{3,1} \\right\\}\n\\]\n\n\\[\n= m_i \\left\\{ -\\frac{b_{1,0}}{a_i} + \\frac{a}{2 a_i^3} \\left( \\sin^2 \\frac{\\theta_i}{2} - \\frac{e^2 + e_i^2}{4} \\right) b_{3,1} \\right\\}\n\\]\nThe succeeding terms admit of similar simplifications, so that\n\n\\[ R = m_i \\left\\{ -\\frac{b_{1,0}}{a_i} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) b_{3,1} - \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,1} \\right\\} \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{a_i^2} \\left( \\cos^2 \\frac{\\theta_i}{2} - \\frac{e^2 + e_i^2}{2} \\right) - \\frac{b_{1,1}}{a_i} + \\frac{a}{a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) \\left( b_{3,0} + \\frac{1}{2} b_{3,2} \\right) \\right\\} \\]\n\n\\[ - \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,2} \\cos(nt - n_i t + \\varepsilon - \\varepsilon_i) \\] [1]\n\n\\[ + m_i \\left\\{ -\\frac{b_{1,2}}{a_i} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) (b_{3,1} + b_{3,3}) + \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,1} \\right\\} \\]\n\n\\[ - \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,3} \\cos(2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i) \\] [2]\n\n\\[ + m_i \\left\\{ -\\frac{b_{1,3}}{a_i} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) (b_{3,2} + b_{3,4}) + \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,2} \\right\\} \\]\n\n\\[ - \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,4} \\cos(3nt - 3n_i t - 3\\varepsilon - 3\\varepsilon_i) \\] [3]\n\n\\[ + m_i \\left\\{ -\\frac{b_{1,4}}{a_i} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} + \\frac{e^2 + e_i^2}{2} \\right) (b_{3,3} + b_{3,5}) + \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,3} \\right\\} \\]\n\n\\[ - \\frac{3}{2.4a_i^3} (e^2 + e_i^2) b_{3,5} \\cos(4nt - 4n_i t + 4\\varepsilon - 4\\varepsilon_i) \\] [4]\n\nThe coefficient of \\( \\cos(\\varpi - \\varpi_i) \\), Argument 41,\n\n\\[ = m_i \\left\\{ -\\frac{9}{4a_i^3} b_{3,0} - \\frac{a}{8a_i^3} b_{3,2} + \\frac{3.6}{2.4a_i^3} \\left( \\frac{a^2 + a_i^2}{a_i^2} \\right) a a_i b_{5,0} - \\frac{3.7}{2.4a_i^3} \\left( \\frac{a^2 + a_i^2}{a_i^2} \\right) a a_i b_{5,1} \\right\\} ee_i \\]\n\n\\[ = m_i \\left\\{ -\\frac{9}{4a_i^3} b_{3,0} - \\frac{a}{8a_i^3} b_{3,2} + \\frac{3.6}{2.4a_i^3} \\left( \\frac{a^2 + a_i^2}{a_i^2} \\right) b_{5,0} - \\frac{a}{a_i} b_{5,1} \\right\\} \\]\n\n\\[ - \\frac{3}{8a_i^3} \\left\\{ \\frac{(a^2 + a_i^2)}{a_i^2} b_{5,2} - \\frac{a}{a_i} b_{5,1} - \\frac{a}{a_i} b_{5,3} \\right\\} + \\frac{9}{16a_i^3} \\left\\{ b_{5,1} - b_{5,3} \\right\\} \\]\n\n\\[ = m_i \\left\\{ -\\frac{9}{4a_i^3} b_{3,0} - \\frac{a}{8a_i^3} b_{3,2} + \\frac{9}{4a_i^3} b_{3,0} - \\frac{3}{8a_i^3} b_{3,2} + \\frac{2.3}{8a_i^3} b_{3,2} \\right\\} ee_i \\]\n\n\\[ = m_i \\frac{a}{4a_i^3} b_{3,2} \\]\n\nSo that the part of \\( R \\) which is independent of \\( nt, n_i t \\)\n\n\\[ = m_i \\left\\{ -\\frac{b_{1,0}}{a_i^2} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} - \\frac{e^2 + e_i^2}{4} \\right) b_{3,1} + \\frac{a}{4a_i^3} b_{3,2} ee_i \\cos(\\varpi - \\varpi_i) \\right\\} \\]\n\n\\[ = m_i \\left\\{ -\\frac{b_{1,0}}{a_i} - \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\theta_i}{2} - \\frac{e^2 + e_i^2}{4} \\right) b_{3,1} - \\left\\{ \\frac{3}{2a_i^3} b_{3,0} - \\frac{a^2 + a_i^2}{2a_i^3} b_{3,1} \\right\\} ee_i \\cos(\\varpi - \\varpi_i) \\right\\} \\]\nIn the general case, when \\( \\iota_1, \\iota_2 \\) are the inclinations of the orbits of the planets \\( P \\) and \\( P_1 \\) to any plane, the direction of which is arbitrary,\n\n\\[\n\\cos \\iota = \\cos \\iota_1 \\cos \\iota_2 + \\sin \\iota_1 \\sin \\iota_2 \\cos (\\nu_1 - \\nu_2)\n\\]\n\n\\[\n\\sin^2 \\frac{\\iota}{2} = \\frac{1 - \\cos \\iota_1 \\cos \\iota_2 - \\sin \\iota_1 \\sin \\iota_2 \\cos (\\nu_1 - \\nu_2)}{2}\n\\]\n\nThe part of \\( R \\) which is independent of \\( n t, n_1 t \\)\n\n\\[\n= m_i \\left[ -\\frac{b_{1,0}}{a_i} + \\frac{a}{4a_i^3} \\left\\{ 1 - \\cos \\iota_1 \\cos \\iota_2 - \\sin \\iota_1 \\sin \\iota_2 \\cos (\\nu_1 - \\nu_2) - \\frac{e^2 + e_i^2}{2} \\right\\} b_{3,1} \\right.\n\\]\n\n\\[\n- \\left\\{ \\frac{3a}{2a_i^3} b_{3,0} - \\frac{(a^2 + a_i^2)}{2a_i^3} b_{3,1} \\right\\} e e_i \\cos (\\varpi - \\varpi_i) \\right]\n\\]\n\n\\[\n\\cos \\iota = \\frac{1}{\\sqrt{1 + \\tan^2 \\iota}} = 1 - \\frac{1}{2} \\tan^2 \\iota, \\quad \\sin \\iota = \\frac{\\tan \\iota}{\\sqrt{1 + \\tan^2 \\iota}} = \\tan \\iota \\quad \\text{nearly}.\n\\]\n\n\\[\n= m_i \\left[ -\\frac{b_{1,0}}{a_i} + \\frac{a}{8a_i^3} \\left\\{ (\\tan^2 \\iota_1 + \\tan^2 \\iota_2 - 2 \\tan \\iota_1 \\tan \\iota_2 \\cos (\\nu_1 - \\nu_2)) - e^2 - e_i^2 \\right\\} b_{3,1} \\right.\n\\]\n\n\\[\n- \\left\\{ \\frac{3a}{2a_i^3} b_{3,0} - \\frac{(a^2 + a_i^2)}{2a_i^3} b_{3,1} \\right\\} e e_i \\cos (\\varpi - \\varpi_i) \\right]\n\\]\n\n\\[\n= m_i \\left[ -\\frac{b_{1,0}}{a_i} + \\frac{a}{8a_i^3} \\left\\{ (\\tan \\iota_1 \\cos \\nu_1 - \\tan \\iota_2 \\cos \\nu_2)^2 + (\\tan \\iota_1 \\sin \\nu_1 - \\tan \\iota_2 \\sin \\nu_2)^2 - e^2 - e_i^2 \\right\\} b_{3,1} \\right.\n\\]\n\n\\[\n- \\left\\{ \\frac{3a}{2a_i^3} b_{3,0} - \\frac{(a^2 + a_i^2)}{2a_i^3} b_{3,1} \\right\\} e e_i \\cos (\\varpi - \\varpi_i) \\right]\n\\]\n\nand if \\( \\tan \\iota \\sin \\nu = p, \\tan \\iota \\cos \\nu = q \\), this quantity\n\n\\[\n= m_i \\left[ -\\frac{b_{1,0}}{a_i} + \\frac{a}{8a_i^3} \\left\\{ (p_1 - p_2)^2 + (q_1 - q_2)^2 - e^2 - e_i^2 \\right\\} b_{3,1} \\right.\n\\]\n\n\\[\n- \\left\\{ \\frac{3a}{2a_i^3} b_{3,0} - \\frac{(a^2 + a_i^2)}{2a_i^3} b_{3,1} \\right\\} e e_i \\cos (\\varpi - \\varpi_i) \\right]\n\\]\n\nwhich evidently agrees with the result given by M. de Pontécoulant, Théor. Anal. vol. i. p. 363. All the other coefficients of terms multiplied by the squares and products of the eccentricities are susceptible of reductions similar to those in the two preceding pages, and finally;\n\n\\[\nR = m_i \\left[ -\\frac{b_{1,0}}{a_i} + \\frac{a}{2a_i^3} \\left( \\sin^2 \\frac{\\iota}{2} - \\frac{e^2 + e_i^2}{4} \\right) b_{3,1} \\right]\n\\]\n\n\\[\n+ m_i \\left[ -\\frac{a}{a_i^2} \\left( \\cos^2 \\frac{\\iota}{2} - \\frac{e^2 + e_i^2}{2} \\right) - \\frac{b_{1,1}}{a_i} + \\frac{a}{a_i^3} \\sin^2 \\frac{\\iota}{2} (b_{3,0} + \\frac{1}{2} b_{3,2}) \\right.\n\\]\n\n\\[\n+ \\frac{a}{a_i^3} \\left( \\frac{e^2 + e_i^2}{8} \\right) (4 b_{3,0} - 4 b_{3,2}) \\right] \\cos (n t - n_1 t + \\varepsilon - \\varepsilon_i)\n\\]\n\\[ + m_i \\left\\{ - \\frac{b_{1,2}}{a_i} + \\frac{a}{2a_i^2} \\sin^2 \\frac{\\theta_i}{2} (b_{3,1} + b_{3,3}) \\right. \\\\\n+ \\left. \\frac{a}{a_i^2} \\left( e^2 + e_i^2 \\right) \\frac{5}{8} (5 b_{3,1} - 7 b_{3,3}) \\right\\} \\cos (2 n t - 2 n_i t + 2 \\varepsilon - 2 \\varepsilon_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{b_{1,3}}{a_i} + \\frac{a}{2a_i^2} \\sin^2 \\frac{\\theta_i}{2} (b_{3,2} + b_{3,4}) \\right. \\\\\n+ \\left. \\frac{a}{a_i^2} \\left( e^2 + e_i^2 \\right) \\frac{8}{8} (8 b_{3,2} - 10 b_{3,4}) \\right\\} \\cos (3 n t - 3 n_i t + 3 \\varepsilon - 3 \\varepsilon_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{b_{1,4}}{a_i} + \\frac{a}{2a_i^2} \\sin^2 \\frac{\\theta_i}{2} (b_{3,3} + b_{3,5}) \\right. \\\\\n+ \\left. \\frac{a}{a_i^2} \\left( e^2 + e_i^2 \\right) \\frac{11}{8} (11 b_{3,3} - 13 b_{3,5}) \\right\\} \\cos (4 n t - 4 n_i t + 4 \\varepsilon - 4 \\varepsilon_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{b_{1,5}}{a_i} + \\frac{a}{2a_i^2} \\sin^2 \\frac{\\theta_i}{2} (b_{3,4} + b_{3,6}) \\right. \\\\\n+ \\left. \\frac{a}{a_i^2} \\left( e^2 + e_i^2 \\right) \\frac{14}{8} (14 b_{3,4} - 16 b_{3,6}) \\right\\} \\cos (5 n t - 5 n_i t + 5 \\varepsilon - 5 \\varepsilon_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{3a}{2a_i^3} + \\frac{3a}{2a_i^3} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} - \\frac{a}{4a_i^2} b_{3,2} \\right\\} e \\cos (n t + \\varepsilon - \\varpi) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a^3}{a_i^3} b_{3,0} + \\frac{a}{2a_i^2} b_{3,1} \\right\\} e \\cos (n t + \\varepsilon - \\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{2a_i^2} - \\frac{a}{2a_i^2} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} + \\frac{3a}{4a_i^2} b_{3,2} \\right\\} e \\cos (2 n t - n_i t + 2 \\varepsilon - \\varepsilon_i - \\varpi) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a}{4a_i^2} b_{3,1} - \\frac{a^2}{2a_i^3} b_{3,2} + \\frac{3a}{4a_i^2} b_{3,3} \\right\\} e \\cos (3 n t - 2 n_i t + 3 \\varepsilon - 2 \\varepsilon_i - \\varpi) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a}{4a_i^2} b_{3,2} - \\frac{a^2}{2a_i^3} b_{3,3} + \\frac{3a}{4a_i^2} b_{3,4} \\right\\} e \\cos (4 n t - 3 n_i t + 4 \\varepsilon - 3 \\varepsilon_i - \\varpi) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a}{4a_i^2} b_{3,3} - \\frac{a}{2a_i^3} b_{3,4} + \\frac{3a}{4a_i^2} b_{3,5} \\right\\} e \\cos (5 n t - 4 n_i t + 5 \\varepsilon - 4 \\varepsilon_i - \\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,1} - \\frac{a^2}{2a_i^3} b_{3,2} - \\frac{a}{4a_i^2} b_{3,3} \\right\\} e \\cos (n t - 2 n_i t + \\varepsilon - 2 \\varepsilon_i + \\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,2} - \\frac{a^2}{2a_i^3} b_{3,3} - \\frac{a}{4a_i^2} b_{3,4} \\right\\} e \\cos (2 n t - 3 n_i t + 2 \\varepsilon - 3 \\varepsilon_i + \\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,3} - \\frac{a^2}{2a_i^3} b_{3,4} - \\frac{a}{4a_i^2} b_{3,5} \\right\\} e_i \\cos (3 n t - 4 n_i t + 3 \\varepsilon - 4 \\varepsilon_i + \\varpi) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a^3}{a_i^3} b_{3,0} + \\frac{a}{2a_i^2} b_{3,1} \\right\\} e_i \\cos (n_i t + \\varepsilon_i - \\varpi_i) \\]\n\n* These numbers indicate the arguments which are symmetrical with regard to \\( n t \\) and \\( n_i t \\).\nDevelopment of $R.$\n\n\\[ + m_i \\left\\{ \\frac{3a}{2a_i^2} b_{3,0} - \\frac{a_i^2}{2a_i^3} b_{3,1} - \\frac{a}{4a_i^2} b_{3,2} \\right\\} e_i \\cos (nt + \\varepsilon - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,1} - \\frac{a_i^2}{2a_i^3} b_{3,2} - \\frac{a}{4a_i^2} b_{3,3} \\right\\} e_i \\cos (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,2} - \\frac{a_i^2}{2a_i^3} b_{3,3} - \\frac{a}{4a_i^2} b_{3,4} \\right\\} e_i \\cos (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{4a_i^2} b_{3,3} - \\frac{a_i^2}{2a_i^3} b_{3,4} - \\frac{a}{4a_i^2} b_{3,5} \\right\\} e_i \\cos (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{2a_i^2} - \\frac{a}{2a_i^2} b_{3,0} - \\frac{a_i^2}{2a_i^3} b_{3,1} + \\frac{3a}{4a_i^2} b_{3,2} \\right\\} e_i \\cos (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ -\\frac{a}{4a_i^2} b_{3,1} - \\frac{a_i^2}{2a_i^3} b_{3,2} + \\frac{3a}{4a_i^2} b_{3,3} \\right\\} e_i \\cos (2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ -\\frac{a}{4a_i^2} b_{3,2} - \\frac{a_i^2}{2a_i^3} b_{3,3} + \\frac{3a}{4a_i^2} b_{3,4} \\right\\} e_i \\cos (3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ -\\frac{a}{4a_i^2} b_{3,3} - \\frac{a_i^2}{2a_i^3} b_{3,4} + \\frac{3a}{4a_i^2} b_{3,5} \\right\\} e_i \\cos (4nt - 5n_i t + 4\\varepsilon - 5\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ -\\frac{5a}{8a_i^2} b_{3,1} + \\frac{a^2}{2a_i^3} b_{3,2} \\right\\} e^2 \\cos (2nt + 2\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{8a_i^2} - \\frac{a}{8a_i^2} b_{3,0} - \\frac{a}{16a_i^2} b_{3,2} \\right\\} e^2 \\cos (nt + n_i t + \\varepsilon + \\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{a^2}{a_i^3} b_{3,0} + \\frac{3a}{8a_i^2} b_{3,1} \\right\\} e^2 \\cos (2nt + 2\\varepsilon - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{3a}{8a_i^2} - \\frac{3a}{8a_i^2} b_{3,0} - \\frac{a^2}{a_i^3} b_{3,1} + \\frac{17a}{16a_i^2} b_{3,2} \\right\\} e^2 \\cos (3nt - n_i t + 3\\varepsilon - \\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{a}{4a_i^2} b_{3,1} - \\frac{3a}{2a_i^2} b_{3,2} + \\frac{13a}{8a_i^2} b_{3,3} \\right\\} e^2 \\cos (4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{5a}{16a_i^2} b_{3,2} - \\frac{2a^2}{a_i^3} b_{3,3} + \\frac{35a}{16a_i^2} b_{3,4} \\right\\} e^2 \\cos (5nt - 3n_i t + 5\\varepsilon - 3\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{3a}{8a_i^2} b_{3,3} - \\frac{5a^2}{a_i^3} b_{3,4} + \\frac{11a}{4a_i^2} b_{3,5} \\right\\} e^2 \\cos (6nt - 4n_i t + 6\\varepsilon - 4\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{7a}{16a_i^2} b_{3,4} - \\frac{3a^2}{a_i^3} b_{3,5} + \\frac{53a}{16a_i^2} b_{3,6} \\right\\} e^2 \\cos (7nt - 5n_i t + 7\\varepsilon - 5\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{19a}{16a_i^2} b_{3,5} + \\frac{a^2}{a_i^3} b_{3,6} + \\frac{a}{16a_i^2} b_{3,4} \\right\\} e^2 \\cos (nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ -\\frac{7a}{4a_i^2} b_{3,3} + \\frac{3a^2}{a_i^3} b_{3,4} + \\frac{1a}{8a_i^2} b_{3,5} \\right\\} e^2 \\cos (2nt - 4n_i t + 2\\varepsilon - 4\\varepsilon_i + 2\\varpi) \\]\n\\[ + m_i \\left\\{ -\\frac{37}{16} \\frac{a}{a_i^2} b_{3,4} + 2 \\frac{a^2}{a_i^3} b_{3,5} + \\frac{3}{16} \\frac{a}{a_i^2} b_{3,6} \\right\\} e^2 \\cos (3nt - 5n_t + 3\\varepsilon - 5\\varepsilon_i + 2\\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{a_i^2} b_{3,1} ee_i \\cos (nt - n_t + \\varepsilon - \\varepsilon_i - \\varpi + \\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{a_i^2} - \\frac{a}{a_i^2} b_{3,0} + \\frac{3}{4} \\frac{a}{a_i^3} b_{3,2} \\right\\} ee_i \\cos (2nt - 2n_t + 2\\varepsilon - 2\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{7}{8} \\frac{a}{a_i^2} b_{3,1} - \\frac{9}{8} \\frac{a}{a_i^2} b_{3,3} \\right\\} ee_i \\cos (3nt - 3n_t + 3\\varepsilon - 3\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{5}{4} \\frac{a}{a_i^2} b_{3,2} - \\frac{3}{2} \\frac{a}{a_i^3} b_{3,4} \\right\\} ee_i \\cos (4nt - 4n_t + 4\\varepsilon - 4\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{13}{8} \\frac{a}{a_i^2} b_{3,3} - \\frac{15}{8} \\frac{a}{a_i^2} b_{3,5} \\right\\} ee_i \\cos (5nt - 5n_t + 5\\varepsilon - 5\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{2}{a_i^2} b_{3,4} - \\frac{9}{4} \\frac{a}{a_i^3} b_{3,6} \\right\\} ee_i \\cos (6nt - 6n_t + 6\\varepsilon - 6\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{4a_i^2} b_{3,2} ee_i \\cos (\\varpi - \\varpi_i) \\right\\} \\]\n\n\\[ - m_i \\left\\{ \\frac{3}{8} \\frac{a}{a_i^2} b_{3,1} - \\frac{5}{8} \\frac{a}{a_i^2} b_{3,3} \\right\\} ee_i \\cos (nt - n_t + \\varepsilon - \\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{3}{4} \\frac{a}{a_i^2} b_{3,2} - \\frac{a}{a_i^3} b_{3,4} \\right\\} ee_i \\cos (2nt - 2n_t + 2\\varepsilon - 2\\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{9}{8} \\frac{a}{a_i^2} b_{3,3} - \\frac{11}{8} \\frac{a}{a_i^2} b_{3,5} \\right\\} ee_i \\cos (3nt - 3n_t + 3\\varepsilon - 3\\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[ - m_i \\left\\{ \\frac{3}{2} \\frac{a}{a_i^2} b_{3,4} - \\frac{7}{4} \\frac{a}{a_i^3} b_{3,6} \\right\\} ee_i \\cos (4nt - 4n_t + 4\\varepsilon - 4\\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[ + m' \\left\\{ -\\frac{3a}{a_i^2} + \\frac{1}{a_i} b_{3,1} - \\frac{3}{4} \\frac{a}{a_i^2} b_{3,2} \\right\\} ee_i \\cos (2nt + 2\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{4a_i^2} b_{3,1} ee_i \\cos (nt + n_t + \\varepsilon + \\varepsilon_i - \\varpi - \\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ \\frac{a^2}{a_i^3} b_{3,1} - \\frac{3}{4} \\frac{a}{a_i^2} b_{3,2} \\right\\} ee_i \\cos (2nt + 2\\varepsilon - \\varpi - \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{21}{8} \\frac{a}{a_i^2} b_{3,1} - \\frac{2}{a_i} b_{3,2} - \\frac{3}{8} \\frac{a}{a_i^2} b_{3,3} \\right\\} ee_i \\cos (3nt - n_t + 3\\varepsilon - \\varepsilon_i - \\varpi - \\varpi) \\]\n\n\\[ + m_i \\left\\{ \\frac{15}{4} \\frac{a}{a_i^2} b_{3,2} - \\frac{3}{a_i} b_{3,3} - \\frac{1}{2} \\frac{a}{a_i^2} b_{3,4} \\right\\} ee_i \\cos (4nt - 2n_t + 4\\varepsilon - 2\\varepsilon_i - \\varpi - \\varpi) \\]\nDevelopment of \\( R \\).\n\n\\[\n+ m_i \\left\\{ \\frac{39}{8} \\frac{a}{a_i^2} b_{3,3} - \\frac{4}{a_i} b_{3,4} - \\frac{5}{8} \\frac{a}{a_i^2} b_{3,5} \\right\\} e_i \\cos (5nt - 3n_i t + 5\\epsilon - 3\\epsilon_i - \\omega - \\omega_i) \\quad [51] \\quad [55]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{6}{a_i^2} b_{3,4} - \\frac{5}{a_i} b_{3,5} - \\frac{3}{4} \\frac{a}{a_i^2} b_{3,6} \\right\\} e_i \\cos (6nt - 4n_i t + 6\\epsilon - 4\\epsilon_i - \\omega - \\omega_i) \\quad [52] \\quad [56]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{21}{8} \\frac{a}{a_i^2} b_{3,1} - \\frac{2}{a_i^2} b_{3,2} - \\frac{3}{8} \\frac{a}{a_i^2} b_{3,3} \\right\\} e_i \\cos (nt - 3n_i t + s - 3\\epsilon_i + \\omega + \\omega_i) \\quad [53] \\quad [49]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{15}{4} \\frac{a}{a_i^2} b_{3,2} - \\frac{3}{a_i^3} b_{3,3} - \\frac{a}{2a_i^2} b_{3,4} \\right\\} e_i \\cos (2nt - 4n_i t + 2\\epsilon - 4\\epsilon_i + \\omega + \\omega_i) \\quad [54] \\quad [50]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{39}{8} \\frac{a}{a_i^3} b_{3,3} - \\frac{4}{a_i^2} b_{3,4} - \\frac{5}{8} \\frac{a}{a_i^2} b_{3,5} \\right\\} e_i \\cos (3nt - 5n_i t + 3\\epsilon - 5\\epsilon_i - \\omega + \\omega_i) \\quad [55] \\quad [51]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{6}{a_i^2} b_{3,4} - \\frac{5}{a_i^2} b_{3,5} - \\frac{3}{4} \\frac{a}{a_i^2} b_{3,6} \\right\\} e_i \\cos (4nt - 6n_i t + 4\\epsilon - 6\\epsilon_i - \\omega + \\omega_i) \\quad [56] \\quad [52]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{1}{a_i} b_{3,0} + \\frac{3}{8} \\frac{a}{a_i^2} b_{3,1} \\right\\} e_i^2 \\cos (2n_i t + 2\\epsilon_i - 2\\omega_i) \\quad [57] \\quad [26]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{a}{8a_i^3} - \\frac{a}{8a_i^2} b_{3,0} - \\frac{a}{16a_i^2} b_{3,2} \\right\\} e_i^2 \\cos (nt + n_i t + \\epsilon + \\epsilon_i - 2\\omega_i) \\quad [58] \\quad [25]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{5}{8} \\frac{a}{a_i^2} b_{3,1} + \\frac{1}{2a_i} b_{3,2} \\right\\} e_i^2 \\cos (2nt + 2\\epsilon - 2\\omega_i) \\quad [59] \\quad [24]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{19}{16} \\frac{a}{a_i^2} b_{3,2} + \\frac{1}{a_i} b_{3,3} + \\frac{a}{16a_i^2} b_{3,4} \\right\\} e_i^2 \\cos (3nt - n_i t + 3\\epsilon - \\epsilon_i - 2\\omega_i) \\quad [60] \\quad [32]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{7}{4} \\frac{a}{a_i^2} b_{3,3} + \\frac{3}{2a_i} b_{3,4} + \\frac{a}{8a_i^2} b_{3,5} \\right\\} e_i \\cos (4nt - 2n_i t + 4\\epsilon - 2\\epsilon_i - 2\\omega_i) \\quad [61] \\quad [33]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{37}{16} \\frac{a}{a_i^2} b_{3,4} + \\frac{2}{a_i} b_{3,5} + \\frac{3}{16} \\frac{a}{a_i^2} b_{3,6} \\right\\} e_i^2 \\cos 5nt - 3n_i t + 5\\epsilon - 3\\epsilon_i - 2\\omega_i) \\quad [62] \\quad [34]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{27}{8} \\frac{a}{a_i^3} - \\frac{3}{8} \\frac{a}{a_i^2} b_{3,0} - \\frac{1}{a_i} b_{3,1} + \\frac{17}{16} \\frac{a}{a_i^2} b_{3,3} \\right\\} e_i^2 \\cos (nt - 3n_i t + \\epsilon - 3\\epsilon_i + 2\\omega_i) \\quad [63] \\quad [27]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{a}{4a_i^2} b_{3,1} - \\frac{3}{2a_i} b_{3,2} + \\frac{13}{8} \\frac{a}{a_i^2} b_{3,3} \\right\\} e_i^2 \\cos (2nt - 4n_i t + 2\\epsilon - 4\\epsilon_i + 2\\omega_i) \\quad [64] \\quad [28]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{5}{16} \\frac{a}{a_i^2} b_{3,2} - \\frac{2}{a_i} b_{3,3} + \\frac{35}{16} \\frac{a}{a_i^2} b_{3,4} \\right\\} e_i^2 \\cos (3nt - 5n_i t + 3\\epsilon - 5\\epsilon_i + 2\\omega_i) \\quad [65] \\quad [29]\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{3}{8} \\frac{a}{a_i^2} b_{3,3} - \\frac{5}{2a_i} b_{3,4} + \\frac{11}{4} \\frac{a}{a_i^2} b_{3,5} \\right\\} e_i^2 \\cos (4nt - 6n_i t + 4\\epsilon - 6\\epsilon_i + 2\\omega) \\quad [66] \\quad [30]\n\\]\nIn the lunar theory, the small value of the quantity \\( \\frac{a}{a_i} \\) makes it desirable to ordain the results according to powers of this quantity. Transforming therefore the preceding expression for \\( R \\) by means of the equations given in the former part of this paper, Phil. Trans. for 1830, p. 346, neglecting terms multiplied by \\( \\frac{a^4}{a_i^3} \\), and supposing \\( i = 0 \\),\n\n\\[\nR = m_i \\left\\{ -\\frac{1}{a_i} \\left( 1 + \\frac{a^2}{4a_i^2} \\right) + \\frac{3}{2} \\left( \\sin^2 \\frac{i}{2} - \\frac{(e^2 + e_i^2)}{4} \\right) \\frac{a^2}{a_i^2} \\right\\}\n\\]\n\n\\[\n+ m_i \\left\\{ \\left( \\cos^2 \\frac{i}{2} - \\frac{(e^2 + e_i^2)}{2} \\right) \\frac{a}{a_i} - \\frac{a}{a_i} \\left( 1 + \\frac{3}{8} \\frac{a}{a_i^2} \\right) + \\sin^2 \\frac{i}{2} \\frac{a}{a_i} \\left( 1 + \\frac{33}{8} \\frac{a^2}{a_i^2} \\right) \\right. \\\\\n\\left. + \\frac{(e^2 + e_i^2)}{2} \\frac{a}{a_i} \\left( 1 - \\frac{3}{2} \\frac{a^2}{a_i^2} \\right) \\right\\} \\cos (nt - n_i t + \\varepsilon - \\varepsilon_i)\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{3}{4} \\frac{a^2}{a_i^3} + \\frac{3}{2} \\sin^2 \\frac{i}{2} \\frac{a^2}{a_i^2} + \\frac{15}{8} (e^2 + e_i^2) \\frac{a^2}{a_i^3} \\right\\} \\cos (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i)\n\\]\n\n\\[\n+ m_i \\left\\{ -\\frac{5}{8} \\frac{a^3}{a_i^4} + \\frac{15}{8} \\sin^2 \\frac{i}{2} \\frac{a^3}{a_i^4} + \\frac{15}{4} (e^2 + e_i^2) \\frac{a^3}{a_i^4} \\right\\} \\cos (3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i)\n\\]\nDevelopment of $R$ according to powers of $\\frac{a}{a_i}$.\n\n\\[\n+ m_i \\frac{15}{16} \\frac{a^3}{a_i^4} e \\cos (n_i t + \\varepsilon_i - \\varpi) \\\\\n+ m_i \\frac{3}{2} \\frac{a^2}{a_i^3} e \\cos (n t + \\varepsilon - \\varpi) \\\\\n+ m_i \\frac{3}{16} \\frac{a^3}{a_i^4} e \\cos (2 n t - n_i t + 2 \\varepsilon - \\varepsilon_i - \\varpi) \\\\\n- m_i \\frac{3}{4} \\frac{a^2}{a_i^3} e \\cos (3 n t - 2 n_i t + 3 \\varepsilon - 2 \\varepsilon_i - \\varpi) \\\\\n- m_i \\frac{15}{16} \\frac{a^3}{a_i^4} e \\cos (4 n t - 3 n_i t + 4 \\varepsilon - 3 \\varepsilon_i - \\varpi) \\\\\n+ m_i \\frac{9}{4} \\frac{a^2}{a_i^3} e \\cos (n t - 2 n_i t + \\varepsilon - 2 \\varepsilon_i + \\varpi) \\\\\n+ m_i \\frac{45}{16} \\frac{a^3}{a_i^4} e \\cos (2 n t - 3 n_i t + 2 \\varepsilon - 3 \\varepsilon_i + \\varpi) \\\\\n- m_i \\left\\{ \\frac{1}{a_i} + \\frac{3}{4} \\frac{a^2}{a_i^3} \\right\\} e_i \\cos (n_i t + \\varepsilon_i - \\varpi_i) \\\\\n- m_i \\frac{3}{8} \\frac{a^3}{a_i^4} e_i \\cos (n t + \\varepsilon_i - \\varpi_i) \\\\\n+ m_i \\frac{3}{8} \\frac{a^2}{a_i^3} e_i \\cos (2 n t - n_i t + 2 \\varepsilon - \\varepsilon_i - \\varpi_i) \\\\\n+ m_i \\frac{5}{8} \\frac{a^3}{a_i^4} e_i \\cos (3 n t - 2 n_i t + 3 \\varepsilon - 2 \\varepsilon_i - \\varpi_i) \\\\\n- m_i \\frac{9}{8} \\frac{a^5}{a_i^4} e_i \\cos (n t - 2 n_i t + \\varepsilon - 2 \\varepsilon_i + \\varpi_i) \\\\\n- m_i \\frac{21}{8} \\frac{a^2}{a_i^3} e_i \\cos (2 n t - 3 n_i t + 2 \\varepsilon - 3 \\varepsilon_i + \\varpi_i) \\\\\n- m_i \\frac{25}{8} \\frac{a^3}{a_i^4} e_i \\cos (3 n t - 4 n_i t + 3 \\varepsilon - 4 \\varepsilon_i + \\varpi_i) \\\\\n- m_i \\frac{15}{8} \\frac{a^2}{a_i^3} e^2 \\cos (2 n_i t + 2 \\varepsilon_i - 2 \\varpi) \\\\\n- m_i \\frac{33}{64} \\frac{a^5}{a_i^4} e^2 \\cos (n t + n_i t + \\varepsilon + \\varepsilon_i - 2 \\varpi) \\\\\n+ m_i \\frac{8}{8} \\frac{a^2}{a_i^3} e^2 \\cos (2 n t + 2 \\varepsilon - 2 \\varpi) \\\\\n+ m_i \\frac{9}{64} \\frac{a^3}{a_i^4} e^2 \\cos (3 n t - n_i t + 3 \\varepsilon - \\varepsilon_i - 2 \\varpi)\n\\]\n\\[- m_i \\frac{3}{4} \\frac{a^2}{a_i^3} e^2 \\cos (4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - 2\\varpi) \\]\n\n\\[- m_i \\frac{75}{64} \\frac{a^3}{a_i^4} e^2 \\cos (5nt - 3n_i t + 5\\varepsilon - 3\\varepsilon_i - 2\\varpi) \\]\n\n\\[- m_i \\frac{285}{64} \\frac{a^3}{a_i^4} e^2 \\cos (nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\varpi) \\]\n\n\\[+ m_i \\frac{3}{4} \\frac{a^2}{a_i^3} ee_i \\cos (nt - n_i t + \\varepsilon - \\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[+ m_i \\frac{9}{16} \\frac{a^3}{a_i^4} ee_i \\cos (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[- m_i \\frac{21}{8} \\frac{a^2}{a_i^3} ee_i \\cos (3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i - \\varpi + \\varpi_i) \\]\n\n\\[+ m_i \\frac{15}{16} \\frac{a^3}{a_i^4} ee_i \\cos (\\varpi - \\varpi_i) \\]\n\n\\[- m_i \\frac{9}{8} \\frac{a^2}{a_i^3} ee_i \\cos (nt - n_i t + \\varepsilon - \\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[- m_i \\frac{45}{16} \\frac{a^3}{a_i^4} ee_i \\cos (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i + \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{45}{16} \\frac{a^3}{a_i^4} ee_i \\cos (2n_i t + 2\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{3}{4} \\frac{a^2}{a_i^3} ee_i \\cos (nt + n_i t + \\varepsilon + \\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{3}{16} \\frac{a^3}{a_i^4} ee_i \\cos (2nt + 2\\varepsilon - \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{3}{8} \\frac{a^2}{a_i^3} ee_i \\cos (3nt - n_i t + 3\\varepsilon - \\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{15}{16} \\frac{a^3}{a_i^4} ee_i \\cos (4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[+ m_i \\frac{63}{8} \\frac{a^2}{a_i^3} ee_i \\cos (nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[+ m_i \\frac{225}{16} \\frac{a^3}{a_i^4} ee_i \\cos (2nt - 4n_i t + 2\\varepsilon - 4\\varepsilon_i + \\varpi + \\varpi_i) \\]\nDevelopment of \\( R \\) according to powers of \\( \\frac{a}{a_i} \\):\n\n\\[\n+ m_i \\left\\{ \\frac{1}{a_i} - \\frac{9}{8} \\frac{a^2}{a_i^3} \\right\\} e_i \\cos (2 n_i t + 2 \\varepsilon_i - 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{33}{64} \\frac{a^3}{a_i^4} e_i^5 \\cos (n t + n_i t + \\varepsilon + \\varepsilon_i - 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{5}{64} \\frac{a^3}{a_i^4} e_i^3 \\cos (3 n_i t - n_i t + 3 \\varepsilon - \\varepsilon_i - 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{159}{64} \\frac{a^3}{a_i^4} e_i^3 \\cos (n t - 3 n_i t + \\varepsilon - 3 \\varepsilon_i + 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{51}{8} \\frac{a^3}{a_i^3} e_i^2 \\cos (2 n t - 4 n_i t + 2 \\varepsilon - 4 \\varepsilon_i + 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{635}{64} \\frac{a^3}{a_i^4} e_i^2 \\cos (3 n t - 5 n_i t + 3 \\varepsilon - 5 \\varepsilon_i + 2 \\varpi_i)\n\\]\n\n\\[\n- m_i \\frac{9}{4} \\frac{a^3}{a_i^3} \\sin^2 \\frac{t}{2} \\cos (n t + n_i t + \\varepsilon + \\varepsilon_i - 2 \\nu)\n\\]\n\n\\[\n- m_i \\frac{3}{2} \\frac{a^2}{a_i^3} \\sin^2 \\frac{t}{2} \\cos (2 n t + 2 \\varepsilon - 2 \\nu)\n\\]\n\n\\[\n- m_i \\frac{3}{2} \\frac{a^2}{a_i^3} \\sin^2 \\frac{t}{2} \\cos (2 n_i t + 2 \\varepsilon_i - 2 \\nu)\n\\]\n\n\\[\n- m_i \\frac{15}{8} \\frac{a^3}{a_i^4} \\sin^3 \\frac{t}{2} \\cos (3 n t - n_i t + 3 \\varepsilon - \\varepsilon_i - 2 \\nu)\n\\]\n\n\\[\n- m_i \\frac{15}{8} \\frac{a^3}{a_i^4} \\sin^3 \\frac{t}{2} \\cos (n t - 3 n_i t + \\varepsilon - 3 \\varepsilon_i + 2 \\nu)\n\\]\n\nIf according to the notation of M. Damoiseau, (Théorie Lunaire, p. 547, Mémoires des Savans Étrangers,) \\( n t - n_i t + \\varepsilon - \\varepsilon_i = t \\), and \\( x \\) and \\( z \\) be put for the mean anomalies of \\( m \\) and \\( m_i \\) respectively and \\( y \\) for the distance of the planet \\( m \\) from its node, or, what is the same in the Lunar Theory, the distance of the moon from her node reckoned on the ecliptic (\\( \\iota_i = 0 \\)),\n\n\\[\nR = m_i \\left\\{ - \\frac{1}{a_i} \\left( 1 + \\frac{a^2}{4 a_i^2} \\right) + \\frac{3}{2} \\left( \\sin^2 \\frac{t}{2} - \\frac{(e^2 + e_i^2)}{4} \\right) \\frac{a^2}{a_i^3} \\right\\}\n\\]\n\n\\[\n+ \\left\\{ \\left( \\cos^2 \\frac{t}{2} - \\frac{(e^2 + e_i^2)}{2} \\right) \\frac{a}{a_i^2} - \\frac{a}{a_i^2} \\left( 1 + \\frac{3}{8} \\frac{a^2}{a_i^2} \\right) + \\sin^2 \\frac{t}{2} \\frac{a}{a_i^2} \\left( 1 + \\frac{33}{8} \\frac{a^2}{a_i^2} \\right) \\right\\}\n\\]\n\n\\[\n+ \\left( \\frac{(e^2 + e_i^2)}{2} \\frac{a}{a_i^2} \\left( 1 - \\frac{3}{2} \\frac{a^2}{a_i^2} \\right) \\right) \\cos t + \\left\\{ - \\frac{3}{4} \\frac{a^2}{a_i^2} + \\frac{3}{2} \\sin^2 \\frac{t}{2} \\frac{a^2}{a_i^2} + \\frac{15}{8} \\frac{(e^2 + e_i^2)}{a_i^2} \\right\\} \\cos 2 t\n\\]\n\\[ + \\left\\{ -\\frac{5}{8} \\frac{a^3}{a_i^4} + \\frac{15}{8} \\sin^2 \\frac{t}{2} \\frac{a^3}{a_i^4} + \\frac{15}{4} (e^2 + e_i^2) \\frac{a^2}{a_i^4} \\right\\} \\cos 3t + \\frac{15}{16} \\frac{a^3}{a_i^4} \\cos (t - x) \\]\n\n\\[ + \\frac{a^3}{2a_i^3} e \\cos x + \\frac{3}{16} \\frac{a^3}{a_i^4} e \\cos t + x) - \\frac{3}{4} \\frac{a^2}{a_i^3} e \\cos (2t + x) - \\frac{15}{16} \\frac{a^3}{a_i^4} e \\cos (3t + x) \\]\n\n\\[ + \\frac{9}{4} \\frac{a^2}{a_i^3} e \\cos (2t - x) + \\frac{45}{16} \\frac{a^3}{a_i^4} e \\cos (3t - x) - \\left\\{ \\frac{1}{a_i} + \\frac{3}{4} \\frac{a^2}{a_i^3} \\right\\} e_i \\cos z \\]\n\n\\[ - \\frac{3}{8} \\frac{a^3}{a_i^4} e_i \\cos (t + z) + \\frac{3}{8} \\frac{a^3}{a_i^4} e_i \\cos (2t + z) + \\frac{5}{8} \\frac{a^3}{a_i^4} e_i \\cos (3t + z) \\]\n\n\\[ - \\frac{9}{8} \\frac{a^3}{a_i^4} e_i \\cos (t - z) - \\frac{21}{8} \\frac{a^3}{a_i^4} e_i \\cos (2t - z) - \\frac{25}{8} \\frac{a^3}{a_i^4} e_i \\cos (3t - z) \\]\n\n\\[ - \\frac{15}{8} \\frac{a^2}{a_i^3} e^2 \\cos (2t - 2x) - \\frac{33}{64} \\frac{a^3}{a_i^4} e^2 \\cos (t - 2x) + \\frac{a^2}{8a_i^3} e^2 \\cos 2x \\]\n\n\\[ + \\frac{9}{64} \\frac{a^3}{a_i^4} e^2 \\cos (t + 2x) - \\frac{3}{4} \\frac{a^2}{a_i^3} e^2 \\cos (2t + 2x) - \\frac{75}{64} \\frac{a^3}{a_i^4} e^2 \\cos (3t + 2x) \\]\n\n\\[ - \\frac{285}{64} \\frac{a^3}{a_i^4} e^2 \\cos (3t - 2x) + \\frac{3}{4} \\frac{a^2}{a_i^3} ee_i \\cos (x - z) + \\frac{9}{16} \\frac{a^3}{a_i^4} ee_i \\cos (t - z + x) \\]\n\n\\[ - \\frac{27}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t - z + x) + \\frac{15}{16} \\frac{a^3}{a_i^4} ee_i \\cos (t + z - x) - \\frac{9}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t + z - x) \\]\n\n\\[ - \\frac{45}{16} \\frac{a^3}{a_i^4} ee_i \\cos (3t + z - x) + \\frac{45}{16} \\frac{a^3}{a_i^4} ee_i \\cos (t - z - x) + \\frac{3}{4} \\frac{a^2}{a_i^3} ee_i \\cos (x + z) \\]\n\n\\[ + \\frac{3}{16} \\frac{a^3}{a_i^4} ee_i \\cos (t + z + x) + \\frac{3}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t + z + x) + \\frac{15}{16} \\frac{a^3}{a_i^4} ee_i \\cos (3t + z + x) \\]\n\n\\[ + \\frac{63}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t - z - x) - \\frac{225}{16} \\frac{a^3}{a_i^4} ee_i \\cos (3t - z - x) + \\left\\{ \\frac{1}{a_i} - \\frac{9}{8} \\frac{a^2}{a_i^3} \\right\\} ee_i^2 \\cos 2z \\]\n\n\\[ - \\frac{33}{64} \\frac{a^3}{a_i^4} ee_i^2 \\cos (t + 2z) - \\frac{5}{64} \\frac{a^3}{a_i^4} ee_i^2 \\cos (3t + 2z) - \\frac{159}{64} \\frac{a^3}{a_i^4} ee_i^2 \\cos (t - 2z) \\]\n\n\\[ - \\frac{51}{8} \\frac{a^2}{a_i^3} ee_i^2 (2t - 2z) - \\frac{635}{64} \\frac{a^3}{a_i^4} ee_i^2 \\cos (3t - 2z) - \\frac{9}{4} \\frac{a^3}{a_i^4} \\sin^2 \\frac{t}{2} \\cos (t - 2y) \\]\n\\[\n-\\frac{3}{2} \\frac{a^3}{a_i^3} \\sin^2 \\frac{t}{2} \\cos 2y - \\frac{3}{2} \\frac{a^3}{a_i^3} \\sin^2 \\frac{t}{2} \\cos (2t - 2y) - \\frac{15}{8} \\frac{a^3}{a_i^4} \\cos (t + 2y)\n\\]\n\nLet \\( \\frac{a}{r} = 1 + e \\cos \\left( n(1 + k)t + \\varepsilon - \\varpi \\right) + e^2 \\left( 1 + r_{48} \\right) \\cos \\left( 2n(1 + k_3)t + 2\\varepsilon - 2\\varpi \\right) \\)\n\n\\[ + e^3 r_{59} \\cos \\left( 2n(1 + k_{10})t + 2\\varepsilon - 2\\varpi \\right) + r_0 + r_1 \\cos (nt + n_t t + \\varepsilon - \\varepsilon_t) \\]\n\n\\[ + r_2 \\cos (2nt - 2n_t t + 2\\varepsilon - 2\\varepsilon_t) + \\&c. + e r_6 \\cos (n_t t + \\varepsilon_t - \\varpi) + \\&c. \\]\n\n\\[\n\\frac{d^2 r^2}{dt^2} - \\frac{\\mu}{r} + \\frac{\\mu}{a} + 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) = 0\n\\]\n\nLet \\( \\delta \\cdot \\frac{1}{r} \\) denote that part of \\( \\frac{1}{r} \\) which depends on the first power of the disturbing force. It is more simple to obtain \\( \\delta \\cdot \\frac{1}{r} \\) from the above differential equation than \\( \\delta r \\); and the circumstance that the elliptic value of \\( r^3 \\) does not contain any term \\( e^2 \\cos (2nt + 2\\varepsilon - 2\\varpi) \\), gives an additional facility.\n\n\\[\n\\frac{d^2 r^3 \\delta \\cdot \\frac{1}{r}}{dt^2} - \\mu \\delta \\cdot \\frac{1}{r} + 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) = 0\n\\]\n\nWhen the disturbing force is neglected\n\n\\[\nr^3 = a^3 \\left\\{ 1 + 3e^2 \\left( 1 + \\frac{e^2}{12} \\right) - 3e \\left( 1 + \\frac{3}{8} e^2 \\right) \\cos (nt + \\varepsilon - \\varpi) - \\frac{3}{8} e^4 \\cos (2nt + 2\\varepsilon - 2\\varpi) \\right. \\]\n\n\\[\n+ \\frac{e^3}{8} \\cos (3nt + 3\\varepsilon - 3\\varpi) + \\frac{e^4}{8} \\cos (4nt + 4\\varepsilon - 4\\varpi) \\}\n\\]\n\nIntegrating the above differential equation by the method of indeterminate co-efficients, \\( q_n \\) being the co-efficient of the \\( n \\)th argument in the development of \\( 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) \\),\n\n\\[- r_0 + \\frac{m_t}{\\mu} a q_0 = 0 \\]\n\n\\[\n\\frac{(n - n_t)^2}{n^2} \\left\\{ (1 + 3e^2) r_1 - \\frac{3}{2} e^2 (r_6 + r_8) \\right\\} - r_1 + \\frac{m_t}{\\mu} a q_1 = 0\n\\]\n\n\\[\n\\frac{(2n - 2n_t)^2}{n^2} \\left\\{ (1 + 3e^2) r_2 - \\frac{3}{2} e^2 (r_9 + r_{13}) \\right\\} - r_2 + \\frac{m_t}{\\mu} a q_2 = 0\n\\]\n\\[\n\\frac{(3n-3n_i)^2}{n^2} \\left\\{ (1 + 3e^2) r_3 - \\frac{3e^2}{2} (r_{10} + r_{13}) \\right\\} - r_3 + \\frac{m_i}{\\mu} a q_3 = 0\n\\]\n\n\\[\n\\frac{(4n-4n_i)^2}{n^2} \\left\\{ (1 + 3e^2) r_4 - \\frac{3e^2}{2} (r_{11} + r_{14}) \\right\\} - r_4 + \\frac{m_i}{\\mu} a q_4 = 0\n\\]\n\n\\[\n\\frac{(5n-5n_i)^2}{n^2} (1 + 3e^2) r_5 - r_5 + \\frac{m_i}{\\mu} a q_5 = 0\n\\]\n\n\\[\n\\frac{n^2}{n_i^2} \\left\\{ r_6 - \\frac{3}{2} r_1 \\right\\} - r_6 + \\frac{m_i}{\\mu} a q_6 = 0\n\\]\n\n\\[\n*(2n-n_i)^2 \\left\\{ r_8 - \\frac{3}{2} r_1 \\right\\} - r_8 + \\frac{m_i}{\\mu} a q_8 = 0\n\\]\n\n\\[\n(3n-2n_i)^2 \\left\\{ r_9 - \\frac{3}{2} r_2 \\right\\} - r_9 + \\frac{m_i}{\\mu} a q_9 = 0\n\\]\n\n\\[\n(4n-3n_i)^2 \\left\\{ r_{10} - \\frac{3}{2} r_3 \\right\\} - r_{10} + \\frac{m_i}{\\mu} a q_{10} = 0\n\\]\n\n\\[\n(5n-4n_i)^2 \\left\\{ r_{11} - \\frac{3}{2} r_4 \\right\\} - r_{11} + \\frac{m_i}{\\mu} a q_{11} = 0\n\\]\n\n\\[\n(n-2n_i)^2 \\left\\{ r_{12} - \\frac{3}{2} r_3 \\right\\} - r_{12} + \\frac{m_i}{\\mu} a q_{12} = 0\n\\]\n\n\\[\n(2n-3n_i)^2 \\left\\{ r_{13} - \\frac{3}{2} r_3 \\right\\} - r_{13} + \\frac{m_i}{\\mu} a q_{13} = 0\n\\]\n\n\\[\n(3n-4n_i)^2 \\left\\{ r_{14} - \\frac{3}{2} r_4 \\right\\} - r_{14} + \\frac{m_i}{\\mu} a q_{14} = 0\n\\]\n\n\\[\n\\frac{n^2}{n_i^2} r_{15} - r_{15} + \\frac{m_i}{\\mu} a q_{15} = 0\n\\]\n\n\\[\n(2n-n_i)^2 r_{17} - r_{17} + \\frac{m_i}{\\mu} a q_{17} = 0\n\\]\n\n\\[\n(3n-2n_i)^2 r_{18} - r_{18} + \\frac{m_i}{\\mu} a q_{18} = 0\n\\]\n\n\\[\n(4n-3n_i)^2 r_{19} - r_{19} + \\frac{m_i}{\\mu} a q_{19} = 0\n\\]\n\n\\[\n(n-3n_i)^2 r_{20} - r_{20} + \\frac{m_i}{\\mu} a q_{20} = 0\n\\]\n\n\\[\n(2n-3n_i)^2 r_{21} - r_{21} + \\frac{m_i}{\\mu} a q_{21} = 0\n\\]\n\n\\[\n(3n-4n_i)^2 r_{22} - r_{22} + \\frac{m_i}{\\mu} a q_{22} = 0\n\\]\n\n* For the determination of the quantity \\( k \\), see p. 50.\n\nMDCCCXXXI.\nEquations which serve to determine the coefficients of the inequalities of the reciprocal of the radius vector.\n\n\\[\n\\frac{(4n - 5n_i)^2}{n^2} r_{23} - r_{23} + \\frac{m_i}{\\mu} a q_{23} = 0\n\\]\n\n\\[\n\\frac{4n_i^2}{n^2} \\left\\{ r_{24} - \\frac{3}{2} r_{12} \\right\\} - r_{24} + \\frac{m_i}{\\mu} a q_{24} = 0\n\\]\n\n\\[\n\\frac{(n + n_i)^2}{n^2} \\left\\{ r_{25} - \\frac{3}{2} r_6 \\right\\} - r_{25} + \\frac{m_i}{\\mu} a q_{25} = 0\n\\]\n\n\\[\n\\frac{(3n - n_i)^2}{n^2} \\left\\{ r_{27} - \\frac{3}{2} r_8 \\right\\} - r_{27} + \\frac{m_i}{\\mu} a q_{27} = 0\n\\]\n\n\\[\n\\frac{(4n - 2n_i)^2}{n^2} \\left\\{ r_{28} - \\frac{3}{2} r_9 \\right\\} - r_{28} + \\frac{m_i}{\\mu} a q_{28} = 0\n\\]\n\n\\[\n\\frac{(5n - 3n_i)^2}{n^2} \\left\\{ r_{29} - \\frac{3}{2} r_{10} \\right\\} - r_{29} + \\frac{m_i}{\\mu} a q_{29} = 0\n\\]\n\n\\[\n\\frac{(6n - 4n_i)^2}{n^2} \\left\\{ r_{30} - \\frac{3}{2} r_{11} \\right\\} - r_{30} + \\frac{m_i}{\\mu} a q_{30} = 0\n\\]\n\n\\[\n\\frac{(7n - 5n_i)^2}{n^2} r_{31} - r_{31} + \\frac{m_i}{\\mu} a q_{31} = 0\n\\]\n\n\\[\n\\frac{(n - 3n_i)^2}{n^2} \\left\\{ r_{32} - \\frac{3}{2} r_{13} \\right\\} - r_{32} + \\frac{m_i}{\\mu} a q_{32} = 0\n\\]\n\n\\[\n\\frac{(2n - 4n_i)^2}{n^2} \\left\\{ r_{33} - \\frac{3}{2} r_{14} \\right\\} - r_{33} + \\frac{m_i}{\\mu} a q_{33} = 0\n\\]\n\n\\[\n\\frac{(3n - 5n_i)^2}{n^2} r_{34} - r_{34} + \\frac{m_i}{\\mu} a q_{34} = 0\n\\]\n\n\\[\n\\frac{(n - n_i)^2}{n^2} \\left\\{ r_{35} - \\frac{3}{2} r_{15} \\right\\} - r_{35} + \\frac{m_i}{\\mu} a q_{35} = 0\n\\]\n\n\\[\n\\frac{(2n - 2n_i)^2}{n^2} \\left\\{ r_{36} - \\frac{3}{2} r_{20} \\right\\} - r_{36} + \\frac{m_i}{\\mu} a q_{36} = 0\n\\]\n\n\\[\n\\frac{(3n - 3n_i)^2}{n^2} \\left\\{ r_{37} - \\frac{3}{2} r_{21} \\right\\} - r_{37} + \\frac{m_i}{\\mu} a q_{37} = 0\n\\]\n\n\\[\n\\frac{(4n - 4n_i)^2}{n^2} \\left\\{ r_{38} - \\frac{3}{2} r_{22} \\right\\} - r_{38} + \\frac{m_i}{\\mu} a q_{38} = 0\n\\]\n\n\\[\n\\frac{(5n - 5n_i)^2}{n^2} \\left\\{ r_{39} - \\frac{3}{2} r_{23} \\right\\} - r_{39} + \\frac{m_i}{\\mu} a q_{39} = 0\n\\]\n\n\\[\n\\frac{(6n - 6n_i)^2}{n^2} r_{40} - r_{40} + \\frac{m_i}{\\mu} a q_{40} = 0\n\\]\n\n\\[\n- r_{41} + \\frac{m_i}{\\mu} a q_{41} = 0\n\\]\n\n\\[\n\\frac{(n - n_i)^2}{n^2} \\left\\{ r_{42} - \\frac{3}{2} r_{17} \\right\\} - r_{42} + \\frac{m_i}{\\mu} a q_{42} = 0\n\\]\n\\[\n\\frac{(2n - 2n_i)^2}{n^2} \\left\\{ r_{43} - \\frac{3}{2} r_{18} \\right\\} - r_{43} + \\frac{m_i}{\\mu} a q_{43} = 0\n\\]\n\n\\[\n\\frac{(3n - 3n_i)^2}{n^2} \\left\\{ r_{44} - \\frac{3}{2} r_{19} \\right\\} - r_{44} + \\frac{m_i}{\\mu} a q_{44} = 0\n\\]\n\n\\[\n\\frac{(4n - 4n_i)^2}{n^2} \\left\\{ r_{45} - r_{45} + \\frac{m_i}{\\mu} a q_{45} = 0\n\\]\n\n\\[\n\\frac{4n_i^2}{n^2} \\left\\{ r_{46} - \\frac{3}{2} r_{20} \\right\\} - r_{46} + \\frac{m_i}{\\mu} a q_{46} = 0\n\\]\n\n\\[\n\\frac{(n + n_i)^2}{n^2} \\left\\{ r_{47} - \\frac{3}{2} r_{15} \\right\\} - r_{47} + \\frac{m_i}{\\mu} a q_{47} = 0\n\\]\n\n\\[\n\\frac{(3n - 3n_i)^2}{n^2} \\left\\{ r_{49} - \\frac{3}{2} r_{17} \\right\\} - r_{49} + \\frac{m_i}{\\mu} a q_{49} = 0\n\\]\n\n\\[\n\\frac{(4n - 2n_i)^2}{n^2} \\left\\{ r_{50} - \\frac{3}{2} r_{18} \\right\\} - r_{50} + \\frac{m_i}{\\mu} a q_{50} = 0\n\\]\n\n\\[\n\\frac{(5n - 5n_i)^2}{n^2} \\left\\{ r_{51} - \\frac{3}{2} r_{19} \\right\\} - r_{51} + \\frac{m_i}{\\mu} a q_{51} = 0\n\\]\n\n\\[\n\\frac{(6n - 4n_i)^2}{n^2} \\left\\{ r_{52} - r_{52} + \\frac{m_i}{\\mu} a q_{52} = 0\n\\]\n\n\\[\n\\frac{(n - n_i)^2}{n^2} \\left\\{ r_{53} - \\frac{3}{2} r_{21} \\right\\} - r_{53} + \\frac{m_i}{\\mu} a q_{53} = 0\n\\]\n\n\\[\n\\frac{(2n - 4n_i)^2}{n^2} \\left\\{ r_{54} - \\frac{3}{2} r_{22} \\right\\} - r_{54} + \\frac{m_i}{\\mu} a q_{54} = 0\n\\]\n\n\\[\n\\frac{(3n - 5n_i)^2}{n^2} \\left\\{ r_{55} - \\frac{3}{2} r_{23} \\right\\} - r_{55} + \\frac{m_i}{\\mu} a q_{55} = 0\n\\]\n\n\\[\n\\frac{(4n - 6n_i)^2}{n^2} r_{56} - r_{56} + \\frac{m_i}{\\mu} a q_{56} = 0\n\\]\n\n\\[\n\\frac{4n_i^2}{n^2} r_{57} - r_{57} + \\frac{m_i}{\\mu} a q_{57} = 0\n\\]\n\n\\[\n\\frac{(m + n_i)^2}{n^2} r_{58} - r_{58} + \\frac{m_i}{\\mu} a q_{58} = 0\n\\]\n\n\\[\n\\frac{(3n - n_i)^2}{n^2} r_{60} - r_{60} + \\frac{m_i}{\\mu} a q_{60} = 0\n\\]\n\n\\[\n\\frac{(4n - 2n_i)^2}{n^2} r_{61} - r_{61} + \\frac{m_i}{\\mu} a q_{61} = 0\n\\]\n\n\\[\n\\frac{(5n - 3n_i)^2}{n^2} r_{62} - r_{62} + \\frac{m_i}{\\mu} a q_{62} = 0\n\\]\n\n\\[\n\\frac{(n - 3n_i)^2}{n^2} r_{63} - r_{63} + \\frac{m_i}{\\mu} a q_{63} = 0\n\\]\nEquations which serve to determine the coefficients of the inequalities of the reciprocal of the radius vector.\n\n\\[\n\\frac{(2n - 4n)}{n^2} r_{64} - r_{64} + \\frac{m_1}{\\mu} a q_{64} = 0\n\\]\n\n\\[\n\\frac{(3n - 5n)}{n^2} r_{65} - r_{65} + \\frac{m_1}{\\mu} a q_{65} = 0\n\\]\n\n\\[\n\\frac{(4n - 6n)}{n^2} r_{66} - r_{66} + \\frac{m_1}{\\mu} a q_{66} = 0\n\\]\n\n\\[\n\\frac{(5n - 7n)}{n^2} r_{67} - r_{67} + \\frac{m_1}{\\mu} a q_{67} = 0\n\\]\n\n\\[\n\\frac{(n + n)}{n^2} r_{68} - r_{68} + \\frac{m_1}{\\mu} a q_{68} = 0\n\\]\n\nIn order to obtain the values of the coefficients of the inequalities from these equations when the cubes of the eccentricities are neglected, as has been the case throughout, the values of \\(r_0, r_1, r_2, r_3, r_4, r_5\\), found from the first six equations by neglecting the terms multiplied by \\(e^2\\), may be substituted in the succeeding equations, which will then serve to determine \\(r_6, r_8, r_9, \\&c.\\) and these values of \\(r_6, r_8, r_9, \\&c.\\) being substituted in the terms multiplied by \\(e^2\\), of the equations which determine \\(r_1, r_2, r_3, r_4, \\&c.\\) more accurate values of those quantities may be obtained. All the other coefficients of which the general symbol is \\(r\\) with a numerical index at foot, may then be obtained in succession without any difficulty.\n\n\\[\n\\frac{d\\lambda}{dt} = \\frac{h}{r^2} - \\frac{1}{r^2} \\int \\frac{dR}{d\\lambda} dt\n\\]\n\nLet \\(r\\) denote the elliptic value of \\(r\\), then\n\n\\[\n\\frac{d\\lambda}{dt} = \\frac{h}{r^2} + \\frac{2h}{r} \\delta \\cdot \\frac{1}{r} - \\frac{1}{r^2} \\int \\frac{dR}{d\\lambda} dt\n\\]\n\n\\[\n\\frac{a^2}{r^2} = 1 + \\frac{e^2}{2} + \\frac{3}{8} e^4 + 2e \\left(1 + \\frac{3e^2}{8}\\right) \\cos (nt - \\varphi) + \\frac{5e^2}{2} \\left(1 + \\frac{2}{15} e^2\\right) \\cos (2nt - 2\\varphi)\n\\]\n\n\\[\n+ \\frac{13}{4} e^3 \\cos (3nt - 3\\varphi) + \\frac{103}{24} e^4 \\cos (4nt - 4\\varphi)\n\\]\n\n\\[\n\\frac{a}{r} = 1 + e \\left(1 - \\frac{e^2}{8}\\right) \\cos (nt - \\varphi) + e^2 \\left(1 - \\frac{e^2}{3}\\right) \\cos (2nt - 2\\varphi)\n\\]\n\n\\[\n+ \\frac{9}{8} e^3 \\cos (3nt - 3\\varphi) + \\frac{4}{3} e^4 \\cos (4nt - 4\\varphi)\n\\]\nLet $R_n$ denote the coefficient of the cosine in the development of $R$ which corresponds to the number $n$ multiplied by $e$ or $e_i$, &c. Thus\n\n$$R_n = \\left\\{ -\\frac{3a}{a_i^2} + \\frac{3a}{2a_i^2} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} - \\frac{a}{4a_i^2} b_{3,2} \\right\\} \\quad \\text{See p.31}$$\n\nSince $\\frac{dR}{d\\lambda} = \\frac{dR}{d\\varepsilon}$, $\\frac{dR}{d\\varepsilon}$ being the differential of $R$ with respect to $\\varepsilon$, considering $\\varepsilon - \\varpi$ and $\\varepsilon - \\nu$ constant\n\n$$\\lambda = n \\left\\{ 1 + 2r_0 \\right\\} t + \\varepsilon$$\n\n$$+ \\left\\{ 2 \\left( r_1 \\left( 1 - \\frac{e^2}{2} \\right) + \\frac{e^2}{2} (r_6 + r_8) \\right) - \\frac{m_i}{\\mu} \\left\\{ \\left( 1 + \\frac{e^2}{2} \\right) \\frac{a R_1 n}{(n-n_i)} - \\frac{e^2}{n_i} a R_6 n + \\frac{e^2 a R_8 n}{(2n-n_i)} \\right\\} \\right\\}$$\n\n$$\\frac{n}{(n-n_i)} \\sin (nt - n_i t + \\varepsilon - \\varepsilon_i)$$\n\n$$+ \\left\\{ 2 \\left( r_2 \\left( 1 - \\frac{e^2}{2} \\right) + \\frac{e^2}{2} (r_9 + r_{12}) \\right) - \\frac{m_i}{\\mu} \\left\\{ \\left( 1 + \\frac{e^2}{2} \\right) \\frac{a R_2 n}{(n-n_i)} + \\frac{2 e^2 a R_9 n}{(3n-2n_i)} + \\frac{2 e^2 a R_{12} n}{(n-2n_i)} \\right\\} \\right\\}$$\n\n$$\\frac{n}{(2n-2n_i)} \\sin (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i)$$\n\n$$+ \\left\\{ 2 \\left( r_3 \\left( 1 - \\frac{e^2}{2} \\right) + \\frac{e^2}{2} (r_{10} + r_{13}) \\right) - \\frac{m_i}{\\mu} \\left\\{ \\left( 1 + \\frac{e^2}{2} \\right) \\frac{a R_3 n}{(n-n_i)} + \\frac{3 e^2 a R_{10} n}{(4n-3n_i)} + \\frac{3 e^2 a R_{13} n}{(2n-3n_i)} \\right\\} \\right\\}$$\n\n$$\\frac{n}{(3n-3n_i)} \\sin (3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i)$$\n\n$$+ \\left\\{ 2 \\left( r_4 \\left( 1 - \\frac{e^2}{2} \\right) + \\frac{e^2}{2} (r_{11} + r_{14}) \\right) - \\frac{m_i}{\\mu} \\left\\{ \\left( 1 + \\frac{e^2}{2} \\right) \\frac{a R_4 n}{(n-n_i)} + \\frac{4 e^2 a R_{11} n}{(5n-4n_i)} + \\frac{4 e^2 a R_{14} n}{(3n-4n_i)} \\right\\} \\right\\}$$\n\n$$\\frac{n}{(4n-4n_i)} \\sin (4nt - 4n_i t + 4\\varepsilon - 4\\varepsilon_i)$$\n\n$$+ \\left\\{ 2 \\left( r_6 + \\frac{r_1}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ -\\frac{a R_6 n}{n_i} + \\frac{a R_1 n}{(n-n_i)} \\right\\} \\right\\} \\frac{e}{n_i} \\sin (nt + \\varepsilon_i - \\varpi)$$\n\n$$+ \\left\\{ 2 \\left( r_8 + \\frac{r_1}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{a R_8 n}{(2n-n_i)} + \\frac{a R_1 n}{(n-n_i)} \\right\\} \\right\\} \\frac{ne}{(n-n_i)} \\sin (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi)$$\n\n$$+ \\left\\{ 2 \\left( r_9 + \\frac{r_2}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{2 a R_9 n}{(3n-2n_i)} + \\frac{a R_2 n}{(n-n_i)} \\right\\} \\right\\} \\frac{ne}{(3n-2n_i)} \\sin (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi)$$\n\n$$+ \\left\\{ 2 \\left( r_{10} + \\frac{r_3}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3 a R_{10} n}{(4n-3n_i)} + \\frac{a R_3 n}{(n-n_i)} \\right\\} \\right\\} \\frac{ne}{(4n-3n_i)} \\sin (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi)$$\nExpression for the longitude.\n\n\\[ + \\left\\{ 2 \\left( r_{11} + \\frac{r_3}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{4aR_{11}n}{(5n-4n_i)} + \\frac{aR_{11}n}{(n-n_i)} \\right] \\right\\} \\frac{n e}{(5n-4n_i)} \\sin (5nt - 4n_i t + 5\\varepsilon - 4\\varepsilon_i - \\varpi) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{12} + \\frac{r_2}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{2aR_{12}n}{(n-2n_i)} + \\frac{aR_{12}n}{(n-n_i)} \\right] \\right\\} \\frac{n e}{(n-2n_i)} \\sin (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{13} + \\frac{r_3}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{3aR_{13}n}{(2n-3n_i)} + \\frac{aR_{13}n}{(n-n_i)} \\right] \\right\\} \\frac{n e}{(2n-3n_i)} \\sin (2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{14} + \\frac{r_4}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{4aR_{14}n}{(3n-4n_i)} + \\frac{aR_{14}n}{(n-n_i)} \\right] \\right\\} \\frac{n e}{(3n-4n_i)} \\sin (3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi) \\]\n\n\\[ + 2r_{15} \\frac{n e_i}{n_i} \\sin (n_i t + \\varepsilon_i - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{17} - \\frac{2m_i aR_{17}n}{\\mu (2n-n_i)} \\right\\} \\frac{n e_i}{(2n-n_i)} \\sin (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{18} - \\frac{3m_i aR_{18}n}{\\mu (3n-2n_i)} \\right\\} \\frac{n e_i}{(3n-2n_i)} \\sin (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{19} - \\frac{4aR_{19}n}{\\mu (4n-3n_i)} \\right\\} \\frac{n e_i}{(4n-3n_i)} \\sin (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{20} - \\frac{m_i aR_{20}n}{\\mu (n-2n_i)} \\right\\} \\frac{n e_i}{(n-2n_i)} \\sin (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{21} - \\frac{2m_i aR_{21}n}{\\mu (2n-3n_i)} \\right\\} \\frac{n e_i}{(2n-3n_i)} \\sin (2nt - 3n_i t + 2\\varepsilon - 2\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{22} - \\frac{3m_i aR_{22}n}{\\mu (3n-4n_i)} \\right\\} \\frac{n e_i}{(3n-4n_i)} \\sin (3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{23} - \\frac{4m_i aR_{23}n}{\\mu (4n-5n_i)} \\right\\} \\frac{n e_i}{(4n-5n_i)} \\sin (4nt - 5n_i t + 4\\varepsilon - 5\\varepsilon_i + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{24} + \\frac{r_{12}}{2} \\right) \\frac{m_i}{\\mu} \\left[ \\frac{aR_{24}n}{m_i} + \\frac{2aR_{12}n}{(n-2n_i)} \\right] \\right\\} \\frac{n e^2}{n_i} \\sin (2n_i t + 2\\varepsilon_i - 2\\varpi) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{25} + \\frac{r_6}{2} + \\frac{r_1}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{aR_{25}n}{(n+n_i)} - \\frac{aR_6n}{n_i} + \\frac{5aR_1n}{4(n-n_i)} \\right] \\right\\} \\frac{n e^2}{n_i} \\sin (nt + n_i t + \\varepsilon + \\varepsilon_i - 2\\varpi) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{27} + \\frac{r_8}{2} + \\frac{r_1}{2} \\right) - \\frac{m_i}{\\mu} \\left[ \\frac{aR_{27}n}{(3n-n_i)} + \\frac{aR_8n}{(2n-n_i)} + \\frac{5aR_1n}{4(n-n_i)} \\right] \\right\\} \\frac{n e^2}{(3n-n_i)} \\sin (3nt - n_i t + 3\\varepsilon - \\varepsilon_i - 2\\varpi) \\]\n\\[\n+ \\left\\{ 2 \\left( r_{28} + \\frac{r_9}{2} + \\frac{r_9}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{2aR_{38}n}{4n-2n_i} + \\frac{2aR_{38}n}{3n-2n_i} + \\frac{5aR_{38}n}{4(n-n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e^2}{(4n-2n_i)} \\sin (4nt-2n_i t + 4\\varepsilon - 2\\varepsilon_i - 2\\varpi)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{30} + \\frac{r_{10}}{2} + \\frac{r_3}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{30}n}{5n-3n_i} + \\frac{3aR_{10}n}{4n-3n_i} + \\frac{5aR_{3}n}{4(n-n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e^2}{(5n-3n_i)} \\sin (5nt-3n_i t + 5\\varepsilon - 3\\varepsilon_i - 2\\varpi)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{30} + \\frac{r_{11}}{2} + \\frac{r_2}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{30}n}{6n-4n_i} + \\frac{4aR_{11}n}{5n-4n_i} + \\frac{5aR_{4}n}{4(n-n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e^2}{(6n-4n_i)} \\sin (6nt-4n_i t + 6\\varepsilon - 4\\varepsilon_i - 2\\varpi)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{32} + \\frac{r_{13}}{2} + \\frac{r_3}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{32}n}{(n-3n_i)} + \\frac{3aR_{13}n}{(2n-3n_i)} + \\frac{5aR_{3}n}{4(n-n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e^2}{(n-3n_i)} \\sin (nt-3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\varpi)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{33} + \\frac{r_{14}}{2} + \\frac{r_2}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{33}n}{(2n-4n_i)} + \\frac{4aR_{14}n}{(3n-4n_i)} + \\frac{5aR_{4}n}{4(n-n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e^2}{(2n-4n_i)} \\sin (2nt-4n_i t + 2\\varepsilon - 4\\varepsilon_i + 2\\varpi)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{36} + \\frac{r_{20}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{aR_{36}n}{(2n-2n_i)} + \\frac{aR_{20}n}{(n-2n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(2n-2n_i)} \\sin (2nt-2n_i t + 2\\varepsilon - 2\\varepsilon_i - \\varpi + \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{37} + \\frac{r_{21}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{2aR_{37}n}{(3n-3n_i)} + \\frac{2aR_{21}n}{(2n-3n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(3n-3n_i)} \\sin (3nt-3n_i t + 3\\varepsilon - 3\\varepsilon_i - \\varpi + \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{38} + \\frac{r_{22}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{38}n}{(4n-4n_i)} + \\frac{3aR_{22}n}{(3n-4n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(4n-4n_i)} \\sin (4nt-4n_i t + 4\\varepsilon - 4\\varepsilon_i - \\varpi + \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{39} + \\frac{r_{23}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{39}n}{(5n-5n_i)} + \\frac{4aR_{23}n}{(4n-5n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(5n-5n_i)} \\sin (5nt-5n_i t + 5\\varepsilon - 5\\varepsilon_i - \\varpi + \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{43} + \\frac{r_{18}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{43}n}{(2n-2n_i)} + \\frac{3aR_{18}n}{(3n-2n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(2n-2n_i)} \\sin (2nt-2n_i t + 2\\varepsilon - 2\\varepsilon_i + \\varpi - \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{44} + \\frac{r_{19}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{44}n}{(3n-3n_i)} + \\frac{4aR_{19}n}{(4n-3n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{(3n-3n_i)} \\sin (3nt-3n_i t + 3\\varepsilon - 3\\varepsilon_i + \\varpi - \\varpi_i)\n\\]\n\n\\[\n+ \\left\\{ 2 \\left( r_{46} + \\frac{r_{20}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{-aR_{46}n}{2n_i} + \\frac{aR_{20}n}{(n-2n_i)} \\right\\} \\right\\}\n\\]\n\n\\[\n\\frac{n e e_i}{2n_i} \\sin (2n_i t + 2\\varepsilon - \\varpi - \\varpi_i)\n\\]\nExpression for the longitude.\n\n\\[ + \\left\\{ 2 \\left( r_{49} + \\frac{r_{17}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{2aR_{49}n}{(3n-n_i)} + \\frac{2aR_{17}n}{(2n-n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(3n-n_i)} \\sin (3nt - n_i t + 3\\varepsilon - \\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{50} + \\frac{r_{18}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{50}n}{(4n-2n_i)} + \\frac{3aR_{18}n}{(3n-2n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(4n-2n_i)} \\sin (4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{51} + \\frac{r_{19}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{51}n}{(5n-3n_i)} + \\frac{4aR_{19}n}{(4n-3n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(5n-3n_i)} \\sin (5nt - 3n_i t + 5\\varepsilon - 3\\varepsilon_i - \\varpi - \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{53} + \\frac{r_{21}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{2aR_{53}n}{(n-3n_i)} + \\frac{2aR_{21}n}{(2n-3n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(n-3n_i)} \\sin (nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{54} + \\frac{r_{22}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{3aR_{54}n}{(2n-4n_i)} + \\frac{3aR_{22}n}{(3n-4n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(2n-4n_i)} \\sin (2nt - 4n_i t + 2\\varepsilon - 4\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2 \\left( r_{55} + \\frac{r_{23}}{2} \\right) - \\frac{m_i}{\\mu} \\left\\{ \\frac{4aR_{55}n}{(3n-5n_i)} + \\frac{4aR_{23}n}{(4n-5n_i)} \\right\\} \\right\\} \\frac{n e e_i}{(3n-5n_i)} \\sin (3nt - 5n_i t + 3\\varepsilon - 5\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{58} - \\frac{m_i}{\\mu} \\frac{aR_{58}n}{(n+n_i)} \\right\\} \\frac{n e e_i^2}{(n+n_i)} \\sin (nt + n_i t + \\varepsilon + \\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{59} - \\frac{m_i}{\\mu} \\frac{aR_{59}n}{\\mu} \\right\\} e_i^2 \\sin (2nt + 2\\varepsilon - 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{60} - \\frac{m_i}{\\mu} \\frac{3aR_{60}n}{(3n-n_i)} \\right\\} \\frac{n e e_i^2}{(3n-n_i)} \\sin (3nt - n_i t + 3\\varepsilon - \\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{61} - \\frac{m_i}{\\mu} \\frac{4aR_{61}n}{(4n-2n_i)} \\right\\} \\frac{n e e_i^2}{(4n-2n_i)} \\sin (4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{62} - \\frac{m_i}{\\mu} \\frac{5aR_{62}n}{(5n-3n_i)} \\right\\} \\frac{n e e_i^2}{(5n-3n_i)} \\sin (5nt - 3n_i t + 5\\varepsilon - 3\\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{63} - \\frac{m_i}{\\mu} \\frac{aR_{63}n}{(n-3n_i)} \\right\\} \\frac{n e e_i^2}{(n-3n_i)} \\sin (nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{64} - \\frac{m_i}{\\mu} \\frac{2aR_{64}n}{(2n-4n_i)} \\right\\} \\frac{n e e_i^2}{(2n-4n_i)} \\sin (2nt - 4n_i t + 2\\varepsilon - 4\\varepsilon_i + 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{65} - \\frac{m_i}{\\mu} \\frac{3aR_{65}n}{(3n-5n_i)} \\right\\} \\frac{n e e_i^2}{(3n-5n_i)} \\sin (3nt - 5n_i t + 3\\varepsilon - 5\\varepsilon_i + 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{66} - \\frac{m_i}{\\mu} \\frac{4aR_{66}n}{(4n-6n_i)} \\right\\} \\frac{n e e_i^2}{(4n-6n_i)} \\sin (4nt - 6n_i t + 4\\varepsilon - 6\\varepsilon_i + 2\\varpi_i) \\]\n\n\\[ + \\left\\{ 2r_{67} - \\frac{m_i}{\\mu} \\frac{5aR_{67}n}{(5n-7n_i)} \\right\\} \\frac{n e e_i^2}{(5n-7n_i)} \\sin (5nt - 7n_i t + 5\\varepsilon - 7\\varepsilon_i + 2\\varpi_i) \\]\nIn order to convert the coefficients of the inequalities of the longitude into sexagesimal seconds, they must be multiplied by \\( \\frac{1296000}{2\\pi} \\), the logarithm of which number is 5.3144251, the corresponding logarithm for centesimal seconds is 5.8038801.\n\n\\[ \\frac{dR}{dz} = m_i \\left[ \\frac{z_i}{r_i^3} + \\frac{z - z_i}{\\{ r^2 - 2rr_i \\cos(\\lambda - \\lambda_i) + r_i^2 \\}^{3/2}} \\right] \\]\n\nIf \\( z = 0 \\), and the products \\( e \\tan \\iota_p \\), \\( e_i \\tan \\iota_i \\) be neglected,\n\n\\[ \\frac{dR}{dz} = \\frac{m_i}{a_i^2} \\tan \\iota_i \\sin(n_i t + \\varepsilon_i - \\nu_i) \\{ 1 - b_{3,0} - b_{3,1} \\cos(nt - n_i t + \\varepsilon - \\varepsilon_i) - b_{3,2} \\cos(2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i) \\} \\]\n\n\\[ \\frac{a^2}{d^2s}{d^2t} + \\mu s + \\frac{m_i}{a_i^2} a^2 \\tan \\iota_i \\sin(n_i t + \\varepsilon_i - \\nu_i) \\{ 1 - b_{3,0} - b_{3,1} \\cos(nt - n_i t + \\varepsilon - \\varepsilon_i) - \\&c. \\} = 0 \\]\n\n\\[ s = -\\frac{m_i}{\\mu} \\frac{n^2}{(n - n_i)(n + n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i \\{ 1 - b_{3,0} \\} \\sin(n_i t + \\varepsilon_i - \\nu_i) \\]\n\n\\[ + i \\sin((1 + l)n_i t + \\varepsilon - \\nu_i) \\]\n\n\\[ -\\frac{m_i}{\\mu} \\frac{n^2}{4n_i(2n - 2n_i)} \\frac{a^2}{a_i^3} \\tan \\iota_i b_{3,1} \\sin(nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\nu_i) \\]\n\n\\[ -\\frac{m_i}{\\mu} \\frac{n^2}{2(n - n_i)(3n - n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,2} \\sin(2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\nu_i) \\]\n\n\\[ +\\frac{m_i}{\\mu} \\frac{n^2}{2(n - 3n_i)(3n - 3n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,2} \\sin(2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\nu_i) \\]\n\n\\[ -\\frac{m_i}{\\mu} \\frac{n^2}{2(2n - 2n_i)(4n - 2n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,3} \\sin(3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\nu_i) \\]\n\n\\[ +\\frac{m_i}{\\mu} \\frac{n^2}{2(2n - 4n_i)(4n - 4n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,4} \\sin(3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\nu_i) \\]\n\n\\[ -\\frac{m_i}{\\mu} \\frac{n^2}{2(3n - 3n_i)(5n - 3n_i)} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,4} \\sin(4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\nu_i) \\]\n\n\\[ l(2 + l)i = \\frac{m}{2\\mu_i} \\frac{a^2}{a_i^2} \\tan \\iota_i b_{3,1} \\]\n\n\\[ R_6 = -\\frac{3a}{2a_i^2} + \\frac{3a}{2a_i^2} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} - \\frac{a}{4a_i^2} b_{3,2} \\quad (\\text{See p. 31.}) \\]\n\n\\[ q_6 = -\\frac{3a}{2a_i^2} + \\frac{3a}{2a_i^2} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} - \\frac{a}{4a_i^2} b_{3,2} - \\frac{3}{2a_i^3} \\left\\{ \\frac{a}{a_i} b_{5,0} - \\frac{1}{2} b_{5,1} \\right\\} \\]\n\n\\[ + \\frac{3}{2} \\frac{a^3}{a_i^4} \\left\\{ \\frac{a}{a_i} b_{5,1} - b_{5,0} - \\frac{1}{2} b_{5,2} \\right\\} + \\frac{3}{4} \\frac{a^2}{a_i^3} \\left\\{ \\frac{a}{a_i} b_{5,2} - \\frac{1}{2} b_{5,1} - \\frac{1}{2} b_{5,3} \\right\\} \\]\n\nMDCCCXXXI.\n\\[ q_6 = -\\frac{3a}{2a_i^2} + \\frac{3}{2} \\frac{a}{a_i^2} b_{s,0} - \\frac{a^2}{2a_i^3} b_{s,1} - \\frac{a}{4a_i^3} b_{s,2} + \\frac{3}{2} \\frac{a^2}{a_i^3} \\left( 1 + \\frac{a^3}{a_i^3} \\right) b_{s,1} - 2 \\frac{a}{a_i} b_{s,0} - \\frac{a}{a_i^2} b_{s,2} \\]\n\n\\[ - \\frac{3a^3}{a_i^4} \\left\\{ b_{s,0} - \\frac{1}{2} b_{s,2} \\right\\} + \\frac{3}{8} \\frac{a^3}{a_i^4} \\left\\{ b_{s,1} - b_{s,3} \\right\\} \\]\n\n\\[ = -\\frac{3a}{2a_i^2} + \\frac{3}{2} \\frac{a}{a_i^2} b_{s,0} - \\frac{a^2}{2a_i^3} b_{s,1} - \\frac{a}{4a_i^3} b_{s,2} + \\frac{3}{2} \\frac{a^2}{a_i^3} b_{s,1} - \\frac{a^2}{a_i^3} b_{s,1} + \\frac{2}{4} \\frac{a}{a_i^2} b_{s,2} \\]\n\n\\[ = -\\frac{3a}{2a_i^2} + \\frac{3}{2} \\frac{a}{a_i^2} b_{s,0} - \\frac{a^2}{2a_i^3} b_{s,1} - \\frac{a}{4a_i^3} b_{s,2} \\]\n\nThe quantities of which the general symbol is \\( q \\), and which refer to the terms in the development of \\( R \\) multiplied by the eccentricities, admit of similar reductions; so that\n\n\\[ q_7 = -\\frac{a^2}{a_i^3} b_{s,0} + \\frac{a}{a_i^2} b_{s,1} \\]\n\n\\[ q_{16} = -\\frac{a}{2a_i^2} b_{s,2} \\]\n\nConsidering only the terms in \\( 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) \\), of which the arguments are \\( nt + \\varepsilon - \\varpi \\), and \\( nt + \\varepsilon - \\varpi_i \\)\n\n\\[ 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) = m_i q_7 e \\cos (nt + \\varepsilon - \\varpi) + m_i q_{16} e_i \\cos (nt + \\varepsilon - \\varpi_i) \\]\n\n\\[ = m_i q \\cos (nt + \\varepsilon - \\varpi_i) \\]\n\nprovided\n\n\\[ q \\cos \\varpi_1 = q_7 e \\cos \\varpi + q_{16} e_i \\cos \\varpi_i \\]\n\n\\[ q \\cos \\varpi_1 = q_7 e \\sin \\varpi + q_{16} e_i \\sin \\varpi_i \\]\n\nAnd if\n\n\\[ \\frac{a}{r} = 1 + r_0 + e \\cos \\left( n(1+k)t + \\varepsilon - \\varpi_1 \\right) + &c. \\]\n\n\\[ (1+k)^2 (1-3r_0) - 1 + \\frac{m_i a}{\\mu e} q = 0 \\]\n\n\\[ k = \\frac{3}{2} r_0 - \\frac{m_i a}{2\\mu e} q \\]\n\n\\[ r_0 = \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{s,0} - \\frac{a^2}{2a_i^2} b_{s,1} \\right\\} \\]\n\n\\[ q^3 = q_7^2 e^2 + 2q_7 q_{16} e e_i \\cos (\\varpi - \\varpi_i) + q_{16}^2 e_i^2 \\]\n\n\\[ q = q_7 e + q_{16} e_i \\cos (\\varpi - \\varpi_i) \\text{ nearly} \\]\n\n\\[ \\frac{q}{e} = -\\frac{a^2}{a_i^3} b_{s,0} + \\frac{a}{a_i^2} b_{s,1} - \\frac{a}{2a_i^2} e_i b_{s,2} \\cos (\\varpi - \\varpi_i) \\]\n\n\\[ k = \\frac{m_i}{\\mu} \\left\\{ \\frac{3a^3}{2a_i^3} b_{s,0} - \\frac{3a^2}{4a_i^2} b_{s,1} + \\frac{a^3}{2a_i^3} b_{s,0} - \\frac{a^2}{2a_i^2} b_{s,1} + \\frac{a^2}{4a_i^2} b_{s,2} e_i \\cos (\\varpi - \\varpi_i) \\right\\} \\]\n\n\\[ = \\frac{m_i}{\\mu} \\left\\{ \\frac{2a^3}{a_i^3} b_{s,0} - \\frac{5}{4} \\frac{a^2}{a_i^2} b_{s,1} + \\frac{a^2}{4a_i^2} b_{s,2} e_i \\cos (\\varpi - \\varpi_i) \\right\\} \\]\n\\[\n\\frac{d\\lambda}{dt} = n(1 + 2r_0)t + \\varepsilon + 2(1 + r_0)e \\cos(n(1 + k)t + \\varepsilon - \\varpi_i) \\\\\n- \\frac{m_i}{\\mu}aR_{16}e_i \\cos(n(1 + k)t + \\varepsilon - \\varpi_i)\n\\]\n\nand neglecting the square of the disturbing force\n\n\\[\n\\frac{d\\lambda}{dt} = n(1 + 2r_0)t + \\varepsilon + 2(1 + r_0)e \\cos(n(1 + k)t + \\varepsilon - \\varpi_i) \\\\\n- \\frac{m_i}{\\mu}aR_{16}e_i \\cos(n(1 + k)t + \\varepsilon - \\varpi_i)\n\\]\n\nIf\n\n\\[e(1 + k)\\cos\\varpi_2 = (1 + r_0)e\\cos\\varpi_1 - \\frac{m_i}{2\\mu}aR_{16}e_i \\cos\\varpi_i\\]\n\n\\[e(1 + k)\\sin\\varpi_2 = (1 + r_0)e\\sin\\varpi_1 - \\frac{m_i}{2\\mu}aR_{16}e_i \\sin\\varpi_i\\]\n\n\\[\\lambda = n(1 + 2r_0)t + \\varepsilon 2e \\sin(n(1 + k)t + \\varepsilon - \\varpi_2)\\]\n\n\\[e(1 + r_0)\\cos\\varpi_1 = e(1 + k)\\cos\\varpi_2 + \\frac{m_i}{2\\mu}aR_{16}e_i \\cos\\varpi_i\\]\n\n\\[e(1 + r_0)\\sin\\varpi_1 = e(1 + k)\\sin\\varpi_2 + \\frac{m_i}{2\\mu}aR_{16}e_i \\sin\\varpi_i\\]\n\n\\[e(1 + r_0) = e(1 + k)\\left\\{1 + \\frac{m_i}{2\\mu}aR_{16}e_i \\cos(\\varpi_2 - \\varpi_i)\\right\\}\\]\n\n\\[\\cos\\varpi_1 = \\cos\\varpi_2 \\left\\{1 - \\frac{m_i}{2\\mu}aR_{16}e_i \\cos(\\varpi_2 - \\varpi_i)\\right\\} + \\frac{m_i}{2\\mu}aR_{16}e_i \\cos\\varpi_i\\]\n\n\\[\\sin\\varpi_1 = \\sin\\varpi_2 \\left\\{1 - \\frac{m_i}{2\\mu}aR_{16}e_i \\cos(\\varpi_2 - \\varpi_i)\\right\\} + \\frac{m_i}{2\\mu}aR_{16}e_i \\sin\\varpi_i\\]\n\n\\[\\sin(\\varpi_2 - \\varpi_i) = \\frac{m_i}{2\\mu}aR_{16}e_i \\cos\\varpi_i\\]\n\ntherefore neglecting the square of the disturbing force\n\n\\[\\frac{e}{a} = \\frac{e}{a} \\left\\{1 + k - r_0 + \\frac{m_i}{2\\mu}aR_{16}e_i \\cos(\\varpi - \\varpi_i)\\right\\}\\]\n\n\\[= \\frac{e}{a} \\left\\{1 + \\frac{m_i}{\\mu} \\left\\{\\frac{2a^3}{a_i^3}b_{3,0} - \\frac{5}{4}\\frac{a^2}{a_i^2}b_{3,1} - \\frac{a^3}{a_i^3}b_{3,0} + \\frac{a^3}{2a_i^3}b_{3,1}\\right\\}\\right.\\]\n\n\\[\\left.+ \\frac{m_i}{\\mu} \\left\\{\\frac{a^2}{4a_i^2}b_{3,2} + \\frac{3}{4}\\frac{a^2}{a_i^2}b_{3,0} - \\frac{a}{4a_i}b_{3,1} - \\frac{a^2}{8a_i^2}b_{3,2}\\right\\}e_i \\cos(\\varpi - \\varpi_i)\\right\\}\\]\n\n\\[= \\frac{e}{a} \\left\\{1 + \\frac{m_i}{\\mu} \\left\\{\\frac{a^3}{a_i^3}b_{3,0} - \\frac{3}{4}\\frac{a^3}{a_i^3}b_{3,1} - \\frac{a^3}{4a_i^3}b_{3,1}e_i \\cos(\\varpi - \\varpi_i)\\right\\}\\right\\}\\]\nLet\n\n\\[ n \\left\\{ 1 + \\frac{2m_i}{\\mu} \\left( \\frac{a^3}{a_i^3} b_{3,o} - \\frac{a^2}{a_i^2} b_{3,1} \\right) \\right\\} = n \\]\n\n\\[ n (1 + k) = n \\left\\{ 1 - \\frac{m_i a^2}{4 \\mu a_i^2} b_{3,1} + \\frac{m_i a^2 e_i}{4 \\mu a_i^2 e} b_{3,2} \\cos (\\varphi - \\varphi_i) \\right\\} \\]\n\nLet \\( \\frac{\\mu}{a^3} = n^2 \\), then \\( a^3 = a^3 \\left\\{ 1 + \\frac{4m_i}{\\mu} \\left( \\frac{a^3}{a_i^3} b_{3,o} - \\frac{a^2}{a_i^2} b_{3,1} \\right) \\right\\} \\)\n\n\\[ e = a \\left\\{ 1 + \\frac{4m_i}{3 \\mu} \\left( \\frac{a^3}{a_i^3} b_{3,o} - \\frac{a^2}{a_i^2} b_{3,1} \\right) \\right\\} \\]\n\n\\[ \\lambda = n t + \\varepsilon + 2e \\sin \\left( n \\left\\{ 1 - \\frac{m_i a^2}{4 \\mu a_i^2} b_{3,1} + \\frac{m_i a^2 e_i}{4 \\mu a_i^2 e} b_{3,2} \\cos (\\varphi - \\varphi_i) \\right\\} t + \\varepsilon - \\varphi_2 \\right) + \\&c. \\]\n\n\\[ \\frac{a}{r} = 1 - \\frac{m_i}{3 \\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,o} - \\frac{a^2}{2 a_i^2} b_{3,1} \\right\\} \\]\n\n\\[ + e \\left\\{ 1 - \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,o} + \\frac{a^2}{12 a_i^2} b_{3,1} - \\frac{a^2}{4 a_i^2} b_{3,2} \\frac{e_i}{e} \\cos (\\varphi - \\varphi_i) \\right\\} \\right\\} \\]\n\n\\[ \\cos \\left( n \\left\\{ 1 - \\frac{m_i a^2}{4 \\mu a_i^2} b_{3,1} + \\frac{m_i a^2 e_i}{4 \\mu a_i^2 e} b_{3,2} \\cos (\\varphi - \\varphi_i) \\right\\} t + \\varepsilon - \\varphi_2 \\right) \\]\n\n\\[ + \\frac{m_i}{\\mu} \\left\\{ \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,o} - \\frac{a}{4 a_i} b_{3,1} - \\frac{a^2}{8 a_i^2} b_{3,2} \\right\\} e_i \\cos (nt + \\varepsilon - \\varphi_i) \\]\n\n\\[ \\frac{r}{a} = 1 + \\frac{m_i}{3 \\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,o} - \\frac{a^2}{2 a_i^2} b_{3,1} \\right\\} \\]\n\n\\[ - e \\left\\{ 1 + \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,o} - \\frac{5}{12} \\frac{a^2}{a_i^2} b_{3,1} + \\frac{a^2}{a_i^2} b_{3,2} \\frac{e_i}{e} \\cos (\\varphi - \\varphi_i) \\right\\} \\right\\} \\]\n\n\\[ \\cos \\left( n \\left\\{ 1 - \\frac{m_i a^2}{4 \\mu a_i^2} b_{3,1} + \\frac{m_i a^2 e_i}{4 \\mu a_i^2 e} b_{3,2} \\cos (\\varphi - \\varphi_i) \\right\\} t + \\varepsilon - \\varphi_2 \\right) \\]\n\n\\[ - \\frac{m_i}{\\mu} \\left\\{ \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,o} - \\frac{a}{4 a_i} b_{3,1} - \\frac{a^2}{8 a_i^2} b_{3,2} \\right\\} e_i \\cos (nt + \\varepsilon - \\varphi_i) + \\&c. \\]\n\nIf\n\n\\[ \\frac{a}{r} = 1 + r_0 + e(1+f) \\cos \\left( n (1+k^*) t + \\varepsilon - \\varphi \\right) + e_i f_i \\cos \\left( n(1+k_i) t + \\varepsilon - \\varphi_i \\right) + \\&c. \\]\n\n\\[ (1+f) \\{(1+k)^3(1-3r_0)-1\\} + \\frac{m_i a q_7}{\\mu} = 0 \\]\n\n* This quantity \\( k \\) must not be confounded with the quantity \\( k \\) above.\nwhence neglecting $k^2$, $kf$, &c.\n\n$$k = \\frac{3}{2} r_0 - \\frac{m_i}{2\\mu} a q_7$$\n\n$$= \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,0} - \\frac{a^2}{2a_i^2} b_{3,1} \\right\\}$$\n\nsimilarly\n\n$$f_i \\left\\{ (1 + k_i)^2 - 1 \\right\\} + \\frac{m_i}{\\mu} a q_{16} = 0$$\n\n$$k_i^2 + 2k_i + \\frac{m_i}{\\mu} \\frac{a}{f_i} q_{16} = 0, \\quad q_{16} = -\\frac{a}{2a_i^2} b_{3,2}$$\n\n$$k_i = -2 - \\frac{m_i}{4\\mu} \\frac{a^2}{f_i a_i^2} b_{3,2} \\quad \\text{or} \\quad + \\frac{m_i}{4\\mu} \\frac{a^2}{f_i a_i^2} b_{3,2}$$\n\n$$\\lambda = n \\left\\{ 1 + 2r_0 \\right\\} t + \\varepsilon + \\frac{2e(1+f)(1+r_0)}{(1+k)} \\sin \\left( n(1+k)t + \\varepsilon - \\varpi \\right)$$\n\n$$+ 2e_i \\left\\{ \\frac{f_i}{1+k_i} - \\frac{m_i}{2\\mu} a R_{16} \\right\\} \\sin \\left( n(1+k_i)t + \\varepsilon - \\varpi_i \\right) + \\text{&c. nearly}$$\n\nIf\n\n$$\\frac{e(1+f)(1+r_0)}{1+k} = e$$\n\n$$(1+f) = e(1+k-r_0)$$\n\nIf $\\frac{f_i}{1+k_i} - \\frac{m_i}{2\\mu} a R_{16} = 0$ and if $k_i^2$ be neglected\n\n$$f_i = \\frac{m_i}{4\\mu} \\frac{a^2}{a_i^2} b_{3,2} + \\frac{m_i}{2\\mu} a R_{16}$$\n\nIf $n(1+2r_0) = n$, and $n^2 = \\frac{\\mu}{a^3}$\n\n$$a = a \\left\\{ 1 + \\frac{4}{3} r_0 \\right\\}$$\n\n$$\\frac{a}{r} = 1 - \\frac{1}{3} r_0 + e \\left\\{ 1 + k - \\frac{7}{3} r_0 \\right\\} \\cos \\left( n(1+k)t + \\varepsilon - \\varpi \\right) + e_i f_i \\cos \\left( n(1+k_i)t + \\varepsilon - \\varpi_i \\right)$$\n\n$$= 1 - \\frac{1}{3} r_0 + e \\left\\{ 1 - \\frac{5}{6} r_0 - \\frac{m_i}{2\\mu} a q_7 \\right\\} \\cos \\left( n(1+\\frac{1}{2} r_0 - \\frac{m_i}{2\\mu} a q_7)t + \\varepsilon - \\varpi \\right)$$\n\n$$- e_i \\left\\{ \\frac{m_i}{2\\mu} a q_{16} - \\frac{m_i}{2\\mu} a R_{16} \\right\\} \\cos \\left( n(1+k_i)t + \\varepsilon - \\varpi_i \\right)$$\n\n$$\\frac{r}{a} = 1 + \\frac{1}{3} r_0 - e \\left\\{ 1 - \\frac{1}{6} r_0 - \\frac{m_i}{2\\mu} a q_7 \\right\\} \\cos \\left( n(1-\\frac{1}{2} r_0 - \\frac{m_i}{2\\mu} a q_7)t + \\varepsilon - \\varpi \\right)$$\n\n$$+ e \\left\\{ \\frac{m_i}{2\\mu} a q_{16} - \\frac{m_i}{2\\mu} a R_{16} \\right\\} \\cos \\left( n(1+k_i)t + \\varepsilon - \\varpi_i \\right)$$\n\\[ r_a = 1 + \\frac{1}{3} r_o - e \\left\\{ 1 - \\frac{1}{6} r_o - \\frac{m_i}{2 \\mu} a q_7 \\right\\} \\cos (nt + \\varepsilon - \\varpi) + e_i \\left\\{ \\frac{m_i}{2 \\mu} a q_{16} - \\frac{m_i}{2 \\mu} a R_{16} \\right\\} \\cos (nt + \\varepsilon - \\varpi_i) \\]\n\n\\[ - e \\left\\{ \\frac{1}{2} r_o + \\frac{m_i}{2 \\mu} a q_7 \\right\\} n t \\sin (nt + \\varepsilon - \\varpi) - e_i \\frac{m_i}{2 \\mu} a q_{16} n t \\sin (nt + \\varepsilon - \\varpi_i) \\]\n\nIn the notation of the Mécanique Céleste\n\n\\[ r_o = \\frac{a^2}{2} \\left( \\frac{d A^{(0)}}{da} \\right) \\]\n\\[ a q_7 = - \\frac{a^3}{2} \\left( \\frac{d^2 A^{(0)}}{da^2} \\right) - \\frac{3 a^2}{2} \\left( \\frac{d A^{(0)}}{da} \\right) \\]\n\n\\[ a q_{16} = - \\frac{1}{2} \\left\\{ a^2 a' \\left( \\frac{d^2 A^{(1)}}{da da'} \\right) + 2 a^3 \\left( \\frac{d A^{(1)}}{da} \\right) + 2 a' a \\left( \\frac{d A^{(1)}}{da'} \\right) + 4 a A^{(1)} \\right\\} \\]\n\n\\[ = - \\frac{1}{2} \\left\\{ 2 a A - 2 a^3 \\left( \\frac{d A^{(1)}}{da} \\right) - a^3 \\left( \\frac{d^2 A^{(1)}}{da^2} \\right) \\right\\} \\]\n\n\\[ R_{16} = - \\frac{1}{2} a' \\left( \\frac{d A^{(1)}}{da'} \\right) - A^{(1)} \\]\n\n\\[ = \\frac{a}{2} \\left( \\frac{d A^{(1)}}{da} \\right) - \\frac{1}{2} A^{(1)} \\]\n\n\\[ \\frac{\\mu}{6 m_i} r_o + \\frac{a}{2} q_7 = \\frac{a^2}{12} \\left( \\frac{d A^{(0)}}{da} \\right) - \\frac{a^3}{4} \\left( \\frac{d^2 A^{(0)}}{da^2} \\right) - \\frac{3 a^2}{4} \\left( \\frac{d A^{(0)}}{da} \\right) \\]\n\n\\[ = - \\frac{2 a^2}{3} \\left( \\frac{d A^{(0)}}{da} \\right) - \\frac{a^3}{4} \\left( \\frac{d^2 A^{(0)}}{da^2} \\right) = - f \\]\n\n\\[ \\frac{a}{2} q_{16} - \\frac{a}{2} R_{16} = - \\frac{a A^{(1)}}{2} + \\frac{a^2}{2} \\left( \\frac{d A^{(1)}}{da} \\right) + \\frac{a^3}{4} \\left( \\frac{d^2 A^{(1)}}{da^2} \\right) + \\frac{a A^{(1)}}{4} - \\frac{a^2}{4} \\left( \\frac{d A^{(1)}}{da} \\right) \\]\n\n\\[ = - \\frac{a A^{(1)}}{4} + \\frac{a^2}{4} \\left( \\frac{d A^{(1)}}{da} \\right) + \\frac{a^3}{4} \\left( \\frac{d^2 A^{(1)}}{da^2} \\right) = - f' \\]\n\n\\[ \\frac{\\mu}{2 m_i} r_o + a q_7 = \\frac{a^3}{4} \\left( \\frac{d A^{(0)}}{da} \\right) - \\frac{a^3}{4} \\left( \\frac{d^2 A^{(0)}}{da^2} \\right) - \\frac{3 a^2}{4} \\left( \\frac{d A^{(0)}}{da} \\right) \\]\n\n\\[ = - \\frac{a^2}{2} \\left( \\frac{d A^{(0)}}{da} \\right) - \\frac{a^3}{4} \\left( \\frac{d^2 A^{(0)}}{da^2} \\right) = - \\frac{C}{2} = \\frac{a^2}{4 a_i^3} b_{3,1} \\]\n\n\\[ - \\frac{1}{2} a q_{16} = \\frac{1}{2} \\left\\{ a A - a^2 \\left( \\frac{d A^{(1)}}{da} \\right) - \\frac{a^3}{2} \\left( \\frac{d^2 A^{(1)}}{da^2} \\right) \\right\\} = - \\frac{D}{2} = \\frac{a^2}{2 a_i^3} b_{3,2} \\]\n\nThese results evidently agree with those given in the Mécanique Céleste, vol. i. p. 279, with the exception of the sign in the value of \\( f' \\) marked with an asterisk, which I think requires alteration in that work.\n\nFinally, neglecting the quantities multiplied by \\( t \\), which may be made to de-\npend upon the secular inequalities of the constants $\\varepsilon$, $\\varpi$, &c., and the squares of the eccentricities.\n\n$$\\frac{r}{a} = 1 + \\frac{m_i}{3\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,0} - \\frac{a^2}{2a_i^2} b_{3,1} \\right\\} - e \\left\\{ 1 + \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,0} - \\frac{5}{12} \\frac{a^2}{a_i^2} b_{3,1} \\right\\} \\right\\} \\cos (nt + \\varepsilon - \\varpi)$$\n\n$$- \\left\\{ \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,0} - \\frac{a}{4a_i} b_{3,1} + \\frac{a^2}{8a_i^2} b_{3,2} \\right\\} e_i \\cos (nt + \\varepsilon - \\varpi_i)$$\n\n$$- r_1 \\cos (nt - n_t t + \\varepsilon - \\varepsilon_i) - r_2 \\cos (2nt - 2n_t t + 2\\varepsilon - 2\\varepsilon_i) - r_3 \\cos (3nt - 3n_t t + 3\\varepsilon - 3\\varepsilon_i) - \\&c.$$ \n\n$$- \\{ r_6 - r_1 \\} e \\cos (nt + \\varepsilon_i - \\varpi) - \\{ r_8 - r_1 \\} e \\cos (2nt - n_t t + 2\\varepsilon - \\varepsilon_i - \\varpi)$$\n\n$$- \\{ r_9 - r_2 \\} e \\cos (3nt - 2n_t t + 3\\varepsilon - 2\\varepsilon_i - \\varpi)$$\n\n$$- \\{ r_{10} - r_3 \\} e \\cos (4nt - 3n_t t + 4\\varepsilon - 3\\varepsilon_i - \\varpi) - \\{ r_{11} - r_4 \\} e \\cos (5nt - 4n_t t + 5\\varepsilon - 4\\varepsilon_i - \\varpi)$$\n\n$$- \\{ r_{12} - r_5 \\} e \\cos (nt - 2n_t t + \\varepsilon - 2\\varepsilon_i + \\varpi) - \\{ r_{13} - r_5 \\} e \\cos (2nt - 3n_t t + 2\\varepsilon - 3\\varepsilon_i + \\varpi)$$\n\n$$- \\{ r_{14} - r_6 \\} e \\cos (3nt - 4n_t t + 3\\varepsilon - 4\\varepsilon_i + \\varpi) - r_{16} e_i \\cos (nt + \\varepsilon_i - \\varpi_i) - r_{17} e_i \\cos (nt + \\varepsilon - \\varpi_i) - \\&c.$$ \n\nThe constant part of $R$\n\n$$= m_i \\left\\{ - \\frac{b_{1,0}}{a_i} + \\frac{a}{2a_i^2} \\left( \\sin^2 \\frac{\\varepsilon}{2} - \\frac{e^2}{4} + \\frac{e_i^2}{4} \\right) b_{3,1} + \\frac{a^2}{4a_i^2} b_{3,2} \\cos (\\varpi - \\varpi_i) \\right\\} \\text{ See p. 29.}$$\n\nIf this quantity $= - F$ according to the notation of the Théor. Anal. vol. i. p. 336.\n\n$$d\\varepsilon = (1 - \\sqrt{1 - e^2}) d\\varpi - \\frac{2a^2 n}{\\mu} \\left( \\frac{dF}{da} \\right) dt$$\n\n$$d\\varpi = an \\frac{\\sqrt{1 - e^2}}{\\mu e} \\left( \\frac{dF}{de} \\right) dt$$\n\n$$\\frac{d\\varepsilon - d\\varpi}{dt} = \\frac{m_i}{\\mu} \\left\\{ \\frac{2a^2}{a_i^3} b_{3,0} - \\frac{5}{4} \\frac{a^2}{a_i^2} b_{3,1} + \\frac{a^2}{4a_i^2} b_{3,2} \\frac{e_i}{e} \\cos (\\varpi - \\varpi_i) \\right\\}$$\n\nLet, as hitherto, (Phil. Trans. for 1830. p. 336.)\n\n$$r = \\frac{h^2 \\sqrt{1 + s^2}}{\\mu \\cos i^2 \\{ \\sqrt{1 + s^2} + e \\cos (\\lambda' - \\varpi) \\}} = a \\{ 1 - e' \\cos (v - a) \\}$$\n\nFig. 1.\n\nFig. 2.\nLet \\( P \\) be the place of the planet, \\( P' \\) its projection on the fixed plane \\( LN \\) (fig. 1 & 2.) \\( SN \\) the line of nodes, \\( SL \\) the line from which longitudes are reckoned. The angle \\( LSP' = \\lambda \\). Let \\( SA \\) be the line of apsides. (fig. 2.)\n\nIn the notation of M. de Pontécoulant, vol. i. p. 316, the angle \\( ASN = g, \\frac{\\varepsilon - \\varpi}{n} = l \\), M. de Pontécoulant has given expressions for the variations of the constants \\( a, g, e', l, \\iota \\) and \\( \\nu \\) in terms of the partial differences of the quantity \\( R \\) with regard to these quantities. It is easy from these to find similar expressions for the variations or differentials with regard to the time of the constants \\( a, \\varpi, e, \\varepsilon, \\iota \\) and \\( \\nu \\).\n\nLet \\( SAB \\) be a plane cutting the plane of the orbit at right angles, so that the angle \\( SAB = 90^\\circ, ANB = \\iota, BSN = \\varpi - \\nu \\)\n\n\\[\n\\frac{dr^2}{dt^2} + \\frac{h^2}{r^3 \\cos \\iota^2} - \\frac{2\\mu}{r} + \\frac{\\mu}{a} = 0\n\\]\n\n\\[\nr = a \\left( 1 - e' \\cos (\\nu - a) \\right)\n\\]\n\nWhen \\( r \\) is a maximum or minimum \\( \\frac{dr}{dt} = 0 \\),\n\n\\[\n\\frac{a h^2}{\\mu \\cos \\iota^2} - 2ar + r^2 = 0, \\quad \\text{whence} \\quad r = a \\pm \\sqrt{a - \\frac{h^2}{\\mu \\cos \\iota^2}}\n\\]\n\n\\[\nr = a \\left( 1 \\pm e' \\right)\n\\]\n\n\\[\n\\frac{h^2}{\\mu \\cos \\iota^2} = a \\left( 1 - e'^2 \\right)\n\\]\n\nBy the equation of p. 336, line 12, (Phil. Trans. 1830.)\n\n\\[\n\\frac{h^2}{\\mu \\cos \\iota^2} = a \\left( 1 - e^2 + e^2 \\sin^2 \\iota \\sin^2 (\\nu - \\varpi) \\right)\n\\]\n\n\\[\ne'^2 = e^2 \\left( 1 - \\sin^2 \\iota \\sin^2 (\\nu - \\varpi) \\right) * = e^2 \\cos^2 ASB\n\\]\n\nConsidering \\( R \\) first as a function of the quantities \\( a, g, e', l, \\iota \\) and \\( \\nu \\), and then of the quantities \\( a, \\varpi, e, \\varepsilon, \\iota \\) and \\( \\nu \\), we have\n\n\\[\n\\left( \\frac{dR}{da} \\right) da + \\left( \\frac{dR}{dg} \\right) dg + \\left( \\frac{dR}{de'} \\right) de' + \\left( \\frac{dR}{dl} \\right) dl + \\left( \\frac{dR}{d\\iota} \\right) d\\iota + \\left( \\frac{dR}{d\\nu} \\right) d\\nu\n\\]\n\n\\[\n= \\left( \\frac{dR}{da} \\right) da + \\left( \\frac{dR}{d\\varpi} \\right) d\\varpi + \\left( \\frac{dR}{de} \\right) de + \\left( \\frac{dR}{d\\varepsilon} \\right) d\\varepsilon + \\left( \\frac{dR}{dl} \\right) dl + \\left( \\frac{dR}{d\\nu} \\right) d\\nu\n\\]\n\n* The equation I gave, Phil. Trans. for 1830, p. 336, line 17, is not correct.\nBy means of this equation, the equations\n\n\\[ \\varepsilon - \\varpi = n l, \\quad e^2 = e^2 \\{ 1 - \\sin^2 \\iota \\sin^2 (\\nu - \\varpi) \\} \\]\n\nand the equations given by M. DE PONTECOULANT, vol. i. p. 328, the values of \\( \\frac{da}{dt}, \\frac{d\\varpi}{dt}, \\frac{de}{dt}, \\frac{d\\varepsilon}{dt}, \\frac{d\\nu}{dt} \\) may be easily obtained in terms of the quantities \\( a, \\varpi, e, \\varepsilon, \\iota, \\nu \\), and the partial differential coefficients \\( \\left( \\frac{dR}{da} \\right), \\left( \\frac{dR}{d\\varpi} \\right), \\left( \\frac{dR}{de} \\right), \\left( \\frac{dR}{d\\varepsilon} \\right), \\left( \\frac{dR}{d\\iota} \\right) \\) and \\( \\left( \\frac{dR}{d\\nu} \\right) \\).\n\nSubstituting in the equations of p. 40, for \\( q_n \\) their values and neglecting \\( e^2, e_1^2, e e_p \\) and \\( \\sin^2 \\frac{\\iota}{2} \\), when \\( a \\) is less than \\( a_i \\);\n\n\\[\n\\frac{a}{r} = 1 - \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,0} - \\frac{a^2}{a_i^2} b_{3,1} \\right\\} + e \\left\\{ 1 - \\frac{m_i}{\\mu} \\left( \\frac{a^3}{a_i^3} b_{3,0} + \\frac{a^2}{12 a_i^2} b_{3,1} \\right) \\right\\} \\cos (nt + \\varepsilon - \\varpi)\n\\]\n\n\\[\n+ \\frac{m_i}{\\mu} \\left( \\frac{n^2}{2n - n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\left( \\frac{a^2}{a_i^2} - \\frac{a}{a_i} b_{1,1} \\right) + \\frac{a^2}{a_i^2} \\right\\}\n\\]\n\n\\[\n+ \\frac{a^3}{a_i^3} \\left( \\frac{a}{a_i} b_{3,1} - b_{3,0} - \\frac{1}{2} b_{3,2} \\right) \\right\\} \\cos (nt - n_i t + \\varepsilon - \\varepsilon_i)\n\\]\n\n\\[\n+ \\frac{m_i}{\\mu} \\left( \\frac{n^2}{3n - 2n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\frac{a}{a_i} b_{1,2}\n\\]\n\n\\[\n- \\frac{a^2}{a_i^2} \\left( \\frac{a}{a_i} b_{3,2} - \\frac{1}{2} b_{3,1} - \\frac{1}{2} b_{3,3} \\right) \\right\\} \\cos (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i)\n\\]\n\n\\[\n+ \\frac{m_i}{\\mu} \\left( \\frac{n^2}{4n - 3n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\frac{a}{a_i} b_{1,3}\n\\]\n\n\\[\n- \\frac{a^2}{a_i^2} \\left( \\frac{a}{a_i} b_{3,3} - \\frac{1}{2} b_{3,2} - b_{3,4} \\right) \\right\\} \\cos (3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i)\n\\]\n\n\\[\n+ \\frac{m_i}{\\mu} \\left( \\frac{n^2}{5n - 4n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\frac{a}{a_i} b_{1,4}\n\\]\n\n\\[\n- \\frac{a^2}{a_i^2} \\left( \\frac{a}{a_i} b_{3,4} - \\frac{1}{2} b_{3,3} - \\frac{1}{2} b_{3,5} \\right) \\right\\} \\cos (4nt - 4n_i t + 4\\varepsilon + 4\\varepsilon_i)\n\\]\n\n\\[\n+ \\frac{m_i}{\\mu} \\left( \\frac{n^2}{6n - 5n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\frac{a}{a_i} b_{1,5}\n\\]\n\n\\[\n- \\frac{a^2}{a_i^2} \\left( \\frac{a}{a_i} b_{3,5} - \\frac{1}{2} b_{3,4} - \\frac{1}{2} b_{3,6} \\right) \\right\\} \\cos (5nt - 5n_i t + 5\\varepsilon - 5\\varepsilon_i)\n\\]\n\n* In the terms multiplied by \\( \\frac{m_i}{\\mu} \\) the quantities \\( a & a_i, e & e \\) may be used indifferently.\n\nMDCCXXXI.\nExpression for the reciprocal of the radius vector, when \\(a < a_r\\):\n\n\\[\n-\\frac{n^2}{(n-n_i)(n+n_i)} \\left\\{ \\frac{3a_i^2}{2} r_1 + \\frac{m_i}{\\mu} \\left[ \\frac{3a_i^2}{2} b_{3,0} + \\frac{a_i^3}{2} b_{3,1} - \\frac{a_i^2}{4} b_{3,2} \\right] \\right\\} e \\cos(n_i t + \\varepsilon_i - \\varphi) [6]\n\\]\n\n\\[\n+\\frac{n^2}{(3n-n_i)(n-n_i)} \\left\\{ \\frac{3(2n-n_i)^2}{2n^2} r_1 - \\frac{m_i}{\\mu} \\left[ \\frac{4n}{(2n-n_i)} \\left\\{ \\frac{a_i^2}{2} b_{3,0} - \\frac{a_i^2}{2} b_{3,1} - \\frac{a_i^3}{2} b_{3,2} + \\frac{3a_i^2}{4} b_{3,3} \\right\\} \\right] \\right\\} e \\cos(2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varphi) [8]\n\\]\n\n\\[\n+\\frac{n^2}{(4n-2n_i)(2n-2n_i)} \\left\\{ \\frac{3(3n-2n_i)^2}{2n^2} r_2 - \\frac{m_i}{\\mu} \\left[ \\frac{6n}{(3n-2n_i)} \\left\\{ -\\frac{a_i^2}{4} b_{3,1} - \\frac{a_i^3}{2} b_{3,2} + \\frac{3a_i^2}{4} b_{3,3} \\right\\} \\right] \\right\\} e \\cos(3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varphi) [9]\n\\]\n\n\\[\n+\\frac{n^2}{(5n-3n_i)(3n-3n_i)} \\left\\{ \\frac{3(4n-3n_i)^2}{2n^2} r_3 - \\frac{m_i}{\\mu} \\left[ \\frac{8n}{(4n-3n_i)} \\left\\{ -\\frac{a_i^2}{4} b_{3,2} - \\frac{a_i^3}{2} b_{3,3} + \\frac{3a_i^2}{4} b_{3,4} \\right\\} \\right] \\right\\} e \\cos(4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varphi) [10]\n\\]\n\n\\[\n+\\frac{n^2}{(6n-4n_i)(4n-4n_i)} \\left\\{ \\frac{3(5n-4n_i)^2}{2n^2} r_4 - \\frac{m_i}{\\mu} \\left[ \\frac{10n}{(5n-4n_i)} \\left\\{ -\\frac{a_i^2}{4} b_{3,3} - \\frac{a_i^3}{2} b_{3,4} + \\frac{3a_i^2}{4} b_{3,5} \\right\\} \\right] \\right\\} e \\cos(5nt - 4n_i t + 5\\varepsilon - 4\\varepsilon_i - \\varphi) [11]\n\\]\n\n\\[\n-\\frac{n^2}{2n_i(2n-2n_i)} \\left\\{ \\frac{3(n-2n_i)^2}{2n^2} r_2 - \\frac{m_i}{\\mu} \\left[ \\frac{2n}{(n-2n_i)} \\left\\{ \\frac{3a_i^2}{4} b_{3,1} - \\frac{a_i^3}{2} b_{3,2} - \\frac{a_i^2}{4} b_{3,3} \\right\\} \\right] \\right\\} e \\cos(nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varphi) [12]\n\\]\n\n\\[\n+\\frac{n^2}{(n-3n_i)(3n-3n_i)} \\left\\{ \\frac{3(2n-3n_i)^2}{2n^2} r_3 - \\frac{m_i}{\\mu} \\left[ \\frac{4n}{(2n-3n_i)} \\left\\{ \\frac{3a_i^2}{4} b_{3,2} - \\frac{a_i^3}{2} b_{3,3} - \\frac{a_i^2}{4} b_{3,4} \\right\\} \\right] \\right\\} e \\cos(2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varphi) [13]\n\\]\n\n\\[\n+\\frac{n^2}{(2n-4n_i)(4n-4n_i)} \\left\\{ \\frac{3(3n-4n_i)^2}{2n^2} r_4 - \\frac{m_i}{\\mu} \\left[ \\frac{6n}{(3n-4n_i)} \\left\\{ \\frac{3a_i^2}{4} b_{3,3} - \\frac{a_i^3}{2} b_{3,4} - \\frac{a_i^2}{4} b_{3,5} \\right\\} \\right] \\right\\} e \\cos(3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varphi) [14]\n\\]\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(n-n_i)(n+n_i)} \\frac{a^2}{2a_i^2} b_{3,1} e_i \\cos(n_i t + \\varepsilon_i - \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(n-n_i)(3n-n_i)} \\left[ \\frac{4n}{(2n-n_i)} \\left\\{ \\frac{3a^2}{4a_i^2} b_{3,1} - \\frac{a}{2a_i} b_{3,2} - \\frac{a^2}{4a_i^2} b_{3,3} \\right\\} \\right]\n\\]\n\n\\[\n-\\frac{9}{4} \\frac{a^2}{a_i^2} b_{3,1} + \\frac{2a}{a_i} b_{3,2} - \\frac{a^2}{4a_i^2} b_{3,3} \\right\\} e_i \\cos(2n_i t - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(2n-n_i)(4n-2n_i)} \\left[ \\frac{6n}{(3n-2n_i)} \\left\\{ \\frac{3a^2}{4a_i^2} b_{3,2} - \\frac{a}{2a_i} b_{3,3} - \\frac{a^2}{4a_i^2} b_{3,4} \\right\\} \\right]\n\\]\n\n\\[\n-\\frac{3a^2}{a_i^2} b_{3,2} + \\frac{3a}{a_i} b_{3,3} - \\frac{a^2}{2a_i^2} b_{3,4} \\right\\} e_i \\cos(3n_i t - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^3}{(3n-2n_i)(5n-3n_i)} \\left[ \\frac{8n}{(4n-3n_i)} \\left\\{ \\frac{3a^2}{4a_i^2} b_{3,3} - \\frac{a}{2a_i} b_{3,4} - \\frac{a^2}{4a_i^2} b_{3,5} \\right\\} \\right]\n\\]\n\n\\[\n-\\frac{15}{4} \\frac{a^2}{a_i^2} b_{3,3} + \\frac{4a}{a_i} b_{3,4} - \\frac{3a^2}{4a_i^2} b_{3,5} \\right\\} e_i \\cos(4n_i t - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi_i)\n\\]\n\n\\[\n+\\frac{m_i}{\\mu} \\frac{n^2}{2n_i(2n-2n_i)} \\left[ \\frac{2n}{(n-2n_i)} \\left\\{ \\frac{2a^2}{a_i^2} - \\frac{a^2}{2a_i^2} b_{3,0} - \\frac{a}{2a_i} b_{3,1} + \\frac{3a^2}{4a_i^2} b_{3,2} \\right\\} \\right]\n\\]\n\n\\[\n+\\frac{2a^2}{a_i^2} b_{3,0} - \\frac{a}{a_i} b_{3,1} \\right\\} e_i \\cos(n_i t - 2n_i t + \\varepsilon - \\varepsilon_i + \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(n-3n_i)(3n-3n_i)} \\left[ \\frac{4n}{(2n-3n_i)} \\left\\{ -\\frac{a^2}{4a_i^2} b_{3,1} - \\frac{a}{2a_i} b_{3,2} + \\frac{3a^2}{4a_i^2} b_{3,3} \\right\\} \\right]\n\\]\n\n\\[\n+\\frac{3a^2}{4a_i^2} b_{3,1} - \\frac{2a}{a_i} b_{3,2} + \\frac{3a^2}{4a_i^2} b_{3,3} \\right\\} e_i \\cos(2n_i t - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(2n-4n_i)(4n-4n_i)} \\left[ \\frac{6n}{(3n-4n_i)} \\left\\{ -\\frac{a^2}{4a_i^2} b_{3,2} - \\frac{a}{2a_i} b_{3,3} + \\frac{3a^2}{4a_i^2} b_{3,4} \\right\\} \\right]\n\\]\n\n\\[\n+\\frac{a^2}{a_i^2} b_{3,2} - \\frac{3a}{a_i} b_{3,3} + \\frac{3a^2}{2a_i^2} b_{3,4} \\right\\} e_i \\cos(3n_i t - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi_i)\n\\]\n\n\\[\n-\\frac{m_i}{\\mu} \\frac{n^2}{(3n-5n_i)(5n-5n_i)} \\left[ \\frac{8n}{(4n-5n_i)} \\left\\{ -\\frac{a^2}{4a_i^2} b_{3,3} - \\frac{a}{2a_i} b_{3,4} + \\frac{3a^2}{4a_i^2} b_{3,5} \\right\\} \\right]\n\\]\n\n\\[\n+\\frac{5}{4} \\frac{a^2}{a_i^2} b_{3,3} - \\frac{4a}{a_i} b_{3,4} + \\frac{9}{4} \\frac{a^2}{a_i^2} b_{3,5} \\right\\} e_i \\cos(4n_i t - 5n_i t + 4\\varepsilon - 5\\varepsilon_i + \\varpi_i)\n\\]\n\nSubstituting in the equations of p. 45, for $R_n$ their values and neglecting $e^2$, $e e_i$, $e_i^2$, and $\\sin^2 \\frac{i}{2}$,\nExpression for the longitude when \\(a < a_i\\).\n\n\\[\n\\lambda = n t + \\varepsilon + e \\sin (nt + \\varepsilon - \\varphi)\n\\]\n\n\\[\n+ \\frac{n}{(n-n_i)} \\left\\{ 2 r_1 - \\frac{m_i n}{\\mu (n-n_i)} \\left( \\frac{a^2}{a_i^2} - \\frac{a}{a_i} b_{1,1} \\right) \\right\\} \\sin (nt - n_i t + \\varepsilon - \\varepsilon_i) \\quad [1]\n\\]\n\n\\[\n+ \\frac{n}{(2n-2n_i)} \\left\\{ 2 r_2 + \\frac{m_i n}{\\mu (n-n_i)} \\frac{a}{a_i} b_{1,2} \\right\\} \\sin (2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i) \\quad [2]\n\\]\n\n\\[\n+ \\frac{n}{(3n-3n_i)} \\left\\{ 2 r_3 + \\frac{m_i n}{\\mu (n-n_i)} a_i b_{1,3} \\right\\} \\sin (3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i) \\quad [3]\n\\]\n\n\\[\n+ \\frac{n}{(4n-4n_i)} \\left\\{ 2 r_4 + \\frac{m_i n}{\\mu (n-n_i)} a_i b_{1,4} \\right\\} \\sin (4nt - 4n_i t + 4\\varepsilon - 4\\varepsilon_i) \\quad [4]\n\\]\n\n\\[\n+ \\frac{n}{(5n-5n_i)} \\left\\{ 2 r_5 + \\frac{m_i n}{\\mu (n-n_i)} a_i b_{1,5} \\right\\} \\sin (5nt - 5n_i t + 5\\varepsilon - 5\\varepsilon_i) \\quad [5]\n\\]\n\n\\[\n+ \\frac{n}{n_i} \\left\\{ 2 \\left( r_6 + \\frac{r_1}{2} \\right) + \\frac{m_i n}{\\mu n_i} \\left( -\\frac{3a^2}{2a_i^2} + \\frac{3a^2}{2a_i^2} b_{3,0} - \\frac{a^3}{2a_i^3} b_{3,1} - \\frac{a^2}{4a_i^2} b_{3,2} \\right) \\right\\} e \\sin (nt + \\varepsilon - \\varphi) \\quad [6]\n\\]\n\n\\[\n+ \\frac{n}{(2n-n_i)} \\left\\{ 2 \\left( r_8 + \\frac{r_1}{2} \\right) - \\frac{m_i n}{\\mu (2n-n_i)} \\left( \\frac{a^2}{2a_i^2} - \\frac{a^2}{2a_i^2} b_{3,0} - \\frac{a^3}{2a_i^3} b_{3,1} + \\frac{3a^2}{4a_i^2} b_{3,2} \\right) \\right\\} e \\sin (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varphi) \\quad [8]\n\\]\n\n\\[\n+ \\frac{n}{(3n-2n_i)} \\left\\{ 2 \\left( r_9 + \\frac{r_2}{2} \\right) - \\frac{2m_i n}{\\mu (3n-2n_i)} \\left( -\\frac{a^3}{4a_i^3} b_{3,1} - \\frac{a^3}{2a_i^3} b_{3,2} + \\frac{3a^2}{4a_i^2} b_{3,3} \\right) \\right\\} e \\sin (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varphi) \\quad [9]\n\\]\n\n\\[\n+ \\frac{n}{(4n-3n_i)} \\left\\{ 2 \\left( r_{10} + \\frac{r_3}{2} \\right) - \\frac{3m_i n}{\\mu (4n-3n_i)} \\left( -\\frac{a^2}{4a_i^2} b_{3,2} - \\frac{a^3}{2a_i^3} b_{3,3} + \\frac{3a^2}{4a_i^2} b_{3,4} \\right) \\right\\} e \\sin (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varphi) \\quad [10]\n\\]\n\n\\[\n+ \\frac{n}{(5n-4n_i)} \\left\\{ 2 \\left( r_{11} + \\frac{r_4}{2} \\right) - \\frac{4m_i n}{\\mu (5n-4n_i)} \\left( -\\frac{a^2}{4a_i^2} b_{3,3} - \\frac{a^3}{2a_i^3} b_{3,4} + \\frac{3a^2}{4a_i^2} b_{3,5} \\right) \\right\\} e \\sin (5nt - 4n_i t + 5\\varepsilon - 4\\varepsilon_i - \\varphi) \\quad [11]\n\\]\n\n\\[\n+ \\frac{n}{(n-2n_i)} \\left\\{ 2 \\left( r_{12} + \\frac{r_5}{2} \\right) - \\frac{2m_i n}{\\mu (n-2n_i)} \\left( \\frac{3a^2}{4a_i^2} b_{3,1} - \\frac{a^3}{2a_i^3} b_{3,2} - \\frac{a^2}{4a_i^2} b_{3,3} \\right) \\right\\} e \\sin (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varphi) \\quad [12]\n\\]\n\\[\n\\begin{align*}\n&+ \\frac{n}{(3n - 3n_i)} \\left\\{ 2 \\left( r_{13} + \\frac{r_3}{2} \\right) - \\frac{3m_n}{\\mu (2n - 3n_i)} \\left( \\frac{3a^2}{4a_i^3} b_{3,2} - \\frac{a^3}{2a_i^3} b_{3,3} - \\frac{a^2}{4a_i^2} b_{3,4} \\right) \\right. \\\\\n&\\quad + \\frac{m_n a}{\\mu (n - n_i)} a_i b_{1,3} \\bigg\\} e \\sin (2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi) \\\\\n&+ \\frac{n}{(3n - 4n_i)} \\left\\{ 2 \\left( r_{14} + \\frac{r_4}{2} \\right) - \\frac{4m_n}{\\mu (3n - 4n_i)} \\left( \\frac{3a^2}{4a_i^2} b_{3,3} - \\frac{a^3}{2a_i^3} b_{3,4} - \\frac{a^2}{4a_i^2} b_{3,5} \\right) \\right. \\\\\n&\\quad + \\frac{m_n a}{\\mu (n - n_i)} a_i b_{1,4} \\bigg\\} e \\sin (3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi) \\\\\n&+ \\frac{2n}{n_i} r_{15} e_i \\sin (n_i t + \\varepsilon_i - \\varpi_i) \\\\\n&+ \\frac{n}{(2n - n_i)} \\left\\{ 2r_{17} - \\frac{2m_n}{\\mu (2n - n_i)} \\left( \\frac{3a^2}{4a_i^2} b_{3,1} - \\frac{a}{2a_i} b_{3,2} \\right. \\right. \\\\\n&\\quad \\left. \\left. - \\frac{a^2}{4a_i^2} b_{3,3} \\right) \\right\\} e_i \\sin (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi_i) \\\\\n&+ \\frac{n}{(3n - 2n_i)} \\left\\{ 2r_{18} - \\frac{3m_n}{\\mu (3n - 2n_i)} \\left( \\frac{3a^2}{4a_i^2} b_{3,2} - \\frac{a}{2a_i} b_{3,3} \\right. \\right. \\\\\n&\\quad \\left. \\left. - \\frac{a^3}{4a_i^3} b_{3,4} \\right) \\right\\} e_i \\sin (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi_i) \\\\\n&+ \\frac{n}{(4n - 3n_i)} \\left\\{ 2r_{19} - \\frac{4m_n}{\\mu (4n - 3n_i)} \\left( \\frac{3a^2}{4a_i^2} b_{3,3} - \\frac{a}{2a_i} b_{3,4} \\right. \\right. \\\\\n&\\quad \\left. \\left. - \\frac{a^3}{4a_i^3} b_{3,3} \\right) \\right\\} e_i \\sin (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi_i) \\\\\n&+ \\frac{n}{(n - 2n_i)} \\left\\{ 2r_{20} - \\frac{m_n}{\\mu (n - 2n_i)} \\left( \\frac{2a^2}{a_i^2} - \\frac{a^3}{2a_i^3} b_{3,0} - \\frac{a}{2a_i} b_{3,1} \\right. \\right. \\\\\n&\\quad \\left. \\left. + \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,2} \\right) \\right\\} e_i \\sin (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi_i) \\\\\n&+ \\frac{n}{(2n - 3n_i)} \\left\\{ 2r_{21} - \\frac{m_n}{\\mu (2n - 3n_i)} \\left( - \\frac{a^2}{4a_i^2} b_{3,1} - \\frac{a}{2a_i} b_{3,2} \\right. \\right. \\\\\n&\\quad \\left. \\left. + \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,3} \\right) \\right\\} e_i \\sin (2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi_i) \\\\\n&+ \\frac{n}{(3n - 4n_i)} \\left\\{ 2r_{22} - \\frac{m_n}{\\mu (3n - 4n_i)} \\left( - \\frac{a^2}{4a_i^2} b_{3,2} - \\frac{a}{2a_i} b_{3,3} \\right. \\right. \\\\\n&\\quad \\left. \\left. + \\frac{3}{4} \\frac{a^2}{a_i^2} b_{3,4} \\right) \\right\\} e_i \\sin (3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi_i)\n\\end{align*}\n\\]\n\nThe expression for the tangent of the latitude has already been given, p. 49.\nWhen the latitude is reckoned from the plane of the orbit of the planet P, the following terms must be added, in the general case where \\( \\iota \\) is not equal to zero, to that expression;\n\n\\[\n-\\frac{n^2}{2(n-n_i)(3n-n_i)} \\frac{a^2}{a_i^2} \\tan \\iota b_{3,1} \\sin (2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\nu)\n\\]\n\n\\[\n+\\frac{n^2}{2(n-n_i)(n+n_i)} \\frac{a^2}{a_i^2} \\tan \\iota b_{3,1} \\sin (n_i t + \\varepsilon_i - \\nu_i)\n\\]\n\n\\[\n-\\frac{n^2}{2(2n-2n_i)(4n-2n_i)} \\frac{a^3}{a_i^2} \\tan \\iota b_{3,2} \\sin (3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\nu)\n\\]\n\n\\[\n-\\frac{n^2}{4n_i(2n-2n_i)} \\frac{a^2}{a_i^2} \\tan \\iota b_{3,2} \\sin (nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\nu)\n\\]\n\n\\[\n-\\frac{n^2}{2(3n-3n_i)(5n-3n_i)} \\frac{a^3}{a_i^2} \\tan \\iota b_{3,3} \\sin (4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\nu)\n\\]\n\n\\[\n+\\frac{n^2}{2(n-3n_i)(3n-3n_i)} \\frac{a^2}{a_i^2} \\tan \\iota b_{3,3} \\sin (2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\nu)\n\\]\n\nAll the equations hitherto given apply to the case of an inferior disturbed by a superior planet, or when \\( a_i > a \\), in order to render them applicable to the case when \\( a_i < a \\) it is necessary to write \\( a \\) instead of \\( a_i \\) in the denominator of the terms multiplied by \\( b_{1,n} \\), and \\( a^3 \\) instead of \\( a_i^3 \\) in the denominator of the terms multiplied by \\( b_{3,n} \\) in the disturbing function \\( R \\), but the expressions for the quantities \\( q \\) are not the same in this case.\n\nIt will I think be admitted that the expressions which occur in the theory of the disturbances of the planets are more simple in terms of the quantities of which the general symbol is \\( b \\), than in terms of the partial differential coefficients of the quantity called \\( A \\) in the notation of the Mécanique Céleste. The development of the disturbing function \\( R \\) in terms of the differential coefficients \\( \\frac{dA}{da}, \\frac{dA}{da_i}, \\&c. \\) admits of reductions, so that it may be expressed in terms of the differential coefficients of \\( A \\) with respect to \\( a \\) only. In this state it has been left by Laplace as may be seen, vol. ii. p. 12, but the coefficients of the terms multiplied by the squares and products of the eccentricities may be expressed very simply in terms of the quantities of which the general symbol is \\( b \\), by means of reductions, of which two exam-\nples are given in the Théor. Anal. vol. i. p. 362. Similar reductions are applicable to the terms in \\( r \\left( \\frac{dR}{dr} \\right) \\) multiplied by the first power of the eccentricities.\n\nIn Laplace's notation\n\n\\[\nb_{1,0} = \\frac{1}{2} b_{\\frac{1}{2}}^{(0)}, \\quad b_{1,1} = b_{\\frac{1}{2}}^{(1)}, \\quad b_{1,2} = b_{\\frac{1}{2}}^{(2)}, \\quad b_{1,3} = b_{\\frac{1}{2}}^{(3)}\n\\]\n\n\\[\nb_{3,0} = \\frac{1}{2} b_{\\frac{3}{2}}^{(0)}, \\quad b_{3,1} = b_{\\frac{3}{2}}^{(1)}, \\quad b_{3,2} = b_{\\frac{3}{2}}^{(2)}, \\quad b_{3,3} = b_{\\frac{3}{2}}^{(3)}\n\\]\n\n\\[\nb_{5,0} = \\frac{1}{2} b_{\\frac{5}{2}}^{(0)}, \\quad b_{5,1} = b_{\\frac{5}{2}}^{(1)}, \\quad b_{5,2} = b_{\\frac{5}{2}}^{(2)}, \\quad b_{5,3} = b_{\\frac{5}{2}}^{(3)}\n\\]\n\n\\[\n\\frac{a}{a_i} b_{3,0} - \\frac{1}{2} b_{3,1} = -\\frac{1}{2} \\frac{db_{\\frac{3}{2}}^{(0)}}{da}, \\quad \\frac{a}{a_i} b_{3,1} - b_{3,0} - \\frac{1}{2} b_{3,2} = -\\frac{db_{\\frac{3}{2}}^{(1)}}{da}\n\\]\n\n\\[\n\\frac{a}{a_i} b_{3,2} - \\frac{1}{2} b_{3,1} - \\frac{1}{2} b_{3,3} = -\\frac{db_{\\frac{3}{2}}^{(2)}}{da}, \\quad \\frac{a}{a_i} b_{3,3} - \\frac{1}{2} b_{3,2} - \\frac{1}{2} b_{3,4} = -\\frac{db_{\\frac{3}{2}}^{(3)}}{da}\n\\]\n\nThe numerical values of these quantities are given for the principal planets in the third volume of the Mécanique Céleste.\n\nThe following numerical examples will serve to explain the expressions given above, and to show their accuracy, the results agreeing exactly with those given in the Mécanique Céleste.\n\n\\[\n\\frac{r}{a} = 1 + \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{a_i^3} b_{3,0} - \\frac{a^2}{2a_i^3} b_{3,1} \\right\\} - e \\left\\{ 1 + \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{3a_i^3} b_{3,0} - \\frac{5}{12} \\frac{a^2}{a_i^3} b_{3,1} \\right\\} \\right\\} \\cos (nt + \\varepsilon - \\varpi)\n\\]\n\n\\[\n- r_1 \\cos (nt - n't + \\varepsilon - \\varepsilon_i) \\quad \\text{See p. 55.}\n\\]\n\n* The coefficient of \\( \\cos (nt + \\varepsilon - \\varpi) \\) p. 52 line 11, and p. 55 line 3, should be\n\n\\[\n- e \\left\\{ 1 + \\frac{m_i}{\\mu} \\left\\{ \\frac{a^3}{3a_i^3} b_{3,0} - \\frac{5}{12} \\frac{a^2}{a_i^3} b_{3,1} \\right\\} \\right\\}.\n\\]\n\\[ r_i = \\frac{m_i}{\\mu} \\left( \\frac{n^2}{(2n - n_i)n_i} \\right) \\left\\{ \\frac{2n}{(n - n_i)} \\left( \\frac{a_i^2}{a_i^2} - \\frac{a_i b_{i,1}}{a_i^2} \\right) + \\frac{a_i^2}{a_i^2} \\left( \\frac{a_i b_{i,1}}{a_i^2} - b_{i,0} - \\frac{1}{2} b_{i,2} \\right) \\right\\} \\] See p. 57.\n\n\\[ \\lambda = n t + \\varepsilon + \\frac{n}{(n - n_i)} \\left\\{ 2 r_i - \\frac{m_i n}{\\mu (n - n_i)} \\left( \\frac{a_i^2}{a_i^2} - \\frac{a_i b_{i,1}}{a_i^2} \\right) \\right\\} \\sin (nt - n_i t + \\varepsilon - \\varepsilon_i) \\]\n\nIn the theory of Jupiter disturbed by Saturn,\n\n\\[ \\frac{a}{a_i} = 0.54531725, \\quad b_{i,1} = 0.6206406, \\quad b_{i,0} = \\frac{4.358387}{2} = 2.179193 \\]\n\n\\[ b_{i,1} = 3.185493 \\quad b_{i,2} = 2.082131 \\quad \\frac{m_i}{\\mu} = \\frac{1}{3359.4} \\]\n\n\\[ n = 337210.78 \\quad n_i = 135792.34 \\]\n\n\\[ a = 5.20116636 \\quad e = 0.0480767 \\]\n\nSee Méc. Cél. vol. iii. p. 61 & 82.\n\nWhence\n\n\\[ \\frac{a^3}{a_i^3} b_{i,0} = 0.353381, \\quad \\frac{a^3}{2 a_i^3} b_{i,1} = 0.473636 \\quad \\frac{a^3}{a_i^3} b_{i,0} - \\frac{a^3}{2 a_i^3} b_{i,1} = -0.120255 \\]\n\n\\[ \\log_9 0.120255 = 9.0801033 \\]\n\n\\[ \\log_5 5.20116636 = 0.7161007 \\]\n\n\\[ \\log_9 9.7962040 \\]\n\n\\[ \\log_3 \\frac{3 \\mu}{m_i} = 4.0033829 \\]\n\n\\[ 5.7928211 = \\log_9 0.000620613 \\text{ minus } \\]\n\nLaplace has . . . . . . . . . . 0000620566 minus\n\nSee Méc. Cél. vol. iii. p. 121, line 5.\n\n\\[ \\log_9 \\left\\{ \\frac{5}{12} \\frac{a^2}{a_i^3} b_{i,1} - \\frac{a^3}{3 a_i^3} b_{i,0} \\right\\} = 9.4423277 \\]\n\n\\[ \\log_9 e = 8.6819347 \\]\n\n\\[ \\log_9 a = 0.7161007 \\]\n\n\\[ \\log_9 8.8403631 \\]\n\n\\[ \\log_9 \\frac{\\mu}{m_i} = 3.5262617 \\]\n\n\\[ 5.3141014 = \\log_9 0.000206111 + \\]\n\nLaplace has . . . . . . . . . . 0000206111 the sign omitted,\n\nMéc. Cél. vol. iii. p. 122, line 28.\nCalculation of \\( r_1 \\), see p. 57.\n\n\\[\n\\log n_i = 5.1328751 \\\\\n\\log n = 5.5279013 \\\\\n\\frac{a^2}{a_i^2} - \\frac{a}{a_i} b_{i,1} = -0.041075 \\\\\n\\log 0.041075 = 8.6135776 \\\\\n\\frac{n - n_i}{n} = 0.7761940 \\\\\n\\frac{a}{a_i} b_{3,1} - b_{3,0} - \\frac{1}{2} b_{3,2} = -1.48315 \\\\\n\\log 1.48315 = 1.711851 \\\\\n\\log \\frac{a^2}{a_i^2} = 9.4732982 \\\\\n\\frac{a^2}{a_i^2} = 2.97371 \\\\\n\\frac{a}{a_i} = 0.281209 \\\\\n\\log 9.6444333 = \\log 0.441045 \\\\\n\\log 0.281209 = 9.4490293 \\\\\n\\log 3.3346239 = 3.3346259 \\\\\n\\log 6.1144054 = \\log 0.0001301383 \\\\\n\\log a = 0.7161007 \\\\\n\\log 6.8305061 = \\log 0.000676871 + \\\\\n\\text{LAPLACE has . . . . . 0.000676876 + p. 120.}\n\\]\n\nSimilarly \\( r_3 = 0.000556924 \\), \\( r_3 = 0.0005809 \\), \\( r_4 = 0.00015045 \\), \\( r_5 = 0.000049785 \\).\n\nCalculation of the coefficient of \\( \\sin (nt - n_i t + \\epsilon - \\epsilon) \\) in the value of \\( \\lambda \\). See p. 59.\n\n\\[\n2r_1 = -0.002602766 \\\\\n\\log 0.0002398064 = 6.3798601 \\\\\n\\log \\frac{n - n_i}{n} = 9.7761940 \\\\\n\\log 6.6036661 \\\\\n\\log 5.8038801 \\text{ see p. 49.} \\\\\n\\log 2.4075461 = \\log 255.591 \\text{ minus.}\n\\]\n\n\\text{LAPLACE has . . . . 255.5917 the sign omitted, p. 120.}\n\nIn the notation of the Méc. Cél. vol. iii. p. 120.\n\n\\[\nn = n^{iv}, \\quad n_i = n^v.\n\\]\n\nMDCCCCXXXI.\nThe reader is requested to make the following corrections.\n\nThe expression for \\( r^{-3} \\) Phil. Trans. 1830, p. 345, should be\n\n\\[\nr^{-3} = a^{-3} \\left\\{ 1 + \\frac{3}{2} e^2 \\left( 1 + \\frac{5}{4} e^2 \\right) + 3e \\left( 1 + \\frac{9}{8} e^2 \\right) \\cos (nt - \\varpi) \\right. \\\\\n+ \\frac{9}{2} e^2 \\left( 1 + \\frac{7}{9} e^2 \\right) \\cos (2nt - 2\\varpi) + \\frac{53}{8} e^3 \\cos (3nt - 3\\varpi) \\\\\n+ \\frac{77}{8} e^4 \\cos (4nt - 4\\varpi)\n\\]\n\nThe third term (multiplied by \\( dt \\)) in the equations of p. 23, line 4 and 6, and the second term in the equation of the same page, line 11, must be suppressed.\n\np. 30, line 21, read \\( + m_i \\left\\{ \\frac{a}{a_i^3} \\left( \\cos^3 \\frac{t_i}{2} - \\frac{e^2 + e_i^2}{2} \\right) \\right\\} \\) for \\( m_i \\left\\{ -\\frac{a}{a_i^3} \\left( \\cos^3 \\frac{t_i}{2} - \\frac{e^2 + e_i^2}{2} \\right) \\right\\} \\)\n\np. 32, line 5, read \\( + m_i \\left\\{ \\frac{2a}{a_i^3} - \\frac{a}{2a_i^3} b_{3,0} - \\frac{a_i^3}{2a_i^3} b_{3,1} \\right\\} \\&c. \\) for \\( + m_i \\left\\{ \\frac{a}{2a_i^3} - \\frac{a}{2a_i^3} b_{3,0} - \\frac{a_i^3}{2a_i^3} b_{3,1} \\right\\} \\&c. \\)\n\np. 39, line 9, read \\( - \\frac{21}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t - z + x) \\) for \\( - \\frac{27}{8} \\frac{a^2}{a_i^3} ee_i \\cos (2t - z + x) \\)\n\np. 41, line 4, read \\( \\frac{n_i^2}{n_i^2} \\) for \\( \\frac{n_i^2}{n_i^2} \\).",
  "source": "olmocr",
  "added": "2026-01-12",
  "created": "2026-01-12",
  "metadata": {
    "Source-File": "/home/jic823/projects/def-jic823/royalsociety/pdfs/107921.pdf",
    "olmocr-version": "0.3.4",
    "pdf-total-pages": 51,
    "total-input-tokens": 87052,
    "total-output-tokens": 61733,
    "total-fallback-pages": 0
  },
  "attributes": {
    "pdf_page_numbers": [
      [
        0,
        0,
        1
      ],
      [
        0,
        1955,
        2
      ],
      [
        1955,
        3844,
        3
      ],
      [
        3844,
        6185,
        4
      ],
      [
        6185,
        8472,
        5
      ],
      [
        8472,
        10568,
        6
      ],
      [
        10568,
        12552,
        7
      ],
      [
        12552,
        14910,
        8
      ],
      [
        14910,
        16880,
        9
      ],
      [
        16880,
        19314,
        10
      ],
      [
        19314,
        21809,
        11
      ],
      [
        21809,
        23725,
        12
      ],
      [
        23725,
        26000,
        13
      ],
      [
        26000,
        28923,
        14
      ],
      [
        28923,
        31663,
        15
      ],
      [
        31663,
        34562,
        16
      ],
      [
        34562,
        37720,
        17
      ],
      [
        37720,
        40558,
        18
      ],
      [
        40558,
        43877,
        19
      ],
      [
        43877,
        45276,
        20
      ],
      [
        45276,
        47157,
        21
      ],
      [
        47157,
        48935,
        22
      ],
      [
        48935,
        51440,
        23
      ],
      [
        51440,
        53870,
        24
      ],
      [
        53870,
        56017,
        25
      ],
      [
        56017,
        57685,
        26
      ],
      [
        57685,
        59792,
        27
      ],
      [
        59792,
        61685,
        28
      ],
      [
        61685,
        63734,
        29
      ],
      [
        63734,
        66573,
        30
      ],
      [
        66573,
        69468,
        31
      ],
      [
        69468,
        72753,
        32
      ],
      [
        72753,
        75981,
        33
      ],
      [
        75981,
        78675,
        34
      ],
      [
        78675,
        81220,
        35
      ],
      [
        81220,
        83300,
        36
      ],
      [
        83300,
        85753,
        37
      ],
      [
        85753,
        87721,
        38
      ],
      [
        87721,
        90480,
        39
      ],
      [
        90480,
        92904,
        40
      ],
      [
        92904,
        95251,
        41
      ],
      [
        95251,
        97923,
        42
      ],
      [
        97923,
        100345,
        43
      ],
      [
        100345,
        103119,
        44
      ],
      [
        103119,
        105801,
        45
      ],
      [
        105801,
        108433,
        46
      ],
      [
        108433,
        110807,
        47
      ],
      [
        110807,
        112845,
        48
      ],
      [
        112845,
        114541,
        49
      ],
      [
        114541,
        115781,
        50
      ],
      [
        115781,
        117092,
        51
      ]
    ],
    "primary_language": [
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "en",
      "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,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      true,
      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,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      0,
      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,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      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,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false,
      false
    ]
  },
  "jstor_metadata": {
    "identifier": "jstor-107921",
    "title": "Researches in Physical Astronomy",
    "authors": "John William Lubbock",
    "year": 1831,
    "volume": "121",
    "journal": "Philosophical Transactions of the Royal Society of London",
    "page_count": 51,
    "jstor_url": "https://www.jstor.org/stable/107921"
  }
}