{
  "id": "15212acb999d65c3231f18f322b76916d0a0c80d",
  "text": "XXII. Researches in physical astronomy. By John William Lubbock, Esq.\nFellow of the Royal Society.\n\nRead April 29, 1830.\n\nIn the first volume of the Mécanique Céleste, Laplace has given expressions for the variations of the elliptic constants, which are true when the square and higher powers of the disturbing force are neglected; and he has proved, upon the supposition that the planets move in the same direction, in orbits nearly circular and little inclined one to another, that the eccentricities and inclinations vary within small limits, thereby demonstrating within these conditions the stability of the planetary system. But these conditions are not necessary to the stability of a system of bodies, subject to the law of attraction, which obtains in our system. I have given in the following investigation the expressions for the variations of the elliptic constants, which are rigorously true whatever power of the disturbing force be retained; and it is easy to conclude from the form of their expressions, 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 no term which varies with the time; their variations are all periodic, and they oscillate therefore within certain limits. This theorem is no longer true if the planet moves in a resisting medium.\n\nI have also given some equations which obtain when an angle is taken for the independent variable, which in the elliptic movement is the eccentric anomaly, which are of remarkable simplicity, and which, as far as I know, have never been noticed, and the development of the disturbing function $R$ to the quantities involving the squares and products of the eccentricities inclusive.\n\nMDCCCXXX.\nLet \\( x, y, z \\) denote the rectangular co-ordinates\n\n\\[\n\\begin{align*}\nr & \\quad \\text{distance from the sun} \\\\\nr' & \\quad \\text{distance from the sun projected upon the plane } x, y \\\\\n*\\lambda & \\quad \\text{longitude reckoned upon the plane of its orbit} \\\\\n\\lambda' & \\quad \\text{longitude reckoned upon the plane } x y \\\\\ns & \\quad \\text{tangent of the latitude} \\\\\nv & \\quad \\text{a variable, which in the elliptic theory is the eccentric anomaly} \\\\\nm & \\quad \\text{the mass} \\\\\na & \\quad \\text{semiaxis major} \\\\\ne & \\quad \\text{eccentricity} \\\\\n\\varpi & \\quad \\text{longitude of the perihelion} \\\\\n\\varepsilon & \\quad \\text{longitude of the epoch} \\\\\n\\nu & \\quad \\text{longitude of the ascending node} \\\\\ni & \\quad \\text{inclination of the orbit to the plane } x y \\\\\n\\alpha & \\quad \\text{a constant quantity which accompanies } v \\\\\nM & \\quad \\text{the mass of the sun.}\n\\end{align*}\n\\]\n\n\\[ M + m = \\mu, \\quad \\sqrt{\\frac{\\mu}{a^3}} = n. \\]\n\n\\[ x = r' \\cos \\lambda', \\quad y = r' \\sin \\lambda', \\quad z = r's, \\quad \\text{and in the elliptic motion } r'^2 d \\lambda' = h dt. \\]\n\n\\[ R = m_i \\left\\{ \\frac{x x_i + y y_i + z z_i}{(x - x_i)^2 + (y - y_i)^2 + (z - z_i)^2} \\right\\}^{1/2} \\]\n\n\\[ = m_i \\left\\{ \\frac{r' \\{ \\cos (\\lambda' - \\lambda_i) + ss_i \\}}{r'^2 (1 + s^2)} \\right\\}^{1/2} \\]\n\n\\[ = \\frac{1}{\\{ r'^2 (1 + s^2) - 2 r' r'_i \\{ \\cos (\\lambda' - \\lambda_i) + ss_i \\} + r'^2 (1 + s^2) \\}^{1/2}} \\]\n\n* Laplace uses the letter \\( v \\) to denote longitude, \\( u \\) the eccentric anomaly, and \\( \\phi \\) the inclination of the orbit to a fixed plane; but as \\( v \\) is very frequently used to signify velocity, and \\( \\phi \\) geographical latitude, and as the letters of the Greek alphabet are generally used for angles, I have taken the letters \\( \\lambda, v, \\) and \\( i \\) for these quantities.\n\\[ \\frac{1}{r_i \\sqrt{1 + s_i^2}} + \\frac{r_i^3}{4 r_i^3} \\left\\{ 1 + 3 \\cos (2 \\lambda' - \\lambda_i) + 12 s_i \\cos (\\lambda' - \\lambda_i) - 2 s_i^3 \\right\\} \\]\n\n\\[ = -m_i \\left\\{ \\frac{r_i^3}{8 r_i^4} \\left\\{ 3 (1 - 4 s_i^2) \\cos (\\lambda' - \\lambda_i) + 5 \\cos (3 \\lambda' - 3 \\lambda_i) \\right\\} \\right. \\]\n\n\\[ + \\frac{r_i^4}{64 r_i^5} \\left\\{ 9 + 20 \\cos (2 \\lambda' - 2 \\lambda_i) + 35 \\cos (4 \\lambda' - 4 \\lambda_i) \\right\\} \\]\n\n\\[ \\frac{dR}{dr} = m_i \\left\\{ \\frac{\\cos (\\lambda' - \\lambda_i) + ss_i}{r_i^2 (1 + s_i^2)^{3/2}} + \\frac{r_i (1 + s_i^2) - r_i \\cos (\\lambda' - \\lambda_i) + ss_i}{r_i^2 (1 + s_i^2)^{3/2}} \\right\\} \\]\n\n\\[ \\frac{dR}{d\\lambda} = -m_i \\left\\{ \\frac{r_i \\sin (\\lambda' - \\lambda_i)}{r_i^2 (1 + s_i^2)^{3/2}} - \\frac{r_i \\sin (\\lambda' - \\lambda_i)}{r_i^2 (1 + s_i^2)^{3/2}} \\right\\} \\]\n\n\\[ \\frac{dR}{ds} = m_i \\left\\{ \\frac{r_i s_i}{r_i^2 (1 + s_i^2)^{3/2}} + \\frac{r_i^2 s_i - r_i r_i s_i}{r_i^2 (1 + s_i^2)^{3/2}} \\right\\} \\]\n\nLet \\( P \\) be the place of the planet \\( m \\), \\( S \\) the place of the sun, \\( SN \\) the intersection of the orbits of \\( m \\) and \\( m_p \\), \\( SL \\) the line from which longitudes are reckoned, \\( P_i \\), the projection of \\( P_i \\) upon the plane of the orbit of \\( P \\); then if the plane \\( xy \\) coincide with the orbit of \\( m \\), \\( SP = r \\), \\( PS_L = \\lambda \\), \\( NS_L = \\nu \\), \\( P_i SN + NS_L = \\lambda_i \\), \\( P_i SL = \\lambda_i \\), \\( SP_i = r_i \\), \\( SP_i = r_i = r_i (1 + s_i^2)^{1/2} \\)\n\n\\[ R = m_i \\left\\{ \\frac{SP \\times SP_i \\cos PSP_i}{SP_i^3} - \\frac{1}{\\{SP^2 - 2 SP \\times SP_i \\cos PSP_i + SP_i^2\\}^{1/2}} \\right\\} \\]\n\n\\[ \\frac{dR}{dr} = m_i \\left\\{ \\frac{SP_i \\cos PSP_i}{SP_i^3} + \\frac{SP - SP_i \\cos PSP_i}{\\{SP^2 - 2 SP \\times SP_i \\cos PSP_i + SP_i^2\\}^{1/2}} \\right\\} \\]\n\n\\[ \\frac{dR}{d\\lambda} = m_i \\left\\{ \\frac{SP \\times SP_i \\cos PSP_i}{SP_i^3} + \\frac{SP^2 - SP \\times SP_i \\cos PSP_i}{\\{SP^2 - 2 SP \\times SP_i \\cos PSP_i + SP_i^2\\}^{1/2}} \\right\\} \\]\n\n\\[ \\frac{dR}{ds} = m_i \\left\\{ \\frac{SP \\times SP_i \\sin PSP_i}{SP_i^3} - \\frac{SP \\times SP_i \\sin PSP_i}{\\{SP^2 - 2 SP \\times SP_i \\cos PSP_i + SP_i^2\\}^{1/2}} \\right\\} \\]\n\n\\[ \\cos (\\lambda - \\lambda_i) = \\cos (\\lambda - \\nu) \\cos (\\lambda_i - \\nu) + \\cos \\nu \\sin (\\lambda - \\nu) \\sin (\\lambda_i - \\nu) \\]\n\n\\[ = \\cos^2 \\frac{i}{2} \\cos (\\lambda - \\lambda_i) + \\sin^2 \\frac{i}{2} \\cos (\\lambda + \\lambda_i - 2 \\nu) \\]\n\n\\[ \\tan (\\lambda_i - \\nu) = \\cos \\nu \\tan (\\lambda_i - \\nu) \\]\n\n\\[ \\sin (\\lambda_i - \\nu) = \\frac{\\cos \\nu \\tan (\\lambda_i - \\nu)}{(1 + \\cos^2 \\nu \\tan^2 (\\lambda_i - \\nu))^{1/2}} \\]\n\\[\n\\cos (\\lambda_i' - \\nu) = \\frac{1}{(1 + \\cos^2 \\nu \\tan^2 (\\lambda_i - \\nu))^{1/2}}\n\\]\n\n\\[\nr r_i' \\sin (\\lambda - \\lambda_i') = r r_i' \\{\\sin \\lambda \\cos \\lambda_i' - \\cos \\lambda \\sin \\lambda_i'\\}\n\\]\n\n\\[\n\\cos \\lambda_i' = \\cos (\\lambda_i' - \\nu) \\cos \\nu - \\sin (\\lambda_i' - \\nu) \\sin \\nu\n\\]\n\n\\[\n= \\frac{\\cos (\\lambda_i - \\nu) \\cos \\nu - \\sin (\\lambda_i - \\nu) \\sin \\nu \\left(1 - 2 \\sin^2 \\frac{\\nu}{2}\\right)}{(1 - \\sin^2 \\nu \\sin^2 (\\lambda_i - \\nu))^{1/2}}\n\\]\n\n\\[\n\\cos \\lambda_i' + 2 \\sin^2 \\frac{\\nu}{2} \\sin (\\lambda_i - \\nu) \\sin \\nu\n\\]\n\n\\[\n= \\frac{\\cos \\lambda_i' + 2 \\sin^2 \\frac{\\nu}{2} \\sin (\\lambda_i - \\nu) \\sin \\nu}{(1 - \\sin^2 \\nu \\sin^2 (\\lambda_i - \\nu))^{1/2}}\n\\]\n\n\\[\n\\sin \\lambda_i' = \\sin (\\lambda_i' - \\nu) \\cos \\nu + \\cos (\\lambda_i' - \\nu) \\sin \\nu\n\\]\n\n\\[\n= \\frac{\\sin (\\lambda_i - \\nu) \\cos \\nu \\left(1 - 2 \\sin^2 \\frac{\\nu}{2}\\right) + \\cos (\\lambda_i - \\nu) \\sin \\nu}{(1 - \\sin^2 \\nu \\sin^2 (\\lambda_i - \\nu))^{1/2}}\n\\]\n\n\\[\n\\sin P_i S P_i' = \\sin \\nu \\sin (\\lambda_i - \\nu), \\quad r_i' = r_i \\cos P_i S P_i' = r_i \\left(1 - \\sin^2 \\nu \\sin^2 (\\lambda_i - \\nu)\\right)^{1/2}\n\\]\n\ntherefore,\n\n\\[\nr r_i' \\sin (\\lambda - \\lambda_i') = r r_i \\left\\{\\sin (\\lambda - \\lambda_i) + 2 \\sin^2 \\frac{\\nu}{2} \\sin (\\lambda_i - \\nu) \\cos (\\lambda - \\nu)\\right\\}\n\\]\n\n\\[\n= r r_i \\left\\{\\cos^2 \\frac{\\nu}{2} \\sin (\\lambda - \\lambda_i) + \\sin^2 \\frac{\\nu}{2} \\sin (\\lambda + \\lambda_i - 2 \\nu)\\right\\}\n\\]\n\nsimilarly it may be shown that\n\n\\[\nr r_i' \\cos (\\lambda - \\lambda_i') = r r_i \\left\\{\\cos^2 \\frac{\\nu}{2} \\cos (\\lambda - \\lambda_i) + \\sin^2 \\frac{\\nu}{2} \\cos (\\lambda + \\lambda_i - 2 \\nu)\\right\\}\n\\]\n\n\\[\nS P_i = \\frac{a_i (1 - e_i^2)}{1 + e_i \\cos (P_i S N + NSL - \\varpi_i)} = \\frac{a_i (1 - e_i^2)}{1 + e_i \\cos (\\lambda_i - \\varpi_i)}\n\\]\n\n\\[\n\\frac{d^2 x}{dt^2} + \\frac{\\mu x}{r^3} + \\left(\\frac{d R}{dx}\\right) = 0, \\quad \\frac{d^2 y}{dt^2} + \\frac{\\mu y}{r^3} + \\left(\\frac{d R}{dy}\\right) = 0, \\quad \\frac{d^2 z}{dt^2} + \\frac{\\mu z}{r^3} + \\left(\\frac{d R}{dz}\\right) = 0\n\\]\n\n\\[\n\\frac{r^2 d \\lambda^2 + dr^2 (1 + s^2) + 2 r' s dr' ds + r'^2 ds^2}{dt^2} - \\frac{2 \\mu}{r'(1 + s^2)^{3/2}} + \\frac{\\mu}{a} + 2 \\int d R = 0\n\\]\n\n\\(d R\\) being the differential of \\(R\\) with regard only to the co-ordinates of the planet \\(m\\).\n\n\\[\n\\frac{d^2 r'}{dt^2} (1 + s^2) - r' d \\lambda^2 + 2 s dr' ds + r'^2 ds^2 + \\frac{\\mu}{r'(1 + s^2)^{3/2}} + \\left(\\frac{d R}{dr'}\\right) = 0\n\\]\n\n\\[\n\\frac{d}{dt} \\left(\\frac{r'^2 d \\lambda'}{dt}\\right) + \\left(\\frac{d R}{d \\lambda'}\\right) = 0\n\\]\n\n\\[\n\\frac{r'^2 d^2 r'}{dt^2} + 2 r' dr' ds + r'^2 ds^2 + \\frac{s}{r'(1 + s^2)^{3/2}} + \\left(\\frac{d R}{ds}\\right) = 0\n\\]\n\\[\n\\frac{r^2 d^2 r}{dt^2} - r^2 \\frac{d^2 \\lambda}{dt^2} + \\frac{\\mu}{r(1+s^2)^{3/2}} + r \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) = 0\n\\]\n\n\\[\n\\frac{r^2 d^2 s}{dt^2} + 2r \\frac{dr}{dt} \\frac{ds}{dt} - r^2 \\frac{d^2 \\lambda}{dt^2} + (1+s^2) \\left( \\frac{dR}{ds} \\right) - r \\left( \\frac{dR}{dr} \\right) = 0\n\\]\n\n\\[\n\\frac{d^2 r}{dt^2} - \\frac{\\mu}{r(1+s^2)^{3/2}} + \\frac{\\mu}{a} + 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) = 0\n\\]\n\nMaking \\( \\lambda \\) the independent variable instead of \\( t \\),\n\n\\[\n\\frac{d^2 r}{dt^2} - 2 \\frac{dr}{dt} \\frac{d^2 t}{dt^2} - r^2 \\frac{d^2 \\lambda}{dt^2} + \\frac{\\mu}{r(1+s^2)^{3/2}} + r \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) = 0\n\\]\n\n\\[\nr^4 \\frac{d^2 \\lambda}{dt^2} = h^2 - 2 \\int r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda, \\quad h \\text{ being a constant},\n\\]\n\n\\[\n4r^3 \\frac{dr}{dt} \\frac{d^2 \\lambda}{dt^2} - 2r^4 \\frac{d\\lambda}{dt} \\frac{d^2 t}{dt^2} = -2r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda\n\\]\n\n\\[\n\\left\\{ \\frac{d^2}{d\\lambda^2} + \\frac{1}{r^2} \\right\\} \\left\\{ 1 - \\frac{2}{h^2} \\int r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda \\right\\} - \\frac{\\mu}{h^2(1+s^2)^{3/2}}\n\\]\n\n\\[\n- \\frac{r}{h^2} \\left\\{ r \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\frac{1}{r} \\left( \\frac{dR}{d\\lambda} \\right) \\frac{dr}{d\\lambda} \\right\\} = 0\n\\]\n\n\\[\n\\left\\{ \\frac{d^2 s}{d\\lambda^2} + s \\right\\} \\left\\{ 1 - \\frac{2}{h^2} \\int r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda \\right\\}\n\\]\n\n\\[\n+ \\frac{r^2}{h^2} \\left\\{ (1+s^2) \\left( \\frac{dR}{ds} \\right) - r \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\frac{ds}{d\\lambda} \\right\\} = 0\n\\]\n\nWhen the disturbing force is neglected\n\n\\[\n\\frac{d^2 r}{dt^2} = 0, \\quad \\frac{d^2 \\lambda}{dt^2} = 0, \\quad \\frac{d^2 s}{d\\lambda^2} + s = 0\n\\]\n\nof which equations the integrals are,\n\n\\[\nr^2 d\\lambda = h dt, \\quad \\frac{1}{r} = \\frac{\\mu \\cos i^2}{h^2} \\left\\{ (1+s^2)^{1/2} + e \\cos (\\lambda - \\omega) \\right\\}\n\\]\n\n\\[\ns = \\tan i \\sin (\\lambda - \\nu).\n\\]\n\nIf \\( dt = \\sqrt{\\frac{a}{\\mu}} r \\frac{dv}{v} \\), and \\( v \\) be taken for the independent variable\n\n\\[\nr^2 \\frac{d^2 r}{dt^2} - r \\frac{dr}{dt} \\frac{d^2 t}{dt^2} - r^2 \\frac{d^2 \\lambda}{dt^2} + \\frac{\\mu}{r(1+s^2)^{3/2}} + r \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) = 0\n\\]\n\\[ \\frac{d^2 t}{d v^2} = \\left( \\frac{a}{\\mu} \\right)^{\\frac{1}{2}} d r' d v \\]\n\n\\[ \\frac{d^2 r'}{d v^2} + \\left( \\frac{d r' s}{d v} \\right)^2 - \\frac{a}{(1 + s^2)^{\\frac{3}{2}}} + r' + \\frac{ar'}{\\mu} \\left\\{ 2 \\int d R + r' \\left( \\frac{d R}{d r'} \\right) - s \\left( \\frac{d R}{d s} \\right) \\right\\} = 0. \\]\n\nIf the orbit of the planet \\( m \\) coincide with the plane \\( x y \\), \\( s \\) is of the order of the disturbing force, of which therefore neglecting the square, \\( r' = r \\), \\( \\lambda' = \\lambda \\)\n\n\\[ \\frac{d^2 r}{d v^2} - a + r + \\frac{ar}{\\mu} \\left\\{ 2 \\int d R + r \\left( \\frac{d R}{d r} \\right) \\right\\} = 0. \\]\n\nWhen the disturbing force is neglected, the integral of this equation is\n\n\\[ r = a \\{ 1 - e \\cos (\\nu - \\alpha) \\}, \\nu \\text{ being the eccentric anomaly.} \\]\n\n\\[ n t + \\varepsilon - \\varpi = \\nu - \\alpha - e \\sin (\\nu - \\alpha) \\]\n\n\\[ \\tan \\frac{\\lambda - \\varpi}{2} = \\left\\{ \\frac{1 + e}{1 - e} \\right\\}^{\\frac{1}{2}} \\tan \\frac{\\nu - \\alpha}{2}. \\]\n\nIf \\( Q \\) be put for the quantity \\( \\frac{ar}{\\mu} \\left\\{ 2 \\int d R + r \\left( \\frac{d R}{d r} \\right) \\right\\} \\) and the constant \\( \\alpha \\), which may afterwards be replaced, be omitted for the present,\n\n\\[ r = a \\{ 1 - e \\cos \\nu \\} - \\sin \\nu \\int Q \\cos \\nu d \\nu + \\cos \\nu \\int Q \\sin \\nu d \\nu \\]\n\n\\[ d t = \\sqrt{\\frac{a}{\\mu}} \\left\\{ a \\{ 1 - e \\cos \\nu \\} - \\sin \\nu \\int Q \\cos \\nu d \\nu + \\cos \\nu \\int Q \\sin \\nu d \\nu \\right\\} d \\nu \\]\n\n\\[ n t + \\varepsilon - \\varpi = \\nu - e \\sin \\nu - \\frac{1}{a} \\left\\{ \\int Q d \\nu - \\cos \\nu \\int Q \\cos \\nu d \\nu - \\sin \\nu \\int Q \\sin \\nu d \\nu \\right\\} \\]\n\nIf \\( \\nu = f (n t + \\varepsilon - \\varpi) \\) in the elliptic theory, then neglecting the square of the disturbing force,\n\n\\[ \\nu = f (n t + \\varepsilon - \\varpi) + \\frac{d.f(n t + \\varepsilon - \\varpi)}{d t} \\frac{1}{a} \\left\\{ \\int Q d \\nu - \\cos \\nu \\int Q \\cos \\nu d \\nu - \\sin \\nu \\int Q \\sin \\nu d \\nu \\right\\} d \\nu \\]\n\nIf \\( \\delta \\nu, \\delta r \\) denote the values of those parts of \\( r \\) and \\( \\nu \\) which are due to the first power of the disturbing force,\n\n\\[ \\delta \\nu = \\frac{d \\nu}{d nd t} \\left\\{ \\int Q d \\nu - \\cos \\nu \\int Q \\cos \\nu d \\nu - \\sin \\nu \\int Q \\sin \\nu d \\nu \\right\\} \\]\n\n\\[ \\delta r = a e \\sin \\nu \\delta \\nu - \\sin \\nu \\int Q \\cos \\nu d \\nu + \\cos \\nu \\int Q \\sin \\nu d \\nu \\]\n\n\\[ = \\frac{\\sin \\nu}{1 - e \\cos \\nu} \\int Q \\{ e - \\cos \\nu \\} d \\nu - \\frac{e - \\cos \\nu}{1 - e \\cos \\nu} \\int Q \\sin \\nu d \\nu. \\]\n\nIn the elliptic theory,\n\n\\[ d \\nu = \\frac{n d t}{1 - e \\cos \\nu}, \\quad \\frac{\\cos \\nu - e}{1 - e \\cos \\nu} = \\cos \\lambda, \\quad \\frac{\\sin \\nu (1 - e^2)^{\\frac{1}{2}}}{1 - e \\cos \\nu} = \\sin \\lambda. \\]\ntherefore,\n\n\\[ \\delta r = \\frac{a \\cos \\lambda \\int r \\sin \\lambda \\left\\{ 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) \\right\\} dt - a \\sin \\lambda \\int r \\cos \\lambda \\left\\{ 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) \\right\\} dt}{\\mu (1 - e^2)^{\\frac{3}{2}}} \\]\n\nwhich is the equation X of the Mécanique Céleste, vol. i. p. 258.\n\nMultiplying the equation of p. 114, l. 2, by \\( \\frac{2r^2 d\\lambda}{dt} \\), and integrating, \\( \\frac{r^2 d\\lambda^2}{dt^2} \\)\n\n\\[ + 2 \\int r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda = h_0^2, \\text{if } h^2 = h_0^2 - 2 \\int r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda, \\text{h being variable, and } h_0 \\text{ the value of } h \\text{ at a given epoch}, \\]\n\n\\[ h d h = - r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda, r^2 d\\lambda = h d t, \\text{and making } \\lambda \\]\n\nthe independent variable instead of \\( t \\),\n\n\\[ 2 r^2 d r^2 d\\lambda = d h d t + h d^2 t \\]\n\n\\[ r^2 \\frac{d^2 r^2}{dt^2} - r^2 \\frac{dr^2 dt}{dt^2} - r^2 \\frac{d\\lambda^2}{dt^2} + \\frac{\\mu}{r^2(1 + s^2)^{\\frac{3}{2}}} + r^2 \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) = 0 \\]\n\n\\[ r^2 \\left( \\frac{d^2 r^2}{dt^2} \\right) + r^2 \\frac{dr^2 dt}{dt^2} - 2 r^2 dr^2 d\\lambda - r^2 \\frac{d\\lambda^2}{dt^2} + \\frac{\\mu}{r^2(1 + s^2)^{\\frac{3}{2}}} + r^2 \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) = 0 \\]\n\n\\[ \\frac{d^2}{d\\lambda^2} \\frac{1}{r^2} + \\frac{1}{r^2} - \\frac{\\mu}{h^2(1 + s^2)^{\\frac{3}{2}}} - \\frac{r^2}{h^2} \\left\\{ r^2 \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\frac{1}{r^2} \\left( \\frac{dR}{d\\lambda} \\right) \\frac{dr^2}{d\\lambda} \\right\\} = 0 \\]\n\n\\[ \\frac{d^2 s}{d\\lambda^2} + s + \\frac{r^2}{h^2} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\frac{ds}{d\\lambda} \\right\\} = 0 \\]\n\nIf all the constants in the elliptic integrals are supposed to vary, subject to the condition that they still satisfy these differential equations, and that the form of the first differential coefficients \\( \\frac{dr^2}{d\\lambda}, \\frac{ds}{d\\lambda} \\) remains unaltered,\n\n\\[ \\frac{d}{d\\lambda} \\frac{1}{r^2} = \\frac{\\mu \\cos i^2}{h^2} \\left\\{ \\frac{s}{(1 + s^2)^{\\frac{3}{2}}} \\frac{ds}{d\\lambda} - e \\sin (\\lambda' - \\varpi) \\right\\}, \\quad \\frac{ds}{d\\lambda} = \\tan i \\cos (\\lambda' - \\nu) \\]\n\n\\[ - 2 \\left\\{ (1 + s^2)^{\\frac{3}{2}} + e \\cos (\\lambda' - \\varpi) \\right\\} \\left\\{ h^2 \\sin i d i + \\cos i h d h \\right\\} \\]\n\n\\[ + h^2 \\cos i \\cos (\\lambda' - \\varpi) d e + h^2 \\cos i e \\sin (\\lambda' - \\varpi) d \\varpi = 0 \\]\n\n\\[ - 2 \\left\\{ \\frac{s \\tan i \\cos (\\lambda' - \\nu)}{(1 + s^2)^{\\frac{3}{2}}} - e \\sin (\\lambda' - \\varpi) \\right\\} \\left\\{ h^2 \\sin i d i + \\cos i h d h \\right\\} \\]\n\n\\[ - h^2 \\cos i \\sin (\\lambda' - \\varpi) d e + h^2 \\cos i e \\cos (\\lambda' - \\varpi) d \\varpi \\]\n\n\\[ - \\frac{\\cos i s r^2}{(1 + s^2)^{\\frac{3}{2}}} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - r^2 \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda' \\]\n\\[ - \\frac{h^2 r}{\\mu \\cos i} \\left\\{ r' \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\frac{1}{r'} \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{dr}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\n\n\\[ \\sin (\\lambda' - \\nu) \\frac{d\\iota}{\\cos i^2} - \\tan \\iota \\cos (\\lambda' - \\nu) d\\nu = 0 \\]\n\n\\[ \\cos (\\lambda' - \\nu) \\frac{d\\iota}{\\cos i^2} + \\tan \\iota \\sin (\\lambda' - \\nu) d\\nu + \\frac{r^2}{h^2} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\n\n\\[ h d h = - r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda \\]\n\nWhence by elimination,\n\n\\[ h^2 \\cos i d e + 2 \\left\\{ (1 + s^2)^{\\frac{1}{2}} \\cos (\\lambda' - \\varphi) + e - \\frac{\\tan \\iota \\cos (\\lambda' - \\nu) \\sin (\\lambda' - \\varphi)}{(1 + s^2)^{\\frac{1}{2}}} \\right\\} \\]\n\n\\[ \\left\\{ \\sin \\iota \\cos^2 i r^2 \\cos (\\lambda' - \\nu) \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} + \\cos i r^2 \\left( \\frac{dR}{d\\lambda} \\right) \\right\\} d\\lambda \\]\n\n\\[ + \\frac{s \\cos i r^2 \\sin (\\lambda' - \\varphi)}{(1 + s^2)^{\\frac{1}{2}}} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - r' \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda \\]\n\n\\[ + \\frac{h^2 r}{\\mu \\cos i} \\sin (\\lambda' - \\varphi) \\left\\{ r' \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\frac{1}{r'} \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{dr}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\n\n\\[ h^2 \\cos i d \\varphi + 2 \\left\\{ (1 + s^2)^{\\frac{1}{2}} \\sin (\\lambda' - \\varphi) + \\frac{s \\tan \\iota \\cos (\\lambda' - \\nu) \\cos (\\lambda' - \\varphi)}{(1 + s^2)^{\\frac{1}{2}}} \\right\\} \\]\n\n\\[ \\left\\{ \\sin \\iota \\cos^2 i r^2 \\cos (\\lambda' - \\nu) \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} + \\cos i r^2 \\left( \\frac{dR}{d\\lambda} \\right) \\right\\} d\\lambda \\]\n\n\\[ - \\frac{s \\cos i r^2 \\cos (\\lambda' - \\varphi)}{(1 + s^2)^{\\frac{1}{2}}} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - r' \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda \\]\n\n\\[ - \\frac{h^2 r}{\\mu \\cos i} \\cos (\\lambda' - \\varphi) \\left\\{ r' \\left( \\frac{dR}{dr} \\right) - s \\left( \\frac{dR}{ds} \\right) - \\frac{1}{r'} \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{dr}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\n\n\\[ d\\iota + \\frac{r^2 \\cos i^2 \\cos (\\lambda' - \\nu)}{h^2} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\n\n\\[ d\\nu + \\frac{r^2 \\sin (\\lambda' - \\nu)}{h^2 \\tan \\iota} \\left\\{ (1 + s^2) \\left( \\frac{dR}{ds} \\right) - s \\left( \\frac{dR}{dr} \\right) - \\left( \\frac{dR}{d\\lambda} \\right) \\left( \\frac{ds}{d\\lambda} \\right) \\right\\} d\\lambda = 0 \\]\nIf \\( x = \\frac{r \\cos \\lambda}{(1 + s^2)^{\\frac{1}{2}}} \\), \\( y = \\frac{r \\sin \\lambda}{(1 + s^2)^{\\frac{1}{2}}} \\), \\( z = \\frac{rs}{(1 + s^2)^{\\frac{1}{2}}} \\),\n\n\\[\n\\frac{d^2 r}{dt^2} + \\frac{r^2}{1 + s^2} \\left( \\frac{ds^2}{dt^2} + d\\lambda^2 \\right) - \\frac{2\\mu}{r} + \\frac{\\mu}{a} + 2 \\int dR = 0\n\\]\n\n\\[\n\\frac{d^2 r}{dt^2} - \\frac{r}{1 + s^2} \\left( \\frac{ds^2}{dt^2} + d\\lambda^2 \\right) + \\frac{\\mu}{r^3} + \\left( \\frac{dR}{dr} \\right) = 0\n\\]\n\n\\[\nd \\cdot \\frac{r^2}{1 + s^2} \\cdot \\frac{d\\lambda}{dt} + \\left( \\frac{dR}{d\\lambda} \\right) = 0\n\\]\n\n\\[\n\\frac{r^3 ds^2 + 2r dr ds}{(1 + s^2)} + (1 + s^2)^2 \\left( \\frac{dR}{ds} \\right) = 0\n\\]\n\nOf which equations the integrals are\n\n\\[\n\\frac{r^3}{1 + s^2} \\cdot d\\lambda = h dt, \\quad \\frac{1}{r} = \\frac{\\mu \\cos i^2}{h^2(1 + s^2)^{\\frac{1}{2}}} \\left\\{ (1 + s^2)^{\\frac{1}{2}} + e \\cos (\\lambda - \\varpi) \\right\\}\n\\]\n\n\\( s = \\tan i \\sin (\\lambda - \\nu) \\)\n\nIf \\( dt = \\sqrt{\\frac{a}{\\mu}} r dv \\), and \\( v \\) be taken for the independent variable\n\n\\[\n\\frac{d^2 r}{dv^2} - a + r + \\frac{ar}{\\mu} \\left\\{ 2 \\int dR + r \\left( \\frac{dR}{dr} \\right) \\right\\} = 0\n\\]\n\n\\( r = a \\{1 - e' \\cos (v - \\alpha)\\} \\) in the elliptic motion.\n\n\\( e' \\) is accented for the present in order to distinguish it from \\( e \\).\n\nIf the constants in the elliptic integrals are supposed to vary, subject to the condition that they still satisfy these differential equations, and that the form of the first differential coefficient \\( \\frac{dr}{dv} \\) remains unaltered,\n\n\\[\nd^2 t = \\frac{rdvda}{2\\sqrt{a\\mu}} + \\sqrt{\\frac{a}{\\mu}} dr dv, \\quad \\frac{r^2}{dt^2} = \\frac{rdr^2}{dt^2} - \\frac{rdr^2}{dt^2}\n\\]\n\n\\[\n\\frac{d^2 r}{dv^2} - a + r + \\frac{ar^2}{\\mu} \\left( \\frac{dR}{dr} \\right) - \\frac{drda}{2adv^2} = 0\n\\]\n\n\\[\n(1 - e' \\cos (v - \\alpha)) da - a \\cos (v - \\alpha) de' - ae' \\sin (v - \\alpha) da = 0\n\\]\n\n\\( e' \\sin (v - \\alpha) da + a \\sin (v - \\alpha) de' - ae' \\cos (v - \\alpha) da \\)\n\nMDCCXXX.\n\\[ \\frac{a r^2}{\\mu} \\left( \\frac{dR}{dr} \\right) d\\nu + \\frac{a^2 e}{\\mu} \\sin (\\nu - \\alpha) dR = 0 \\]\n\n\\[ da = -\\frac{2a^2}{\\mu} dR \\]\n\n\\[ de' = \\frac{a}{\\mu} \\{ 2e' - 2\\cos(\\nu - \\alpha) - e'\\sin(\\nu - \\alpha)^2 \\} dR - \\frac{r^3}{\\mu} \\sin(\\nu - \\alpha) \\left( \\frac{dR}{dr} \\right) d\\nu \\]\n\n\\[ e'da = \\frac{a}{\\mu} \\sin(\\nu - \\alpha) \\{ e'\\cos(\\nu - \\alpha) - 2 \\} dR + \\frac{r^3}{\\mu} \\cos(\\nu - \\alpha) \\left( \\frac{dR}{dr} \\right) d\\nu \\]\n\n\\[ \\int n dt + \\varepsilon - \\varpi = \\nu - \\alpha - e'\\sin(\\nu - \\alpha) \\]\n\n\\[ d\\varepsilon - d\\varpi = -da - \\sin(\\nu - \\alpha) de' + e'\\cos(\\nu - \\alpha) da. \\]\n\nThe equation of condition which obtains between the constants \\(a, e, \\varpi, \\nu, \\iota\\) and \\(h\\) may be found from the equation\n\n\\[ \\frac{dr^3}{dt^2} + \\frac{r^3}{1+s^3} \\left( \\frac{ds^2}{1+s^3} + d\\lambda^3 \\right) - \\frac{2\\mu}{r} + \\frac{\\mu}{a} = 0 \\]\n\nwhich gives\n\n\\[ \\frac{\\mu \\cos \\iota^2}{h^3} \\left\\{ e^2 \\cos^2 \\iota \\sin^2(\\nu - \\varpi) - \\{ 1 + e \\cos(\\nu - \\varpi) \\} \\{ 1 - e \\cos(\\nu - \\varpi) \\} \\right\\} + \\frac{1}{a} = 0 \\]\n\nEquating the values of \\(r\\) which have been found,\n\n\\[ \\frac{h^3 \\sqrt{1+s^3}}{\\mu \\cos \\iota^3 \\{ \\sqrt{1+s^3} + e \\cos(\\lambda - \\varpi) \\}} = a \\{ 1 - e'\\cos(\\nu - \\alpha) \\} \\]\n\nsince the origin of the angle \\(\\nu\\) is arbitrary, we may suppose \\(\\lambda - \\varpi = 0\\), and \\(\\nu - \\alpha = 0\\) at the same time,\n\nso that\n\n\\[ \\frac{h^3 \\sqrt{1+\\tan^2 \\sin(\\varpi - \\nu)^3}}{\\mu \\cos \\iota^3 \\{ \\sqrt{1+\\tan^2 \\sin(\\varpi - \\nu)^3} + e \\}} = a (1 - e') \\]\n\nAll the equations which have hitherto been proved are rigorously true, whatever powers of the disturbing function be retained. They are susceptible of simplification when the square and higher powers of the disturbing function are neglected: in this case, if the orbit be supposed to coincide with the plane \\(x y\\), \\(\\tan \\iota = 0\\), and if the longitude be reckoned from the perihelion of the planet P,\n\n\\[ h^2 de + \\{ 2 \\cos \\lambda + e + e \\cos \\lambda^2 \\} r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda + \\frac{h^3 r^3 \\sin \\lambda}{\\mu} \\left( \\frac{dR}{dr} \\right) d\\lambda = 0 \\]\n\\[ h^2 e \\sin \\varpi + \\{2 + \\cos \\lambda\\} \\sin \\lambda r^2 \\left( \\frac{dR}{d\\lambda} \\right) d\\lambda - h^2 r^2 \\cos \\lambda \\left( \\frac{dR}{dr} \\right) d\\lambda = 0 \\]\n\n\\[ d\\iota + r^2 \\frac{\\cos (\\lambda - \\nu)}{\\mu} \\left( \\frac{dR}{ds} \\right) d\\lambda = 0 \\]\n\n\\[ d\\nu + r^2 \\frac{\\sin (\\lambda - \\nu)}{\\mu} \\left( \\frac{dR}{ds} \\right) d\\lambda = 0 \\]\n\nby equation, p. 336, line 12, \\( \\mu a (1 - e^2) = h^2 \\), and by equation, line 17, \\( e = e' \\).\n\n\\[ d e = \\frac{a}{\\mu} \\left\\{ 2e - 2 \\cos (\\nu - \\alpha) - e \\sin (\\nu - \\alpha)^2 \\right\\} dR - \\frac{r^2}{\\mu} \\sin (\\nu - \\alpha) \\left( \\frac{dR}{dr} \\right) d\\nu \\]\n\n\\[ e d\\alpha = \\frac{a}{\\mu} \\left\\{ e \\cos (\\nu - \\alpha) - 2 \\right\\} \\sin (\\nu - \\alpha) dR + \\frac{r^2}{\\mu} \\cos (\\nu - \\alpha) \\left( \\frac{dR}{dr} \\right) d\\nu \\]\n\n\\[ \\int n dt + \\varepsilon - \\varpi = \\nu - \\alpha - e \\sin (\\nu - \\alpha) \\]\n\n\\[ d\\varepsilon - d\\varpi = - \\{1 - e \\cos (\\nu - \\alpha)\\} d\\alpha - \\sin (\\nu - \\alpha) d e \\]\n\n\\[ d\\varepsilon - d\\varpi = - \\frac{r^2}{a^2 \\sqrt{1 - e^2}} d\\varpi - \\left\\{ \\frac{r}{a(1 - e^2)} + 1 \\right\\} \\sin (\\nu - \\alpha) d e \\]\n\n\\[ \\frac{h^2}{\\mu} \\frac{1 + e \\cos (\\lambda - \\varpi)}{1 + e \\cos (\\lambda - \\varpi)} = a \\{1 - e \\cos (\\nu - \\alpha)\\} \\]\n\n\\[ d e = - \\frac{a ndt}{\\mu \\sqrt{1 - e^2}} \\left\\{ 2 \\cos (\\lambda - \\varpi) + e + e \\cos (\\lambda - \\varpi)^2 \\right\\} \\left( \\frac{dR}{d\\lambda} \\right) \\]\n\n\\[ - \\frac{a^2 ndt \\sqrt{1 - e^2}}{\\mu} \\sin (\\lambda - \\varpi) \\left( \\frac{dR}{dr} \\right) \\]\n\n\\[ e d\\varpi = - \\frac{a ndt}{\\mu \\sqrt{1 - e^2}} \\left\\{ 2 + e \\cos (\\lambda - \\varpi) \\right\\} \\sin (\\lambda - \\varpi) \\left( \\frac{dR}{d\\lambda} \\right) \\]\n\n\\[ + \\frac{a^2 ndt \\sqrt{1 - e^2}}{\\mu} \\cos (\\lambda - \\varpi) \\left( \\frac{dR}{dr} \\right) \\]\n\n\\[ d\\varepsilon - d\\varpi = - \\sqrt{1 - e^2} d\\varpi + \\frac{2 a^2 n (1 - e^2) dt}{\\mu \\{1 + e \\cos (\\lambda - \\varpi)\\}} \\left( \\frac{dR}{dr} \\right) \\]\n\nIf the longitudes be reckoned from the perihelion of the planet P,\n\n\\[ \\frac{de}{dt} = - \\frac{an(\\cos \\lambda + e)}{\\mu \\sqrt{1 - e^2}} \\left( \\frac{dR}{d\\lambda} \\right) - \\frac{m_i a^2 n \\sqrt{1 - e^2}}{\\mu} \\left\\{ \\frac{r'_i \\sin \\lambda'_i}{r'_i^3} + \\frac{r \\sin \\lambda - r'_i \\sin \\lambda'_i}{\\{r^2 - 2 rr'_i \\cos (\\lambda - \\lambda'_i) + r'_i^2\\}^{3/2}} \\right\\} \\]\n\\[\n\\frac{d\\varpi}{dt} = -\\frac{an \\sin \\lambda}{\\mu \\sqrt{1-e^2}} \\left( \\frac{dR}{d\\lambda} \\right) - \\frac{m_1 a^2 n \\sqrt{1-e^2}}{\\mu} \\left\\{ \\frac{r_i' \\cos \\lambda_i'}{r_i'^3} + \\frac{r \\cos \\lambda - r_i' \\cos \\lambda_i'}{\\{r^2 - 2rr_i' \\cos (\\lambda - \\lambda_i') + r_i'^2\\}^{3/2}} \\right\\}\n\\]\n\nand since \\( ndt = \\frac{r}{a} dv \\),\n\n\\[\n\\frac{de}{dv} = -\\frac{r(\\cos \\lambda + e)}{\\mu \\sqrt{1-e^2}} \\left( \\frac{dR}{d\\lambda} \\right) - \\frac{m_1 ar \\sqrt{1-e^2}}{\\mu} \\left\\{ \\frac{r_i' \\sin \\lambda_i'}{r_i'^3} + \\frac{r \\sin \\lambda - r_i' \\sin \\lambda_i'}{\\{r^2 - 2rr_i' \\cos (\\lambda - \\lambda_i') + r_i'^2\\}^{3/2}} \\right\\}\n\\]\n\n\\[\n\\frac{ed\\varpi}{dv} = -\\frac{r \\sin \\lambda}{\\mu \\sqrt{1-e^2}} \\left( \\frac{dR}{d\\lambda} \\right) + \\frac{m_1 ar \\sqrt{1-e^2}}{\\mu} \\left\\{ \\frac{r_i' \\cos \\lambda_i'}{r_i'^3} + \\frac{r \\cos \\lambda - r_i' \\cos \\lambda_i'}{\\{r^2 - 2rr_i' \\cos (\\lambda - \\lambda_i') + r_i'^2\\}^{3/2}} \\right\\}\n\\]\n\n\\[\n\\frac{d\\varepsilon - d\\varpi}{dv} = -\\frac{r^2}{a^2 \\sqrt{1-e^2}} \\frac{d\\varpi}{dv} - \\left\\{ \\frac{r}{a(1-e^2)} + 1 \\right\\} \\sin v \\frac{de}{dv}\n\\]\n\n\\[\nda = -\\frac{2a^2}{\\mu} dR\n\\]\n\nand since \\( r^2 d\\lambda = \\sqrt{\\mu a (1-e^2)} dt = a \\sqrt{1-e^2} r dv \\)\n\n\\[\n\\frac{d\\lambda}{dv} + \\frac{a \\sqrt{1-e^2} r \\cos \\lambda}{\\mu} \\left( \\frac{dR}{ds} \\right) = 0\n\\]\n\n\\[\n\\frac{dv}{dv} + \\frac{a \\sqrt{1-e^2} r \\sin \\lambda}{\\mu} \\left( \\frac{dR}{ds} \\right) = 0.\n\\]\n\nThe last six equations serve to determine the perturbations of a comet.\n\nLet \\( (\\Delta e)_n \\) be the variation of any element \\( e \\) during the variation \\( \\Delta v \\) of \\( v \\) at any given epoch \\( n \\), neglecting the square and higher powers of \\( \\Delta v \\),\n\n\\[\n(\\Delta e)_n = \\left( \\frac{de}{dv} \\right)_n \\Delta v\n\\]\n\nIf the values of \\( (\\Delta e)_n \\) be calculated for the epochs 0, 1, 2 ... \\( m \\) corresponding to the values \\( v, v + i \\Delta v, v + 2i \\Delta v, \\) &c., differing from each other by \\( i \\Delta v \\), then the whole variation \\( (\\delta e) \\) of \\( v \\) corresponding to the variation \\( i \\Delta v \\) of \\( v \\)\n\n\\[\n= i \\left\\{ (\\Delta e)_0 + (\\Delta e)_1 + \\ldots + (\\Delta e)_m \\right\\}\n\\]\n\n\\[\n+ \\frac{i-1}{2} \\left\\{ (\\Delta e)_m - (\\Delta e)_0 \\right\\} - \\frac{(i-1)(i+1)}{12i} \\left\\{ (\\Delta^2 e)_m - (\\Delta^2 e)_0 \\right\\} + \\text{&c.}\n\\]\n\n\\[\n= i \\left\\{ \\left( \\frac{de}{dv} \\right)_0 + \\left( \\frac{de}{dv} \\right)_1 + \\ldots + \\left( \\frac{de}{dv} \\right)_m \\right\\} \\Delta v\n\\]\n\n\\[\n+ \\frac{i-1}{2} \\left\\{ \\left( \\frac{de}{dv} \\right)_m - \\left( \\frac{de}{dv} \\right)_0 \\right\\} \\Delta v - \\frac{(i-1)(i+1)}{12i} \\left\\{ \\Delta \\left( \\frac{de}{dv} \\right)_m - \\Delta \\left( \\frac{de}{dv} \\right)_0 \\right\\} \\Delta v.\n\\]\nWhen the interval $\\Delta \\nu$ is indefinitely diminished, $i \\Delta \\nu$ is still equal to the variation of $\\nu$ between the epochs for which the quantities $\\left(\\frac{d e}{d \\nu}\\right)_0$, $\\left(\\frac{d e}{d \\nu}\\right)_1$, &c. are calculated, and\n\n$$\\delta e = i \\Delta \\nu \\left\\{ \\left(\\frac{d e}{d \\nu}\\right)_0 + \\left(\\frac{d e}{d \\nu}\\right)_1 + \\ldots + \\left(\\frac{d e}{d \\nu}\\right)_{m-1} + \\frac{1}{2} \\left\\{ \\left(\\frac{d e}{d \\nu}\\right)_m - \\left(\\frac{d e}{d \\nu}\\right)_0 \\right\\} \\right\\}$$\n\n$$- \\frac{1}{12} \\left\\{ \\Delta \\left(\\frac{d e}{d \\nu}\\right)_m - \\Delta \\left(\\frac{d e}{d \\nu}\\right)_0 \\right\\} + \\text{&c.} \\right\\}$$\n\nIf the radius be taken for unity, and $i \\Delta \\nu$ is the $m$th part of the circumference; $i \\Delta \\nu = \\frac{2 \\times 3.14159}{m}$, or, in other words, the resulting values of $\\delta e$ and $\\delta a$ in the equations given above must be multiplied by $2 \\times 3.14159$, and divided by $360^\\circ$ expressed in the same unit as $i \\Delta \\nu$.\n\n$n$ is equal to the angular circumference divided by the periodic time expressed in the same unit as $t$; so that if a degree be taken as the unity of angular circumference, $n = 360^\\circ$ divided by the periodic time expressed in the same unit as $t$.\n\nIn the elliptic movement or first approximation\n\n$$\\lambda = n t + \\varepsilon + \\text{a series of sines of arcs multiples of } n t \\text{ &c.}$$\n\n$$\\lambda_i = n_i t + \\varepsilon_i + \\text{a series of sines of arcs multiples of } n_i t \\text{ &c.}$$\n\n$$\\frac{a}{r} = \\text{constant} + \\text{a series of cosines.}$$\n\n$s = \\text{a series of sines.}$\n\n$s_i = \\text{a series of sines.}$\n\nThese values being substituted in the equations of p. 334 give $\\frac{d e}{d \\lambda}$, $\\frac{d h}{d \\lambda}$, and $\\frac{d i}{d \\lambda}$ each equal to a series of sines without any constant quantity, and $\\frac{d \\varphi}{d \\lambda}$, $\\frac{d \\nu}{d \\lambda}$, and $\\frac{d \\varepsilon}{d \\lambda}$ each equal to a series of cosines + a constant quantity.\n\nIn the second approximation the values of $\\lambda$, $r$ and $s$ retain the same form; and it is easy to see from the form of the expressions for $\\frac{d e}{d \\lambda}$, $\\frac{d h}{d \\lambda}$, &c. p. 334, that the form of the values of these quantities is not altered however far the approximation be carried.\nIf the sun or primary be a spheroid, \\( \\omega \\) the angle which the plane of the sun's equator makes with the plane of the orbit of the planet; and if the longitude be reckoned from the line of intersection of the sun's equator with the orbit of the planet; \\( R \\) is increased by the quantity \\( c \\left\\{ \\frac{3 \\sin^2 \\omega \\sin^2 \\lambda}{r^3} - 1 \\right\\} \\), \\( c \\) being a constant dependent upon the figure of the sun; but the partial differential coefficients of this quantity, which are introduced into the values of \\( d e, d \\varpi, \\&c. \\) do not change the form of the expressions for those quantities.\n\nIf the planet move in a medium which resists according to any power \\( n \\) of the velocity, if \\( c \\) be a constant and \\( v \\) the velocity, the term \\( 2c \\int v^{n+1} dt \\) must be added to \\( 2 \\int dR \\),\n\n\\[ cv^{n-1} \\left\\{ (1 + s^2) \\frac{dr'}{dt} + r' s \\frac{ds}{dt} \\right\\} \\text{ to } \\left( \\frac{dR}{dr'} \\right), \\]\n\n\\[ cv^{n-1} r'^2 \\frac{d\\lambda'}{dt} \\text{ to } \\left( \\frac{dR}{d\\lambda'} \\right), \\quad \\text{and } cv^{n-1} r' \\frac{d.r's}{dt} \\text{ to } \\left( \\frac{dR}{ds} \\right) \\]\n\nin the equations of p. 330.\n\nIf the orbit of the planet be supposed to coincide with the plane \\( xy \\), so that \\( s = 0 \\), then by the equations of p. 337 after reductions\n\n\\[ da = -2c a^2 \\left( \\frac{\\mu}{a} \\right)^{\\frac{n+1}{2}} \\left\\{ \\frac{1 + e \\cos v}{1 - e \\cos v} \\right\\}^{\\frac{n+1}{2}} \\frac{(1 - e \\cos v)}{n} dv \\]\n\n\\[ de = -2c \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ \\frac{1 + e \\cos v}{1 - e \\cos v} \\right\\}^{\\frac{n-1}{2}} \\left( \\frac{1 - e^2}{n} \\right) \\cos v dv \\]\n\n\\[ e d \\varpi = -2c \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ \\frac{1 + e \\cos v}{1 - e \\cos v} \\right\\}^{\\frac{n-1}{2}} \\sqrt{1 - e^2} \\sin v dv \\]\n\n\\[ d\\varepsilon - d\\varpi = 2c \\left( \\frac{\\mu}{a} \\right)^{\\frac{n-1}{2}} \\left\\{ \\frac{1 + e \\cos v}{1 - e \\cos v} \\right\\}^{\\frac{n-1}{2}} \\frac{\\sin v}{n} \\left\\{ \\frac{1 - e^3 \\cos v}{e} \\right\\} dv \\]\n\nThe form of these equations differs from that which obtained before, now the variations of \\( e, \\varpi \\) and \\( \\varepsilon \\) are periodical, while that of \\( a \\) has a term which varies with the time. \\( \\frac{de}{dv} \\) contains only odd powers of \\( \\cos v \\) and for that reason has no constant term. The periods of the periodic inequalities of all the elliptic constants due to the action of the resisting medium are fractional parts of the periodic time of the planet.\nIf the origin of \\( t \\) coincides with the instant of the perihelion passage, by Lagrange's theorem\n\n\\[\n\\cos v = \\cos n t - e \\sin n t^2 - \\frac{e^2}{2} \\cdot \\frac{d \\sin n t^3}{d.n t} - \\frac{e^3}{2 \\cdot 3} \\cdot \\frac{d^2 \\sin n t^4}{(d.n t)^2} - \\frac{e^4}{2 \\cdot 3 \\cdot 4} \\cdot \\frac{d^3 \\sin n t^5}{(d.n t)^3} - &c.\n\\]\n\n\\[\n\\sin v = \\sin n t + e \\sin n t \\cos n t + \\frac{e^2}{2} \\cdot \\frac{d \\sin n t^3 \\cos n t}{d.n t} + \\frac{e^3}{2 \\cdot 3} \\cdot \\frac{d^2 \\sin n t^4 \\cos n t}{(d.n t)^2} + \\frac{e^4}{2 \\cdot 3 \\cdot 4} \\cdot \\frac{d^3 \\sin n t^5 \\cos n t}{(d.n t)^3}\n\\]\n\n\\[\n\\sin n t^2 = \\frac{1 - \\cos 2 n t}{2}, \\quad \\sin n t^3 = \\frac{3 \\sin n t - \\sin 3 n t}{4}\n\\]\n\n\\[\n\\sin n t^4 = \\frac{3 - 4 \\cos 2 n t + \\cos 4 n t}{8}, \\quad \\sin n t^5 = \\frac{10 \\sin n t - 5 \\sin 3 n t + \\sin 5 n t}{16}\n\\]\n\n\\[\n\\frac{d \\sin n t^3}{d.n t} = \\frac{3 \\cos n t - 3 \\cos 3 n t}{4}\n\\]\n\n\\[\n\\frac{d \\sin n t^4}{d.n t} = \\frac{2 \\sin 2 n t - \\sin 4 n t}{2}, \\quad \\frac{d^2 \\sin n t^4}{(d.n t)^2} = 2 \\cos 2 n t - 2 \\cos 4 n t\n\\]\n\n\\[\n\\frac{d \\sin n t^5}{d.n t} = \\frac{10 \\cos n t - 15 \\cos 3 n t + 5 \\cos 5 n t}{16}\n\\]\n\n\\[\n\\frac{d^2 \\sin n t^5}{(d.n t)^2} = \\frac{-10 \\sin n t + 45 \\sin 3 n t - 25 \\sin 5 n t}{16}\n\\]\n\n\\[\n\\frac{d^3 \\sin n t^5}{(d.n t)^3} = \\frac{-10 \\cos n t + 135 \\cos 3 n t - 125 \\cos 5 n t}{16}\n\\]\n\n\\[\n\\cos v = \\cos n t - \\frac{e}{2} + \\frac{e}{2} \\cos 2 n t - \\frac{e^2}{2} \\left\\{ \\frac{3 \\cos n t - 3 \\cos 3 n t}{4} \\right\\} - \\frac{e^3}{2 \\cdot 3} \\left\\{ 2 \\cos 2 n t - 2 \\cos 4 n t \\right\\}\n\\]\n\n\\[\n- \\frac{e^4}{2 \\cdot 3 \\cdot 4} \\left\\{ \\frac{-10 \\cos n t + 135 \\cos 3 n t - 125 \\cos 5 n t}{16} \\right\\}\n\\]\n\nIf the origin of the time does not coincide with the perihelion passage, \\( n t + \\epsilon - \\varpi \\) must be substituted for \\( n t \\), but as \\( \\epsilon \\) always accompanies \\( n t \\), it may be suppressed at present for convenience, and afterwards replaced.\n\n\\[\n\\cos v = \\left\\{ 1 - \\frac{3}{8} e^2 + \\frac{5 e^4}{192} \\right\\} \\cos (n t - \\varpi) - \\frac{e}{2} + \\frac{e}{2} \\left\\{ 1 - \\frac{2 e^2}{3} \\right\\} \\cos (2nt - 2 \\varpi)\n\\]\n\n\\[\n+ \\frac{3}{8} e^2 \\left\\{ 1 - \\frac{15}{16} e^2 \\right\\} \\cos (3nt - 3\\varpi) + \\frac{e^3}{3} \\cos (4nt - 4\\varpi) + \\frac{125}{384} e^4 \\cos (5nt - 5\\varpi)\n\\]\n\n\\[\n\\sin n t^2 \\cos n t = \\frac{d \\sin n t^3}{3 d.n t} = \\frac{\\cos n t - \\cos 3 n t}{4}, \\quad \\frac{d \\sin n t^3 \\cos n t}{d.n t} = \\frac{-\\sin n t + 3 \\sin 3 n t}{4}\n\\]\n\\[\n\\sin n t^3 \\cos n t = \\frac{d \\cdot \\sin n t^4}{4 d \\cdot n t} = \\frac{2 \\sin 2 n t - \\sin 4 n t}{8}, \\quad \\frac{d \\cdot \\sin n t^3 \\cos n t}{d \\cdot n t} = \\frac{\\cos 2 n t - \\cos 4 n t}{2}\n\\]\n\n\\[\n\\frac{d^2 \\cdot \\sin n t^3 \\cos n t}{(d \\cdot n t)^2} = -\\sin 2 n t + 2 \\sin 4 n t\n\\]\n\n\\[\n\\sin n t^4 \\cos n t = \\frac{d \\cdot \\sin n t^5}{5 d \\cdot n t} = \\frac{2 \\cos n t - 3 \\cos 3 n t + \\cos 5 n t}{16}\n\\]\n\n\\[\n\\frac{d \\cdot \\sin n t^4 \\cos n t}{d \\cdot n t} = -2 \\sin n t + 9 \\sin 3 n t - 5 \\sin 5 n t\n\\]\n\n\\[\n\\frac{d^2 \\cdot \\sin n t^4 \\cos n t}{(d \\cdot n t)^2} = -2 \\cos n t + 27 \\cos 3 n t - 25 \\cos 5 n t\n\\]\n\n\\[\n\\frac{d^3 \\cdot \\sin n t^4 \\cos n t}{(d \\cdot n t)^3} = \\frac{2 \\sin n t - 81 \\sin 3 n t + 125 \\sin 5 n t}{16}\n\\]\n\n\\[\n\\sin v = \\sin n t + \\frac{e^2}{2} \\sin 2 n t + \\frac{e^3}{2} \\left\\{ \\frac{-\\sin n t + 3 \\sin 3 n t}{4} \\right\\} + \\frac{e^3}{2 \\cdot 3} \\left\\{ -\\sin 2 n t + 2 \\sin 4 n t \\right\\}\n\\]\n\n\\[\n+ \\frac{e^4}{2 \\cdot 3 \\cdot 4} \\left\\{ \\frac{2 \\sin n t - 81 \\sin 3 n t + 125 \\sin 5 n t}{16} \\right\\}\n\\]\n\nand replacing \\( n t \\) by \\( n t - \\varpi \\)\n\n\\[\n\\sin v = \\left\\{ 1 - \\frac{3 e^2}{8} + \\frac{5}{192} e^4 \\right\\} \\sin (n t - \\varpi) + \\frac{e}{2} \\left\\{ 1 - 2 \\frac{e^2}{3} \\right\\} \\sin (2 n t - 2 \\varpi)\n\\]\n\n\\[\n+ \\left\\{ \\frac{2 e^3}{8} - \\frac{4}{192} e^4 \\right\\} \\sin (n t - \\varpi) + \\frac{e^3}{6} \\sin (2 n t - 2 \\varpi)\n\\]\n\n\\[\n+ \\frac{3 e^2}{8} \\left\\{ 1 - \\frac{15}{16} e^2 \\right\\} \\sin (3 n t - 3 \\varpi) + \\frac{e^3}{3} \\sin (4 n t - 4 \\varpi) + \\frac{125}{384} e^4 \\sin (5 n t - 5 \\varpi)\n\\]\n\n\\[\n+ \\frac{3}{8} \\cdot \\frac{6}{16} e^4 \\sin (3 n t - 3 \\varpi)\n\\]\n\n\\[\n\\cos (v + \\varpi) = \\left\\{ 1 - \\frac{e^2}{4} + \\frac{e^4}{64} \\right\\} \\cos n t - \\frac{e}{2} \\cos \\varpi + \\frac{e}{2} \\left\\{ 1 - \\frac{e^2}{2} \\right\\} \\cos (2 n t - \\varpi)\n\\]\n\n\\[\n+ \\frac{3}{8} e^2 \\left\\{ 1 - \\frac{3}{4} e^2 \\right\\} \\cos (3 n t - 2 \\varpi) + \\frac{e^3}{3} \\cos (4 n t - 3 \\varpi) + \\frac{125}{384} e^4 \\cos (5 n t - 4 \\varpi)\n\\]\n\n\\[\n- \\frac{e^3}{8} \\left\\{ 1 - \\frac{e^2}{12} \\right\\} \\cos (n t - 2 \\varpi) - \\frac{e^3}{12} \\cos (2 n t - 3 \\varpi) - \\frac{9}{128} e^4 \\cos (3 n t - 4 \\varpi)\n\\]\n\n\\[\n\\sin (v + \\varpi) = \\left\\{ 1 - \\frac{e^2}{4} + \\frac{e^4}{64} \\right\\} \\sin n t - \\frac{e}{2} \\sin \\varpi + \\frac{e}{2} \\left\\{ 1 - \\frac{e^2}{2} \\right\\} \\sin (2 n t - \\varpi)\n\\]\n\n\\[\n+ \\frac{3}{8} e^2 \\left\\{ 1 - \\frac{3}{4} e^2 \\right\\} \\sin (3 n t - 2 \\varpi) + \\frac{e^3}{3} \\sin (4 n t - 3 \\varpi) + \\frac{125}{384} e^4 \\sin (5 n t - 4 \\varpi)\n\\]\n\\[\n\\cos(\\nu - \\varpi) = \\left(1 - \\frac{e^2}{4}\\right)\\cos(nt - 2\\varpi) - \\frac{e}{2}\\cos\\varpi + \\frac{e}{2}\\cos(2nt - 3\\varpi) + \\frac{3}{8}e^2\\cos(3nt - 4\\varpi) - \\frac{e^2}{8}\\cos nt\n\\]\n\n\\[\n\\sin(\\nu - \\varpi) = \\left(1 - \\frac{e^2}{4}\\right)\\sin(nt - 2\\varpi) + \\frac{e}{2}\\sin\\varpi + \\frac{e}{2}\\sin(2nt - 3\\varpi) + \\frac{3}{8}e^2\\sin(3nt - 4\\varpi) + \\frac{e^2}{8}\\sin nt\n\\]\n\n\\[\nr\\cos\\lambda = r\\cos(\\lambda - \\varpi)\\cos\\varpi - r\\sin(\\lambda - \\varpi)\\sin\\varpi = a\\left\\{\\cos\\nu - e\\cos\\varpi - (1 - e^2)\\frac{1}{2}\\sin\\nu\\sin\\varpi\\right\\}\n\\]\n\n\\[\nr\\sin\\lambda = r\\sin(\\lambda - \\varpi)\\cos\\varpi + r\\cos(\\lambda - \\varpi)\\sin\\varpi = a\\left\\{(1 - e^2)\\frac{1}{2}\\sin\\nu\\cos\\varpi + (\\cos\\nu - e)\\sin\\varpi\\right\\}\n\\]\n\n\\[\nr\\cos\\lambda = a\\left\\{\\left(1 - \\frac{e^2}{4} - \\frac{e^4}{16}\\right)\\cos(\\nu + \\varpi) - e\\cos\\varpi + \\frac{e^2}{4}\\left(1 + \\frac{e^2}{4}\\right)\\cos(\\nu - \\varpi)\\right\\}\n\\]\n\n\\[\nr\\sin\\lambda = a\\left\\{\\left(1 - \\frac{e^2}{4} - \\frac{e^4}{16}\\right)\\sin(\\nu + \\varpi) - e\\sin\\varpi - \\frac{e^2}{4}\\left(1 + \\frac{e^2}{4}\\right)\\sin(\\nu - \\varpi)\\right\\}\n\\]\n\n\\[\nr\\cos\\lambda = a\\left\\{\\left(1 - \\frac{e^2}{2} - \\frac{e^4}{64}\\right)\\cos nt - \\frac{3}{2}e\\cos\\varpi + \\frac{e}{2}\\left(1 - \\frac{3}{4}e^2\\right)\\cos(2nt - \\varpi)\\right.\n\\]\n\n\\[\n+ \\frac{3}{8}e^2\\left(1 - e^2\\right)\\cos(3nt - 2\\varpi) + \\frac{e^3}{3}\\cos(4nt - 3\\varpi) + \\frac{125}{384}e^4\\cos(5nt - 4\\varpi)\n\\]\n\n\\[\n+ \\frac{e^2}{8}\\left(1 + \\frac{e^2}{3}\\right)\\cos(nt - 2\\varpi) + \\frac{e^3}{24}\\cos(2nt - 3\\varpi) + \\frac{3}{128}e^4\\cos(3nt - 4\\varpi) &c.\n\\]\n\n\\[\nrr_i\\cos\\{\\lambda - \\lambda_i\\} = aa_i\\left\\{\\left(1 - \\frac{e^2}{2} - \\frac{e^4}{64}\\right)\\left(1 - \\frac{e_i^2}{2} - \\frac{e_i^4}{64}\\right)\\cos(nt - n_i t)\\right.\n\\]\n\n\\[\n- \\frac{3}{2}e\\left(1 - \\frac{e^2}{2}\\right)\\cos(n_i t - \\varpi) + \\frac{e}{2}\\left(1 - \\frac{3}{4}e^2\\right)\\left(1 - \\frac{e^2}{2}\\right)\\cos(2nt - n_i t - \\varpi)\n\\]\n\n\\[\n+ \\frac{3}{8}e^2\\left(1 - e^2\\right)\\left(1 - \\frac{e_i^2}{2}\\right)\\cos(3nt - n_i t - 2\\varpi) + \\frac{e^3}{3}\\cos(4nt - n_i t - 3\\varpi)\n\\]\n\n\\[\n+ \\frac{125}{384}e^4\\cos(5nt - n_i t - 4\\varpi) + \\frac{e^3}{8}\\left(1 + \\frac{e^2}{3}\\right)\\left(1 - \\frac{e_i^2}{2}\\right)\\cos(nt + n_i t - 2\\varpi)\n\\]\n\n\\[\n+ \\frac{e^3}{24}\\cos(2nt + n_i t - 3\\varpi) + \\frac{3}{128}e^4\\cos(3nt + n_i t - 4\\varpi)\n\\]\n\n\\[\n- \\frac{3}{2}e_i\\left(1 - \\frac{e^2}{2}\\right)\\cos(nt - \\varpi_i) + \\frac{9}{4}ee_i\\cos(\\varpi - \\varpi_i)\n\\]\n\n\\[\n- \\frac{3}{4}e_e_i\\left(1 - \\frac{3}{4}e^2\\right)\\cos(2nt - \\varpi - \\varpi_i) - \\frac{9}{16}e^2e_i\\cos(3nt - 2\\varpi - \\varpi_i)\n\\]\n\n\\[\n- \\frac{e^2e_i}{2}\\cos(4nt - 3\\varpi - \\varpi_i) + \\frac{3}{16}e^2e_i\\cos(nt - 2\\varpi + \\varpi_i) - \\frac{e^3e_i}{16}\\cos(2nt - 3\\varpi + \\varpi_i)\n\\]\n\\[\n\\begin{align*}\n&+ \\frac{e_i}{2} \\left(1 - \\frac{3}{4} e_i^2\\right) \\left(1 - \\frac{e_i^2}{2}\\right) \\cos(nt - 2n_i t + \\varpi_i) - \\frac{3}{4} ee_i \\left(1 - \\frac{3}{4} e_i^2\\right) \\cos(2n_i t - \\varpi - \\varpi_i) \\\\\n&+ \\frac{ee_i}{4} \\left(1 - \\frac{3}{4} e_i^2\\right) \\left(1 - \\frac{3}{4} e_i^2\\right) \\cos(2nt - 2n_i t - \\varpi + \\varpi_i) + \\frac{3}{16} e^2 e_i \\cos(3nt - 2n_i t - 2\\varpi + \\varpi_i) \\\\\n&+ \\frac{e^3 e_i}{6} \\cos(4nt - 2n_i t - 3\\varpi + \\varpi_i) + \\frac{e^2 e_i}{16} \\cos(nt + 2n_i t - 2\\varpi - \\varpi_i) \\\\\n&- \\frac{e^3 e_i}{48} \\sin(2nt + 2n_i t - 3\\varpi - \\varpi_i) + \\frac{3}{8} e_i^2 \\left(1 - e_i^2\\right) \\left(1 - \\frac{e_i^2}{2}\\right) \\cos(nt - 3n_i t + 2\\varpi_i) \\\\\n&- \\frac{9}{16} ee_i^2 \\cos(3n_i t - \\varpi - 2\\varpi_i) + \\frac{3}{16} ee_i^2 \\cos(2nt - 3n_i t - \\varpi + 2\\varpi_i) \\\\\n&+ \\frac{9}{64} e^2 e_i^2 \\cos(3nt - 3n_i t - 2\\varpi + 2\\varpi_i) + \\frac{3}{64} e^2 e_i^2 \\cos(nt + 3n_i t - 2\\varpi - 2\\varpi_i) \\\\\n&+ \\frac{e^3}{3} \\cos(nt - 4n_i t + 3\\varpi_i) + \\frac{ee_i^3}{2} \\cos(4n_i t - \\varpi - 3\\varpi_i) \\\\\n&+ \\frac{ee_i^3}{6} \\cos(2nt - 4n_i t - \\varpi + 3\\varpi_i) + \\frac{125}{384} e_i^4 \\cos(nt - 5n_i t + 4\\varpi_i) \\\\\n&+ \\frac{e^2}{8} \\left(1 + \\frac{e^2}{3} e_i^2\\right) \\left(1 - \\frac{e_i^2}{2}\\right) \\cos(nt + n_i t - 2\\varpi_i) - \\frac{3}{16} ee_i^2 \\cos(n_i t + \\varpi - 2\\varpi_i) \\\\\n&+ \\frac{ee_i^3}{16} \\cos(2nt + n_i t - \\varpi - 2\\varpi_i) + \\frac{3}{64} e^2 e_i^2 \\cos(3nt + n_i t - 2\\varpi - 2\\varpi_i) \\\\\n&+ \\frac{e^3 e_i^2}{64} \\cos(nt - n_i t - 2\\varpi + 2\\varpi_i) + \\frac{e^3}{24} \\cos(nt + 2n_i t - 3\\varpi_i) \\\\\n&- \\frac{ee_i^3}{16} \\cos(2n_i t + \\varpi - 3\\varpi_i) + \\frac{ee_i^3}{48} \\cos(2nt + 2n_i t - \\varpi - 3\\varpi_i) \\\\\n&+ \\frac{1}{128} e_i^4 \\cos(nt + 3n_i t - 4\\varpi_i) + \\text{&c.}\n\\end{align*}\n\\]\n\n\\[rr_i \\cos\\{\\lambda + \\lambda_i\\} = aa_i \\left\\{ \\left(1 - \\frac{e^2 + e_i^2}{2}\\right) \\cos(nt + n_i t)\\right.\\]\n\\[- \\frac{3}{2} e \\cos(n_i t + \\varpi) + \\frac{e}{2} \\cos(2nt + n_i t - \\varpi)\\]\n\\[- \\frac{3}{8} e^2 \\cos(3nt + n_i t - 2\\varpi) + \\frac{e^2}{8} \\cos(nt - n_i t - 2\\varpi)\\]\n\\[- \\frac{3}{2} e_i \\cos(nt + \\varpi) + \\frac{9}{4} ee_i \\cos(\\varpi + \\varpi_i) - \\frac{3}{4} ee_i \\cos(2nt - \\varpi + \\varpi_i)\\]\n\\[ r = a \\left\\{ 1 + \\frac{e^2}{2} - e \\left( 1 - \\frac{3}{8} e^2 \\right) \\cos(nt - \\varpi) - \\frac{e^2}{2} \\left( 1 - \\frac{2e^2}{3} \\right) \\cos(2nt - 2\\varpi) \\right. \\\\\n- \\frac{3}{8} e^3 \\cos(3nt - 3\\varpi) - \\frac{e^4}{3} \\cos(4nt - 4\\varpi) + &c. \\]\n\n\\[ r^2 = a^2 \\left\\{ 1 + \\frac{3e^2}{2} - 2e \\left( 1 - \\frac{e^2}{8} \\right) \\cos(nt - \\varpi) - \\frac{e^2}{2} \\left( 1 - \\frac{e^2}{3} \\right) \\cos(2nt - 2\\varpi) \\right. \\\\\n- \\frac{e^2}{4} \\cos(3nt - 3\\varpi) - \\frac{e^4}{6} \\cos(4nt - 4\\varpi) + &c. \\]\n\n\\[ r^{-3} = a^{-3} \\left\\{ 1 + \\frac{3}{2} e^2 \\left( 1 + \\frac{e^2}{2} \\right) + 3e \\left( 1 + \\frac{9}{8} e^2 \\right) \\cos(nt - \\varpi) \\right. \\\\\n+ \\frac{9}{2} e^2 \\left( 1 - \\frac{e^2}{18} \\right) \\cos(2nt - 2\\varpi) + \\frac{53}{8} e^3 \\cos(3nt - 3\\varpi) \\\\\n+ \\frac{31}{4} e^4 \\cos(4nt - 4\\varpi) + &c. \\]\n\nIf the coefficient of \\( \\cos n \\theta \\) in the development of \\( \\left\\{ 1 - \\frac{2a}{a_i} \\cos \\theta + \\frac{a^2}{a_i^2} \\right\\}^{-m} \\)\naccording to cosines of arcs multiples of \\( \\theta \\) is called \\( b_{2m,n} \\)\n\n\\[ b_{1,0} = 1 + \\left( \\frac{1}{2} \\right)^2 \\frac{a^2}{a_i^2} + \\left( \\frac{1 \\cdot 3}{2 \\cdot 4} \\right)^2 \\frac{a^4}{a_i^4} + \\left( \\frac{1 \\cdot 3 \\cdot 5}{2 \\cdot 4 \\cdot 6} \\right)^2 \\frac{a^6}{a_i^6} + &c. \\]\n\\[ b_{1,1} = \\frac{a}{a_i} + \\frac{1.3}{2.4} \\frac{a^3}{a_i^3} + \\frac{1.3.3.5}{2.4.4.6} \\frac{a^5}{a_i^5} + &c. \\]\n\n\\[ b_{1,2} = \\frac{3}{4} \\frac{a^2}{a_i^2} + \\frac{1.3.5}{2.4.6} \\frac{a^4}{a_i^4} + \\frac{1.3.3.5.7}{2.4.4.6.8} \\frac{a^6}{a_i^6} + &c. \\]\n\n\\[ b_{1,3} = \\frac{3.5}{4.6} \\frac{a^3}{a_i^3} + \\frac{1.3.5.7}{2.4.6.8} \\frac{a^5}{a_i^5} + \\frac{1.3.3.5.7.9}{2.4.4.6.8.10} \\frac{a^7}{a_i^7} + &c. \\]\n\n\\[ b_{3,0} = 1 + \\left( \\frac{3}{2} \\right)^2 \\frac{a^2}{a_i^2} + \\left( \\frac{3.5}{2.4} \\right)^2 \\frac{a^4}{a_i^4} + \\left( \\frac{3.5.7}{2.4.6} \\right)^2 \\frac{a^6}{a_i^6} + &c. \\]\n\n\\[ b_{3,1} = \\frac{3}{a_i} + \\frac{3.3.5}{2.4} \\frac{a^3}{a_i^3} + \\frac{3.5.3.5.7}{2.4.4.6} \\frac{a^5}{a_i^5} + &c. \\]\n\n\\[ b_{3,2} = \\frac{3.5}{4} \\frac{a^2}{a_i^2} + \\frac{3.3.5.7}{2.4.6} \\frac{a^4}{a_i^4} + \\frac{3.5.3.5.7}{2.4.4.6.8} \\frac{a^6}{a_i^6} + &c. \\]\n\n\\[ b_{3,3} = \\frac{3.5.7}{4.6} \\frac{a^3}{a_i^3} + \\frac{3.3.5.7.9}{2.4.6.8} \\frac{a^5}{a_i^5} + \\frac{3.5.3.5.7.9.11}{2.4.4.6.8.10} \\frac{a^7}{a_i^7} + &c. \\]\n\n\\[ b_{5,0} = 1 + \\left( \\frac{5}{2} \\right)^2 \\frac{a^2}{a_i^2} + \\left( \\frac{5.7}{2.4} \\right)^2 \\frac{a^4}{a_i^4} + \\left( \\frac{5.7.9}{2.4.6} \\right)^2 \\frac{a^6}{a_i^6} + &c. \\]\n\n\\[ b_{5,1} = 5 \\frac{a}{a_i} + \\frac{5.5.7}{2.4} \\frac{a^3}{a_i^3} + \\frac{5.7.5.7.9}{2.4.4.6} \\frac{a^5}{a_i^5} + &c. \\]\n\n\\[ b_{5,2} = \\frac{5.7}{4} \\frac{a^2}{a_i^2} + \\frac{5.5.7.9}{2.4.6} \\frac{a^4}{a_i^4} + \\frac{5.7.5.7.9.11}{2.4.4.6.8} \\frac{a^6}{a_i^6} + &c. \\]\n\n\\[ b_{5,3} = \\frac{5.7.9}{4.6} \\frac{a^3}{a_i^3} + \\frac{5.5.7.9.11}{2.4.6.8} \\frac{a^5}{a_i^5} + \\frac{5.7.5.7.9.11.13}{2.4.4.6.8.10} \\frac{a^7}{a_i^7} \\]\n\n\\[ b_{3,1} = 3 \\frac{a}{a_i} \\left\\{ b_{5,0} - \\frac{1}{2} b_{5,2} \\right\\} \\quad b_{1,1} = \\frac{a}{a_i} \\left\\{ b_{3,0} - \\frac{1}{2} b_{3,2} \\right\\} \\]\n\n\\[ 2 b_{3,2} = \\frac{3}{2} \\frac{a}{a_i} \\left\\{ b_{5,1} - b_{5,3} \\right\\} \\quad 2 b_{1,2} = \\frac{1}{2} \\frac{a}{a_i} \\left\\{ b_{3,1} - b_{3,3} \\right\\} \\]\n\n\\[ 3 b_{3,3} = \\frac{3}{2} \\frac{a}{a_i} \\left\\{ b_{5,2} - b_{5,4} \\right\\} \\quad 3 b_{1,3} = \\frac{1}{2} \\frac{a}{a_i} \\left\\{ b_{3,2} - b_{3,4} \\right\\} \\]\n\n\\[ R = m_r \\left\\{ \\frac{r r'_i \\cos (\\lambda - \\lambda')}{r_i^3} - \\frac{1}{\\{ r^3 - 2 r r'_i \\cos (\\lambda - \\lambda') + r_i^3 \\}^{1/2}} \\right\\} \\]\n\nand since\n\n\\[ rr'_i \\cos (\\lambda - \\lambda') = rr'_i \\left\\{ \\cos^2 \\frac{\\lambda'}{2} \\cos (\\lambda - \\lambda') + \\sin^2 \\frac{\\lambda'}{2} \\cos (\\lambda + \\lambda' - 2 \\nu) \\right\\}, \\text{ p. 330; } \\]\n\\[ R = m_i \\left\\{ rr_i \\left[ \\cos^2 \\frac{\\lambda}{2} \\cos (\\lambda - \\lambda_i) + \\sin^2 \\frac{\\lambda}{2} \\cos (\\lambda + \\lambda_i - 2\\nu_i) \\right] \\right\\} \\]\n\n\\[ - \\left\\{ r^2 - 2rr_i \\left[ \\cos^2 \\frac{\\lambda}{2} \\cos (\\lambda - \\lambda_i) + \\sin^2 \\frac{\\lambda}{2} \\cos (\\lambda + \\lambda_i - 2\\nu_i) \\right] + r_i^2 \\right\\}^{\\frac{1}{2}} \\]\n\nNeglecting the cubes and higher powers of \\( e \\)\n\n\\[ R = m_i \\frac{a}{a_i^2} \\left\\{ \\left( 1 - \\sin^2 \\frac{\\lambda}{2} - \\frac{e^2 + e_i^2}{2} \\right) \\cos (nt - n_i t) - \\frac{3}{2} e \\cos (n_i t - \\varpi) - \\frac{3}{2} e_i \\cos (nt - \\varpi_i) \\right. \\]\n\n\\[ + \\frac{e}{2} \\cos (2nt - n_i t - \\varpi) + \\frac{e_i}{2} \\cos (nt - 2n_i t + \\varpi_i) + \\frac{3}{8} e^2 \\cos (3nt - n_i t - 2\\varpi) \\]\n\n\\[ + \\frac{3}{8} e_i^2 \\cos (nt - 3n_i t + 2\\varpi_i) + \\frac{e^2}{8} \\cos (nt + n_i t - 2\\varpi) + \\frac{e_i^2}{8} \\cos (nt + n_i t - 2\\varpi_i) \\]\n\n\\[ + \\frac{9}{4} ee_i \\cos (\\varpi - \\varpi_i) - \\frac{3}{4} ee_i \\cos (2nt - \\varpi - \\varpi_i) - \\frac{3}{4} ee_i \\cos (2n_i t - \\varpi - \\varpi_i) \\]\n\n\\[ + \\frac{ee_i}{4} \\cos (2nt - 2n_i t - \\varpi + \\varpi_i) + \\sin^2 \\frac{\\lambda}{2} \\cos (nt + n_i t - 2\\nu) \\right\\} \\left\\{ 1 + \\frac{3}{2} e^2, \\right. \\]\n\n\\[ + 3e \\cos (nt - \\varpi) + \\frac{9}{2} e^2 \\cos (2nt - 2\\varpi) + \\text{&c.} \\}\n\n\\[ - m_i \\left\\{ a^2 \\left[ 1 + \\frac{3}{2} e^2 - 2e \\cos (nt - \\varpi) - \\frac{e^2}{2} \\cos (2nt - 2\\varpi) \\right] \\right\\} \\]\n\n\\[ - 2aa_i \\left\\{ \\left( 1 - \\sin^2 \\frac{\\lambda}{2} - \\frac{e^2 + e_i^2}{2} \\right) \\cos (nt - n_i t) - \\frac{3}{2} e \\cos (n_i t - \\varpi) \\right. \\]\n\n\\[ - \\frac{3}{2} e_i \\cos (nt - \\varpi) + \\frac{e}{2} \\cos (2nt - n_i t - \\varpi) + \\frac{e_i}{2} \\cos (nt - 2n_i t + \\varpi_i) \\]\n\n\\[ + \\frac{3}{8} e^2 \\cos (3nt - n_i t - 2\\varpi) + \\frac{3}{8} e_i^2 \\cos (nt - 3n_i t + 2\\varpi_i) + \\frac{e^2}{8} \\cos (nt + n_i t - 2\\varpi) \\]\n\n\\[ + \\frac{e_i^2}{8} \\cos (nt + n_i t - 2\\varpi_i) + \\frac{9}{4} ee_i \\cos (\\varpi - \\varpi_i) - \\frac{3}{4} ee_i \\cos (2nt - \\varpi - \\varpi_i) \\]\n\n\\[ - \\frac{3}{4} ee_i \\cos (2n_i t - \\varpi - \\varpi_i) + \\frac{ee_i}{4} \\cos (2nt - 2n_i t - \\varpi + \\varpi_i) \\]\n\n\\[ + \\sin^2 \\frac{\\lambda}{2} \\cos (nt + n_i t - 2\\nu) \\right\\} + a_i^2 \\left\\{ 1 + \\frac{3}{2} e_i^2 \\right. \\]\n\n\\[ - 2e_i \\cos (n_i t - \\varpi) - \\frac{e_i^2}{2} \\cos (2n_i t - 2\\varpi_i) + \\text{&c.} \\} \\right\\}^{\\frac{1}{2}} \\]\nterms independent of \\( b \\)\n\n\\[\n-\\frac{m_i}{a_i} \\left\\{ b_{1,0} + b_{1,1} \\cos(nt - n_t) + b_{1,2} \\cos(2nt - 2n_t) + \\text{etc.} \\right\\}\n\\]\n\n\\[\n+\\frac{m_i}{2a_i^3} \\left\\{ b_{3,0} + b_{3,1} \\cos(nt - n_t) + b_{3,2} \\cos(2nt - 2n_t) \\text{etc.} \\right\\} \\left\\{ a^2 \\left\\{ \\frac{3}{2} e^2 \\right. \\right.\n\\]\n\n\\[\n-2e \\cos(nt - \\varphi) - \\frac{e^2}{2} \\cos(2nt - 2\\varphi) \\right\\} - 2aa_i \\left\\{ \\left(-\\sin^2 \\frac{n_t}{2} - \\frac{e^2 + e_i^2}{2}\\right) \\cos(nt - n_t)\n\\]\n\n\\[\n-\\frac{3}{2} e \\cos(n_t - \\varphi) - \\frac{3}{2} e_i \\cos(nt - \\varphi) + \\frac{e}{2} \\cos(2nt - n_t - \\varphi)\n\\]\n\n\\[\n+\\frac{e_i}{2} \\cos(nt - 2n_t + \\varphi_i) + \\frac{3}{8} e^2 \\cos(3nt - n_t - 2\\varphi) + \\frac{3}{8} e_i^2 \\cos(nt - 3n_t + 2\\varphi_i)\n\\]\n\n\\[\n+\\frac{e^2}{8} \\cos(nt + n_t - 2\\varphi) + \\frac{e_i^2}{8} \\cos(nt + n_t - 2\\varphi_i) + \\frac{9}{4} ee_i \\cos(\\varphi - \\varphi_i)\n\\]\n\n\\[\n-\\frac{3}{4} ee_i \\cos(2nt - \\varphi - \\varphi_i) - \\frac{3}{4} ee_i \\cos(2nt - \\varphi - \\varphi_i) + \\frac{ee_i}{4} \\cos(2nt - 2n_t - \\varphi + \\varphi_i)\n\\]\n\n\\[\n+\\sin^2 \\frac{n_t}{2} \\cos(nt + n_t - 2\\varphi_i) \\right\\} + a_i^2 \\left\\{ \\frac{3}{2} e_i^2 - 2e_i \\cos(n_t - \\varphi) - \\frac{e_i^2}{2} \\cos(2n_t - 2\\varphi_i) \\text{etc.} \\right\\}\n\\]\n\n\\[\n-\\frac{1}{2} \\cdot \\frac{3m_i}{4a_i^5} \\left\\{ b_{5,0} + b_{5,1} \\cos(nt - n_t) + b_{5,2} \\cos(2nt - 2n_t) + \\text{etc.} \\right\\} \\left\\{ 2a^2 e \\cos(nt - \\varphi)\n\\]\n\n\\[\n-3aa_i e \\cos(nt - \\varphi) - 3aa_i e \\cos(nt - \\varphi_i) + aa_i e \\cos(2nt - n_t - \\varphi)\n\\]\n\n\\[\n+a_a e_i \\cos(nt - 2n_t + \\varphi_i) + 2a_i^2 e_i \\cos(n_t - \\varphi_i) \\right\\}^2 + \\text{etc.}\n\\]\n\n\\[\n\\left\\{ 2a^2 e \\cos(nt - \\varphi) - 3aa_i e \\cos(nt - \\varphi) - 3aa_i e \\cos(nt - \\varphi_i)\n\\]\n\n\\[\n+aa_i e \\cos(2nt - n_t - \\varphi) + aa_i e \\cos(nt - 2n_t + \\varphi_i) + 2a_i^2 e_i \\cos(n_t - \\varphi_i) \\right\\}^2\n\\]\n\n\\[\n= 2a^4 e^2 + 5a^2 a_i^2 (e^2 + e_i^2) + 2a_i^4 e_i^2 + a^2(2a^2 - 3a_i^2)e^2 \\cos(2nt - 2\\varphi)\n\\]\n\n\\[\n+\\frac{9}{2} a^2 a_i^2 e^2 \\cos(2nt - 2\\varphi) + \\frac{9}{2} a^2 a_i^2 e_i^2 \\cos(2nt - 2\\varphi_i) + \\frac{a^2 a_i^2}{2} e^2 \\cos(4nt - 2n_t - 2\\varphi)\n\\]\n\n\\[\n+\\frac{a^2 a_i^2}{2} e_i^2 \\cos(2nt - 4n_t + 2\\varphi) + a_i^2(2a^2 - 3a_i^2)e_i^2 \\cos(2nt - 2\\varphi_i)\n\\]\n\n\\[\n-4(a^2 e^2 + a_i^2 e_i^2)aa_i \\cos(nt - n_t) - 6a^3 a_i e^2 \\cos(nt + n_t - 2\\varphi)\n\\]\n\\[-6aa^3e^2\\cos(nt+n_t-2\\varpi) - 6(a^2+a_i^2)aa_iee_i\\cos(\\varpi-\\varpi_i)\\]\n\n\\[-2(3a^2-a_i^2)aa_iee_i\\cos(2nt-\\varpi-\\varpi_i) - 2(3a_i^2-a^2)aa_iee_i\\cos(2n_t-\\varpi-\\varpi_i)\\]\n\n\\[+2a^3a_ie^2\\cos(3nt-n_t-2\\varpi) + 2a^3a_ie^2\\cos(nt-3n_t+2\\varpi)\\]\n\n\\[+2(a^2+a_i^2)aa_iee_i\\cos(2nt-2n_t-\\varpi+\\varpi_i) + 9a^2a_i^2ee_i\\cos(nt-n_t+\\varpi-\\varpi_i)\\]\n\n\\[+14a^2a_i^2ee_i\\cos(nt+n_t-\\varpi-\\varpi_i) - 3a^2a_i^2(e^2+e_i^2)\\cos(2nt-2n_t)\\]\n\n\\[-3a^2a_i^2ee_i\\cos(nt-3n_t+\\varpi+\\varpi_i) - 3a^2a_i^2ee_i\\cos(3nt-n_t-\\varpi-\\varpi_i)\\]\n\n\\[-2a^2a_i^2ee_i\\cos(nt-n_t-\\varpi+\\varpi_i) + a^2a_i^2ee_i\\cos(3nt-3n_t-\\varpi+\\varpi_i)\\]\n\nReplacing \\(e\\) which accompanies \\(nt\\),\n\n\\[R = m_i\\left\\{-\\frac{b_{1,0}}{a_i} + \\frac{3(a^2e^2+a_i^2e_i^2)}{2.2a_i^3}b_{3,0} + \\frac{aa_i}{2a_i^3}\\left(\\sin^2\\frac{\\varpi_i}{2} + \\frac{e^2+e_i^2}{2}\\right)b_{3,1}\\right.\\]\n\n\\[-\\frac{1.3}{2.4a_i^5}\\left(2a^4e^2+5a^2a_i^2(e^2+e_i^2)+2a_i^4e_i^2\\right)b_{5,0} + \\frac{1.3.2}{2.4a_i^5}(a^2e^2+a_i^2e_i^2)aa_ib_{5,1}\\]\n\n\\[\\left.+\\frac{1.3.3a^2a_i^2}{2.4.2a_i^5}(e^2+e_i^2)b_{5,2}\\right\\}\\]\n\n\\[+m_i\\left\\{\\frac{a}{a_i^2}\\left\\{\\cos^2\\frac{\\varpi_i}{2} - \\frac{e^2+e_i^2}{2}\\right\\} - \\frac{b_{1,1}}{a_i} + \\frac{aa_i}{a_i^3}\\left\\{\\sin^2\\frac{\\varpi_i}{2} + \\frac{e^2+e_i^2}{2}\\right\\}\\left\\{b_{3,0} + \\frac{1}{2}b_{3,2}\\right\\}\\right.\\]\n\n\\[\\left.+\\frac{3(a^2e^2+a_i^2e_i^2)}{2.2a_i^3}b_{3,1} + \\frac{1.3.4}{2.4a_i^5}(a^2e^2+a_i^2e_i^2)aa_ib_{5,0}\\right.\\]\n\n\\[-\\frac{1.3}{2.4a_i^5}\\left\\{2a^4e^2+\\frac{7}{2}a^2a_i^2(e^2+e_i^2)+2a_i^4e_i^2\\right\\}b_{5,1} + \\frac{1.3.2}{2.4a_i^5}(a^2e^2+a_i^2e_i^2)aa_ib_{5,2}\\]\n\n\\[\\left.+\\frac{1.3.3a^2a_i^2}{2.4.2a_i^5}(e^2+e_i^2)b_{5,3}\\right\\}\\cos(nt-n_t+\\varepsilon-\\varepsilon_i)\\]\n\n\\[+m_i\\left\\{-\\frac{b_{1,2}}{a_i} + \\frac{aa_i}{2a_i^3}\\left\\{\\sin^2\\frac{\\varpi_i}{2} + \\frac{e^2+e_i^2}{2}\\right\\}\\left\\{b_{3,1} + b_{3,3}\\right\\} + \\frac{3(a^2e^2+a_i^2e_i^2)}{2.2a_i^3}b_{3,2}\\right.\\]\n\n\\[\\left.+\\frac{1.3.3a^2a_i^2}{2.4a_i^5}(e^2+e_i^2)b_{5,0} + \\frac{1.3.2}{2.4}(a^2e^2+a_i^2e_i^2)aa_ib_{5,1}\\right\\}\\]\n\n*The numbers at the side serve to distinguish the arguments.*\n\\[\n-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^4 e^2 + 5a^2 a_i^2 (e^2 + e_i^2) + 2a_i^4 e_i^2}{a_i^5} b_{5,2}\n\\]\n\\[\n+\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2 e^2 + a_i^2 e_i^2}{a_i^5} b_{5,3} \\cos(2nt - 2n_i t + 2\\varepsilon - 2\\varepsilon_i)\n\\]\n\\[\n+ m_1 \\left\\{ -\\frac{b_{1,3}}{a_i} + \\frac{a a_i}{2 a_i^3} \\left\\{ \\sin^2 \\frac{t}{2} + \\frac{e^2 + e_i^2}{2} \\right\\} b_{3,2} + \\frac{3(a^2 e^2 + a_i^2 e_i^2)}{2 \\cdot 2 a_i^3} b_{3,3} \\right\\}\n\\]\n\\[\n+\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2 a_i^3 (e^2 + e_i^2)}{a_i^5} b_{5,1} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(a^2 e^2 + a_i^2 e_i^2)}{4 a_i^5} a a_i b_{5,2}\n\\]\n\\[\n-\\frac{1}{2} \\cdot \\frac{3}{4} \\left\\{ 2a^4 e^2 + 5a^2 a_i^2 (e^2 + e_i^2) + 2a_i^4 e_i^2 \\right\\} b_{5,3} \\cos(3nt - 3n_i t + 3\\varepsilon - 3\\varepsilon_i)\n\\]\n\\[\n+ m_2 \\left\\{ -\\frac{b_{1,4}}{a_i} + \\frac{a a_i}{2 a_i^3} \\left\\{ \\sin^2 \\frac{t}{2} + \\frac{e^2 + e_i^2}{2} \\right\\} b_{3,3} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2 a_i^3 (e^2 + e_i^2)}{a_i^5} b_{5,2} \\right\\}\n\\]\n\\[\n+\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(a^2 e^2 + a_i^2 e_i^2)}{a_i^5} a a_i b_{5,3} \\cos(4nt - 4n_i t + 4\\varepsilon - 4\\varepsilon_i)\n\\]\n\\[\n+ m_3 \\left\\{ -\\frac{b_{1,5}}{a_i} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2 a_i^3 (e^2 + e_i^2)}{a_i^5} b_{5,3} \\right\\} \\cos(5nt - 5n_i t + 5\\varepsilon - 5\\varepsilon_i)\n\\]\n\\[\n+ m_4 \\left\\{ -\\frac{3a}{2 a_i^2} + \\frac{3a a_i}{2 a_i^3} b_{3,0} - \\frac{a^2}{2 a_i^3} b_{3,1} - \\frac{a a_i}{2 a_i^3} b_{3,2} \\right\\} e \\cos(nt + \\varepsilon - \\varpi)\n\\]\n\\[\n+ m_5 \\left\\{ -\\frac{a^2}{a_i^3} b_{3,0} + \\frac{a a_i}{2 a_i^3} b_{3,1} \\right\\} e \\cos(nt + \\varepsilon - \\varpi)\n\\]\n\\[\n+ m_6 \\left\\{ \\frac{1}{2} \\cdot \\frac{a}{a_i^2} - \\frac{a a_i}{2 a_i^3} b_{3,0} - \\frac{a^2}{2 a_i^3} b_{3,1} + \\frac{3}{2} \\cdot \\frac{a a_i}{a_i^3} b_{3,2} \\right\\} e \\cos(2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi)\n\\]\n\\[\n+ m_7 \\left\\{ -\\frac{a a_i}{2 a_i^3} b_{3,1} - \\frac{a^2}{2 a_i^3} b_{3,2} + \\frac{3}{2} \\cdot \\frac{a a_i}{a_i^3} b_{3,3} \\right\\} e \\cos(3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi)\n\\]\n\\[\n+ m_8 \\left\\{ -\\frac{a a_i}{2 a_i^3} b_{3,2} - \\frac{a^2}{2 a_i^3} b_{3,3} \\right\\} e \\cos(4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi)\n\\]\n\\[\n- m_9 \\frac{a a_i}{2 a_i^3} b_{3,3} e \\cos(5nt - 4n_i t + 5\\varepsilon - 4\\varepsilon_i - \\varpi)\n\\]\n\\[\n+ m_{10} \\left\\{ \\frac{3}{2} \\cdot \\frac{a a_i}{a_i^3} b_{3,1} - \\frac{a^2}{2 a_i^3} b_{3,2} - \\frac{a a_i}{2 a_i^3} b_{3,3} \\right\\} e \\cos(nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi)\n\\]\n\\[\n+ m_{11} \\left\\{ \\frac{3}{2} \\cdot \\frac{a a_i}{a_i^3} b_{3,2} - \\frac{a^2}{2 a_i^3} b_{3,3} \\right\\} e \\cos(2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi)\n\\]\n\\[\n+ m_{12} \\frac{3a a_i}{2 a_i^3} b_{3,3} e \\cos(3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi)\n\\]\n\\[\n+ m_i \\left\\{ - \\frac{a_i^2}{a_i^3} b_{3,0} - \\frac{a_i a_i}{2 a_i^3} b_{3,1} \\right\\} e_i \\cos(n_i t + \\varepsilon_i - \\varpi_i) \\quad [15]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{3 a_i a_i}{2 a_i^3} b_{3,1} - \\frac{a_i^3}{2 a_i^3} b_{3,2} - \\frac{a_i a_i}{2 a_i^3} b_{3,3} \\right\\} e_i \\cos(2nt - n_i t + 2\\varepsilon - \\varepsilon_i - \\varpi_i) \\quad [17]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{3 a_i a_i}{2 a_i^3} b_{3,2} - \\frac{a_i^3}{2 a_i^3} b_{3,3} \\right\\} e_i \\cos(3nt - 2n_i t + 3\\varepsilon - 2\\varepsilon_i - \\varpi_i) \\quad [18]\n\\]\n\n\\[\n+ m_i \\frac{3 a_i a_i}{2 a_i^3} b_{3,3} e_i \\cos(4nt - 3n_i t + 4\\varepsilon - 3\\varepsilon_i - \\varpi_i) \\quad [19]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{2 a_i}{a_i^3} - \\frac{a_i a_i}{2 a_i^3} b_{3,0} - \\frac{a_i^3}{2 a_i^3} b_{3,1} + \\frac{3 a_i a_i}{2 a_i^3} b_{3,2} \\right\\} e_i \\cos(nt - 2n_i t + \\varepsilon - 2\\varepsilon_i + \\varpi_i) \\quad [20]\n\\]\n\n\\[\n+ m_i \\left\\{ - \\frac{a_i a_i}{2 a_i^3} b_{3,1} - \\frac{a_i^3}{2 a_i^3} b_{3,2} + \\frac{3 a_i a_i}{2 a_i^3} b_{3,3} \\right\\} e_i \\cos(2nt - 3n_i t + 2\\varepsilon - 3\\varepsilon_i + \\varpi_i) \\quad [21]\n\\]\n\n\\[\n+ m_i \\left\\{ - \\frac{a_i a_i}{2 a_i^3} b_{3,2} - \\frac{a_i^3}{2 a_i^3} b_{3,3} \\right\\} e_i \\cos(3nt - 4n_i t + 3\\varepsilon - 4\\varepsilon_i + \\varpi_i) \\quad [22]\n\\]\n\n\\[\n- m_i \\frac{a_i a_i}{2 a_i^3} b_{3,3} e_i \\cos(4nt - 5n_i t + 4\\varepsilon - 5\\varepsilon_i + \\varpi_i) \\quad [23]\n\\]\n\n\\[\n+ m_i \\left\\{ - \\frac{a_i}{2 a_i^3} b_{3,1} - \\frac{a_i^3}{2 a_i^3} b_{3,2} - \\frac{3 a_i a_i}{2 a_i^3} b_{3,3} - \\frac{1.3 \\cdot 9 a_i^3 a_i^2}{2 a_i^5} b_{5,0} + \\frac{1.3 \\cdot 3 a_i^3 a_i}{2 a_i^5} b_{5,1} \\right\\} e^2 \\cos(2n_i t + 2\\varepsilon_i - 2\\varpi_i) \\quad [24]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{a_i}{8 a_i^3} - \\frac{a_i a_i}{2 a_i^3} b_{3,0} - \\frac{a_i^3}{2 a_i^3} b_{3,1} - \\frac{3 a_i a_i}{2 a_i^3} b_{3,2} + \\frac{1.3 \\cdot 6 a_i^3 a_i}{2 a_i^5} b_{5,0} \\right\\} e^2 \\cos(nt + n_i t + \\varepsilon_i - \\varpi_i) \\quad [25]\n\\]\n\n\\[\n+ m_i \\left\\{ - \\frac{a_i^2}{2 a_i^3} b_{3,0} - \\frac{a_i a_i}{2 a_i^3} b_{3,1} - \\frac{1.3 a_i^2 (2 a_i^2 - 3 a_i^3)}{2 a_i^5} b_{5,0} \\right\\} e^2 \\cos(2nt + 2\\varepsilon_i - 2\\varpi_i) \\quad [26]\n\\]\n\n\\[\n+ m_i \\left\\{ \\frac{3 a_i}{8 a_i^3} - \\frac{3 a_i a_i}{8 a_i^3} b_{3,0} - \\frac{a_i a_i}{2 a_i^3} b_{3,1} - \\frac{a_i^3}{2 a_i^3} b_{3,2} - \\frac{1.3 a_i^2 (4 a_i^3 - 5 a_i^3)}{2 a_i^5} b_{5,0} \\right\\} e^2 \\cos(2nt + 2\\varepsilon_i - 2\\varpi_i)\n\\]\n\\[\n-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{9}{a_i^5} b_{5,3} e^2 \\cos (3nt - n_t t + 3\\varepsilon - \\varepsilon_i - 2\\varphi)\n\\]\n\n\\[+ m_i \\left\\{ -\\frac{3}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,1} - \\frac{a_i^2}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,2} - \\frac{a_i}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,0}\n\\right.\\]\n\n\\[-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(a_i^2 - 3a_i^3)}{a_i^5} b_{5,2}\n\\]\n\n\\[\\left. + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} \\right\\} e^2 \\cos (4nt - 2n_t t + 4\\varepsilon - 2\\varepsilon_i - 2\\varphi)\n\\]\n\n\\[+ m_i \\left\\{ -\\frac{3}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,2} - \\frac{a_i^2}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,2}\n\\right.\\]\n\n\\[-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} \\right\\} e^2 \\cos (5nt - 3n_t t + 5\\varepsilon - 3\\varepsilon_i - 2\\varphi)\n\\]\n\n\\[+ m_i \\left\\{ -\\frac{3}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,2}\n\\right.\\]\n\n\\[-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} \\right\\} e^2 \\cos (6nt - 4n_t t + 6\\varepsilon - 4\\varepsilon_i - 2\\varphi)\n\\]\n\n\\[-m_i \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} e^2 \\cos (7nt - 5n_t t + 7\\varepsilon - 5\\varepsilon_i - 2\\varphi)\n\\]\n\n\\[+ m_i \\left\\{ -\\frac{aa_i}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,2} - \\frac{a_i^2}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,1} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,2}\n\\right.\\]\n\n\\[-\\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} \\right\\} e^2 \\cos (nt - 3n_t t + \\varepsilon - 3\\varepsilon_i + 2\\varphi)\n\\]\n\n\\[-m_i \\left\\{ -\\frac{aa_i}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,2}\n\\right.\\]\n\n\\[\\left. + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} \\right\\} e^2 \\cos (2nt - 4n_t t + 2\\varepsilon - 4\\varepsilon_i + 2\\varphi)\n\\]\n\n\\[-m_i \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a_i^3}{a_i^5} b_{5,3} e^2 \\cos (3nt - 5n_t t + 3\\varepsilon - 5\\varepsilon_i + 2\\varphi)\n\\]\n\n\\[+ m_i \\left\\{ -\\frac{5}{2} \\cdot \\frac{aa_i}{a_i^3} b_{3,1} + \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a_i^2}{a_i^5} b_{5,0} + \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{(a_i^3 + a_i^2)}{a_i^5} aa_i b_{5,1}\n\\right.\\]\n\n\\[-\\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a_i^3}{a_i^5} b_{5,2} \\right\\} ee_i \\cos (nt - n_t t + \\varepsilon - \\varepsilon_i - \\varphi + \\varphi_i)\n\\]\n\n\\[+ m_i \\left\\{ \\frac{a_i}{a_i^2} - \\frac{a_i}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,0} - \\frac{9}{2} \\cdot \\frac{a_i}{a_i^3} b_{3,2} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{(a_i^3 + a_i^2)}{a_i^5} aa_i b_{5,0}\n\\right.\\]\n\\[\n\\begin{align*}\n&+ \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,1} + \\frac{1.3 \\cdot 3}{2.4} \\frac{a^2 + a_i^2}{a_i^5} aa_i b_{5,2} \\\\\n&- \\frac{1.3 \\cdot 9}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,3} ee_i \\cos(2nt - 2n_t + 2\\varepsilon - 2\\varepsilon - \\varpi + \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{aa_i}{2.4} b_{3,1} - \\frac{9}{2.4} \\frac{aa_i}{a_i^3} b_{3,3} - \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,0} - \\frac{1.3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,1} + \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,2} \\right. \\\\\n&+ \\frac{1.3 \\cdot 3}{2.4} \\frac{a^2 + a_i^2}{a_i^5} aa_i b_{5,3} \\left\\} ee_i \\cos(3nt - 3n_t + 3\\varepsilon - 3\\varepsilon - \\varpi + \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{aa_i}{2.4} b_{3,2} - \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,1} - \\frac{1.3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,2} \\right. \\\\\n&+ \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,3} \\left\\} ee_i \\cos(4nt - 4n_t + 4\\varepsilon - 4\\varepsilon - \\varpi + \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{aa_i}{2.4} b_{3,3} - \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,2} \\right. \\\\\n&- \\frac{1.3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,3} \\left\\} ee_i \\cos(5nt - 5n_t + 5\\varepsilon - 5\\varepsilon - \\varpi + \\varpi) \\\\\n&- m_i \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,3} \\cos(6nt - 6n_t - \\varpi + \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{9}{2.4} \\frac{aa_i}{a_i^3} b_{3,0} - \\frac{aa_i}{2.4} b_{3,2} + \\frac{1.3 \\cdot 6}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,0} - \\frac{1.3 \\cdot 7}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,1} \\right. \\\\\n&- \\frac{1.3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,2} - \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,3} \\left\\} ee_i \\cos(\\varpi - \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{9}{2.4} \\frac{aa_i}{a_i^3} b_{3,1} - \\frac{aa_i}{2.4} b_{3,3} - \\frac{1.3 \\cdot 9}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,0} + \\frac{1.3 \\cdot 3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,1} \\right. \\\\\n&+ \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,2} - \\frac{1.3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,3} \\left\\} ee_i \\cos(nt - n_t + \\varepsilon - \\varepsilon + \\varpi - \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{9}{2.4} \\frac{aa_i}{a_i^3} b_{3,2} - \\frac{1.3 \\cdot 9}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,1} + \\frac{1.3 \\cdot 3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,2} \\right. \\\\\n&+ \\frac{1.3}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,3} \\left\\} ee_i \\cos(2nt - 2n_t + 2\\varepsilon - 2\\varepsilon + \\varpi - \\varpi) \\\\\n&+ m_i \\left\\{ - \\frac{9}{2.4} \\frac{aa_i}{a_i^3} b_{3,3} - \\frac{1.3 \\cdot 9}{2.4} \\frac{a^2 a_i^2}{a_i^5} b_{5,2} \\right. \\\\\n&+ \\frac{1.3 \\cdot 3}{2.4} \\frac{(a^2 + a_i^2)}{a_i^5} aa_i b_{5,3} \\left\\} ee_i \\cos(3nt - 3n_t + 3\\varepsilon - 3\\varepsilon + \\varpi - \\varpi)\n\\end{align*}\n\\]\n\\[ - m_1 \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{9}{2} \\cdot \\frac{a^2}{a_i^5} b_{5,3} ee_t \\cos(4nt - 4n_t + 4\\varepsilon - 4\\varepsilon_t + \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ - \\frac{3}{2} \\cdot \\frac{a}{a_i^3} + \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,0} + \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,2} + \\frac{1}{1} \\cdot \\frac{3}{2} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,0} - \\frac{1}{2} \\cdot \\frac{3}{1} \\cdot \\frac{a^6}{a_i^5} b_{5,1} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,2} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,3} \\right\\} ee_t \\cos(2nt + 2\\varepsilon - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,1} - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{14}{a_i^5} b_{5,0} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(a^2 + a^2)}{a_i^5} aa_i b_{5,1} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,2} \\right\\} ee_t \\cos(nt + n_t + \\varepsilon + \\varepsilon_t - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,0} + \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,2} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,0} - \\frac{1}{2} \\cdot \\frac{3}{1} \\cdot \\frac{a^6}{a_i^5} b_{5,1} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,2} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,3} \\right\\} ee_t \\cos(2nt + 2\\varepsilon - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,1} + \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,3} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{7} \\cdot \\frac{a^2}{a_i^5} b_{5,2} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,3} \\right\\} ee_t \\cos(3nt - n_t + 3\\varepsilon - \\varepsilon_t - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,2} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,1} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,2} \\right\\} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{7} \\cdot \\frac{a^2}{a_i^5} b_{5,3} \\right\\} ee_t \\cos(4nt - 2n_t + 4\\varepsilon - 2\\varepsilon_t - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,3} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} aa_i b_{5,2} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,3} \\right\\} ee_t \\cos(5nt - 3n_t + 5\\varepsilon - 3\\varepsilon_t - \\varphi - \\varphi_t) \\]\n\n\\[ + m_1 \\left\\{ \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,3} ee_t \\cos(6nt - 4n_t + 6\\varepsilon - 4\\varepsilon_t - \\varphi - \\varphi_t) \\right\\} \\]\n\n\\[ + m_1 \\left\\{ \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,1} + \\frac{3}{2} \\cdot \\frac{a}{a_i^3} b_{3,3} + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^2}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(3a^2 - a^2)}{a_i^5} aa_i b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{7} \\cdot \\frac{a^2}{a_i^5} b_{5,2} \\]\n\\[ \\frac{1}{2} \\left( \\frac{3a_i^2 - a_i^3}{a_i^5} \\right) aa_i b_{5,3} \\right\\} ee_i \\cos(nt - 3nt + \\varepsilon - 3\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{3}{2.4} \\frac{aa_i}{a_i^3} b_{5,1} + \\frac{3}{2.4} \\frac{aa_i}{a_i^3} b_{5,2} + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ + \\frac{1}{2.4} \\left( \\frac{3a_i^2 - a_i^3}{a_i^5} \\right) aa_i b_{5,1} + \\frac{1}{2.4} \\left( \\frac{3a_i^2 - a_i^3}{a_i^5} \\right) aa_i b_{5,2} \\]\n\n\\[ - \\frac{1}{2.4} \\left( \\frac{3a_i^2 - a_i^3}{a_i^5} \\right) b_{5,3} \\right\\} ee_i \\cos(2nt - 4nt + 2\\varepsilon - 4\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{3}{2.4} \\frac{aa_i}{a_i^3} b_{5,3} + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,2} \\right\\} \\]\n\n\\[ + \\frac{1}{2.4} \\left( \\frac{3a_i^2 - a_i^3}{a_i^5} \\right) aa_i b_{5,3} \\right\\} ee_i \\cos(3nt - 5nt + 3\\varepsilon - 5\\varepsilon_i + \\varpi + \\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,3} ee_i \\cos(4nt - 6nt + 4\\varepsilon - 6\\varepsilon_i + \\varpi + \\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ - \\frac{a_i^3}{2.2a_i^3} b_{5,0} - \\frac{aa_i}{2.2a_i^3} b_{5,1} - \\frac{1}{2.4} \\frac{(2a_i^3 - 3a_i^2)}{a_i^5} a_i^2 b_{5,0} \\right\\} \\]\n\n\\[ + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,1} - \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,2} \\right\\} e_i^2 \\cos(2nt + 2\\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{a_i}{8a_i^3} - \\frac{aa_i}{2.4a_i^3} b_{5,0} - \\frac{a_i^3}{2.4a_i^3} b_{5,1} - \\frac{3}{2.8a_i^3} b_{5,2} + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ - \\frac{1}{2.4} \\frac{(4a_i^3 + 3a_i^2)}{a_i^5} a_i^2 b_{5,1} - \\frac{1}{2.4} \\frac{aa_i^3}{a_i^5} b_{5,2} \\]\n\n\\[ - \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,3} \\right\\} e_i^2 \\cos(nt + n_i t + \\varepsilon + \\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{a_i}{2.8a_i^3} b_{5,1} - \\frac{a_i^3}{2.4a_i^3} b_{5,2} - \\frac{3}{2.8a_i^3} b_{5,3} - \\frac{1}{2.4} \\frac{9a_i^2}{a_i^5} b_{5,0} + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,1} \\right\\} \\]\n\n\\[ - \\frac{1}{2.4} \\frac{(2a_i^3 - 3a_i^2)}{a_i^5} a_i^2 b_{5,2} - \\frac{1}{2.4} \\frac{aa_i^3}{a_i^5} b_{5,3} \\right\\} e_i^2 \\cos(2nt + 2\\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{aa_i}{2.8a_i^3} b_{5,2} - \\frac{a_i^3}{2.4a_i^3} b_{5,3} - \\frac{1}{2.4} \\frac{9a_i^2}{a_i^5} b_{5,1} + \\frac{1}{2.4} \\frac{3a_i^2}{a_i^5} b_{5,2} \\right\\} \\]\n\n\\[ - \\frac{1}{2.4} \\frac{(2a_i^3 - 3a_i^2)}{a_i^5} a_i^2 b_{5,3} \\right\\} e_i^2 \\cos(3nt - n_i t + 3\\varepsilon - \\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ - \\frac{aa_i}{2.8a_i^3} b_{5,3} - \\frac{1}{2.4} \\frac{9a_i^2}{a_i^5} b_{5,2} \\right\\} \\]\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^3}{a_i^5} b_{5,3} \\left\\{ e_i^2 \\cos(4nt - 2n_i t + 4\\varepsilon - 2\\varepsilon_i - 2\\varpi_i) \\right\\} \\]\n\n\\[ - m_i \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^3}{a_i^5} b_{5,3} e_i^2 \\cos(5nt - 3n_i t + 5\\varepsilon - 3\\varepsilon_i - 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{27}{8} \\frac{a}{a_i^3} - \\frac{3}{2} \\frac{a a_i}{a_i^3} b_{3,0} - \\frac{a_i^3}{2 \\cdot 4 \\cdot a_i^3} b_{3,1} - \\frac{aa_i}{2 \\cdot 8 \\cdot a_i^3} b_{3,2} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a^3}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(4a_i^2 - 5a^2)}{a_i^5} a_i^2 b_{5,1} + \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{aa_i^3}{a_i^5} b_{5,2} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^3}{a_i^5} b_{5,3} \\left\\{ e_i^2 \\cos(nt - 3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ - \\frac{3}{2} \\frac{aa_i}{a_i^3} b_{3,1} - \\frac{a_i^3}{2 \\cdot 4 \\cdot a_i^3} b_{3,2} - \\frac{aa_i}{2 \\cdot 8 \\cdot a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a^3}{a_i^5} b_{5,0} \\right\\} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{aa_i^3}{a_i^5} b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{(2a_i^2 - 3a^2)}{a_i^5} a_i^2 b_{5,2} \\]\n\n\\[ + \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{aa_i^3}{a_i^5} b_{5,3} \\left\\{ e_i^2 \\cos(2nt - 4n_i t + 2\\varepsilon - 4\\varepsilon_i + 2\\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ - \\frac{3}{2} \\frac{aa_i}{a_i^3} b_{3,2} - \\frac{a_i^3}{2 \\cdot 4 \\cdot a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a^3}{a_i^5} b_{5,1} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{aa_i^3}{a_i^5} b_{5,2} \\right\\} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{(2a_i^2 - 3a^2)}{a_i^5} a_i^2 b_{5,3} \\left\\{ e_i^2 \\cos(3nt - 5n_i t + 3\\varepsilon - 3\\varepsilon_i + 2\\varpi_i) \\right\\} \\]\n\n\\[ + m_i \\left\\{ - \\frac{3}{2} \\frac{aa_i}{a_i^3} b_{3,3} - \\frac{1}{2} \\cdot \\frac{3}{2} \\cdot \\frac{a^3}{a_i^5} b_{5,2} \\right\\} \\]\n\n\\[ - \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{aa_i^3}{a_i^5} b_{5,3} \\left\\{ e_i^2 \\cos(4nt - 6n_i t + 4\\varepsilon - 6\\varepsilon_i + 2\\varpi_i) \\right\\} \\]\n\n\\[ - m_i \\frac{1}{2} \\cdot \\frac{3}{4} \\cdot \\frac{a^3}{a_i^5} b_{5,3} e_i^2 \\cos(5nt - 7n_i t + 5\\varepsilon - 7\\varepsilon_i + 2\\varpi_i) \\]\n\n\\[ + m_i \\left\\{ \\frac{a}{a_i^3} \\left( 1 + \\frac{3}{2} e_i^2 \\right) - \\frac{aa_i}{a_i^3} b_{3,0} \\right\\} \\sin^2 \\frac{n_i}{2} \\cos(nt + n_i t + \\varepsilon + \\varepsilon_i - 2\\nu_i) \\]\n\n\\[ - m_i \\frac{a a_i}{2 a_i^3} b_{3,1} \\sin^2 \\frac{n_i}{2} \\cos(2n_i t + 2\\varepsilon_i - 2\\nu_i) \\]\n\n\\[ m_i \\frac{a a_i}{2 a_i^3} b_{3,1} \\sin^2 \\frac{n_i}{2} \\cos(2nt + 2\\varepsilon_i - 2\\nu_i) \\]\n\n\\[ - m_i \\frac{a a_i}{2 a_i^3} b_{3,2} \\sin^2 \\frac{n_i}{2} \\cos(nt + 3n_i t + \\varepsilon - 3\\varepsilon_i + 2\\nu_i) \\]\n\\[-m_i \\frac{a_i}{2} a_i^3 b_{3i} \\sin^2 \\frac{\\theta_i}{2} \\cos(3nt - n_t + 3\\varepsilon_i - \\varepsilon_i - 2\\nu_i)\\]\n\n\\[-m_i \\frac{a_i}{2} a_i^3 b_{3i} \\sin^2 \\frac{\\theta_i}{2} \\cos(2nt - 4n_t + 2\\varepsilon_i - 4\\varepsilon_i + 2\\nu_i)\\]\n\n\\[-m_i \\frac{a_i}{2} a_i^3 b_{3i} \\sin^2 \\frac{\\theta_i}{2} \\cos(4nt - 2n_t + 4\\varepsilon_i - 2\\varepsilon_i - 2\\nu_i)\\]\n\nIn the development of \\(R\\), I have supposed \\(\\iota = 0\\), so that \\(\\iota_i\\) is the angle contained between the orbits of the planets \\(P\\) and \\(P_i\\), or \\(P_1\\) and \\(P_2\\); in the general case, when \\(\\iota_1\\) and \\(\\iota_2\\) are retained, \\(\\cos \\iota = \\cos \\iota_1 \\cos \\iota_2 + \\sin \\iota_1 \\sin \\iota_2 \\cos (\\iota_1 - \\iota_2)\\), \\(\\iota_1\\) and \\(\\iota_2\\) being the inclinations of the orbits of the planets \\(P_1, P_2\\) to any plane \\(xy\\), of which the direction is arbitrary,\n\n\\[r' r_i \\sin (\\lambda' - \\lambda_i) = r_1 r_2 \\left\\{ \\cos^2 \\frac{\\theta_1}{2} \\cos^2 \\frac{\\theta_2}{2} \\sin (\\lambda_1 - \\lambda_2) \\right.\\]\n\n\\[-\\sin^2 \\frac{\\theta_1}{2} \\cos^2 \\frac{\\theta_2}{2} \\sin (\\lambda_1 + \\lambda_2 - 2\\nu_1) + \\sin^2 \\frac{\\theta_2}{2} \\cos^2 \\frac{\\theta_1}{2} \\sin (\\lambda_1 + \\lambda_2 - 2\\nu_2)\\]\n\n\\[-\\sin^2 \\frac{\\theta_1}{2} \\cos^2 \\frac{\\theta_2}{2} \\sin (\\lambda_1 - \\lambda_2 - 2\\nu_1 + 2\\nu_2)\\}\n\n\\[r' r_i \\left\\{ \\cos (\\lambda' - \\lambda_i) + s_i s_i \\right\\} = r_1 r_2 \\left\\{ \\cos^2 \\frac{\\theta_1}{2} \\cos^2 \\frac{\\theta_2}{2} \\cos (\\lambda_1 - \\lambda_2)\\right.\\]\n\n\\[\\sin^2 \\frac{\\theta_1}{2} \\cos^2 \\frac{\\theta_2}{2} \\cos (\\lambda_1 + \\lambda_2 - 2\\nu_1) + \\sin^2 \\frac{\\theta_2}{2} \\cos^2 \\frac{\\theta_1}{2} \\cos (\\lambda_1 + \\lambda_2 - 2\\nu_2)\\]\n\n\\[\\sin^2 \\frac{\\theta_1}{2} \\sin^2 \\frac{\\theta_2}{2} \\cos (\\lambda_1 - \\lambda_2 - 2\\nu_1 + 2\\nu_2)\\]\n\n\\[\\sin \\iota_1 \\sin \\iota_2 \\cos (\\lambda_1 - \\lambda_2 - \\nu_1 + \\nu_2) - \\sin \\iota_1 \\sin \\iota_2 \\cos (\\lambda_1 + \\lambda_2 - \\nu_1 - \\nu_2)\\}\n\nErrata.\n\nIn page 330, for \\(\\iota\\) and \\(\\nu\\), read \\(\\iota_i\\) and \\(\\nu_i\\).",
  "source": "olmocr",
  "added": "2026-01-12",
  "created": "2026-01-12",
  "metadata": {
    "Source-File": "/home/jic823/projects/def-jic823/royalsociety/pdfs/107911.pdf",
    "olmocr-version": "0.3.4",
    "pdf-total-pages": 32,
    "total-input-tokens": 54638,
    "total-output-tokens": 38461,
    "total-fallback-pages": 0
  },
  "attributes": {
    "pdf_page_numbers": [
      [
        0,
        0,
        1
      ],
      [
        0,
        1767,
        2
      ],
      [
        1767,
        3580,
        3
      ],
      [
        3580,
        6145,
        4
      ],
      [
        6145,
        8762,
        5
      ],
      [
        8762,
        11083,
        6
      ],
      [
        11083,
        13750,
        7
      ],
      [
        13750,
        16801,
        8
      ],
      [
        16801,
        19901,
        9
      ],
      [
        19901,
        21837,
        10
      ],
      [
        21837,
        23954,
        11
      ],
      [
        23954,
        26247,
        12
      ],
      [
        26247,
        28930,
        13
      ],
      [
        28930,
        31257,
        14
      ],
      [
        31257,
        33751,
        15
      ],
      [
        33751,
        36168,
        16
      ],
      [
        36168,
        38660,
        17
      ],
      [
        38660,
        41378,
        18
      ],
      [
        41378,
        43603,
        19
      ],
      [
        43603,
        44873,
        20
      ],
      [
        44873,
        47341,
        21
      ],
      [
        47341,
        49714,
        22
      ],
      [
        49714,
        52022,
        23
      ],
      [
        52022,
        54104,
        24
      ],
      [
        54104,
        56927,
        25
      ],
      [
        56927,
        59308,
        26
      ],
      [
        59308,
        62241,
        27
      ],
      [
        62241,
        64935,
        28
      ],
      [
        64935,
        68116,
        29
      ],
      [
        68116,
        70718,
        30
      ],
      [
        70718,
        73501,
        31
      ],
      [
        73501,
        75498,
        32
      ]
    ],
    "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"
    ],
    "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
    ],
    "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
    ],
    "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
    ],
    "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
    ]
  },
  "jstor_metadata": {
    "identifier": "jstor-107911",
    "title": "Researches in Physical Astronomy",
    "authors": "John William Lubbock",
    "year": 1830,
    "volume": "120",
    "journal": "Philosophical Transactions of the Royal Society of London",
    "page_count": 32,
    "jstor_url": "https://www.jstor.org/stable/107911"
  }
}