# Scientific Papers of Josiah Willard Gibbs, Volume 2/Chapter V

V.
ON THE DETERMINATION OF ELLIPTIC ORBITS FROM THREE COMPLETE OBSERVATIONS.
[Memoirs of the National Academy of Sciences, vol. iv. part ii. pp. 79–104, 1889.]
The determination of an orbit from three complete observations by the solution of the equations which represent elliptic motion presents so great difficulties in the general case, that in the first solution of the problem we must generally limit ourselves to the case in which the intervals between the observations are not very long. In this case we substitute some comparatively simple relations between the unknown quantities of the problem, which have an approximate validity for short intervals, for the less manageable relations which rigorously subsist between these quantities. A comparison of the approximate solution thus obtained with the exact laws of elliptic motion will always afford the means of a closer approximation, and by a repetition of this process we may arrive at any required degree of accuracy.
It is therefore a problem not without interest—it is, in fact, the natural point of departure in the study of the determination of orbits—to express in a manner combining as far as possible simplicity and accuracy the relations between three positions in an orbit separated by small or moderate intervals. The problem is not entirely determinate, for we may lay the greater stress upon simplicity or upon accuracy; we may seek the most simple relations which are sufficiently accurate to give us any approximation to an orbit, or we may seek the most exact expression of the real relations, which shall not be too complex to be serviceable.
Derivation of the Fundamental Equation.
​If we adopt such a unit of time that the acceleration due to the sun's action is unity at a unit's distance, and denote the vectors drawn from the sun to the body in its three positions by ${\mathfrak {R}}_{1},{\mathfrak {R}}_{2},{\mathfrak {R}}_{3},$ and the lengths of these vectors (the heliocentric distances) by $r_{1},r_{2},r_{3},$ the accelerations corresponding to the three positions will be represented by $-{\frac {{\mathfrak {R}}_{1}}{r_{1}^{3}}},-{\frac {{\mathfrak {R}}_{2}}{r_{2}^{3}}},-{\frac {{\mathfrak {R}}_{3}}{r_{3}^{3}}}\cdot$ Now the motion between the positions considered may be expressed with a high degree of accuracy by an equation of the form
* ${\mathfrak {R}}={\mathfrak {A}}+t{\mathfrak {B}}+t^{2}{\mathfrak {C}}+t^{3}{\mathfrak {D}}+t^{4}{\mathfrak {E}},$
$\tau _{e}$
$\tau _{1}$
$t=0$
for $t=-\tau _{3},$
* ${\mathfrak {R}}={\mathfrak {R}}_{1},$${\frac {d^{2}{\mathfrak {R}}}{dt^{2}}}=-{\frac {{\mathfrak {R}}_{1}}{r_{1}^{3}}};$
$t=0,$
* ${\mathfrak {R}}={\mathfrak {R}}_{2},$${\frac {d^{2}{\mathfrak {R}}}{dt^{2}}}=-{\frac {{\mathfrak {R}}_{2}}{r_{2}^{3}}};$
$t=\tau _{1},$
* ${\mathfrak {R}}={\mathfrak {R}}_{3},$${\frac {d^{2}{\mathfrak {R}}}{dt^{2}}}=-{\frac {{\mathfrak {R}}_{3}}{r_{3}^{3}}}\cdot$
* ${\begin{aligned}{\mathfrak {R}}_{1}&={\mathfrak {A}}-\tau _{3}{\mathfrak {B}}+\tau _{3}^{2}{\mathfrak {C}}-\tau _{3}^{3}{\mathfrak {D}}+\tau _{3}^{4}{\mathfrak {E}}\\{\mathfrak {R}}_{2}&={\mathfrak {A}}\\{\mathfrak {R}}_{3}&={\mathfrak {A}}+\tau _{3}{\mathfrak {B}}+\tau _{1}^{2}{\mathfrak {C}}+\tau _{1}^{3}{\mathfrak {D}}-\tau _{1}^{4}{\mathfrak {E}}\\-{\frac {{\mathfrak {R}}_{1}}{r_{1}}}&=2{\mathfrak {C}}-6\tau _{3}{\mathfrak {D}}+12\tau _{3}^{2}{\mathfrak {E}}\\-{\frac {{\mathfrak {R}}_{2}}{r_{2}^{3}}}&=2{\mathfrak {E}}\\-{\frac {{\mathfrak {R}}_{3}}{r_{3}^{3}}}&=2{\mathfrak {E}}+6\tau _{1}{\mathfrak {D}}+12\tau _{1}^{2}{\mathfrak {E}}\end{aligned}}$
${\mathfrak {A,B,C,D,E}}$
| $A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right){\mathfrak {R}}_{1}-\left(1-{\frac {B_{2}}{r_{2}^{3}}}\right){\mathfrak {R}}_{2}+A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right){\mathfrak {R}}_{3}=0,$ | (1) |
* $A_{1}={\frac {\tau _{1}}{\tau _{1}+\tau _{3}}},$$A_{3}={\frac {\tau _{3}}{\tau _{1}+\tau _{2}}},$
* $B_{1}={\tfrac {1}{12}}(-\tau _{1}^{2}+\tau _{1}\tau _{3}+\tau _{3}^{2}),$$B_{2}={\tfrac {1}{12}}(\tau _{1}^{2}+3\tau _{1}\tau _{3}+\tau _{3}^{2})$
* $B_{3}={\tfrac {1}{12}}(\tau _{1}^{2}+\tau _{1}\tau _{3}-\tau _{3}^{2}).$
| $n_{1}=A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right),$$n_{2}=\left(1-{\frac {B_{2}}{r_{2}^{3}}}\right),$$n_{3}=A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right)\cdot$ | (2) |
| $n_{1}{\mathfrak {R}}_{2}-n_{2}{\mathfrak {R}}_{2}+n_{3}{\mathfrak {R}}_{2}=0.$ | (3) |
$n_{2}{\mathfrak {R}}_{2}$
$n_{1}{\mathfrak {R}}_{1}$
$n_{3}{\mathfrak {R}}_{3}$
${\mathfrak {R}}_{3}$
${\mathfrak {R}}_{1},$
| $n_{1}{\mathfrak {R}}_{1}\times {\mathfrak {R}}_{3}-n_{2}{\mathfrak {R}}_{2}\times {\mathfrak {R}}_{3}=0,$$-n_{2}{\mathfrak {R}}_{1}\times {\mathfrak {R}}_{2}+n_{3}{\mathfrak {R}}_{1}\times {\mathfrak {R}}_{3}=0,$ | (4) |
| ${\frac {{\mathfrak {R}}_{2}\times {\mathfrak {R}}_{3}}{n_{1}}}={\frac {{\mathfrak {R}}_{1}\times {\mathfrak {R}}_{3}}{n_{2}}}={\frac {{\mathfrak {R}}_{1}\times {\mathfrak {R}}_{2}}{n_{3}}}\cdot$ | (5) |
${\mathfrak {R}}_{1},{\mathfrak {R}}_{2},{\mathfrak {R}}_{3}$
$[r_{2}r_{3}],[r_{1}r_{3}],[r_{1}r_{2}],$
* $A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right)$$\left(1-{\frac {B_{2}}{r_{2}^{3}}}\right)$$A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right)\cdot$
${\mathfrak {E}}_{1},{\mathfrak {E}}_{2},{\mathfrak {E}}_{3}$
${\mathfrak {F}}_{1},{\mathfrak {F}}_{2},{\mathfrak {F}}_{3}$
$\rho _{1},\rho _{2},\rho _{3}$
| ${\mathfrak {R}}_{1}={\mathfrak {E}}_{1}+\rho _{1}{\mathfrak {F}}_{1},$${\mathfrak {R}}_{2}={\mathfrak {E}}_{2}+\rho _{2}{\mathfrak {F}}_{2},$${\mathfrak {R}}_{3}={\mathfrak {E}}_{3}+\rho _{3}{\mathfrak {F}}_{3}.$ | (6) |
| $A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right)({\mathfrak {E}}_{1}+\rho _{1}{\mathfrak {F}}_{1})-\left(1-{\frac {B_{2}}{r_{2}^{3}}}\right)({\mathfrak {E}}_{2}+\rho _{2}{\mathfrak {F}}_{2})+A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right)({\mathfrak {E}}_{3}+\rho _{3}{\mathfrak {F}}_{3})=0,$ | (7) |
$\rho _{1},\rho _{2},\rho _{3},r_{1},r_{2},r_{3}$
| $r_{1}^{2}={\mathfrak {E}}_{3}^{2}+2({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})\rho _{3}+\rho _{3}^{2},\,\,\,r_{2}^{2}={\mathfrak {E}}_{2}^{2}+2({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})\rho _{2}+\rho _{2}^{2},$ | $\scriptstyle {\left.{\begin{matrix}\ \\\ \end{matrix}}\right\}\,$ (8) |
| $\scriptstyle {\left.{\begin{matrix}\ \\\ \end{matrix}}\right\}\,$ (8) ||
$r_{1},r_{2},r_{3}$
$\rho _{1},\rho _{2},\rho _{3},$
$r_{1},r_{2},r_{3},\rho _{1},\rho _{2},\rho _{3};$
$r_{1},r_{2},r_{3}$
$\rho _{1},\rho _{2},\rho _{3},$
$\rho _{1},\rho _{2},\rho _{3}.$
It remains to be shown, first, how the numerical solution of the equation may be performed, and secondly, how such an approximate solution of the actual problem may furnish the basis of a closer approximation.
Solution of the Fundamental Equation.
The relations with which we have to do will be rendered a little more simple if instead of each geocentric distance we introduce the ​distance of the body from the foot of the perpendicular from the sun upon the line of sight. If we set
| $q_{1}=\rho _{1}+({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1}),$$q_{2}=\rho _{2}+({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2}),$$q_{3}=\rho _{3}+({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3}),$ | (9) |
| $p_{1}^{2}={\mathfrak {E}}_{1}^{2}-({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})^{2},$$p_{2}^{2}={\mathfrak {E}}_{2}^{2}-({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})^{2},$$p_{3}^{2}={\mathfrak {E}}_{3}^{2}-({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})^{2},$ | (10) |
| $r_{1}^{2}=q_{1}^{2}+p_{1}^{2},$$r_{2}^{2}=q_{2}^{2}+p_{2}^{2},$$r_{3}^{2}=q_{3}^{2}+p_{3}^{2}.$ | (11) |
| ${\mathfrak {S}}_{1}=A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right)({\mathfrak {E}}_{1}+\rho _{1}{\mathfrak {F}}_{1}),$${\mathfrak {E}}_{2}=-\left(1-{\frac {B_{2}}{r_{2}^{3}}}\right)({\mathfrak {E}}_{2}+\rho _{2}{\mathfrak {F}}_{2}),$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (12) |
| ${\mathfrak {E}}_{3}=A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right)({\mathfrak {E}}_{3}+\rho _{3}{\mathfrak {F}}_{3}).$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (12) |
${\mathfrak {S}}_{1},{\mathfrak {S}}_{2},{\mathfrak {S}}_{3}$
$\rho _{1},\rho _{2},\rho _{3},$
$q_{1},q_{2},q_{3},$
| ${\mathfrak {S}}'={\frac {d{\mathfrak {S}}_{1}}{dq_{1}}},$${\mathfrak {S}}''={\frac {d{\mathfrak {S}}_{2}}{dq_{2}}},$${\mathfrak {S}}'''={\frac {d{\mathfrak {S}}_{3}}{dq_{3}}},$ | (13) |
| ${\mathfrak {S}}={\mathfrak {S}}_{1}+{\mathfrak {S}}_{2}+{\mathfrak {S}}_{3},$ | (14) |
| $d{\mathfrak {S}}={\mathfrak {S}}'dq_{1}+{\mathfrak {S}}''dq_{2}+{\mathfrak {S}}'''dq_{3}.$ | (15) |
${\mathfrak {S}}',$
| ${\mathfrak {S}}'=A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right){\mathfrak {F}}_{1}-A_{1}{\frac {3B_{1}}{r_{1}^{4}}}{\frac {dr_{1}}{dq_{1}}}({\mathfrak {E}}_{1}+\rho _{1}{\mathfrak {F}}_{1}).$ | (16) |
| ${\frac {dr_{1}}{dq_{1}}}={\frac {q_{1}}{r_{1}}}\cdot$ | (17) |
| ${\mathfrak {S}}'=A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right){\mathfrak {F}}_{1}-{\frac {3B_{1}q_{1}}{r_{1}^{5}(1+B_{1}r_{1}^{-3})}}{\mathfrak {S}}_{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (18) |
| ${\mathfrak {S}}''=-\left(1+{\frac {B_{2}}{r_{2}^{3}}}\right){\mathfrak {F}}_{2}+{\frac {3B_{2}q_{2}}{r_{2}^{5}(1+B_{2}r_{2}^{-3})}}{\mathfrak {S}}_{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (18) |
| ${\mathfrak {S}}'''=A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right){\mathfrak {F}}_{3}-{\frac {3B_{3}q_{3}}{r_{3}^{5}(1+B_{3}r_{3}^{-3})}}{\mathfrak {S}}_{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (18) |
Now if any values of $q_{1},q_{2},q_{3}$ (either assumed or obtained by a previous approximation) give a certain residual ${\mathfrak {S}}$ (which would be zero if the values of $q_{1},q_{2},q_{3}$ satisfied the fundamental equation), and we wish to find the corrections $\Delta q_{1},\Delta q_{2},\Delta q_{3}$ which must be added to $q_{1},q_{2},q_{3}$ to reduce the residual to zero, we may apply equation (15) to these finite differences, and will have approximately, when these differences are not very large,
| $-{\mathfrak {S}}={\mathfrak {S}}'\Delta q_{1}+{\mathfrak {S}}''\Delta q_{2}+{\mathfrak {S}}'''\Delta q_{3}.$ | (19) |
| $\Delta q_{1}=-{\frac {({\mathfrak {SS}}''{\mathfrak {S}}''')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}}$$\Delta q_{2}=-{\frac {({\mathfrak {SS}}'''{\mathfrak {S}}')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}}$$\Delta q_{3}=-{\frac {({\mathfrak {S}}{\mathfrak {S}}'{\mathfrak {S}}'')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}}\cdot$ | (20) |
$q_{1},q_{2},q_{3}$
${\mathfrak {S}},$
$q_{1},q_{2},q_{3}.$
It will sometimes be worth while to use formulæ a little less simple for the sake of a more rapid approximation. Instead of equation (19) we may write, with a higher degree of accuracy,
| $-{\mathfrak {S}}={\mathfrak {S}}'\Delta q_{1}+{\mathfrak {S}}''\Delta q_{2}+{\mathfrak {S}}'''\Delta q_{3}+{\tfrac {1}{2}}{\mathfrak {T}}'(\Delta q_{1})^{2}+{\tfrac {1}{2}}{\mathfrak {T}}''(\Delta q_{2})^{2}+{\tfrac {1}{2}}{\mathfrak {T}}'''(\Delta q_{3})^{2},$ | (21) |
| ${\mathfrak {T}}'=$ | ${\frac {d^{2}{\mathfrak {S}}_{1}}{dq_{1}^{2}}}=$ | $2A_{1}B_{1}{\frac {d(r_{1}^{-3})}{dq_{1}}}{\mathfrak {F}}_{1}+$ | ${\frac {B_{1}}{1+B_{1}r_{1}^{-3}}}{\frac {d^{2}(r_{1}^{-3})}{dq_{1}^{2}}}{\mathfrak {S}}_{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (22) |
| ${\mathfrak {T}}''=$ | ${\frac {d^{2}{\mathfrak {S}}_{2}}{dq_{2}^{2}}}=$ | $2B_{2}{\frac {d(r_{2}^{-3})}{dq_{2}}}{\mathfrak {F}}_{2}-$ | ${\frac {B_{2}}{1-B_{2}r_{2}^{-3}}}{\frac {d^{2}(r_{2}^{-3})}{dq_{2}^{2}}}{\mathfrak {S}}_{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (22) |
| ${\mathfrak {T}}'''=$ | ${\frac {d^{2}{\mathfrak {S}}_{3}}{dq_{3}^{2}}}=$ | $2A_{3}B_{3}{\frac {d(r_{3}^{-3})}{dq_{3}}}{\mathfrak {F}}_{3}+$ | ${\frac {B_{3}}{1+B_{3}r_{3}^{-3}}}{\frac {d^{2}(r_{3}^{-3})}{dq_{3}^{2}}}{\mathfrak {S}}_{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ (22) |
It is evident that ${\mathfrak {T}}''$ is generally many times greater than ${\mathfrak {T}}'$ or ${\mathfrak {T}}''',$ the factor $B_{2},$ in the case of equal intervals, being exactly ten times as great as $A_{1}B_{1}$ or $A_{3}B_{3}.$ This shows, in the first place, that the accurate determination of $\Delta q_{2}$ is of the most importance for the subsequent approximations. It also shows that we may attain nearly the same accuracy in writing
| $-{\mathfrak {S}}={\mathfrak {S}}'\Delta q_{1}+{\mathfrak {S}}''\Delta q_{2}+{\mathfrak {S}}'''\Delta q_{3}+{\tfrac {1}{2}}{\mathfrak {T}}''\Delta q_{2}^{2}.$ | (23) |
${\mathfrak {T}}'+{\mathfrak {T}}'''$
${\tfrac {1}{2}}{\mathfrak {T}}''.$
$\Delta q_{1}$
$\Delta q_{3}$
$\Delta q_{2},$
${\tfrac {6}{10}}$
${\tfrac {1}{2}}$
For $\Delta q_{2},$ then, we have the equation
| $-({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')=({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')\Delta q_{2}+{\tfrac {6}{10}}({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')\Delta q_{2}^{2}.$ | (24) |
$({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')$
| $({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')={\frac {1}{q_{2}}}\left(1-5{\frac {q_{2}^{2}}{r_{2}^{2}}}\right)\left(({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')+({\mathfrak {F}}_{2}{\mathfrak {S}}'''{\mathfrak {S}}')\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (25) |
| $-{\frac {B_{2}}{q_{2}r_{2}^{3}}}\left(1+{\frac {q_{2}^{2}}{r_{2}^{2}}}\right)({\mathfrak {F}}_{2}{\mathfrak {S}}'''{\mathfrak {S}}'),$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (25) |
​The quadratic equation (24) gives two values of the correction to be applied to the position of the body. When they are not too large, they will belong to two different solutions of the problem, generally to the two least removed from the values assumed. But a very large value of $\Delta q_{2}$ must not be regarded as affording any trustworthy indication of a solution of the problem. In the majority of cases we only care for one of the roots of the equation, which is distinguished by being very small, and which will be most easily calculated by a small correction to the value which we get by neglecting the quadratic term.
When a comet is somewhat near the earth we may make use of the fact that the earth's orbit is one solution of the problem, i.e., that $-\rho _{2}$ is one value of $\Delta q_{2},$ to save the trifling labor of computing the value of $({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}').$ For it is evident from the theory of equations that if $-\rho _{2}$ and $z$ are the two roots,
* $\rho _{2}-z={\frac {({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}{{\tfrac {3}{5}}({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')}}$$-\rho _{2}z={\frac {({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')}{{\tfrac {3}{5}}({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')}}\cdot$
$({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}'),$
* $(\rho _{2}-z)({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')=-\rho _{2}z({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}'''),$
* ${\frac {1}{z}}={\frac {1}{\rho _{2}}}-{\frac {({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}{({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')}}\cdot$
$-{\frac {({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}}$
$\Delta q_{2},$
$[\Delta q_{2}],$
| $\Delta q_{2}={\frac {[\Delta q_{2}]}{1+{\frac {[\Delta q_{2}]}{\rho _{2}}}}}\cdot$ | (26) |
$\Delta q_{1}$
$\Delta _{2}$
| $-({\mathfrak {S}}{\mathfrak {S}}''{\mathfrak {S}}''')=$ | $({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')\Delta q_{1}+$ | ${\tfrac {6}{10}}({\mathfrak {T}}''{\mathfrak {S}}''{\mathfrak {S}}''')\Delta q_{2}^{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (27) |
| $-({\mathfrak {S}}{\mathfrak {S}}'{\mathfrak {S}}'')=$ | $({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')\Delta q_{3}+$ | ${\tfrac {6}{10}}({\mathfrak {T}}''{\mathfrak {S}}'{\mathfrak {S}}'')\Delta q_{2}^{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (27) |
| $\Delta q_{1}=[\Delta q_{1}]+\Delta q_{2}-[\Delta q_{2}],$$\Delta q_{3}=[\Delta q_{3}]+\Delta q_{2}-[\Delta q_{2}],$ | (28) |
$[\Delta q_{1}],[\Delta q_{2}],[\Delta q_{3}]$
In making successive corrections of the distances $q_{1},q_{2},q_{3},$ it will not be necessary to recalculate the values of ${\mathfrak {S}}',{\mathfrak {S}}'',{\mathfrak {S}}''',$ when these have been calculated from fairly good values of $q_{1},q_{2},q_{3}.$ But when, as is generally the case, the first assumption is only a rude guess, the values of ${\mathfrak {S}}',{\mathfrak {S}}'',{\mathfrak {S}}'''$ should be recalculated after one or two corrections of $q_{1},q_{2},q_{3}.$ To get the best results when we do not recalculate ${\mathfrak {S}}',{\mathfrak {S}}'',{\mathfrak {S}}''',$ we may proceed as follows: Let ${\mathfrak {S}}',{\mathfrak {S}}'',{\mathfrak {S}}'''$ denote the values which have been calculated; $Dq_{1},Dq_{2},Dq_{3},$ respectively, the sum of the corrections of each of the quantities $q_{1},q_{2},q_{3},$ which have been made since the calculation of ${\mathfrak {S}}',{\mathfrak {S}}'',{\mathfrak {S}}''';\,{\mathfrak {S}}$ the residual after all the corrections of $q_{1},q_{2},q_{3},$ which have been made; and $\Delta q_{1},\Delta q_{2},\Delta q_{3}$ the remaining corrections which we are seeking. We have, then, very nearly
| $-{\mathfrak {S}}=\{{\mathfrak {S}}'+{\mathfrak {T}}'(Dq_{1}+{\tfrac {1}{2}}\Delta q_{1})\}\Delta q_{1}+\{{\mathfrak {S}}''+{\mathfrak {T}}''(Dq_{2}+{\tfrac {1}{2}}\Delta q_{2})\}\Delta q_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (29) |
| $+\{{\mathfrak {S}}'''+{\mathfrak {T}}'''(Dq_{3}+{\tfrac {1}{2}}\Delta q_{3})\}\Delta q_{3}.$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \ \end{matrix}}\right\}\,$ (29) |
The same considerations which we applied to equation (21) enable us to simplify this equation also, and to write with a fair degree of accuracy
| $-({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')=\{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')+{\tfrac {6}{5}}({\mathfrak {T}}''{\mathfrak {S}}'''{\mathfrak {S}}')(Dq_{2}+{\tfrac {1}{2}}\Delta q_{2})\}\Delta q_{2},$ | (30) |
| $\Delta q_{1}=[\Delta q_{1}]+\Delta q_{2}-[\Delta q_{2}],$$\Delta q_{3}=[\Delta q_{3}]+\Delta q_{2}-[\Delta q_{2}],$ | (31) |
| $[\Delta q_{1}]=-{\frac {({\mathfrak {S}}{\mathfrak {S}}''{\mathfrak {S}}''')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}},$$[\Delta q_{2}]=-{\frac {({\mathfrak {S}}{\mathfrak {S}}'''{\mathfrak {S}}')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}},$$[\Delta q_{3}]=-{\frac {({\mathfrak {S}}{\mathfrak {S}}'{\mathfrak {S}}'')}{({\mathfrak {S}}'{\mathfrak {S}}''{\mathfrak {S}}''')}}\cdot$ | (32) |
Correction of the Fundamental Equation.
When we have thus determined, by the numerical solution of our fundamental equation, approximate values of the three positions of the body, it will always be possible to apply a small numerical correction to the equation, so as to make it agree exactly with the laws of elliptic motion in a fictitious case differing but little from the actual. After such a correction the equation will evidently apply to the actual case with a much higher degree of approximation.
There is room for great diversity in the application of this principle. The method which appears to the writer the most simple and direct is the following, in which the correction of the intervals for aberration ​is combined with the correction required by the approximate nature of the equation.
The solution of the fundamental equation gives us three points, which must necessarily lie in one plane with the sun, and in the lines of sight of the several observations. Through these points we may pass an ellipse, and calculate the intervals of time required by the exact laws of elliptic motion for the passage of the body between them. If these calculated intervals should be identical with the given intervals, corrected for aberration, we would evidently have the true solution of the problem. But suppose, to fix our ideas, that the calculated intervals are a little too long. It is evident that if we repeat our calculations, using in our fundamental equation intervals shortened in the same ratio as the calculated intervals have come out too long, the intervals calculated from the second solution of the fundamental equation must agree almost exactly with the desired values. If necessary, this process may be repeated, and thus any required degree of accuracy may be obtained, whenever the solution of the uncorrected equation gives an approximation to the true positions. For this it is necessary that the intervals should not be too great. It appears, however, from the results of the example of Ceres, given hereafter, in which the heliocentric motion exceeds 62° but the calculated values of the intervals of time differ from the given values by little more than one part in two thousand, that we have here not approached the limit of the application of our formula.
In the usual terminology of the subject, the fundamental equation with intervals uncorrected for aberration represents the first hypothesis; the same equation with the intervals affected by certain numerical coefficients (differing little from unity) represents the second hypothesis; the third hypothesis, should such be necessary, is represented by a similar equation with corrected coefficients, etc.
In the process indicated there are certain economies of labor which should not be left unmentioned, and certain precautions to be observed in order that the neglected figures in our computations may not unduly infiuence the result.
It is evident, in the first place, that for the correction of our fundamental equation we need not trouble ourselves with the position of the orbit in the solar system. The intervals of time, which determine this correction, depend only on the three heliocentric distances $r_{1},r_{2},r_{3}$ and the two heliocentric angles, which will be represented by $v_{2}-v_{1}$ and $v_{3}-v_{2},$ if we write $v_{1},v_{2},v_{3}$ for the true anomalies. These angles ($v_{2}-v_{1}$ and $v_{3}-v_{2}$) may be determined from $r_{1},r_{2},r_{3}$ and $n_{1},n_{2},n_{3},$ ​and therefore from $r_{1},r_{2},r_{3}$ and the given intervals. For our fundamental equation, which may be written
| $n_{1}{\mathfrak {R}}_{1}-n_{2}{\mathfrak {R}}_{2}+n_{3}{\mathfrak {R}}_{3}=0,$ | (33) |
$n_{1}r_{1},n_{2}r_{2},$
$n_{3}r_{3}$
$s_{1},s_{2},s_{3}$
$s_{1}$
$s_{3}$
$v_{3}-v_{2}$
$v_{2}-v_{1}.$
| $\tan {\frac {v_{3}-v_{2}}{2}}=$ | ${\sqrt {\frac {(s_{1}-s_{2}+s_{3})(s_{1}+s_{2}-s_{3})}{(s_{1}+s_{2}+s_{3})(-s_{1}+s_{2}+s_{3})}}}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$. (34) |
| $\tan {\frac {v_{2}-v_{1}}{2}}=$ | ${\sqrt {\frac {(-s_{1}+s_{2}+s_{3})(s_{1}-s_{2}-s_{3})}{(s_{1}+s_{2}+s_{3})(s_{1}+s_{2}-s_{3})}}}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$. (34) |
As soon, therefore, as the solution of our fundamental equation has given a sufficient approximation to the values of $r_{1},r_{2},r_{3}$ (say five- or six-figure values, if our final result is to be as exact as seven-figure logarithms can make it), we calculate $n_{1},n_{2},n_{3}$ with seven-figure logarithms by equations (2), and the heliocentric angles by equations (34).
The semi-parameter corresponding to these values of the heliocentric distances and angles is given by the equation
| $p={\frac {n_{1}r_{1}-n_{2}r_{2}+n_{3}r_{3}}{n_{1}-n_{2}+n_{3}}}\cdot$ | (35) |
$n_{1}-n_{2}+n_{3},$
$n_{1}r_{1}-n_{2}r_{2}+n_{3}r_{3},$
$s_{1}-s_{2}+s_{3},$
| $n_{1}-n_{2}+n_{3}={\frac {A_{1}B_{1}}{r_{1}^{2}}}+{\frac {B_{2}}{r_{2}^{3}}}+{\frac {A_{3}B_{3}}{r_{3}^{3}}}\cdot$ | (36) |
$s_{1}-s_{2}+s_{3},$
$s_{1}-s_{2}+s_{3}$
$r_{1},r_{2},r_{3}$
$r_{2}$
After the determination of the heliocentric angles and the semi-parameter, the eccentricity and the true anomalies of the three positions may next be determined, and from these the intervals of time. These processes require no especial notice. The appropriate formulæ will be given in the Summary of Formulæ.
Determination of the Orbit from the Three Positions and the Intervals of Time.
The values of the semi-parameter and the heliocentric angles as given in the preceding paragraphs depend upon the quantity $s_{1}-s_{2}+s_{3},$ the numerical determination of which from $s_{1},s_{2},$ and $s_{3},$ a critical to the second degree when the heliocentric angles are small. This was of no consequence in the process which we have called the correction of the fundamental equation. But for the actual determination of the orbit from the positions given by the corrected equation— or by the uncorrected equation, when we judge that to be sufficient—a more accurate determination of this quantity will generally be necessary. This may be obtained in different ways, of which the following is perhaps the most simple. Let us set
| ${\mathfrak {S}}_{4}={\mathfrak {S}}_{3}-{\mathfrak {S}}_{1},$ | (37) |
$s_{4}$
${\mathfrak {S}}_{4},$
$s_{2}$
$s_{4}$
$s_{1}$
$s_{3}.$
$s_{1},s_{2},s_{3}$
$s_{1},s_{3},s_{4},$
| ${\begin{aligned}(s_{1}+s_{2}+s_{3})(-s_{1}+s_{2}+s_{3})(s_{1}-s_{2}+s_{3})(s_{1}+s_{2}-s_{3})\\=(s_{1}+s_{4}+s_{3})(-s_{1}+s_{4}+s_{3})(s_{1}-s_{4}+s_{3})(s_{1}+s_{4}-s_{3}),\end{aligned}}$ | (38) |
| $s_{1}-s_{2}+s_{3}={\frac {(s_{1}+s_{4}+s_{3})(-s_{1}+s_{4}+s_{3})(s_{1}-s_{4}+s_{3})(s_{1}+s_{4}-s_{3})}{(s_{1}+s_{2}+s_{3})(-s_{1}+s_{2}+s_{3})(s_{1}+s_{2}-s_{3})}}\cdot$ | (39) |
$s_{1}-s_{2}+s_{3}$
The eccentricity and the true anomalies may be determined in the same way as in the correction of the formula. The position of the orbit in space may be derived from the following considerations. The vector $-{\mathfrak {S}}_{2}$ is directed from the sun toward the second position of the body; the vector ${\mathfrak {S}}_{4}$ from the first to the third position. If we set
| ${\mathfrak {S}}_{5}={\mathfrak {S}}_{4}-{\frac {{\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}}{s_{2}^{2}}}{\mathfrak {S}}_{2},$ | (40) |
${\mathfrak {S}}_{5}$
$-{\mathfrak {S}}_{2}$
$s_{5}$
${\mathfrak {S}}_{5},$
* $-{\frac {{\mathfrak {S}}_{2}}{s_{2}}}$${\frac {{\mathfrak {S}}_{5}}{s_{5}}}$
${\mathfrak {J}}$
${\mathfrak {J}}'$
${\mathfrak {J}}$
${\mathfrak {J}}'$
${\mathfrak {J}},$
| ${\mathfrak {J}}=-\cos v_{2}{\frac {{\mathfrak {S}}_{2}}{s_{2}}}-\sin v_{2}{\frac {{\mathfrak {S}}_{5}}{s_{5}}},$ | (41) |
| ${\mathfrak {J}}'=-\sin v_{2}{\frac {{\mathfrak {S}}_{2}}{s_{2}}}+\cos v_{2}{\frac {{\mathfrak {S}}_{5}}{s_{5}}}\cdot$ | (42) |
${\text{T}}$
| ${\frac {k}{a^{\frac {3}{2}}}}(t-T)=E-e\sin E,$ | (43) |
$E$
$v.$
$t-T$
$T$
$T$
For the calculation of an ephemeris we have
| ${\mathfrak {R}}=-ae{\mathfrak {J}}+\cos E\,a{\mathfrak {J}}+\sin E\,b{\mathfrak {J}}'$ | (44) |
Sometimes it may be worth while to make the calculations for the correction of the formula in the slightly longer form indicated for ​the determination of the orbit. This will be the case when we wish simultaneously to correct the formula for its theoretical imperfection, and to correct the observations by comparison with others not too remote. The rough approximation to the orbit given by the uncorrected formula may be sufficient for this purpose. In fact, for observations separated by very small intervals, the imperfection of the uncorrected formula will be likely to affect the orbit less than the errors of the observations.
The computer may prefer to determine the orbit from the first and third heliocentric positions with their times. This process, which has certain advantages, is perhaps a little longer than that here given, and does not lend itself quite so readily to successive improvements of the hypothesis. When it is desired to derive an improved hypothesis from an orbit thus determined, the formulæ in § XII of the summary may be used.
SUMMARY OF FORMULÆ
WITH DIRECTIONS FOR USE.
(For the case in which an approximate orbit is known in advance, see XII.)
I.
Preliminary computations relating to the intervals of time.
$t_{1},t_{2},t_{3}=$ times of the observations in days,
$\log k=8.2355814$ (after Gauss)
$\tau _{1}=k(t_{3}-t_{2})$ $\tau _{3}=k(t_{2}-t_{1})$
$A_{1}={\frac {t_{3}-t_{2}}{t_{3}-t_{1}}}$$A_{3}={\frac {t_{2}-t_{1}}{t_{3}-t_{1}}}$
$B_{1}={\frac {\tau _{1}^{2}+\tau _{1}\tau _{3}+\tau _{3}^{2}}{12}}$$B_{2}={\frac {\tau _{1}^{2}+3\tau _{1}\tau _{3}+\tau _{3}^{2}}{12}}$$B_{3}={\frac {\tau _{1}^{2}+\tau _{1}\tau _{3}-\tau _{3}^{2}}{12}}$
For control: $A_{1}B_{1}+B_{2}+A_{3}B_{3}={\tfrac {1}{2}}\tau _{1}\tau _{3}.$
II.
Preliminary computations relating to the first observation.
$X_{1},Y_{1},Z_{1}$ (components of ${\mathfrak {E}}_{1}$) $=$ the heliocentric coordinates of the earth, increased by the geocentric coordinates of the observatory.
$\xi _{1},\eta _{1},\zeta _{1}$ (components of ${\mathfrak {F}}_{1}$) $=$ the direction-cosines of the observed position, corrected for the aberration of the fixed stars.
${\mathfrak {E}}_{1}^{2}=X_{1}^{2}+Y_{1}^{2}+Z_{1}^{2}$$({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})=X_{1}\xi _{1}+Y_{1}\eta _{1}+Z_{1}\zeta _{1}$$p_{1}^{2}={\mathfrak {E}}_{1}^{2}-({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})^{2}$
Preliminary computations relating to the second and third observations.
The formulæ are entirely analogous to those relating to the first observation, the quantities being distinguished by the proper suffixes.
III
Equations of the first hypothesis.
When the preceding quantities have been computed, their numerical values (or their logarithms, when more convenient for computation,) are to be substituted in the following equations:
Components of ${\mathfrak {S}}_{1}$ :
| $q_{1}=\rho _{1}+({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})$ | $\alpha _{1}=A_{1}\xi _{1}(1+R_{1})\left(q_{1}+{\frac {X_{1}}{\xi _{1}}}-({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
| $r_{1}^{2}=q_{1}^{2}+p_{1}^{2}$ | $\beta _{1}=A_{1}\eta _{1}(1+R_{1})\left(q_{1}+{\frac {Y_{1}}{\eta _{1}}}-({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
| $R_{1}={\frac {B_{1}}{r_{1}^{3}}}$ | $\gamma _{1}=A_{1}\zeta _{1}(1+R_{1})\left(q_{1}+{\frac {Z_{1}}{\zeta _{1}}}-({\mathfrak {E}}_{1}.{\mathfrak {F}}_{1})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
$s_{1}^{2}=\alpha _{1}^{2}+\beta _{1}^{2}+\gamma _{1}^{2}=A_{1}^{2}(1+R_{1})^{2}r_{1}^{2}$
Components of ${\mathfrak {S}}'$ :
| $P'={\frac {3R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\alpha '=A_{1}\xi _{1}+A_{1}\xi _{1}R_{1}-P'\alpha _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'$ |
| $P'={\frac {3R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\beta '=A_{1}\eta _{1}+A_{1}\eta _{1}R_{1}-P'\beta _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'$ |
| $P'={\frac {3R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\gamma '=A_{1}\zeta _{1}+A_{1}\zeta _{1}R_{1}-P'\gamma _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'$ |
Components of ${\mathfrak {S}}_{2}$ :
| $q_{2}=\rho _{2}+({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})$ | $\alpha _{2}=-\xi _{2}(1-R_{2})\left(q_{2}+{\frac {X_{2}}{\xi _{2}}}-({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
| $r_{2}^{2}=q_{2}^{2}+p_{2}^{2}$ | $\beta _{2}=-\eta _{2}(1-R_{2})\left(q_{2}+{\frac {Y_{2}}{\eta _{2}}}-({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
| $R_{2}={\frac {B_{2}}{r_{2}^{3}}}$ | $\gamma _{2}=-\zeta _{2}(1+R_{2})\left(q_{2}+{\frac {Z_{2}}{\zeta _{2}}}-({\mathfrak {E}}_{2}.{\mathfrak {F}}_{2})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
$s_{2}^{2}=\alpha _{2}^{2}+\beta _{2}^{2}+\gamma _{2}^{2}=(1-R_{2})^{2}r_{2}^{2}$
Components of ${\mathfrak {S}}''$ :
| $P''={\frac {3R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\alpha ''=-\xi _{2}+\xi _{2}R_{2}+P''\alpha _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}''$ |
| $P''={\frac {3R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\beta ''=-\eta _{2}+A_{1}\eta _{2}R_{2}+P''\beta _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}''$ |
| $P''={\frac {3R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\gamma ''=-\zeta _{2}+\zeta _{2}R_{2}+P''\gamma _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}''$ |
Components of ${\mathfrak {S}}_{3}$ :
| $q_{3}=\rho _{3}+({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})$ | $\alpha _{2}=A_{3}\xi _{3}(1+R_{3})\left(q_{3}+{\frac {X_{3}}{\xi _{3}}}-({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
| $r_{3}^{2}=q_{3}^{2}+p_{3}^{2}$ | $\beta _{3}=A_{3}\eta _{3}(1+R_{3})\left(q_{3}+{\frac {Y_{3}}{\eta _{3}}}-({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
| $R_{3}={\frac {B_{3}}{r_{3}^{3}}}$ | $\gamma _{3}=A_{3}\zeta _{3}(1+R_{3})\left(q_{3}+{\frac {Z_{3}}{\zeta _{3}}}-({\mathfrak {E}}_{3}.{\mathfrak {F}}_{3})\right)$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
$s_{3}^{2}=\alpha _{3}^{2}+\beta _{3}^{2}+\gamma _{3}^{2}=A_{3}(1+R_{3})^{2}r_{3}^{2}$
Components of ${\mathfrak {S}}'''$ :
| $P'''={\frac {3R_{3}q_{3}}{(1+R_{3})r_{3}^{2}}}$ | $\alpha '''=A_{3}\xi _{3}+A_{3}\xi _{3}R_{3}-P'''\alpha _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'''$ |
| $P'''={\frac {3R_{3}q_{3}}{(1+R_{3})r_{3}^{2}}}$ | $\beta '''=A_{3}\eta _{3}+A_{3}\eta _R_{3}-P'''\beta _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'''$ |
| $P'''={\frac {3R_{3}q_{3}}{(1+R_{3})r_{3}^{2}}}$ | $\gamma '''=A_{3}\zeta _{3}+A_{3}\zeta _{3}R_{3}-P'''\gamma _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}'''$ |
The computer is now to assume any reasonable values either of the geocentric distances, $\rho _{1},\rho _{2},\rho _{3},$ or of the heliocentric distances, $r_{1},r_{2},r_{3}$ (the former in the case of a comet, the latter in the case of an asteroid), and from these assumed values to compute the rest of the following quantities:
| By equations ${\text{III}}_{1},{\text{III}}'.$ | By equations ${\text{III}}_{2},{\text{III}}''.$ | By equations ${\text{III}}_{3},{\text{III}}'''.$ |
| $q_{1}$ | $q_{2}$ | $q_{3}$ |
| $\log r_{1}$ | $\log r_{2}$ | $\log r_{3}$ |
| $\log R_{1}$ | $\log R_{2}$ | $\log R_{3}$ |
| $\log(1+R_{1})$ | $\log(1+R_{2})$ | $\log(1+R_{3})$ |
| $\log P'$ | $\log P''$ | $\log P'''$ |
| $\alpha _{1}$ | $\alpha _{2}$ | $\alpha _{3}$ |
| $\beta _{1}$ | $\beta _{2}$ | $\beta _{3}$ |
| $\gamma _{1}$ | $\gamma _{2}$ | $\gamma _{3}$ |
| $\alpha '$ | $\alpha ''$ | $\alpha '''$ |
| $\beta '$ | $\beta ''$ | $\beta '''$ |
| $\gamma '$ | $\gamma ''$ | $\gamma '''$ |
IV.
Calculations relating to differential coefficients.
| Components of ${\mathfrak {S}}''\times {\mathfrak {S}}'''$ | Components of ${\mathfrak {S}}'''\times {\mathfrak {S}}'$ | Components of ${\mathfrak {S}}'\times {\mathfrak {S}}''$ |
| $a_{1}=\beta ''\gamma '''-\gamma ''\beta '''$ | $a_{1}=\beta '''\gamma '-\gamma '''\beta '$ | $a_{3}=\beta '\gamma ''-\gamma '\beta ''$ |
| $b_{1}=\gamma ''\alpha '''-\alpha ''\gamma '''$ | $b_{2}=\gamma '''\alpha '-\alpha '''\gamma '$ | $b_{3}=\gamma '\alpha ''-\alpha '\gamma ''$ |
| $c_{1}=\alpha ''\beta '''-\beta ''\alpha '''$ | $c_{2}=\alpha '''\beta '-\beta '''\alpha '$ | $c_{3}=\alpha '\beta ''-\beta '\alpha ''$ |
These computations are controlled by the agreement of the three values of $G.$
The following are not necessary except when the corrections to be made are large:
* ${\begin{aligned}H&=({\mathfrak {F}}_{2}{\mathfrak {S}}'''{\mathfrak {S}}')=a_{2}\xi _{2}+b_{2}\eta _{2}+c_{2}\zeta _{2}\\L&={\frac {1}{q_{2}}}\left(1+{\frac {H}{G}}\right)\left(1-5{\frac {q_{2}^{2}}{r_{2}^{2}}}\right)-{\frac {R_{2}H}{q_{2}G}}\left(1+{\frac {q_{2}^{2}}{r_{2}^{2}}}\right)\end{aligned}}$
V.
Corrections of the geocentric distances.
| Components of ${\mathfrak {S}}.$ |
| $\alpha =\alpha _{1}+\alpha _{2}+\alpha _{3}$ | $C_{1}=-{\frac {a_{1}\alpha +b_{1}\beta +c_{1}\gamma }{G}}$ |
| $\beta =\beta _{1}+\beta _{2}+\beta _{3}$ | $C_{2}=-{\frac {a_{2}\alpha +b_{2}\beta +c_{2}\gamma }{G}}$ |
| $\gamma =\gamma _{1}+\gamma _{2}+\gamma _{3}$ | $C_{3}=-{\frac {a_{3}\alpha +b_{3}\beta +c_{3}\gamma }{G}}$ |
| $\Delta q_{2}=C_{2}-{\tfrac {6}{10}}L(\Delta q_{2})^{2}.$ |
(This equation will generally be most easily solved by repeated substitutions.)
* $\Delta q_{1}=C_{1}-{\tfrac {6}{10}}L(\Delta q_{2})^{2}$$\Delta q_{3}=C_{3}-{\tfrac {6}{10}}L(\Delta q_{2})^{2}.$
VI.
Successive corrections.
$\Delta q_{1},\Delta q_{2},\Delta q_{3}$ are to be added as corrections to $q_{1},q_{2},q_{3}.$ With the new values thus obtained the computation by equations ${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3}$ are to be recommenced. Two courses are now open:
(a) The work may be carried on exactly as before to the determination of new corrections for $q_{1},q_{2},q_{3}.$
(b) The computations by equations ${\text{III}}',{\text{III}}'',{\text{III}}''',$ and ${\text{IV}}$ may be omitted, and the old values of $a_{1},b_{1},c_{1},a_{2},$ etc., $G,$ and $L$ may be used with the new residuals $\alpha ,\beta ,\gamma$ to get new corrections for $q_{1},q_{2},q_{3}$ by the equations
* $\Delta q_{2}={\frac {C_{2}}{1+{\tfrac {6}{5}}L(Dq_{2}+{\tfrac {1}{2}}C_{2})}},$
* $\Delta q_{1}=C_{1}+\Delta q_{2}-C_{2},$$\Delta q_{3}=C_{3}+\Delta q_{2}-C_{2},$
$Dq_{2}$
$q_{2}.$
$Dq_{2}$
$q_{2}$
$a_{1},b_{1},$
$a_{1},b_{1},$
$G$
$q_{1},q_{2},q_{3},$
$L$
$q_{2},$
$L$
The corrections of $q_{1},q_{2},q_{3}$ should be repeated until the equations
* $\alpha =0$$\beta =0$$\gamma =0$
$r_{1},r_{2},r_{3}$
VII.
Test of the first hypothesis.
$\log r_{1},\log r_{2},\log r_{3}$ (approximate values from the preceding computations).
* ${\begin{aligned}N&=A_{1}B_{1}r_{1}^{-3}+B_{2}r_{2}^{-3}+B_{2}r_{2}^{-3}+A_{3}B_{3}r_{3}^{-3}\\s_{1}&=A_{1}r_{1}+A_{1}B_{1}r_{1}^{-2}\\s_{2}&=r_{2}-B_{2}r_{2}^{-2}\\s_{3}&=A_{3}r_{3}+A_{3}B_{3}r_{3}^{-2}\\s={\tfrac {1}{2}}(s_{1}+s_{2}+s_{3})\\s&-s_{1},s-s_,s-s_{3}.\end{aligned}}$
$s-s_{2}$
| $_{\text{R}}$ | $={\sqrt {\frac {(s-s_{1})(s-s_{2})(s-s_{3})}{s}}}$ | $\tan {\tfrac {1}{2}}(v_{2}-v_{1})$ | = ${\frac {_{\text{R}}}{s-s_{2}}}$ |
| $p$ | $={\frac {2(s-s_{2})}{N}}$ | $\tan {\tfrac {1}{2}}(v_{3}-v_{2})$ | $={\frac {_{\text{R}}}{s-s_{1}}}$ |
| | | $\tan {\tfrac {1}{2}}(v_{3}-v_{1})$ | $={\frac {s-s_{2}}{_{\text{R}}}}$ |
For adjustment of values: ${\tfrac {1}{2}}(v_{3}-v_{1})={\tfrac {1}{2}}(v_{2}-v_{1})+{\tfrac {1}{2}}(v_{3}-v_{2})$}}
* ${\begin{aligned}e\sin {\tfrac {1}{2}}(v_{3}+v_{1})&={\frac {{\frac {p}{r_{1}}}-{\frac {p}{r_{3}}}}{2\sin {\tfrac {1}{2}}(v_{3}-v_{1})}}\\e\cos {\tfrac {1}{2}}(v_{3}+v_{1})&={\frac {{\frac {p}{r_{1}}}+{\frac {p}{r_{3}}}-2}{2\cos {\tfrac {1}{2}}(v_{3}-v_{1})}}\\\tan {\tfrac {1}{2}}(v_{3}+v_{1})&e^{2}\end{aligned}}$
$e\cos v_{2}={\frac {p}{r_{2}}}-1$
* $\epsilon ={\sqrt {\frac {1-e}{1+e}}}$$a={\frac {p}{1-e^{2}}}$
* $\tan {\tfrac {1}{2}}E_{1}=\epsilon \tan {\tfrac {1}{2}}v_{1}$$\tan {\tfrac {1}{2}}E_{2}=\epsilon \tan {\tfrac {1}{2}}v_{2}$$\tan {\tfrac {1}{2}}E_{3}=\epsilon \tan {\tfrac {1}{2}}v_{3}$
* ${\begin{aligned}\tau _{1{\text{ calc.}}}&=a^{\frac {3}{2}}(E_{3}-E_{2})+ea^{\frac {3}{2}}\sin E_{2}-ea^{\frac {3}{2}}\sin E_{3}\\\tau _{2{\text{ calc.}}}&=a^{\frac {3}{2}}(E_{2}-E_{1})+ea^{\frac {3}{2}}\sin E_{1}-ea^{\frac {3}{2}}\sin E_{2}\end{aligned}}$
VIII.
For the second hypothesis.
| $\delta \tau _{1}$ | $=.0057613k(\rho _{2}-\rho _{3})$ | (aberration-constant after Struve.) |
| $\delta \tau _{2}$ | $=.0057613k(\rho _{1}-\rho _{2})$ | $\log(.0057613k)=5.99610$ |
| $\Delta \log \tau _{1}$ | $=\log \tau _{1}-\log(\tau _{1{\text{ calc.}}}-\delta \tau _{1})$ |
| $\Delta \log \tau _{3}$ | $=\log \tau _{3}-\log(\tau _{3{\text{ calc.}}}-\delta \tau _{3})$ |
| $\Delta \log(\tau _{1}\tau _{3})$ | $=\Delta \log \tau _{1}+\Delta \log \tau _{3}$ |
| $\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ | $=\Delta \log \tau _{1}-\Delta \log \tau _{3}$ |
| $\Delta \log A_{1}$ | $=-A_{3}\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ |
| $\Delta \log A_{3}$ | $=-A_{1}\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ |
| $\Delta \log B_{1}$ | $=\Delta \log(\tau _{1}\tau _{3})-{\frac {\tau _{1}^{2}+\tau _{3}^{2}}{12B_{1}}}\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ |
| $\Delta \log B_{2}$ | $=\Delta \log(\tau _{1}\tau _{3})+{\frac {\tau _{1}^{2}-\tau _{3}^{2}}{12B_{2}}}\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ |
| $\Delta \log B_{3}$ | $=\Delta \log(\tau _{1}\tau _{3})+{\frac {\tau _{1}^{2}+\tau _{3}^{2}}{12B_{3}}}\Delta \log {\frac {\tau _{1}}{\tau _{3}}}$ |
These corrections are to be added to the logarithms of $A_{1},A_{3},B_{1},B_{2},B_{3}$ in equations ${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3}$ and the corrected equations used to correct the values of $q_{1},q_{2},q_{3}$ until the residuals $\alpha ,\beta ,\gamma$ vanish. The new values of $A_{1},A_{3}$ must satisfy the relation $A_{1}+A_{3}=1,$ and the corrections $\Delta \log A_{1},\Delta \log A_{3}$ must be adjusted, if necessary, for this end.
Third hypothesis.
A second correction of equations ${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3}$ may be obtained in the same manner as the first, but this will rarely be necessary.
IX.
Determination of the ellipse.
It is supposed that the values of
| $\alpha _{1},\beta _{1},\gamma _{1},$ | $\alpha _{2},\beta _{2},\gamma _{2},$ | $\alpha _{3},\beta _{3},\gamma _{3},$ |
| $r_{1},r_{2},r_{3},$ | $R_{1},R_{2},R_{3},$ | $s_{1},s_{2},s_{3},$ |
${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3}$
$\alpha ,\beta ,\gamma$
$s_{1},s_{2},s_{3}$
| Components of ${\mathfrak {S}}_{4}$ | Components of ${\mathfrak {S}}_{4}$ | Components of ${\mathfrak {S}}_{5}$ | Components of ${\mathfrak {S}}_{5}$ |
| $\alpha _{4}$ | $=\alpha _{3}-\alpha _{1}$ | $\alpha _{5}$ | $=\alpha _{4}-{\frac {\alpha _{4}\alpha _{2}+\beta _{4}\beta _{2}+\gamma _{4}\gamma _{2}}{s_{2}^{2}}}\alpha _{2}$ |
| $\beta _{4}$ | $=\beta _{3}-\beta _{1}$ | $\beta _{5}$ | $=\beta _{4}-{\frac {\alpha _{4}\alpha _{2}+\beta _{4}\beta _{2}+\gamma _{4}\gamma _{2}}{s_{2}^{2}}}\beta _{2}$ |
| $\gamma _{4}$ | $=\gamma _{3}-\gamma _{1}$ | $\gamma _{5}$ | $=\gamma _{4}-{\frac {\alpha _{4}\alpha _{2}+\beta _{4}\beta _{2}+\gamma _{4}\gamma _{2}}{s_{2}^{2}}}\gamma _{2}$ |
| $s_{4}^{2}$ | $=\alpha _{4}^{2}+\beta _{4}^{2}+\gamma _{4}^{2}$ | $s_{5}^{2}$ | $=\alpha _{5}^{2}+\beta _{5}^{2}+\gamma _{5}^{2}$ |
| $s$ | $={\tfrac {1}{2}}(s_{1}+s_{2}+s_{3})$ | $S$ | $={\tfrac {1}{2}}(s_{1}+s_{4}+s_{3})$ |
For control only:
* $s-s_{2}={\frac {S(S-s_{1})(S-s_{4})(S-s_{3})}{s(s-s_{1})(s-s_{3})}}$
| $_{{\text{R}}^{2}}$ | $={\frac {S(S-s_{1})(S-s_{4})(S-s_{3})}{s_{2}^{2}}}$ | $\tan {\tfrac {1}{2}}(v_{2}-v_{1})$ | $={\frac {_{\text{R}}}{s-s_{3}}}$ |
| $N$ | $=A_{1}R_{1}+R_{2}+A_{3}R_{3}$ | $\tan {\tfrac {1}{2}}(v_{3}-v_{2})$ | $={\frac {_{\text{R}}}{s-s_{1}}}$ |
| $p$ | $={\frac {2_{{\text{R}}^{2}}s}{N(s-s_{1})(s-s_{3})}}$ | $\tan {\tfrac {1}{2}}(v_{3}-v_{1})$ | $={\frac {_{\text{R}}s}{(s-s_{1})(s-s_{3})}}$ |
The computer should be careful to use the corrected values of $A_{1},A_{3}.$ (See VIII.) Trifling errors in the angles should be distributed.
* ${\begin{aligned}e\sin {\tfrac {1}{2}}(v_{3}+v_{1})&={\frac {{\frac {p}{r_{1}}}-{\frac {p}{r_{3}}}}{2\sin {\tfrac {1}{2}}(v_{3}-v_{1})}}\\e\cos {\tfrac {1}{2}}(v_{3}+v_{1})&={\frac {{\frac {p}{r_{1}}}+{\frac {p}{r_{3}}}-2}{2\cos {\tfrac {1}{2}}(v_{3}-v_{1})}}\\\tan {\tfrac {1}{2}}(v_{3}+v_{1})&e^{2}\end{aligned}}$
$e\cos v_{2}={\frac {p}{r_{2}}}-1$
* $\epsilon ={\sqrt {\frac {1-e}{1+e}}}$$a={\frac {p}{1-e^{2}}}$$b={\sqrt {a/p}}$
Direction-cosines of semi-major axis. :
| $l$ | $=-{\frac {\cos v_{2}}{s_{2}}}\alpha _{2}-{\frac {\sin v_{2}}{s_{5}}}\alpha _{5}$ |
| $m$ | $=-{\frac {\cos v_{2}}{s_{2}}}\beta _{2}-{\frac {\sin v_{2}}{s_{5}}}\beta _{5}$ |
| $n$ | $=-{\frac {\cos v_{2}}{s_{2}}}\gamma _{2}-{\frac {\sin v_{2}}{s_{5}}}\gamma _{5}$ |
Direction-cosines of semi-minor axis. :
| $\lambda$ | $=-{\frac {\sin v_{2}}{s_{2}}}\alpha _{2}+{\frac {\cos v_{2}}{s_{5}}}\alpha _{5}$ |
| $\mu$ | $=-{\frac {\sin v_{2}}{s_{2}}}\beta _{2}+{\frac {\cos v_{2}}{s_{5}}}\beta _{5}$ |
| $\nu$ | $=-{\frac {\sin v_{2}}{s_{2}}}\gamma _{2}+{\frac {\cos v_{2}}{s_{5}}}\gamma _{5}$ |
Components of the semi-axes. :
| $a_{x}=al$ | $a_{y}=am$ | $a_{z}=an$ |
| $b_{x}=b\lambda$ | $b_{y}=b\mu$ | $b_{z}=b\nu$ |
X.
Time of perihelion passage.
| | Corrections for aberration. |
| $\tan {\tfrac {1}{2}}E_{1}=\epsilon \tan {\tfrac {1}{2}}v_{1}$ | $\delta t_{1}=-.0057613\rho _{1}$ |
| $\tan {\tfrac {1}{2}}E_{2}=\epsilon \tan {\tfrac {1}{2}}v_{2}$ | $\delta t_{2}=-.0057613\rho _{2}$ |
| $\tan {\tfrac {1}{2}}E_{3}=\epsilon \tan {\tfrac {1}{2}}v_{3}$ | $\delta t_{3}=-.0057613\rho _{3}$ |
| $\log .0057613=7.76052$ ||
| $t_{1}+\delta t_{1}-T=k^{-1}a^{\frac {3}{2}}(E_{1}-e\sin E_{1})$ ||
| $t_{2}+\delta t_{2}-T=k^{-1}a^{\frac {3}{2}}(E_{2}-e\sin E_{2})$ ||
| $t_{3}+\delta t_{3}-T=k^{-1}a^{\frac {3}{2}}(E_{3}-e\sin E_{3})$ ||
The threefold determination of $T$ affords a control of the exactness of the solution of the problem. If the discrepancies in the values of $T$ are such as to require another correction of the formulæ (a third hypothesis), this may be based on the equations
* $\Delta \log \tau _{1}=M{\frac {T_{(3)}-T_{(2)}}{t_{3}-t_{2}}}$$\Delta \log \tau _{2}=M{\frac {T_{(2)}-T_{(1)}}{t_{2}-t_{1}}}$
$T_{(1)},T_{(2)},T_{(3)}$
$M$
XI.
For an ephemeris.
${\frac {k}{a^{\frac {3}{2}}}}(t-T)=E-e\sin E$ :
* Heliocentric coordinates. (Components of ${\mathfrak {R}}$.) |
* $x=-ea_{x}+a_{x}\cos E+b_{x}\sin E$
* $y=-ea_{y}+a_{y}\cos E+b_{y}\sin E$
* $z=-ea_{z}+a_{z}\cos E+b_{z}\sin E$
These equations are completely controlled by the agreement of the computed and observed positions and the following relations between the constants:
* $a_{x}b_{x}+a_{y}b_{y}+a_{z}b_{z}=0$$a_{x}^{2}+a_{y}^{2}+a_{z}^{2}=a^{2}$$b_{x}^{2}+b_{y}^{2}+b_{z}^{2}=(1-e)^{2}a^{2}$
XII.
When an approximate orbit is known in advance, we may use it to improve our fundamental equation. The following appears to be the most simple method:
Find the excentric anomalies $E_{1},E_{2},E_{3},$ and the heliocentric distances $r_{1},r_{2},r_{3},$ which belong in the approximate orbit to the times of observation corrected for aberration.
Calculate $B_{1},B_{3},$ as in § I, using these corrected times.
Determine $A_{1},A_{3}$ by the equation
* ${\frac {A_{1}\left(1+{\frac {B_{1}}{r_{1}^{3}}}\right)}{\sin(E_{3}-E_{2})-e\sin E_{3}+e\sin E_{2}}}={\frac {A_{3}\left(1+{\frac {B_{3}}{r_{3}^{3}}}\right)}{\sin(E_{2}-E_{1})-e\sin E_{2}+e\sin E_{1}}}$
$A_{1}+A_{3}=1.$
Determine $B_{2}$ so as to make
* ${\frac {A_{1}{\frac {B_{1}}{r_{1}^{3}}}+{\frac {B_{2}}{r_{2}^{3}}}+A_{3}{\frac {B_{3}}{r_{3}^{3}}}}{4\sin {\tfrac {1}{2}}(E_{2}-E_{1})\sin {\tfrac {1}{2}}(E_{3}-E_{2})\sin {\tfrac {1}{2}}(E_{3}-E_{1})}}$
It is not necessary that the times for which $E_{1},E_{2},E_{3},r_{1},r_{2},r_{3},$ are calculated should precisely agree with the times of observation corrected for aberration. Let the former be represented by $t_{1}',t_{2}',t_{3}'$ and the latter by $t_{1}'',t_{2}'',t_{3}'';$ and let
* ${\begin{aligned}\Delta \log \tau _{1}&=\log(t_{3}''-t_{2}'')-\log(t_{3}'-t_{2}'),\\\Delta \log \tau _{2}&=\log(t_{2}''-t_{1}'')-\log(t_{2}'-t_{1}').\end{aligned}}$
$B_{1},B_{3},A_{1},A_{3},B_{2},$
$t_{1}',t_{2}',t_{3}',$
$\Delta \log \tau _{1},\Delta \log \tau _{2}$
Numerical Example.
To illustrate the numerical computations we have chosen the following example, both on account of the large heliocentric motion, and because Gauss and Oppolzer have treated the same data by their different methods.
The data are taken from the Theoria Motus, § 169, viz.,
| Times, 1805, September | 5.51336 | 139.42711 | 265.39813 |
| Longitudes of Ceres | 95° 32' 18''.56 | 99° 49' 5''.87 | 118° 5' 28''.85 |
| Latitudes of Ceres | -0° 59' 34''.06 | +7° 16' 36''.80 | +7° 38' 49''.39 |
| Longitudes of the Earth | 342° 54' 56''.00 | 117° 12' 43''.25 | 241° 58' 50''.71 |
| Logs of the Sun's distance | 0.0031514 | 9.9929861 | 0.0056974 |
I.
From the given times we obtain the following values:
| | Numbers. | Logarithms. |
| $t_{2}-t_{1}$ | 133.91375 | 2.1268252 |
| $t_{3}-t_{2}$ | 125.97102 | 2.1002706 |
| $t_{3}-t_{1}$ | 259.88477 | 2.4147809 |
| $A_{1}$ | .4847187 | 9.6854897 |
| $A_{3}$ | .5152812 | 9.7120443 |
| $\tau _{1}$ | | .3358520 |
| $\tau _{3}$ | | .3624066 |
| $B_{1}$ | | 9.6692113 |
| $B_{2}$ | | .3183722 |
| $B_{3}$ | | 9.5623916 |
Control:
* ${\begin{aligned}A_{1}B_{1}+B_{2}&+A_{3}B_{3}&=2.4959086\\&{\tfrac {1}{2}}\tau _{1}\tau _{3}&=2.4959081\end{aligned}}$
II.
From the given positions we get:
| $\log X_{1}$ | 9.9835515 | $+$ | $\log X_{2}$ | 9.6531725 | $-$ | $\log X_{3}$ | 9.6775810 | $-$ |
| $\log Y_{1}$ | 9.4711748 | $-$ | $\log Y_{2}$ | 9.9420444 | $+$ | $\log Y_{3}$ | 9.515547 | $-$ |
| $Z_{1}$ | 0 | | $Z_{2}$ | 0 | | $Z_{3}$ | 0 | |
| $\log \xi _{1}$ | 8.9845270 | $-$ | $\log \xi _{2}$ | 9.2282738 | $-$ | $\log \xi _{3}$ | 9.6690294 | $-$ |
| $\log \eta _{1}$ | 9.9979027 | $+$ | $\log \eta _{2}$ | 9.9900800 | $+$ | $\log \eta _{3}$ | 9.9416855 | $+$ |
| $\log \zeta _{1}$ | 8.2387150 | $-$ | $\log \zeta _{2}$ | 9.1026549 | $+$ | $\log \zeta _{3}$ | 9.1240813 | $+$ |
| ${\mathfrak {E}}_{1}.{\mathfrak {F}}_{1}$ | .3874081 | $-$ | ${\mathfrak {E}}_{2}.{\mathfrak {F}}_{2}$ | .9314223 | $+$ | ${\mathfrak {E}}_{3}.{\mathfrak {F}}_{3}$ | .5599304 | $-$ |
| $p_{1}^{2}$ | .8645336 | $+$ | $p_{2}^{2}$ | .10006681 | $+$ | $p_{3}^{2}$ | .7130624 | $+$ |
III.
The preceding computations furnish the numerical values for the equations ${\text{III}}_{1},{\text{III'}},{\text{III}}_{2},{\text{III}}'',{\text{III}}_{3},{\text{III}}''',$ which follow. Brackets indicate that logarithms have been substituted for numbers.
$r_{1},r_{2},r_{3}.$
$r_{1},r_{2},r_{3}.$
| $q_{1}$ | $=\rho _{1}-.3874081$ | $\alpha _{1}$ | $=-[8.6700167](q_{1}-9.5901555)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
| $r_{1}^{2}$ | $=q_{1}^{2}+.8645336$ | $\beta _{1}$ | $=[9.6833924](q_{1}+.090052)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
| $R_{1}$ | $=[9.6692113]r_{1}^{-3}$ | $\gamma _{1}$ | $=-[7.9242047](q_{1}+.3874081)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ |
| $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\alpha '$ | $=-.46775-[8.67002]R_{1}-P'\alpha _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'$ |
| $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\beta '$ | $=.482383+[9.68339]R_{1}-P'\beta _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'$ |
| $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $P'={\frac {[.47712]R_{1}q_{1}}{(1+R_{1})r_{1}^{2}}}$ | $\gamma '$ | $=-.008399-[7.92420]R_{1}-P'\gamma _{1}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'$ |
| $\Delta q_{1}$ | | | –.66731 | –.04558 | –.0010434 | + .0000006 |
| $q_{1}$ | $+$ | 3.22606 | 2.55875 | 2.51317 | 2.5142134 | 2.5142140 |
| $\log r_{1}$ | $+$ | .52600 | .434960 | .4280791 | .4282376 | .4282377 |
| $\log R_{1}$ | $+$ | 9.09121 | 8.364331 | 8.3849740 | 8.3844985 | |
| $\log(1+R_{1})$ | $+$ | .00533 | .009934 | .0104122 | .0104010 | |
| $\log P'$ | $+$ | 8.01967 | 8.369626 | 8.3957468 | 8.3951457 | |
| $\alpha _{1}$ | $+$ | .30136 | .336506 | .3390605 | .3390018 | |
| $\beta _{1}$ | $+$ | 1.61938 | 1.307304 | 1.286223 | 1.2867056 | |
| $\gamma _{1}$ | $-$ | .03072 | .025316 | .0249518 | .0249601 | |
| $\alpha '$ | $-$ | .050505 | .025316 | .0249518 | .0249601 | |
| $\beta '$ | $+$ | .47139 | | .0563438 | | |
| $\gamma '$ | $-$ | .00818 | | .0079821 | | |
| $q_{2}$ | $=\rho _{2}-.9314223$ | $\alpha _{2}$ | $=+[9.2282738](q_{2}+1.7286820)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
| $r_{2}^{2}$ | $=q_{2}^{2}+.100681$ | $\beta _{2}$ | $=-[9.9900800](q_{2}-.0361309)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
| $R_{2}$ | $=[0.38183722]r_{2}^{-3}$ | $\gamma _{2}$ | $=-[9.1026549](q_{2}-.9314223)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ |
| $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\alpha ''$ | $=.169151-[9.22827]R_{2}-P''\alpha _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}''$ |
| $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\beta ''$ | $=-.977417+[9.99008]R_{2}-P''\beta _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}''$ |
| $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $P''={\frac {[.47712]R_{2}q_{2}}{(1-R_{2})r_{2}^{2}}}$ | $\gamma ''$ | $=-.126664+[9.10265]R_{2}-P''\gamma _{2}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}''$ |
| $\Delta q_{2}$ | | | –.77826 | + .005042 | +.0013222 | + .0000021 |
| $q_{2}$ | $+$ | 3.34235 | 2.56409 | 2.569132 | 2.5704542 | 2.5704563 |
| $\log r_{2}$ | $+$ | .52600 | .412233 | .4130733 | .4132934 | .4132937 |
| $\log R_{2}$ | $+$ | 8.74037 | 9.081673 | 9.0791524 | 9.0784920 | |
| $\log(1-R_{2})$ | $+$ | 9.97543 | 9.944142 | 9.9444866 | 9.9445766 | |
| $\log P''$ | $+$ | 8.71411 | 9.199120 | 9.1954270 | 9.1944598 | |
| $\alpha _{2}$ | $+$ | .81059 | .638489 | .6397466 | .6400760 | |
| $\beta _{2}$ | $-$ | 3.05379 | 2.172660 | 2.1787230 | 2.1803116 | |
| $\gamma _{2}$ | $-$ | .28858 | .181843 | .1825486 | .1827338 | |
| $\alpha ''$ | $+$ | .20182 | | .2491854 | | |
| $\beta ''$ | $-$ | 1.08177 | | 1.2018221 | | |
| $\gamma ''$ | $+$ | .13464 | | .1400944 | | |
| $q_{3}$ | $=\rho _{3}-.5599304$ | $\alpha _{3}$ | $=-[9.3810737](q_{3}+1.5798163)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
| $r_{3}^{2}$ | $=q_{3}^{2}+.7130624$ | $\beta _{3}$ | $=-[9.6537308](q_{3}-.4630521)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
| $R_{3}$ | $=[9.5623916]r_{3}^{-3}$ | $\gamma _{3}$ | $=[8.8361256](q_{3}+.5599304)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ |
| $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $\alpha '''$ | $=-.240477-[9.38107]R_{3}-P'''\alpha _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'''$ |
| $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $\beta '''$ | $=+.450537+[9.65373]R_{3}-P'''\beta _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'''$ |
| $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $P'''={\frac {[.47712]R_{3}q_{3}}{(1-R_{3})r_{3}^{2}}}$ | $\gamma '''$ | $=+.068569+[8.83613]R_{3}-P'''\gamma _{3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \ \end{matrix}}\right\}\,$${\text{III}}'''$ |
| $\Delta q_{3}$ | | | –.80780 | –.04055 | +.0025316 | + .0000031 |
| $q_{3}$ | $+$ | 3.24945 | 2.44165 | 2.40110 | 2.4036316 | 2.4036347 |
| $\log r_{3}$ | $+$ | 0.52600 | .412217 | .4057319 | .4061394 | .4061399 |
| $\log R_{3}$ | $+$ | 7.98439 | 8.325742 | 8.3451948 | 8.3439733 | |
| $\log(1+R_{3})$ | $+$ | .00417 | .009099 | .0095108 | .0094843 | |
| $\log P'''$ | $+$ | 7.91715 | 8.357016 | 8.3817516 | 8.3801993 | |
| $\alpha _{3}$ | $-$ | 1.17253 | .987590 | .9785152 | .9790776 | |
| $\beta _{3}$ | $+$ | 1.26749 | .910305 | .8924956 | .8936069 | |
| $\gamma _{3}$ | $-$ | .26373 | .210171 | .2075292 | .2076940 | |
| $\alpha '''$ | $-$ | .22847 | | .2222335 | | |
| $\beta '''$ | $+$ | .44441 | | .4390163 | | |
| $\gamma '''$ | $+$ | .06690 | | .0650888 | | |
IV.
The values of $\alpha ',\beta ',$ etc., furnish the basis for the computation of the following quantities:
| $a_{1}=-.01254$ | $a_{2}=-.03517$ | $a_{3}=-.07232$ |
| $b_{1}=+.01726$ | $b_{2}=-.00525$ | $b_{3}=-.00845$ |
| $c_{1}=-.15746$ | $c_{2}=-.08526$ | $c_{3}=-.04050$ |
For we get three values sensibly identical. Adopting the mean, we set
* $G=.01006.$
* $H=-.00998,$$L=.02322$
V.
Taking the values of $\alpha _{1},\alpha _{2},$ etc., from the columns under ${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3},$ we form the residuals
* $\alpha =-.06058,$$\beta =-.16692,$$\gamma =-.05557.$
* $C_{1}=-.65888,$$C_{2}=-.76983,$$C_{3}=-.79939,$
$q_{1},q_{2},q_{3}.$
* $\Delta q_{2}=C_{2}-{\tfrac {6}{10}}L(\Delta q_{2})^{2},$or$\Delta q_{2}=-.76983-.01393(\Delta q_{2})^{2},$
* $\Delta q_{2}=-.77826.$
$\Delta q_{2}$
$C_{1}$
$C_{2}$
* $\Delta q_{1}=-.66731,$$\Delta q_{3}=-.80780.$
VI.
Applying these corrections to the values of $q_{1},q_{2},q_{3}$ compute the second numerical columns under equations ${\text{III}}_{1},{\text{III}}_{2},$ and ${\text{III}}_{3}.$ We do not go on to the computations by equations ${\text{III}}',$ etc., but content ourselves with the old values of $a_{1},b_{1},$ etc., $G,$ and $L,$ which with the new residuals
* $\alpha =-.012595,$$\beta =.044949,$$\gamma =.003012,$
* $C_{1}=-.04567,$$C_{2}=.004952,$$C_{3}=-.04064.$
* $\Delta q_{2}=C_{2}-L(Dq_{2}+{\tfrac {1}{2}}C_{2})\Delta q_{2}=.004952-.02322(-.77826+.00247)\Delta q_{2}.$
* $\Delta q_{2}=.005042.$
$L$
$\Delta q_{2}$
$C_{1}$
$C_{3},$
* $\Delta q_{1}=-.04558,$$\Delta q_{3}=-.04055.$
${\text{III}}_{1},$
$\alpha ',$
| $a_{1}=-.01671$ | $a_{2}=-.03351815$ | $a_{3}=-.0743299$ |
| $b_{1}=+.0149145$ | $b_{2}=-.0054413$ | $b_{3}=-.0098825$ |
| $c_{1}=-.1576886$ | $c_{2}=-.0779570$ | $c_{3}=-.0474318$ |
| | $G=.0090929.$ | |
The quantities $H$ and $L$ we neglect as of no consequence at this stage of the approximation.
With these values the new residuals
* $\alpha =+.0002919,$$\beta =-.0000044,$$\gamma =+.0000288,$
* $\Delta q_{1}=C_{1}=+.0010434,$$\Delta q_{2}=C_{2}=+.0013222,$ $\Delta q_{3}=C_{3}=+.0025316.$
${\text{III}}_{1},$
* $\alpha =+.0000002,$$\beta =+.0000009,$$\gamma =+.0000001,$
* $\Delta q_{1}=+.0000006,$$\Delta q_{2}=+.0000021,$$\Delta q_{3}=+.0000031.$
$q_{1},q_{2},q_{3}$
* $\log r_{1}=0.428377,$$\log r_{2}=0.4132937,$$\log r_{3}=0.4061399.$
VII.
The computations for the test of the uncorrected formula (the first hypothesis) are as follows:
| | | Number of arc. | Logarithm. |
| $r_{1}$ | | | 0.4282377 |
| $r_{2}$ | | | 0.4132937 |
| $r_{3}$ | | | 0.4061399 |
| $A_{1}B_{1}r_{1}^{-3}$ | $+$ | .01174865 | 8.0699879 |
| $B_{2}r_{2}^{-3}$ | $+$ | .11980944 | 9.0784911 |
| $A_{3}B_{3}r_{3}^{-3}$ | $+$ | .01137670 | 8.0560162 |
| $N$ | $+$ | .14293479 | 9.1551380 |
| $s_{1}$ | $+$ | 1.3308476 | 0.1241283 |
| $s_{2}$ | $+$ | 2.2796616 | 0.3578704 |
| $s_{3}$ | $+$ | 1.3417404 | 0.1276685 |
| $s$ | $+$ | 2.4761248 | 0.3937725 |
| $s-s_{1}$ | $+$ | 1.1452772 | 0.0589106 |
| $s-s_{2}$ | $+$ | 0.1964632 | 9.2932812 |
| $s-s_{3}$ | $+$ | 1.1343844 | 0.0547602 |
| ${R}$ | $+$ | | 9.5065898 |
| $p$ | $+$ | | 0.4391732 |
| $\tan {\tfrac {1}{2}}(v_{2}-v_{1})$ | $+$ | 15° 48' 10''.82 | 9.4518296 |
| $\tan {\tfrac {1}{2}}(v_{3}-v_{2})$ | $+$ | 15° 39' 36''.38 | 9.4476792 |
| $\tan {\tfrac {1}{2}}(v_{3}-v_{1})$ | $+$ | 31° 27' 47''.20 | 9.7866915 |
| $e\sin {\tfrac {1}{2}}(v_{3}+v_{1})$ | $-$ | | 8.7099387 |
| $e\cos {\tfrac {1}{2}}(v_{3}+v_{1})$ | $+$ | | 8.7872701 |
| $\tan {\tfrac {1}{2}}(v_{3}+v_{1})$ | $-$ | –39° 55' 32''.32 | 9.9226686 |
| $e$ | $+$ | | 8.9025438 |
| $\epsilon$ | $+$ | | 9.9652259 |
| $a$ | $+$ | | 0.4419546 |
| $\tan {\tfrac {1}{2}}v_{1}$ | $-$ | –35° 41' 39''.75 | 9.8563809 |
| $\tan {\tfrac {1}{2}}v_{2}$ | $-$ | –19° 53' 28''.93 | 9.5384981 |
| $\tan {\tfrac {1}{2}}v_{3}$ | $-$ | –4° 13' 52''.55 | 8.8691380 |
| $\tan {\tfrac {1}{2}}E_{1}$ | $-$ | –33° 33' 0''.17 | 9.8216068 |
| $\tan {\tfrac {1}{2}}E_{2}$ | $-$ | –18° 28' 6''.35 | 9.5237240 |
| $\tan {\tfrac {1}{2}}E_{3}$ | $-$ | –3° 54' 24''.21 | 8.8343639 |
| $\sin {\tfrac {1}{2}}E_{1}$ | $-$ | –67° 6' 0''.34 | 9.9643473 |
| $\sin {\tfrac {1}{2}}E_{2}$ | $-$ | –36° 56' 12''.70 | 9.7788272 |
| $\sin {\tfrac {1}{2}}E_{3}$ | $-$ | –7° 48' 48''.42 | 9.1333734 |
| $ea^{\frac {3}{2}}\sin {\tfrac {1}{2}}E_{1}$ | $-$ | .3387061 | 9.5298230 |
| $ea^{\frac {3}{2}}\sin {\tfrac {1}{2}}E_{2}$ | $-$ | .2209545 | 9.3343029 |
| $ea^{\frac {3}{2}}\sin {\tfrac {1}{2}}E_{3}$ | $-$ | .0499681 | 8.6988491 |
| $a^{\frac {3}{2}}(E_{2}-E_{1})$ | $+$ | 2.4226307 | 0.3842872 |
| $a^{\frac {3}{2}}(E_{3}-E_{2})$ | $+$ | 2.3391145 | 0.3690515 |
| $\tau _{3{\text{ calc.}}}$ | $+$ | 2.3048791 | 0.3626482 |
| $\tau _{1{\text{ calc.}}}$ | $+$ | 2.1681461 | 0.3360885 |
{{c|VIII. The logarithms of the calculated values of the intervals of time exceed those of the given values by .0002416 for the first interval ($\tau _{3}$) and .0002365 for the second ($\tau _{1}$). Therefore, since the corrections for aberration have been incorporated in the data, we set for the correction of the formula (for the second hypothesis)
* $\Delta \log \tau _{1}=-.0002365$$\Delta \log \tau _{2}=-.0002416$
* $\Delta \log A_{3}=.0000026$$\Delta \log A_{2}=-.0000025$ $\Delta \log B_{1}=-.004872$$\Delta \log B_{2}=-.004782$$\Delta \log B_{3}=-.0004665$
$A_{1},A_{3}$
* $\log A_{1}=9.6854923$$\log A_{3}=9.7120418$
${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3}$
| $r_{1}^{2}$ | $=q_{1}^{2}+.8645336$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ corrected. |
| | $\alpha _{1}=-[8.6700193](q_{1}-9.5901555)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ corrected. |
| $R_{1}$ | $=[9.6687241]r_{1}^{-3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ corrected. |
| | $\beta _{1}=+[9.6833950](q_{1}+.0900552)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ corrected. |
| | $\gamma _{1}=-[7.9242073](q_{1}+.3874081)(1+R_{1})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{1}$ corrected. |
| $\Delta q_{1}$ | | | + .002887 | –.0000217 |
| $q_{1}$ | $+$ | 2.5142140 | 2.5145027 | 2.5144810 |
| $\log r_{1}$ | $+$ | .4282377 | .4282816 | .4282782 |
| $\log R_{1}$ | $+$ | 8.3838110 | 8.3838793 | 8.3838894 |
| $\log(1+R_{1})$ | $+$ | .0103847 | .0103863 | .0103865 |
| $\alpha _{1}$ | $+$ | .3389910 | .3389784 | .3389796 |
| $\beta _{1}$ | $+$ | 1.2866654 | 1.2868124 | 1.2868024 |
| $\gamma _{1}$ | $-$ | .0249593 | .0249619 | .0249617 |
| $\log s_{1}$ | $+$ | | | .1241571 |
| $r_{2}^{2}$ | $=q_{2}^{2}+.1006681$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ corrected. |
| | $\alpha _{2}=+[9.2282738](q_{2}+1.7286820)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ corrected. |
| $R_{2}$ | $=-[9.9900800]r_{2}^{-3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ corrected. |
| | $\beta _{2}=-[9.9900800](q_{2}-.0361309)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ corrected. |
| | $\gamma _{2}=-[9.1026549](q_{2}-.9314223)(1-R_{2})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{2}$ corrected. |
| $\Delta q_{2}$ | | | – .0000955 | +.0000187 |
| $q_{2}$ | $+$ | 2.5704563 | 2.5703608 | 2.5703795 |
| $\log r_{2}$ | $+$ | 0.4132937 | .4132778 | .4132809 |
| $\log R_{2}$ | $+$ | 9.0780129 | 9.0780605 | 9.0780513 |
| $\log(1-R_{2})$ | $+$ | 9.9446418 | 9.9446353 | 9.9446365 |
| $\alpha _{2}$ | $+$ | .6401725 | .6401487 | .6401532 |
| $\beta _{2}$ | $-$ | 2.1806412 | 2.1805261 | 2.1805482 |
| $\gamma _{2}$ | $-$ | .1827615 | .1827481 | .1827507 |
| $\log s_{2}$ | $+$ | | | .3579174 |
| $r_{3}^{2}$ | $=q_{3}^{2}+.7130624$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ corrected. |
| | $\alpha _{3}=-[9.3810712](q_{3}+1.5798163)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ corrected. |
| $R_{3}$ | $=[9.5619251]r_{3}^{-3}$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ corrected. |
| | $\beta _{3}=+[9.6537283](q_{3}-.4630521)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ corrected. |
| | $\gamma _{3}=+[8.8361231](q_{3}+.5599304)(1+R_{3})$ | $\scriptstyle {\left.{\begin{matrix}\ \\\\\ \\\ \\\ \\\ \\\ \\\ \ \end{matrix}}\right\}\,$ ${\text{III}}_{3}$ corrected. |
| $\Delta q_{3}$ | | | + .0003302 | +.0000424 |
| $q_{3}$ | $+$ | 2.4036347 | 2.4039649 | 2.4040073 |
| $\log r_{3}$ | $+$ | .4061399 | .4061929 | .4061998 |
| $\log R_{3}$ | $+$ | 8.3435055 | 8.3433463 | 8.3433257 |
| $\log(1+R_{3})$ | $+$ | .094742 | .0094708 | .0094704 |
| $\alpha _{3}$ | $-$ | .9790500 | .9791236 | .9791329 |
| $\beta _{3}$ | $+$ | .8935824 | .8937277 | .8937461 |
| $\gamma _{3}$ | $+$ | .2076882 | .2077097 | .2077124 |
| $\log s_{3}$ | $+$ | | | .1277120 |
With these corrected equations the last values of $q_{1},q_{2},q_{3}$ give the residuals
* $\alpha =.0001135$$\beta =-.0003934$$\gamma =-.0000326$
* $\Delta q_{1}=.0002887$$\Delta q_{2}=-.0000955$$\Delta q_{3}=.0003302$
* $\alpha =.0000035$$\beta =.0000140$$\gamma =-.0000003$
* $\Delta q_{1}=.-.0000217$$\Delta q_{2}=.0000187$$\Delta q_{3}=.0000424$
* $\alpha =-.0000001$$\beta =.0000003$$\gamma =.0000000$
IX, X.
${\text{III}}_{1},{\text{III}}_{2},{\text{III}}_{3},$
| | | Number of arc. | Logarithm. | | | Number of arc. | Logarithm. |
| $\alpha _{4}$ | $-$ | 1.3181125 | 0.1199525 | $e$ | $+$ | 0.8076104 | 8.9072019 |
| $\beta _{4}$ | $-$ | 3.930563 | 9.5944547 | $\tan {\tfrac {1}{2}}v_{1}$ | $-$ | –35° 26' 25''.87 | 9.8523140 |
| $\gamma _{4}$ | $+$ | .2326741 | 9.3667481 | $\tan {\tfrac {1}{2}}v_{2}$ | $-$ | –19° 38' 18''.20 | 9.5524725 |
| $s_{4}^{2}$ | $+$ | 1.9460513 | 0.2891543 | $\tan {\tfrac {1}{2}}v_{3}$ | $-$ | –5° 58' 47''.06 | 8.8424313 |
| $s_{4}$ | $+$ | 1.3950100 | 0.1445772 | $\epsilon$ | $+$ | | 9.9648493 |
| $s_{1}$ | $+$ | 1.3309358 | | $a$ | $+$ | | 0.4424586 |
| $s_{3}$ | $+$ | 1.3418749 | | $b$ | $+$ | | 0.4410376 |
| $s_{2}$ | $+$ | 2.2799084 | 0.3579174 | $l$ | $-$ | .8139130 | 9.9105780 |
| $s$ | $+$ | 2.4763596 | 0.3938137 | $m$ | $+$ | .5564152 | 9.7453990 |
| $s-s_{1}$ | $+$ | 1.1454238 | 0.0589662 | $n$ | $+$ | .1671768 | 9.2231760 |
| $s-s_{3}$ | $+$ | 1.1344847 | 0.0547986 | $\lambda$ | $-$ | .5517065 | 9.7417081 |
| $S$ | $+$ | 2.0339104 | 0.3083318 | $\mu$ | $-$ | .8304016 | 9.9192881 |
| $S-s_{1}$ | $+$ | .7029746 | 9.8469397 | $\nu$ | $+$ | .07780215 | 8.8909916 |
| $S-s_{4}$ | $+$ | .6389004 | 9.8054332 | $a_{x}$ | $-$ | | 0.3530366 |
| $S-s_{3}$ | $+$ | .6920355 | 9.8401284 | $a_{y}$ | $+$ | | 0.1878576 |
| ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}$ | $-$ | .0292373 | 8.4659373 | $a_{z}$ | $+$ | | 9.6656346 |
| ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}s_{2}^{2}{-2}$ | $-$ | | 7.7501025 | $b_{x}$ | $-$ | | .01827457 |
| ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}s_{2}^{2}{-2}\alpha _{2}$ | $-$ | .003600696 | 7.5563864 | $b_{y}$ | $-$ | | 0.3603257 |
| ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}s_{2}^{2}{-2}\beta _{2}$ | $+$ | .012265017 | 8.0886682 | $\tan {\tfrac {1}{2}}E_{1}$ | $-$ | –33° 16' 50''.08 | 9.8171633 |
| ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}s_{2}^{2}{-2}\gamma _{2}$ | $+$ | .001027926 | 7.0119616 | $\tan {\tfrac {1}{2}}E_{2}$ | $-$ | –18° 12' 58''.10 | 9.5173218 |
| $\alpha _{5}$ | $-$ | 1.3145118 | 0.1187645 | $\tan {\tfrac {1}{2}}E_{3}$ | $-$ | –3° 40' 16''.33 | 8.8072806 |
| $\beta _{5}$ | $-$ | .4053213 | 9.6077994 | $k^{-1}a^{\frac {3}{2}}e\sin E_{1}$ | $-$ | 19.85671 | 1.2979074 |
| $\gamma _{5}$ | $+$ | .2316462 | 9.3648252 | $k^{-1}a^{\frac {3}{2}}e\sin E_{2}$ | $-$ | 12.85291 | 1.1090014 |
| $s_{5}^{2}$ | $+$ | 1.9458866 | 0.2891176 | $k^{-1}a^{\frac {3}{2}}E_{1}$ | $-$ | 311.31839 | 2.4932048 |
| $_{{\text{R}}^{2}}$ | $+$ | | 9.0132057 | $k^{-1}a^{\frac {3}{2}}E_{2}$ | $-$ | 170.40027 | 2.2314703 |
| $_{\text{R}}$ | $+$ | | 9.5066028 | $k^{-1}a^{\frac {3}{2}}E_{3}$ | $-$ | 34.34174 | 1.538223 |
| $N$ | $+$ | .14278017 | 9.1546679 | $T_{(1)}$ | $+$ | 296.97504 | |
| $p$ | $+$ | 2.7517987 | 0.4396167 | $T_{(2)}$ | $+$ | 296.97447 | |
| $\tan {\tfrac {1}{2}}(v_{2}-v_{1})$ | $+$ | 15° 48' 7''.67 | 9.4518042 | $T_{(3)}$ | $+$ | 296.97398 | |
| $\tan {\tfrac {1}{2}}(v_{3}-v_{2})$ | $+$ | 15° 39' 31''.14 | 9.4476366 | $T$ | $+$ | 296.97450 | |
| $\tan {\tfrac {1}{2}}(v_{3}-v_{1})$ | $+$ | 31° 27' 38''.81 | 9.7866517 | | | |
| $e\sin {\tfrac {1}{2}}(v_{3}+v_{1})$ | $-$ | | 8.7099782 | | | |
| $e\cos {\tfrac {1}{2}}(v_{3}+v_{1})$ | $+$ | | 8.7951055 | | | |
| $\tan {\tfrac {1}{2}}(v_{3}+v_{1})$ | $-$ | 39° 25' 12''.94 | 9.9148727 | | | |
XI.
This gives the following equations for an ephemeris:
* $T=1806,{\text{ June, }}23.97450,{\text{ Paris mean time}}$ $[2.8863186](t-T)=E_{\text{in seconds}}-[4.2216270]\sin E$
Heliocentric coordinates relating to the ecliptic. :
* $x=+.1820700-[0.3530366]\cos E-[0.1827457]\sin E$
* $y=-.1244685+[0.1878576]\cos E-[0.3603257]\sin E$
* $z=-.03739370+[9.6656346]\cos E+[9.3320292]\sin E$
The differences of the values of $T_{(1)},T_{(2)},T_{(3)}$ from their mean $T,$ indicate the residual errors of this hypothesis. They indicate differences in the calculated and the observed geocentric positions which are represented by the geocentric angles subtended by the path described by the planet in the following fractions of a day: .00054, .00003, .00052. Since the heliocentric motion of the planet is about one-fourth of a degree per day, and the planet is considerably farther from the earth than from the sun at the times of the first and third observations, the errors will be less than half a second in arc.
If we desire all the accuracy possible with seven-figure logarithms, we may form a third hypothesis based on the following corrections:
* ${\begin{aligned}\Delta \log \tau _{1}&=M{\frac {T_{(3)-T_{(2)}}}{t_{3}-t_{2}}}&=-.0000017,\\\Delta \log \tau _{2}&=M{\frac {T_{(3)-T_{(1)}}}{t_{2}-t_{1}}}&=-.0000018.\end{aligned}}$
* $T=1806,{\text{ June, }}23.96378,{\text{ Paris mean time}}$ $[2.8863140](t-T)=E_{\text{in seconds}}-[4.2216530]\sin E$
Heliocentric coordinates relating to the ecliptic. :
* $x=+.1820765-[0.3530261]\cos E-[0.1827783]\sin E$
* $y=-.1244853+[0.1878904]\cos E-[0.3603153]\sin E$
* $z=-.03739987+[9.6656285]\cos E+[9.3320758]\sin E$
The agreement of the calculated geocentric positions with the data is shown in the following table:
| Times, 1805, September | 2.51336 | 139.42711 | 265.39813 |
| Second hypothesis: | | | |
| longitudes | 95° 32' 18''.88 | 99° 49' 5''.87 | 118° 5' 28''.52 |
| errors | 0''.32 | 0''.00 | –0''.33 |
| latitudes | –0° 59' 34''.01 | 7° 16' 36''.82 | 7° 38' 49''.34 |
| errors | 0''.05 | 0''.02 | –0''.05 |
| Third hypothesis: | | | |
| longitudes | 95° 32' 18''.65 | 99° 49' 5''.82 | 118° 5' 28''.79 |
| errors | 0''.09 | –0''.05 | –0''.06 |
| latitudes | –0° 59' 34''.04 | 7° 16' 36''.78 | 7° 38' 49''.38 |
| errors | 0''.02 | –0''.02 | –0''.01 |
| | $\log r_{1}.$ | $\log r_{2}.$ | $\log r_{3}.$ |
| First hypothesis | .4282377 | .4132937 | .4061399 |
| Second hypothesis | .4282782 | .4132809 | .4061998 |
| Third hypothesis | .4282786 | .4132808 | .4062003 |
| Gauss: | | |
| First hypothesis | .4323934 | .4114726 | .4064712 |
| Second hypothesis | .4291773 | .4129371 | .4071975 |
| Third hypothesis | .4284841 | .4132107 | .4064697 |
| Fourth hypothesis | .4282792 | .4132817 | .4062033 |
| Oppolzer: | | |
| First hypothesis | .4281340 | .413330 | .4061699 |
| Second hypothesis | .4282794 | .4132801 | .4061976 |
| Third hypothesis | .4282787 | | .4062009 |
In comparing the different methods, it should be observed that the determination of the positions in any hypothesis by Gauss's method requires successive corrections of a single independent variable, a corresponding determination by Oppolzer's method requires the successive corrections of two independent variables, while the corresponding determination by the method of the present paper requires the successive corrections of three independent variables.
* ↑ Vectors, or directed quantities, will be represented in this paper by German capitals. The following notations wiU be used in connection with them: The sign $=$ denotes identity in direction as well as length. The sign $+$ denotes geometrical addition, or what is called composition in mechanics. The sign $-$ denotes reversal of direction, or composition after reversal. The notation ${\mathfrak {A}}.{\mathfrak {B}}$ denotes the product of the lengths of the vectors and the cosine of the angle which they include. It will be called the direct product of ${\mathfrak {A}}$ and ${\mathfrak {B}}.$ If $x,y,z$ are the rectangular components of ${\mathfrak {A}}$ and $x',y',z'$ those of ${\mathfrak {B}},$ ${\mathfrak {A}}.{\mathfrak {B}}=xx'+yy'+zz'.$ ${\mathfrak {A.A}}$ may be written ${\mathfrak {A}}^{2}$ and called the square of ${\mathfrak {A}}.$ The notation ${\mathfrak {A}}\times {\mathfrak {B}}$ will be used to denote a vector of which the length is the product of the lengths of ${\mathfrak {A}}$ and ${\mathfrak {B}}$ and the sine of the angle which they include. Its direction is perpendicular to ${\mathfrak {A}}$ and ${\mathfrak {B}}$ and on that side on which a rotation from ${\mathfrak {A}}$ to ${\mathfrak {B}}$ appears counter-clockwise. It will be called the skew product of ${\mathfrak {A}}$ and ${\mathfrak {B}}.$ If the rectangular components of ${\mathfrak {A}}$ and ${\mathfrak {B}}$ are $x,y,z,$ and $x',y',z',$ those of ${\mathfrak {A}}\times {\mathfrak {B}}$ will be $yz'-zy',$$zx'-xz',$$xy'-yx'.$ The notation $({\mathfrak {ABC}})$ denotes the volume of the parallelepiped of which three edges are obtained by laying off the vectors ${\mathfrak {A}},{\mathfrak {B}},$ and ${\mathfrak {C}}$ from any same point, which volume is to be taken positively or negatively, according as the vector ${\mathfrak {C}}$ falls on the side of the plane containing ${\mathfrak {A}}$ and ${\mathfrak {B}},$ on which a rotation from ${\mathfrak {A}}$ to ${\mathfrak {B}}$ appears counter-clockwise, or on the other side. If the rectangular components of ${\mathfrak {A}},{\mathfrak {B}},$ and ${\mathfrak {C}}$ are $x,y,z;x',y',z';$ and $x'',y'',z'',$ $({\mathfrak {ABC}})={\begin{vmatrix}x&y&z\\x'&y'&z'\\x''&y''&z''\\\end{vmatrix}}.$ It follows, from the above definitions, that for any vectors ${\mathfrak {A}},{\mathfrak {B}},$ and ${\mathfrak {C}}$ ${\mathfrak {A}}.{\mathfrak {B}}={\mathfrak {B}}.{\mathfrak {A}},$${\mathfrak {A}}\times {\mathfrak {B}}=-{\mathfrak {B}}\times {\mathfrak {A}},$ $({\mathfrak {ABC}})=({\mathfrak {BCA}})=({\mathfrak {CAB}})=-({\mathfrak {ACB}})=-({\mathfrak {CBA}})=-({\mathfrak {BAC}}),$ and $({\mathfrak {ABC}})={\mathfrak {A}}.({\mathfrak {B}}\times {\mathfrak {C}})={\mathfrak {B}}.({\mathfrak {C}}\times {\mathfrak {A}})={\mathfrak {C}}.({\mathfrak {A}}\times {\mathfrak {B}});$ also that ${\mathfrak {A}}.{\mathfrak {B}},{\mathfrak {A}}\times {\mathfrak {B}}$ are distributive functions of ${\mathfrak {A}}$ and ${\mathfrak {B}},$ and $({\mathfrak {ABC}})$ a distributive function of ${\mathfrak {A}},{\mathfrak {B}},$ and ${\mathfrak {C}},$ for example, that if ${\mathfrak {A}}={\mathfrak {L}}+{\mathfrak {M}},$ ${\mathfrak {A}}.{\mathfrak {B}}={\mathfrak {L}}.{\mathfrak {B}}+{\mathfrak {M}}.{\mathfrak {B}},$${\mathfrak {A}}\times {\mathfrak {B}}={\mathfrak {L}}\times {\mathfrak {B}}+{\mathfrak {M}}\times {\mathfrak {B}},$$({\mathfrak {ABC}})=({\mathfrak {LBC}})+({\mathfrak {MBC}}),$ and so for ${\mathfrak {B}}$ and ${\mathfrak {C}}.$ The notation $({\mathfrak {ABC}})$ is identical with that of Lagrange in the Mécanique Analytique, except that there its use is limited to unit vectors. The signification of ${\mathfrak {A}}\times {\mathfrak {B}}$ is closely related to, but not identical with, that of the notation $[r_{1}r_{2}]$ commonly used to denote the double area of a triangle determined by two positions in an orbit.
* ↑ These equations are obtained by taking the direct products of both members of the preceding equation with ${\mathfrak {S}}''\times {\mathfrak {S}}''',$ ${\mathfrak {S}}'''\times {\mathfrak {S}}',$ and ${\mathfrak {S}}'\times {\mathfrak {S}}'',$ respectively. See footnote on page 119.
* ↑ In the case of Swift's comet (V, 1880), the writer found by the quadratic equation –.247 and –.116 for corrections of the assumed geocentric distance 1250. The first of these numbers gives an approximation to the position of the earth; the second to that of the comet, viz., the geocentric distance .134 instead of the true value .1333. The coefficient ${\tfrac {6}{10}}$ was used in the quadratic equation; with the coefficient the approximations would not be quite so good. The value of the correction obtained by neglecting the quadratic term was .079, which indicates that the approximations (in this very critical case) would be quite tedious without the use of the quadratic term.
* ↑ In the case mentioned in the preceding footnote, from $[\Delta q_{2}]=.079$ and $\rho _{2}=.25,$ we get $\Delta q_{2}=-.1155,$ which is sensibly the same value as that obtained by calculating the quadratic term.
* ↑ When an approximate orbit is known in advance, we may correct the fundamental equation at onoe. The formulæ will be given in the Summary, § xii.
* ↑ It would have been better to omit altogether the oaloulation of $H$ and $L,$ if the small value of the latter coidd have been foreseen. In fact, it will be found that the terms containing $L$ hardly improve the convergence, being smaller than quantities which have been neglected. Nevertheless, the use of these terms in this example will illustrate a process which in other cases may be beneficial.
* ↑ The corrections may be made without rewriting the equation.
* ↑ ${\mathfrak {S}}_{4}.{\mathfrak {S}}_{2}=\alpha _{4}\alpha _{2}+\beta _{4}\beta _{2}+\gamma _{4}\gamma _{2}$
* ↑ Theoria motus, § 159.
* ↑ Lehrbuch zur Bahnbestimmung der Kometen und Planeten, 2nd ed., vol. i, p. 394.