John S. Prell . THEORETICAL ASTRONOMY RELATING TO THE MOTIONS OF THE HEAVENLY BODIES REVOLVING AROUND THE SUN IN ACCORDANCE WITH THE LAW OF UNIVERSAL GRAVITATION EMBRACING A SYSTEMATIC DERIVATION OF THE FORMULA FOR THE CALCULATION OF THE GEOCENTRIC AND UELIO- CENTRIC PLACES, FOR THE DETERMINATION OF THE ORBITS OF PLANETS ANi COMETS, FOB THE CORRECTION OF APPROXIMATE ELEMENTS, AND FOR THE COMPUTATION OF SPECIAL PERTURBATIONS; TOGETHER WITH THE THEORY Of THE COMBI- NATION OF OBSERVATIONS AND THE METHOD OF LEAST SQUARES. Utiih Jtummral femgte and JOHN S. PRELL Civil & Mechanical Engineer. SAN FRANCISCO, CAL. BY JAMES C. WATSON DIRECTOR OF THE OBSERVATORY AT ANN ARBOR, AND PROFESSOR OF ASTRONOMY IN THE UNIVERSITY OF MICHIGAN PHILADELPHIA: J. B. LIPPINCOTT COMPANY. LONDON: 36 SOUTHAMPTON STREET, COVENT GARDEN. 1900. Entered, according to Act of Congress, in the year 1868 '-v J. B. LIPPINCOTT & CO., in the Clerk's Office of the District Court of the United States for the Eastern District of Pennsylvania. Copyright, 1896, by ANNETTE H. WATSON. GIFT 555 JOHN S. PRELL Civil & Mechanical Engineer. SAN FRANCISCO, CAL. PREFACE. THE discovery of the great law of nature, the law of gravitation, by NEWTON, prepared the way for the brilliant achievements which have distinguished the history of astronomical science. A first essential, how- ever, to the solution of those recondite problems which were to exhibit the effect of the mutual attraction of the bodies of our system, was the development of the infinitesimal calculus; and the labors of those who devoted themselves to pure analysis have contributed a most important part in the attainment of the high degree of perfection which character- izes the results of astronomical investigations. Of the earlier efforts to develop the great results following from the law of gravitation, those of EULER stand pre-eminent, and the memoirs which he published have, in reality, furnished the germ of all subsequent investigations in celestial mechanics. In this connection also the names of BERNOTJILLI, CLAIRATJT, and D'ALEMBERT deserve the most honorable mention as having contributed also, in a high degree, to give direction to the inves- tigations which were to unfold so many mysteries of nature. By means of the researches thus inaugurated, the great problems of mechanics were successfully solved, many beautiful theorems relating to the planet- ary motions demonstrated, and many useful formulae developed. It is true, however, that in the early stage of the science methods were developed which have since been found to be impracticable, even if not erroneous; still, enough was effected to direct attention in the proper channel, and to prepare the way for the more complete labors of LAGE-O.NGE and LAPLACE. The genius and the analytical skill of these extraordinary men gave to the progress of Theoretical Astronomy the most rapid strides ; and the intricate investigations which they success- fully performed, served constantly to educe new discoveries, so that of all the problems relating to the mutual attraction of the several planets 3 M718385 4 PREFACE. but little more remained to be accomplished by their successors than to develop and simplify the methods which they made known, and to intro- duce such modifications as should be indicated by experience or rendered possible by the latest discoveries in the domain of pure analysis. The problem of determining the elements of the orbit of a comet moving in a parabola, by means of observed places, which had been considered by NEWTON, EULER, BOSCOVICH, LAMBERT, and others, received from LAGRANGE and LAPLACE the most careful consideration in the light of all that had been previously done. The solution given by the former is analytically complete, but far from being practically complete ; that given by the latter is especially simple and practical so far as regards the labor of computation ; but the results obtained by it are so affected by the unavoidable errors of observation as to be often little more than rude approximations. The method which was found to answer best in actual practice, was that proposed by OLBERS in his work entitled Leichteste und bequemste Methode die Bahn eines Cometen zu berechnen, in which, by making use of a beautiful theorem of para- bolic motion demonstrated by EULER and also by LAMBERT, and by adopting a method of trial and error in the numerical solution of certain equations, he was enabled to effect a solution which could be performed with remarkable ease. The accuracy of the results obtained by OLBERS'S method, and the facility of its application, directed the attention of LEGENDRE, IVORY, GAUSS, and ENCKE to this subject, and by them the method was extended and generalized, and rendered appli cable in the exceptional cases in which the other methods failed. It should be observed, however, that the knowledge of one element, the eccentricity, greatly facilitated the solution ; and, although elliptic elements had been computed for some of the comets, the first hypothesis was that of parabolic motion, so that the subsequent process required simply the determination of the corrections tc be applied to these ele- ments in order to satisfy the observations. The more difficult problem of determining all the elements of planetary motion directly from three observed places, remained unsolved until the discovery of Ceres by PIAZZI in 1801, by which the attention of GAUSS was directed to this subject, the result of which was the subsequent publication of his Theoria Motus Corporum Cwlestium, a most able work, in which he gave to the world, in a finished form, the results of many years of attention PRPJFACE. 5 to the subject of which it treats. His method for .determining all the elements directly from given observed places, as given in the Theoria Motus, and as subsequently given in a revised form by ENCKE, leaves scarcely any thing to be desired on this topic. In the same work he gave the first explanation of the method of least squares, a method which has been of inestimable service in investigations depending on observed data. The discovery of the minor planets directed attention also to the methods of determining their perturbations, since those applied in the case of the major planets were found to be inapplicable. For a long time astronomers were content simply to compute the special perturba- tions of these bodies from epoch to epoch, and it was not until the com- mencement of the brilliant researches by HANSEN that serious hopes were entertained of being able to compute successfully the general per- turbations of these bodies. By devising an entirely new mode of con- sidering the perturbations, namely, by determining what may be called the perturbations of the time, and thus passing from the undisturbed place to the disturbed place, and by other ingenious analytical and mechanical devices, he succeeded in effecting a solution of this most difficult problem, and his latest works contain all the formulse which are required for the cases actually occurring. The refined and difficult analysis and the laborious calculations involved were such that, even after HANSEN'S methods were made known, astronomers still adhered to the method of special perturbations by the variation of constants as developed by LAGRANGE. The discovery of Astrcea by HENCKE was speedily followed by the discovery of other planets, and fortunately indeed it so happened that the subject of special perturbations was to receive a new improvement. The discovery by BOND and ENCKE of a method by which we determine at once the variations of the rectangular co-ordinates of the disturbed body by integrating the fundamental equations of motion by means of mechanical quadrature, directed the attention of HANSEN to this phase of the problem, and soon after he gave formulse for the determination of the perturbations of the latitude, the mean anomaly, and the loga- rithm of the radius-vector, which are exceedingly convenient in the process of integration, and which have been found to give the most satisfactory results. The formulse for the perturbations of the latitude, 6 PREFACE. true longitude, and radius-vector, to be integrated in the same manner, were afterwards given by BRUNNOW. Having thus stated briefly a few historical facts relating to the problems of theoretical astronomy, I proceed to a statement of the object of this work. The discovery of so many planets and comets has furnished a wide field for exercise in the calculations relating to their motions, and it has occurred to me that a work which should contain a development of all the formulae required in determining the orbits of the heavenly bodies directly from given observed places, and in correcting these orbits by means of more extended discussions of series of observa- tions, including also the determination of the perturbations, together with a complete collection of auxiliary tables, and also such practical directions as might guide the inexperienced computer, might add very materially to the progress of the science by attracting the attention of a greater number of competent computers. Having carefully read the works of the great masters, my plan was to prepare a complete work on this subject, commencing with the fundamental principles of dynamics, and systematically treating, from one point of view, all the problems presented. The scope and the arrangement of the work will be best understood after an examination of its contents ; and let it suffice to add that I have endeavored to keep constantly in view the wants of the computer, providing for the exceptional cases as they occur, and giving all the formulae which appeared to me to be best adapted to the problems under consideration. I have not thought it worth while to trace out the geometrical signification of many of the auxiliary quantities introduced. Those who are curious in such matters may readily derive many beau- tiful theorems from a consideration of the relations of some of these auxiliaries. For convenience, the formulae are numbered consecutively through each chapter, and the references to those of a preceding chapter are defined by adding a subscript figure denoting the number of the chapter. Besides having read the works of those who have given special atten tion to these problems, I have consulted the Astronomische Nachrichten, the Astronomical Journal, and other astronomical periodicals, in which 'is to be found much valuable information resulting from the experi- ence of those w r ho have been or are now actively engaged in astro- nomical pursuits. I must also express my obligations to the publishers, PREFACE. 7 Messrs. J. B. LIPPINCOTI & Co., for the generous interest which they have manifested in the publication of the work, and also to Dr. B. A. GOULD, of Cambridge, Mass., and to Dr. OPPOLZER, of Vienna, for valuable suggestions. For the determination of the time from the perihelion and of the true anomaly in very eccentric orbits I have given the method proposed by BESSEL in the Monatliche Correspondent, vol. xii., the tables for which were subsequently given by BRUNNOW in his Astronomical Notices, and also the method proposed by GAUSS, but in a more convenient form. For obvious reasons, I have given the solution for the special case of parabolic motion before completing the solution of the general problem of finding all of the elements of the orbit by means of three observed places. The differential formulae and the other formulae for correcting approximate elements are given in a form convenient for application, and the formulas for finding the chord or the time of describing the subtended arc of the orbit, in the case of very eccentric orbits, will be found very convenient in practice. I have given a pretty full development of the application of the theory of probabilities to the combination of observations, endeavoring to direct the attention of the reader, as far as possible, to the sources of error to be apprehended and to the most advantageous method of treat- ing the problem so as to eliminate the effects of these errors. For the rejection of doubtful observations, according to theoretical considerations, I have given the simple formula, suggested by CHAUVENET, which fol lows directly from the fundamental equations for the probability of errors, and which will answer for the purposes here required as well as the more complete criterion proposed by PEIRCE. In the chapter devoted to the theory of special perturbations I have taken particular pains to develop the whole subject in a complete and practical form, keeping constantly in view the requirements for accurate and convenient numerical application. The time is adopted as the independent variable in the determination of the perturbations of the elements directly, since experience has established the convenience of this form ; and should it be desired to change the independent variable and to use the differential coefficients with respect to the eccentric anomaly, the equations between this function and the mean motion will enable us to effect readily the required transformation. 8 PREFACE. The numerical examples involve data derived from actual observa- tions, and care has been taken to make them complete in every respect, so as to serve as a guide to the efforts of those not familiar with these calculations ; and when different fundamental planes are spoken of, it is presumed that the reader is familiar with the elements of spherical astronomy, so that it is unnecessary to state, in all cases, whether the centre of the sphere is taken at the centre of the earth, or at any other point in space. The preparation of the Tables has cost me a great amount of labor, logarithms of ten decimals being employed in order to be sure of the last decimal given. Several of those in previous use have been recom- puted and extended, and others here given for the first time have been prepared with special care. The adopted value of the constant of the solar attraction is that given by GAUSS, which, as will appear, is not accurately in accordance with the adoption of the mean distance of the earth from the sun as the unit of space; but until the absolute value of the earth's mean motion is known, it is best, for the sake of uniformity and accuracy, to retain GAUSS'S constant. The preparation of this work has been effected amid many intern; p- tions, and with other labors constantly pressing me, by which the progress of its publication has been somewhat delayed, even since the stereo- typing was commenced, so that in some cases I have been anticipated in the publication of formulae which would have here appeared for the first time. I have, however, endeavored to perform conscientiously the self-imposed task, seeking always to secure a logical sequence in the de- velopment of the formulae, to preserve uniformity and elegance in the notation, and to elucidate the successive steps in the analysis, so that the work may be read by those who, possessing a respectable mathematical education, desire to be informed of the means by which astronomers are enabled to arrive at so many grand results connected with the motions of the heavenly bodies, and by which the grandeur and sublimity of creation are unveiled. The labor of the preparation of the work will have been fully repaid if it shall be the means of directing a more general attention to the study of the wonderful mechanism of the hea- vens, the contemplation of which must ever serve to impress upon the mind the reality of the perfection of the OMNIPOTENT, the LIVING GOD ! OBSEBVATOKY, ANN ARBOR, June, 1867. CONTENTS. THEORETICAL ASTRONOMY. CHAPTER I. INVESTIGATION OP THE FUNDAMENTAL EQUATIONS OF MOTION, AND OF THE FOR- MULA FOR DETERMINING, FROM KNOWN ELEMENTS, THE HELIOCENTRIC AND GEOCENTRIC PLACES OF A HEAVENLY BODY, ADAPTED TO NUMERICAL COM- PUTATION FOR CASES OF ANY ECCENTRICITY WHATEVER. PAOB Fundamental Principles 15 Attraction of Spheres 19 Motions of a System of Bodies '. 23 Invariable Plane of the System 29 Motion of a Solid Body 31 The Units of Space, Time, and Mass 36 Motion of a Body relative to the Sun 38 Equations for Undisturbed Motion 42 Determination of the Attractive Force of the Sun 49 Determination of the Place in an Elliptic Orbit 53 Determination of the Place in a Parabolic Orbit 59 Determination of the Place in a Hyperbolic Orbit 65 Methods for finding the True Anomaly and the Time from the Perihelion in the case of Orbits of Great Eccentricity 70 Determination of the Position in Space 81 Heliocentric Longitude and Latitude 83 Reduction to the Ecliptic 85 Geocentric Longitude and Latitude 86 Transformation of Spherical Co-ordinates 87 Direct Determination of the Geocentric Eight Ascension and Declination 90 .Reduction of the Elements from one Epoch to another 99 Numerical Examples 103 Interpolation < 112 Time of Opposition 114 10 CONTENTS. CHAPTER II. INVESTIGATION OF THE DIFFERENTIAL FORMULAE WHICH EXPRESS THE RELATION BETWEEN THE GEOCENTRIC OR HELIOCENTRIC PLACES OF A HEAVENLY BODY AND THE VARIATIONS OF THE ELEMENTS OF ITS ORBIT. PAGB Variation of the Eight Ascension and Declination 118 Case of Parabolic Motion 125 Case of Hyperbolic Motion 128 Case of Orbits differing but little from the Parabola 130 Numerical Examples 135 Variation of the Longitude and Latitude 143 The Elements referred to the same Fundamental Plane as the Geocentric Places 149 Numerical Example 150 Plane of the Orbit taken as the Fundamental Plane to which the Geocentric Places are referred 153 Numerical Example 159 Variation of the Auxiliaries for the Equator 163 CHAPTER III. INVESTIGATION OF FORMULA FOR COMPUTING THE ORBIT OF A COMET MOVING IN A PARABOLA, AND FOR CORRECTING APPROXIMATE ELEMENTS BY THE VARIATION OF THE GEOCENTRIC DISTANCE. Correction of the Observations for Parallax 167 Fundamental Equations 169 Particular Cases 172 Ratio of Two Curtate Distances 178 Determination of the Curtate Distances 181 Relation between Two Radii- Vectores, the Chord joining their Extremities, and the Time of describing the Parabolic Arc 184 Determination of the Node and Inclination 192 Perihelion Distance and Longitude of the Perihelion 194 Time of Perihelion Passage 195 Numerical Example 199 Correction of Approximate Elements by varying the Geocentric Distance 208 Numerical Example 213 CHAPTER IV. DETERMINATION, FROM THREE COMPLETE OBSERVATIONS, OF THE ELEMENTS OF THE ORBIT OF A HEAVENLY BODY, INCLUDING THE ECCENTRICITY OR FORM OF THE CONIC SECTION. Reduction of the Data 220 Corrections for Parallax 223 CONTENTS. ll PAOK Fundamental Equations 225 Formulae for the Curtate Distances 228 Modification of the Formulae in Particular Cases 231 Determination of the Curtate Distance for the Middle Observation 236 Case of a Double Solution 239 Posjtion indicated by the Curvature of the Observed Path of the Body 242 Formulae for a Second Approximation 243 Formulas for finding the Katio of the Sector to the Triangle 247 Final Correction for Aberration 257 Determination of the Elements of the Orbit 259 Numerical Example 264 Correction of the First Hypothesis 278 Approximate Method of finding the Katio of the Sector to the Triangle 279 CHAPTER V. DETERMINATION OF THE ORBIT OP A HEAVENLY BODY FROM FOUR OBSERVA- TIONS, OF WHICH THE SECOND AND THIRD MUST BE COMPLETE. Fundamental Equations 282 Determination of the Curtate Distances , 289 Successive Approximations 293 Determination of the Elements of the Orbit 294 Numerical Example 294 Method for the Final Approximation 307 CHAPTER VI. INVESTIGATION OF VARIOUS FORMULAE FOR THE CORRECTION OF THE APPROXI- MATE ELEMENTS OF THE ORBIT OF A HEAVENLY BODY. Determination of the Elements of a Circular Orbit 311 Variation of Two Geocentric Distances 313 Differential Formulae 318 Plane of the Orbit taken as the Fundamental Plane 320 Variation of the Node and Inclination 324 Variation of One Geocentric Distance 328 Determination of the Elements of the Orbit by means of the Co-ordinates and Velocities 332 Correction of the Ephemeris 335 Final Correction of the Elements 338 Kelation between Two Places in the Orbit 339 Modification when the Semi-Transverse Axis is very large 341 Modification for Hyperbolic Motion 346 Variation of the Semi-Transverse Axis and Eatio of Two Curtate Distances 349 12 CONTENTS. PAQB Variation of the Geocentric Distance and of the Reciprocal of the Semi-Trans- verse Axis 352 Equations of Condition 353 Orbit of a Comet 355 Variation of Two Eadii-Vectores 357 CHAPTER VII. METHOD OF LEAST SQUARES, THEORY OP THE COMBINATION OF OBSERVATIONS, AND DETERMINATION OF THE MOST PROBABLE SYSTEM OF ELEMENTS FROM A SERIES OF OBSERVATIONS. Statement of the Problem 360 Fundamental Equations for the Probability of Errors 362 Determination of the Form of the Function which expresses the Probability ... 363 The Measure of Precision, and the Probable Error 366 Distribution of the Errors 367 The Mean Error, and the Mean of the Errors 368 The Probable Error of the Arithmetical Mean 370 Determination of the Mean and Probable Errors of Observations 371 Weights of Observed Values 372 Equations of Condition 376 Normal Equations 378 Method of Elimination 380 Determination of the Weights of the Resulting Values of the Unknown Quanti- ties 386 Separate Determination of the Unknown Quantities and of their Weights 392 Eelation between the Weights and the Determinants 396 Case in which the Problem is nearly Indeterminate 398 Mean and Probable Errors of the Results 399 Combination of Observations 401 Errors peculiar to certain Observations 408 Rejection of Doubtful Observations 410 Correction of the Elements 412 Arrangement of the Numerical Operations 415 Numerical Example 418 Case of very Eccentric Orbits 423 CHAPTER VIII. INVESTIGATION OF VARIOUS FORMULA FOR THE DETERMINATION OF THE SPECIAL PERTURBATIONS OF A HEAVENLY BODY. Fundamental Equations 426 Statement of the Problem '. 428 Variation of Co-ordinates 429 CONTENTS. 13 PAQB Mechanical Quadrature 433 The Interval for Quadrature 443 Mode of effecting the Integration 445 Perturbations depending on the Squares and Higher Powers of the Masses 446 Numerical Example 448 Change of the Equinox and Ecliptic 455 Determination of New Osculating Elements 459 Variation of Polar Co-ordinates 462 Determination of the Components of the Disturbing Force 467 Determination of the Heliocentric or Geocentric Place 471 Numerical Example 474 Change of the Osculating Elements 477 Variation of the Mean Anomaly, the Radius- Vector, and the Co-ordinate 2...... 480 Fundamental Equations 483 Determination of the Components of the Disturbing Force 489 Case of very Eccentric Orbits 493 Determination of the Place of the Disturbed Body 495 Variation of the Node and Inclination 502 Numerical Example 505 Change of the Osculating Elements 510 Variation of Constants 516 Case of very Eccentric Orbits 523 Variation of the Periodic Time 526 Numerical Example 529 Formulae to be used when the Eccentricity or the Inclination is small 533 Correction of the Assumed Value of the Disturbing Mass 535 Perturbations of Comets 536 Motion about the Common Centre of gravity of the Sun and Planet 537 Keduction of the Elements to the Common Centre of Gravity of the Sun and Planet 538 Keduction by means of Differential Formulae 540 Near Approach of a Comet to a Planet 546 The Sun may be regarded as the Disturbing Body 548 Determination of the Elements of the Orbit about the Planet 550 Subsequent Motion of the Comet 551 Effect of a Resisting Medium in Space 552 Variation of the Elements on account of the Resisting Medium 554 Method to be applied when no Assumption is made in regard to the Density o* the Ether.... 5*6 14 CONTENTS. TABLES. PAOl I. Angle of the Vertical and Logarithm of the Earth's Radius ~ 561 H. For converting Intervals of Mean Solar Time into Equivalent Intervals of Sidereal Time 563 III. For converting Intervals of Sidereal Time into Equivalent Intervals of Mean Solar Time 564 IV. For converting Hours, Minutes, and Seconds into Decimals of a Day... 565 V. For finding the Number of Days from the Beginning of the Year 565 VI. For finding the True Anomaly or the Time from the Perihelion in a Parabolic Orbit 566 VII. For finding the True Anomaly in a Parabolic Orbit when v is nearly 180 611 VIII. For finding the Time from the Perihelion in a Parabolic Orbit 612 IX. For finding the True Anomaly or the Time from the Perihelion in Orbits of Great Eccentricity 614 X. For finding the True Anomaly or the Time from the Perihelion in El- liptic and Hyperbolic Orbits 618 XI. For the Motion in a Parabolic Orbit 619 XII. For the Limits of the Boots of the Equation sin (z 1 )=m sin 4 y! ... 622 XIII. For finding the Eatio of the Sector to the Triangle 624 XIV. For finding the Eatio of the Sector to the Triangle 629 XV. For Elliptic Orbits of Great Eccentricity 632 XVI. For Hyperbolic Orbits 632 XVII. For Special Perturbations 633 XVIII. Elements of the Orbits of the Comets which have been observed 638 XIX. Elements of the Orbits of the Minor Planets 646 XX. Elements of the Orbits of the Major Planets 648 XXI. Constants, &c 649 EXPLANATION OP THE TABLES 651 APPENDIX. Precession 657 Nutation 658 Aberration 659 Intensity of Light 660 Numerical Calculations .. 662 THEORETICAL ASTRONOMY. CHAPTER I. INVESTIGATION OF THE FUNDAMENTAL EQUATIONS OF MOTION, AND OF THE FOR- MULA FOB DETERMINING, FROM KNOWN ELEMENTS, THE HELIOCENTRIC AND GEOCENTRIC PLACES OF A HEAVENLY BODY, ADAPTED TO NUMERICAL COMPUTA- TION FOR CASES OF ANY ECCENTRICITY WHATEVER. 1. THE study of the motions of the heavenly bodies does not re- quire that we should know the ultimate limit of divisibility of the matter of which they are composed, whether it may be subdivided indefinitely, or whether the limit is an indivisible, impenetrable atom. Nor are we concerned with the relations which exist between the separate atoms or molecules, except so far as they form, in the aggre- gate, a definite body whose relation to other bodies of the system it is required to investigate. On the contrary, in considering the ope- ration of the laws in obedience to which matter is aggregated into single bodies and systems of bodies, it is sufficient to conceive simply of its divisibility to a limit which may be regarded as infinitesimal compared with the finite volume of the body, and to regard the mag- nitude of the element of matter thus arrived at as a mathematical point. An element of matter, or a material body, cannot give itself motion; neither can it alter, in any manner whatever, any motion which may have been communicated to it. This tendency of matter to resist all changes of its existing state of rest or motion is known as inertia, and is the fundamental law of the motion of bodies. Ex- perience invariably confirms it as a law of nature; the continuance of motion as resistances are removed, as well as the sensibly unchanged motion of the heavenly bodies during many centuries, affording the 16 16 THEORETICAL ASTRONOMY. most convincing proof of its universality. Whenever, therefore, a material point experiences any change of its state as respects rest or motion, the cause must be attributed to the operation of something external to the element itself, and which we designate by the word force. The nature of forces is generally unknown, and we estimate them by the effects which they produce. They are thus rendered com- parable with some unit, and may be expressed by abstract numbers. 2. If a material point, free to move, receives an impulse by virtue of the action of any force, or if, at any instant, the force by which motion is communicated shall cease to act, the subsequent motion of the point, according to the law of inertia, must be rectilinear and uniform, equal spaces being described in equal times. Thus, if s, v, and t represent, respectively, the space, the velocity, and the time, the measure of v being the space described in a unit of time, we shall have, in this case, s = vt. It is evident, however, that the space described in a unit of time will vary with the intensity of the force to which the motion is due, and, the nature of the force being unknown, we must necessarily compare the velocities communicated to the point by different forces, in order to arrive at the relation of their effects. We are thus led to regard the force as proportional to the velocity; and this also has received the most indubitable proof as being a law of nature. Hence, the principles of the composition and resolution of forces may be applied also to the composition and resolution of velocities. If the force acts incessantly, the velocity will be accelerated, and the force which produces this motion is called an accelerating force. In regard to the mode of operation of the force, however, we may consider it as acting absolutely without cessation, or we may regard it as acting instantaneously at successive infinitesimal intervals repre- sented by dt, and hence the motion as uniform during each of these intervals. The latter supposition is that which is best adapted to the requirements of the infinitesimal calculus; and, according to the fundamental principles of this calculus, the finite result will be the same as in the case of a force whose action is absolutely incessant. Therefore, if we represent the element of space by ds, and the ele- ment of time by dt, the instantaneous velocity will be ds - V = ~dt> which will vary from one instant to another. FUNDAMENTAL PRINCIPLES. 17 3. Since the force is proportional to the velocity, its measure at any instant will be determined by the corresponding velocity. If the accelerating force is constant, the motion will be uniformly accele- rated; and if we designate the acceleration due to the force by/, the unit of/ being the velocity generated in a unit of time, we shall have If, however, the force be variable, we shall have, at any instant, the relation efc / ~~ dt' the force being regarded as constant in its action during the element of time dt. The instantaneous value of v gives, by differentiation, dv _ d*s ~dt~~W and hence we derive so that, in varied motion, the acceleration due to the force is mea- sured by the second differential of the space divided by the square of the element of time. 4. By the mass of the body we mean its absolute quantity of mat- ter. The density is the mass of a unit of volume, and hence the entire mass is equal to the volume multiplied by the density. If it is required to compare the forces which act upon different bodies, it is evident that the masses must be considered. If equal masses receive impulses by the action of instantaneous forces, the forces acting on each will be to each other as the velocities imparted; and if we consider as the unit of force that which gives to a unit of mass the unit of velocity, we have for the measure of a force F, denoting the mass by Jf, F = Mv. This is called the quantity of motion of the body, and expresses its capacity to overcome inertia. By virtue of the inert state of matter, there can be no action of a force without an equal and contrary re- action ; for, if the body to which the force is applied is fixed, the equilibrium between the resistance and the force necessarily implies the development of an equal and contrary force ; and, if the body be free to move, in the change of state, its inertia will oppose equal and 18 THEORETICAL ASTRONOMY. contrary resistance. Hence, as a necessary consequence of inertia, it follows that action and reaction are simultaneous, equal, and contrary. If the body is acted upon by a force such that the motion is varied, the accelerating force upon each element of its mass is represented by "7- , and the entire motive force F is expressed by M being the sum of all the elements, or the mass of the body. Since _ds_ this gives which ^is the expression for the intensity of the motive force, or of the force of inertia developed. For the unit of mass, the measure of the force is d?s and this, therefore, expresses that part of the intensity of the motive force which is impressed upon the unit of mass, and is what is usually called the accelerating force. 5. The force in obedience to which the heavenly bodies perform their journey through space, is known as the attraction of gravitation; and the law of the operation of this force, in itself simple and unique, has been confirmed and generalized by the accumulated researches of modern science. Not only do we find that it controls the motions of the bodies of our own solar system, but that the revolutions of binary systems of stars in the remotest regions of space proclaim the uni- versality of its operation. It unfailingly explains all the phenomena observed, and, outstripping observation, it has furnished the means of predicting many phenomena subsequently observed. The law of this force is that every particle of matter is attracted by every other particle by a force which varies directly as the mass and inversely as the square of the distance of the attracting particle. This reciprocal action is instantaneous, and is not modified, in any degree, by the interposition of other particles or bodies of matter. It is also absolutely independent of the nature of the molecules them- selves, and of their aggregation. ATTRACTION OF SPHERES. 19 If we consider two bodies the masses of which are m and m', and whose magnitudes are so small, relatively to their mutual distance />, that we may regard them as material points, according to the law of gravitation, the action of m on each molecule or unit of m' will be , and the total force on m! will be m'. />' The action of m f on each molecule of m will be expressed by , and its total action by m' -r- The absolute or moving force with which the masses m and m f tend toward each other is, therefore, the same on each body, which result is a necessary consequence of the equality of action and reaction. The velocities, however, with which these bodies would approach each other must be different, the velocity of the smaller mass exceed- ing that of the greater, and in the ratio of the masses moved. The expression for the velocity of m', which would be generated in a unit of time if the force remained constant, is obtained by dividing the absolute force exerted by m by the mass moved, which gives . 7 and this is, therefore, the measure of the acceleration due to the action of m at the distance p. For the acceleration due to the action of m r we derive, in a similar manner, 6. Observation shows that the heavenly bodies are nearly spherical in form, and we shall therefore, preparatory to finding the equations which express the relative motions of the bodies of the system, de- termine the attraction of a spherical mass of uniform density, or varying from the centre to the surface according to any law, for a point exterior to it. If we suppose a straight line to be drawn through the centre of the sphere and the point attracted, the total action of the sphere on the point will be a force acting along this line, since the mass of the sphere is symmetrical with respect to it. Let dm denote an element 20 THEORETICAL ASTRONOMY. of the mass of the sphere, and p its distance from the point attracted ; then will dm express the action of this element on the point attracted. If we sup- pose the density of the sphere to be constant, and equal to unity, the element dm becomes an element of volume, and will be expressed by dm = dx dy dz ; x y y. and z being the co-ordinates of the element referred to a system of rectangular co-ordinates. If we take the origin of co-ordinates at the centre of the sphere, and introduce polar co-ordinates, so that x = r cos
*, with respect to a, gives dp a r sin dr d , we have
7 P 7
r cos = e;
and, since a(l e 2 ) is half the parameter of the transverse axis,
which we have designated by p, we have .
r
1 -}- sin (f> cos v
The angle ) tan |v. (45)
Again.
V\-\- e = 1/1 -f sin y = 1/1 -|- 2 sin p cos fop,
which may be written
1/1 -}- e I/sin 2 -|^ -f c s 2 i^ + 2 sin Jp cos Jp,
or _
1/1 -|- e = sin ^ -f cos ?
In a similar manner we find
1/1 e = sin \ , and 6 = a cos ^, hence
r sin v cos ^ sin v
. _,
sin E =
1 -f- e GOBV
Equation (41) gives
r cos v + ae p cos v
or
_, ^) cos v + ae ag2 cos v
and, putting a cos 2 (p instead of _p, and sin (sin 8 Jf(l e cos IT) )
dM v
Expanding these terms, and performing the operations indicated, we
get
^ = 1 + 2e cos M + - (6 xjos 2 M 4 sin 2 Jf)
-f | (16 cos 8 M 36 sin 2 Jf cos Jf) -f- . . . ,
PLACE IN THE ORBIT. 57
which reduces to
, (51)
i
Equation (22) gives
. 2fdt
=
and, since /= }&|/p(l + m), we have
dv = ^j- dt, (52)
or ___
, &V 1 -f- Wl Cl 2 /
But g = fa and therefore
a?
2 2 CX
By expanding the factor j/1 e 2 , we obtain
and hence
Substituting for its value from equation (51), and integrating, we
get, since v = when M=0,
vM=2e sin M +je 2 sin2Jf + -^- (13 sin 33f 3 sin Jf ) +. . . (53)
which is the expression for the equation of the centre to terms involving
e 3 . In the same manner, this series may be extended to higher powers
of e.
When the eccentricity is very small, this series converges very
rapidly ; and the value of v M for any planet may be arranged in
a table with the argument M.
For the purpose, however, of computing the places of a heavenly
body from the elements of its orbit, it is preferable to solve the
equations which give v and E directly ; and when the eccentricity is
58 THEOEETICAL ASTRONOMY.
very great, this mode is indispensable, since the series will not in
that case be sufficiently convergent.
It will be observed that the formula which must be used in obtain-
ing the eccentric anomaly from the mean anomaly is transcendental,
and hence it can only be solved either by series or by trial. But
fortunately, indeed, it so happens that the circumstances of the celes-
tial motions render these approximations very rapid, the orbits being
usually either nearly circular, or else very eccentric.
If, in equation (50), we put F(E] = E, and consequently F(M]
= M, we shall have, performing the operations indicated and reducing,
E = M -f e sin M + ^ sin 2 M + &c. (54)
Let us now denote the approximate value of E computed from this
equation by E ot then will
in which &E Q is the correction to be applied to the assumed value of E.
Substituting this in equation (39), we get
M= E -f &E Q e sin E Q e cos E *E ;
and, denoting by M Q the value of M corresponding to E w we shall
also have
M Q = E Q esmE .
Subtracting this equation from the preceding one, we obtain
M-M,
- -- = = A A..
1 e cos -c/
It remains, therefore, only to add the value of &E Q found from this
formula to the first assumed value of E, or to E w and then, using
this for a new value of E , to proceed in precisely the same manner
for a second approximation, and so on, until the correct value of E is
obtained. When the values of E for a succession of dates, at equal
intervals, are to be computed, the assumed values of E may be ob-
tained so closely by interpolation that the first approximation, in the
manner just explained, will give the correct value; and in nearly
every case two or three approximations in this manner will suffice.
Having thus obtained the value of E corresponding to M for any
instant of time, we may readily deduce from it, by the formulse
already investigated, the corresponding values of r and v.
In the case of an ellipse of very great eccentricity, corresponding
to the orbits of many of the comets, the most convenient method of
PLACE IN THE CEBIT. 59
computing r and v, for any instant, is somewhat different. The
manner of proceeding in the computation in such cases we shall con-
sider hereafter; and we will now proceed to investigate the formulae
for determining r and v, when the orbit is a parabola, the formulse
for elliptic motion not being applicable, since, in the parabola, a = ' sin B' = A y ,
a sin b f sin _B' = ^t y ,
a tan 4> sin b' cos B f = v y ;
ae sin c' sin C" = >*
a sin c' sin G r = /*
a tan 4 sin c' cos C" = v z .
Then we shall have
x = A x -f fj. K sec .F + v x tan jP,
v = A y 4- ^ sec jF + v y tan F, (106;
2 = A s -j- /7. z sec F -f- v z tan jP.
In a similar manner we may derive expressions for the co-ordinates,
in the case of a hyperbolic orbit, when the auxiliary quantity a is
used instead of F.
39. If we denote by ;:', ', and V the elements which determine
the position of the orbit in space when referred to the equator as the
96 THEORETICAL ASTRONOMY.
fundamental plane, and by o) the angular distance between the
ascending node of the orbit on the ecliptic and its ascending node on
the equator, being measured positively from the equator in the
direction of the motion, we shall have
To find Q,' and i', we have, from the spherical triangle formed by
the intersection of the planes of the orbit, ecliptic, and equator with
the celestial vault,
cos if = cos i cos e sin i sin e cos & ,
sin i r sin &' = sin i sin ^ ,
sin i' cos &' = cos t sin e -f- sin i cos e cos SI
Let us now put
n sin JV = cos i t
n cos JV= sin i cos &,
and these equations reduce to
cos i' = n sin ( JV e),
sin i' sin &' = sin i sin & ,
sin i' cos &' = n cos (JV e) ;
from which we find
HT cot i cos JV
tan JV= , tan ' = ^ ^ tan &,
cos 1 cos (JV e)
cot i' = tan (JV e) cos &'. ( 107 )
Since sin i is always positive, cos JV and cos & must have the same
signs. To prove the numerical calculation, we have
sin i cos Q _ cos JV
sin i' cos ' cos (JV e)'
the value of the second member of which must agree with that used
in computing & '.
In order to find CD O , we have, from the same triangle,
sin CM O sin i' = sin & sin e,
cos a> n sin i' = cos e sin i -{- sin e cos i cos & .
Liet us now take
m sin Jlf = cos e,
m cos Jlf = sin e cos & ;
and we obtain
POSITION IN SPACE. 97
cot M = tan e cos &,
and, also, to check the calculation,
sin e cos & cos Jf
sin i' cos % cos (M i)'
If we apply Gauss's analogies to the same spherical triangle, we
get
cos^' cos^ (' + > ) == cosi cosi (*' + s \ MAQ x
sin It sin I (&' ) sin J& sin i (t e),
The quadrant in which \ (&' + ^ ) or K^ ^o) ^ s situated, must be
so taken that sin \i f and cos \V shall be positive ; and the agreement
of the values of the latter two quantities, computed by means of the
value of \i f derived from tan Ji', will serve to check the accuracy of
the numerical calculation.
For the case in which the motion is regarded as retrograde, we
must use 180 i instead of i in these equations, and we have, also,
We may thus find the elements TT', & ', and i f y in reference to the
equator, from the elements referred to the ecliptic; and using the
elements so found instead of it, &, and i, and using also the places
of the sun referred to the equator, we may derive the heliocentric
and geocentric places with respect to the equator by means of the
formulae already given for the ecliptic as the fundamental plane.
If the position of the orbit with respect to the equator is given,
and its position in reference to the ecliptic is required, it is only
necessary to interchange & and &', as well as i and 180 i', e
remaining unchanged, in these equations. These formulae may
also be used to determine the position of the orbit in reference to
any plane in space; but the longitude & must then be measured
from the place of the descending node of this plane on the ecliptic.
The value of &, therefore, which must be used in the solution of the
equations is, in this case, equal to the longitude of the ascending
node of the orbit on the ecliptic diminished by the longitude of the
descending node of the new plane of reference on the ecliptic. The
quantities & ', i', and w will have the same signification in reference
7
98 THEORETICAL ASTRONOMY.
to this plane that they have in reference to the equator, with this dis-
tinction, however, that & ' is measured from the descending node of
this new plane of reference on the ecliptic ; and e will in this case
denote the inclination of the ecliptic to this plane.
40. We have now derived all the formulae which can be required
in the case of undisturbed motion, for the computation of the helio-
centric or geocentric place of a heavenly body, referred either to the
ecliptic or equator, or to any other known plane, when the elements
of its orbit are known ; and the formulae which have been derived
are applicable to every variety of conic section, thus including all
possible forms of undisturbed orbits consistent with the law of uni-
versal gravitation. The circle is an ellipse of which the eccentricity
is zero, and, consequently, M=v = u, and r = a, for every point of
the orbit. There is no instance of a circular orbit yet known ; but
in the case of the discovery of the asteroid planets between Mars
and Jupiter it is sometimes thought advisable, in order to facilitate
the identification of comparison stars for a few days succeeding the
discovery, to compute circular elements, and from these an ephemeris.
The elements which determine the form of the orbit remain con-
stant so long as the system of elements is regarded as unchanged ;
but those which determine the position of the orbit in space, TT, & ,
and i, vary from one epoch to another on account of the change of
the relative position of the planes to which they are referred. Thus
the inclination of the orbit will vary slowly, on account of the change
of the position of the ecliptic in space, arising from the perturbations
of the earth by the other planets ; while the longitude of the peri-
helion and the longitude of the ascending node will vary, both on
account of this change of the position of the plane of the ecliptic,
and also on account of precession and nutation. If TT, &, and i are
referred to the true equinox and ecliptic of any date, the resulting
heliocentric places will be referred to the same equinox and ecliptic ;
and, further, in the computation of the geocentric places, the longi-
tudes of the sun must be referred to the same equinox, so that the
resulting geocentric longitudes or right ascensions will also be re-
ferred to that equinox. It will appear, therefore, that, on account
of these changes in the values of n, &, and i, the auxiliaries sin a,
sin 6, sin c, A, J5, and C, introduced into the formulae for the co-
ordinates, will not be constants in the computation of the places for
a series of dates, unless the elements are referred constantly, in the
calculation, to a fixed equinox and ecliptic. It is customary, there-
POSITION IN SPACE. 99
fore, to reduce the elements to the ecliptic and mean equinox of the
beginning of the year for which the ephemeris is required, and then
to compute the places of the planet or comet referred to this equinox,
using, in the case of the right ascension and declination, the mean
obliquity of the ecliptic for the date of the fixed equinox adopted, in
the computation of the auxiliary constants and of the co-ordinatey
of the sun. The places thus found may be reduced to the true
equinox of the date by the well-known formulae for precession and
nutation. Thus, for the reduction of the right ascension and declina-
tion from the mean equinox and equator of the beginning of the
year to the apparent or true equinox and equator of any date, usually
the date to which the co-ordinates of the body belong, we have
A<* =/+ g sin (G + a) tan d,
for which the quantities /, g, and G are derived from the data given
either in the solar and lunar tables, or in astronomical ephemerides,
such as have already been mentioned.
The problem of reducing the elements from the ecliptic of one
date t to that of another date t f may be solved by means of equations
(109), making, however, the necessary distinction in regard to the
point from which Q> and &>' are measured. Let d denote the longi-
tude of the descending node of the ecliptic of t' on that of t y and
let T} denote the angle which the planes of the two ecliptics make
with each other, then, in the equations (109), instead of & we must
write & 0, and, in order that Q, ' shall be measured from the
vernal equinox, we must also write & ' 6 in place of Q, ' . Finally,
we must write T] instead of e, and A = 37 35' 0".0, or log 6 = 0.1010188; and lug a
= 0.6020600, to find r and v. First, we compute ^Vfrom
" at
in which log A = 9.6377843, and we obtain
logN= 8.7859356; N= 0.06108514.
The value of F must now be found from the equation
N= el tan F log tan (45 + F).
NUMERICAL EXAMPLES. 109
If we assume F= 30, a more approximate value may be derived
from
which gives F, = 28 40' 23", and hence N, = 0.072678. Then we
compute the correction to be applied to this value of F, by means of
the equation
(N-N,^o S *F
' l(e cosF,-) *'
wherein s = 206264".8; and the result is
*F, = 4.6097 (N N f ~) s = 3 3' 43".0.
Hence, for a second approximation to the value of F } we have
.F, = 25 36'40".0.
The corresponding value of Nis N, 0.0617653, and hence
AJ F, = 5.199 (N NJ s = 12' 9".4.
The third approximation, therefore, gives F, = 25 24' 30".6, and,
repeating the operation, we get
J F T =2524'27".74.
which requires no further correction.
To find r, we have
which gives
log r = 0.2008544.
Then, v is derived from
tan v cot ^4 tan F,
and we find
v = 67 3' 0".0.
When several places are required, it is convenient to compute
and r by means of the equations
110 THEORETICAL ASTRONOMY.
For the given values of a and e we have log V a(e + 1) = 0.4782649,
logl / a(e 1) = 0.0100829, and hence we derive
v = 67 2' 59".92, log r = 0.2008545.
It remains yet to illustrate the calculation of v and r for elliptic
and hyperbolic orbits in which the eccentricity differs but little from
unity. First, in the case of elliptic motion, let t T= 68.25 days;
e = 0.9675212; and log q = 9.7668134. We compute M from
wherein log C = 9.9601277, which gives
log 3f= 2.1404550.
\Yith this as argument we get, from Table VI.,
V= 101 38' 3".74,
and then with this value of V as argument we find, from Table IX.,
A --= 1540".08, B = 9".506, C= 0".062.
^ _ e
Then we have log i = log = 8.217680, and from the equation
1 -f- e
v = V+ -4(1000
we get
v = V+ 42' 22".28 -f 25".90 -f 0".28 = 102 20' 52".20.
The value of r is then found from
-f e cos v '
namely,
log r = 0.1614051.
We may also determine r and v by means of Table X. Thus, we
first compute M from
Assuming B = 1, we get log M= 2.13757, and, entering Table VI.
with this as argument, we find w = 101 25'. Then we compute A
from
5(l--e
A ~ * 1 9e
NUMERICAL EXAMPLES. Ill
which gives A = 0.024985. With this value of A as argument, we
find, from Table X.,
log B = 0.0000047.
The exact value of M is then found to be
log M= 2.1375635,
which, by means of Table VI., gives
w = 101 24' 36".26.
By means of this we derive
A = 0.02497944,
and hence, from Table X.,
log C =0.0043771.
Then we have
which gives
v = 102 20' 52".20,
agreeing exactly with the value already found. Finally, r is given by
" (1 -{-AC 2 ) cos 1 }*'
from which we get
log r = 0.1614052.
Before the time of perihelion passage, t T is negative ; but the
value of v is computed as if this were positive, and is then considered
as negative.
In the case of hyperbolic motion, i is negative, and, with this dis-
tinction, the process when Table IX. is used is precisely the same
as for elliptic motion; but when table X. is used, the value of A
must be found from
5(e 1)
^(f+
and that of r from
(1 .AC 2 ) cos 8 in'
the values of log B and log C being taken from the columns of the
table which belong to hyperbolic motion.
In the calculation of the position of a comet in space, if the motion
112 THEORETICAL ASTRONOMY.
is retrograde and the inclination is regarded as less than 90, the dis-
tinctions indicated in the formulae must be carefully noted.
42. When we have thus computed the places of a planet or comet
for a series of dates equidistant, we may readily interpolate the places
for intermediate dates by the usual formulae for interpolation. The
interval between the dates for which the direct computation is made
should also be small enough to permit us to neglect the effect of the
fourth differences in the process of interpolation. This, however, is
not absolutely necessary, provided that a very extended series of
places is to be computed, so that the higher orders of differences may
be taken into account. To find a convenient formula for this inter-
polation, let us denote any date, or argument of the function, by
a + na>> and the corresponding value of the co-ordinate, or of the
function, for which the interpolation is to be made, by / (a -f nco).
If we have computed the values of the function for the dates, or
arguments, a to, a, a -f a>, a -f- 2w, &c., we may assume that an
expression for the function which exactly satisfies these values will
also give the exact values corresponding to any intermediate value
of the argument. If we regard n as variable, we mav expand the
function into the series
f(a + no>) =/() + An + JSn 9 + On 9 -f &c. (116)
and if we regard the fourth differences as vanishing, it is only neces-
sary to consider terms involving n s in the determination of the
unknown coefficients A, B, and C. If we put n successively equal
to 1, 0, 1, and 2, and then take the successive differences of these
values, we get
I. Diff. II. Diff. III. Diff.
/(<,-) =f(a)-A + B -C
f(a + 2>) =/(o) + 2A + 4B + 8 C
If we symbolize, generally, the difference f(a + nco) f(a -\-(n 1 ) , Jf , and //, are
dx
d and have been
aM d/j. a/j.
found, the partial differential coefficients of the heliocentric co-ordj
nates with respect to the elements tp, M w and // will be completely
determined, and hence, by means of (2), making the necessary
changes, the differential coefficients of a and d with respect to these
elements.
46. If we differentiate the equation
M=E
DIFFERENTIAL FORMULA. 123
we shall have
dM= dE(l e cos E) cos we get
J a cosy
a? cos y j , , . sin v /
dv = - - ---
sn v a cos
2 ^ . \ ,
- + \\d sin v dM a cos . (14)
Further, we have
T being the epoch for which the mean anomaly is M w and
A; 1/1 + m
^= - 3 --
a^
Differentiating these expressions, we get
T / a sin v tan 4 sin v \ ,
ai? = ^,dJ\ n I 1 h : I 4'
r tan F \ r cos ^ sin ^ /
DIFFERENTIAL FORMULA. 129
Bufc, since r sin v = a tan $ tan F, and p a tan 2 ^ this reduces to
* = ^v5i-(| + 1 )^d*. (25)
If we differentiate the equation
we get
r , , T, d-F 7 a tan 4* T
dr = - da -f ae tan* F ^ H -- ^ -- cfy.
a smF cosF cos4-
Trrr
Substituting in this equation the value of ^ we obtain
7 r , , a?e tan F 7 . _ / a 2 e tan 2 jP a \ tan 4*
dr = - da -\ -- dN tt I ---- =f )
a r \ r cos jP/ cos 4
which is easily reduced to
r 7 sinv 7 , T . pi r ae
dr = a da+a ^ dN + -r(^-^
But, since
r ae a
cos F ~~ cos 2 F ~~ cos F'
this reduces to
r
or
7 r , sin v , , r . cos v , , n
-7-, and -7-
dr dv dr dv
in connection with the numerical values last found, according to the
NUMERICAL, EXAMPLES. 137
equations which result from the analytical substitution of the expres-
sions for -j-, -j > -y-, &c., in equation (2), writing successively
and fJL in place of TT. Thus, we have
. da, .da,dr. .da, dv
cos <5 - = cos d . -- k cos <5 -r- -r ,
a?' ar* a^ av a?>
d<5 __ d dr dd dv
d
+ 1.13004A3f -}- 507.264A/.,
A5 = 0.48900 ATT 0.09533 A & 0.78993Ai 0.65307 A
0.38023A^f 179.315A/Z.
To prove the calculation of the coefficients in these equations, we
assign to the elements the increments
Ajf = + 10", ATT == - 20", A ft = - 10", Ai = + 10",
A^ = + 10", A/ = + 0".01,
BO that they become
Epoch = 1864 Jan. 1.0 Greenwich mean time.
M = 1 29' 50".21
TT = 44 20 13 .09 ^
ft = 206 42 30 .13 > Mean Equinox 1864.0
t = 4 37 .51 J
V = 11 16 1 .02
log a = 0.3881288
L = 928".56745
With these elements we compute the geocentric place for 1865 Feb-
ruary 24.5 mean time at Washington ; and the result is
a = 181 8' 34".81, d = - 4 42' 30".58, log A = 0.2450284,
138 THEORETICAL ASTRONOMY.
which are referred to the mean equinox and equator of 1865.0. The
difference between these values of a and d and those already given, as
derived from the unchanged elements, gives
Aa = + 5".52, COS * Aa = + 5".50, A sin v,
dv ~ a 2 cos , we must put
=r-,, (63)
-- -j-, - 5 -
aq aq
and the equations become
= sn
S 9 (64)
ooi^-^iiaU + Jr+ii),
4^ = \ B sin (jB H- H+ u).
In the numerical application of these formulae, the values of the
second members of the equations (63) are found as already exem-
plified for the cases of parabolic orbits and of elliptic and hyperbolic
orbits in which the eccentricity differs but little from unity. In the
same manner, the differential coefficients of A and ft with respect to
any other elements which determine the form of the orbit may be
computed.
NUMERICAL EXAMPLES. 149
In the case of a parabolic orbit, if the parabolic eccentricity is
supposed to be invariable, the terms involving e vanish. Further,
in the case of parabolic elements, we have
. ~ dr ksmv dv
= --
which give
tan G = tan
V2
-, which is the
expression for the linear velocity of a comet moving in a parabola.
Therefore,
[ ' dv dy
df) _ df) dr df) dv
d ? -- = ?cos(0 u G),
,, , i I (79)
dO h , dr) h . f _.
cos r) = cos (0 u L ). ^ = -j- sm TJ sm (0 u M ).
d/j- A dp A
If we express r and v in terms of the elements T, g, and e, the
values of the auxiliaries /, g, A, F, &c. must be found by means of
(64) ; and, in the same manner, any other elements which determine
the form of the orbit and the position of the body in its orbit, may
be introduced.
The partial differential coefficients with respect to the elements
having been found, we have
do do do dO
COS JJ A0 = COS r) -- Ay 4- COS f] -7 ACP -f- COS r] -rTT A M n + COS Tt -7- &fJ.
d% dy dM Q dp
DIFFERENTIAL FORMULA. 157
from which it appears that, by the introduction of as one of the
elements of the orbit, when the geocentric places are referred directly
to the plane of the unchanged orbit as the fundamental plane, the
variation of the geocentric longitude in reference to this plane depends
on only four elements.
57. It remains now to derive the formulae for finding the values
of T) and 6 from those of X and /?. Let X Q , y , z be the geocentric co-
ordinates of the body referred to a system in which the ecliptic is
the plane of xy, the positive axis of x being directed to the point
whose longitude is & ; and let #/, y ', z ' be the geocentric co-ordi-
nates of the body referred to a system in which the axis of x remains
the same, but in which the plane of the unchanged orbit is the plane
of xy; then we shall have
x A cos /? cos (A ft ), x ' = A cos i? cos 0,
y o = J cos /? sin (A Q ), y ' = A cos T? sin 0,
Z Q = J sin , z ' = A sin iy,
and also
< = *0
yo = y cos * -f z o sin *
z ' = y Q sin i + z cos i.
Hence we obtain
COS r) COS = COS /? COS (A & ),
cos T) sin = cos /? sin (A & ) cos i + sin sin t, (80 )
sin 17 = cos sin (A & ) sin i -\- sin /? cos i.
These equations correspond to the relations between the parts of a
spherical triangle of which the sides are i, 90 57, and 90 ft
the angles opposite to 90 ^ and 90 ft being respectively
90 -f- (X ^) and 90 0. Let the other angle of the triangle be
denoted by f, and we have
cos 17 sin Y = sin * cos (A & ),
cos T? cos r = sin i sin (A & ) sin ft -f- cos i cos /?.
(81)
The equations thus obtained enable us to determine y, 0, and f from
I and /9. Their numerical application is facilitated by the intro-
duction of auxiliary angles. Thus, if we put
, ,
n cos ^= cos j9 sin (A ft),
158 THEORETICAL ASTRONOMY.
in which n is always positive, we get
cos TI cos 6 = cos /? cos (A & ),
cos rj sin 6 = n cos (N t), (83)
sin 7] = ra sin (JV i\
from which ^ and may be readily found. If we also put
n' sin N' = cos i,
n' cos JV = sin i sin (A ft ),
we shall have
cot N r - tan i sin (A ^ ),
If f is small, it may be found from the equation
sin* cos (A -S
cos??
The quadrants in which the angles sought must be taken, are easily
determined by the relations of the quantities involved ; and the
accuracy of the numerical calculation may be checked as already
illustrated for similar cases.
If we apply Gauss's analogies to the same spherical triangle, we get
fiin (45 - j,) sin (45 - J (0 + r )) =
cos (45 -f i (A - )) sin (45 J ( + i)),
sin (45 - j?) cos (45 - J (* + r )) =
sin (45 + i (A - &)) sin (45 - J (0 0),
cos (45 J?) sin (45 j (0 r)) = (87)
cos (45 + i (A - 8)) cos (45 - i (p + t)),
cos (45 - ^) cos (45 - J ( - r )) =
sin (45 + C; - )) cos (45 - J (/? - i)),
from which we may derive ^, 0, and y.
When the problem is to determine the corrections to be applied to
the elements of the orbit of a heavenly body, in order to satisfy
given observed places, it is necessary to find the expressions for
cos 7] A# and A^ in terms of cos /9 AA and A/3. If we differentiate the
first and second of equations (80), regarding & and i (which here
determine the position of the fundamental plane adopted) as con-
stant, eliminate the terms containing dy from the resulting equations,
and reduce by means of the relations of the parts of the spheric!
triangle, we get
NUMERICAL EXAMPLE. 159
cos i) do = cos Y cos p dl -f si n F dp.
Differentiating the last of equations (80), and reducing, we find
df] = sin Y cos p ctt -f cos Y dp.
The equations thus derived give the values of the differential co-
efficients of 6 and y with respect to A and /9 ; and if the differences
A^ and A/9 are small, we shall have
cos rj A0 = cos Y cos p A/* -f sin y A/9,
A>? = sin Y cos /5 AA -f cos y A/1
The value of p required in the application of numbers to these
equations may generally be derived with sufficient accuracy from
(86), the algebraic sign of cosf being indicated by the second of
equations (81) ; and the values of rj and 6 required in the calculation
of the differential coefficients of these quantities with respect to the
elements of the orbit, need not be determined with extreme accuracy.
58. EXAMPLE. Since the spherical co-ordinates which are fur-
nished directly by observation are the right ascension and declina-
tion, the formulae will be most frequently required in the form for
finding rj and from a and d. For this purpose, it is only necessary
to write a and d in place of A and /9, respectively, and also & ', V,
co f , %', and u' in place of &, i, a), %, and u, in the equations which
have been derived for the determination of rj and 0, and for the
differential coefficients of these quantities with respect to the elements
of the orbit.
To illustrate this clearly, let it be required to find the expressions
for cos rj A# and &rj in terms of the variations of the elements in the
case of the example already given ; for which we have
/ = 50 10' 7".29, &' = 353 45' 35".87, i' = 19 26' 25".76.
These are the elements which determine the position of the orbit of
Evsrynome @, referred to the mean equinox and equator of 1865.0.
We have, further,
log/= 0.62946, log# = 0.34593, log h = 2.97759,
F= 339 14' 0", G = 350 11' 16", H= 14 30' 48' ; ,
u' = 179 13' 58".
In the first place, we compute ^, 0, and 7- by means of the formulae
160 THEORETICAL ASTRONOMY.
(83) and (85), or by means of (87), writing a, d, &', and V instead
of ^, ft, &, and i 9 respectively. Hence we obtain
e = 188 31' 9", T] = 1 59' 28", r = 19 17' 7".
Since the equator is here considered as the fundamental plane, the
longitude 6 is measured on the equator from the place of the ascend-
ing node of the orbit on this plane. The values of the differential
Coefficients are then found by means of the formulae
dO di) r
cos f) ^7 = 0, jr = - cos T) sm j cos tr,
de dy , r
cos TI -^r = 0, -jj- = 4 - cos rj sm u ,
cos 7] --r = T - cos (0 u') t -- = sin T? sin (0 if),
cos 7] ,,, = cos (e u' G\ ,,, = ^- sin -n sin (0 u' GO,
dM A dM A
Mh c ^(0-u'-H\ ^L = -^ s [n7 i sm(e-u'-H) t
dfj. A ^ dfj. A ' '*
which give
008,^ = 0, ^ = + 0.5072,
^ = 0, A. = + 0.0204,
di! di'
-J- = + 1.5051, -|r = + -
-^- = + 2.0978, -^- = + 0.0422.
dff =-1-1.1922, -^ = + 0.0143,
ClM n
~ = + 538.00, -^?-= -1.71.
Therefore, the equations for cos 7] &6 and &y become
cos ^ A0 = 4 1.5051 A/ 4 2.0978 A? 4 1.1922 Ajtf 4 538.00
AT? = 4 0.0086 A/ 4 0.0422 A?. 4 0.0143 A Jf 1.71 A/*
4 0.5072 A ^'4 0.0204 Ai'.
If we assign to the elements of the orbit the variations
DIFFERENTIAL FORMULA. 161
AO/= -6".64, A&' = 14".12, AI* = 8".86,
A? == -f 10", Aj/ = -|- 10", A// = -f 0".01,
we have
A/ = A/ + cos i' A&' = 19".96 ;
and the preceding equations give
cos y A0 = + 8".24, AT; = 6".96.
With the same values of AO/, A & ', &c., we have already found
cos 8 Aa = + 5".47, Ad == 9".29,
which, by means of the equations (88), writing a and d in place of
A and /9, give
cos ti A0 = -f 8".23, AI? = 6".96.
59. In special cases, in which the differences between the calcu-
lated and the observed values of two spherical co-ordinates are given,
and the corrections to be applied to the assumed elements are sought,
it may become necessary, on account of difficulties to be encountered
in the solution of the equations of condition, to introduce other ele-
ments of the orbit of the body. The relation of the elements chosen
to those commonly used will serve, without presenting any difficulty,
for the transformation of the equations into a form adapted to the
special case. Thus, in the case of the elements which determine the
form of the orbit, we may use a or log a instead of ft, and the
equation
kVl -f m
"= -IT
gives
dft = j-efo = jyefloga, (89)
a /
in which ^ is the modulus of the system of logarithms. Therefore,
the coefficient of A// is transformed into that of A log a by multiply-
ing it by f ; and if the unit of the rath decimal place of the loga-
A
rithms is taken as the unit of A loga, the coefficient must be also
multiplied by 10~ m . The homogeneity of the equation is not disturbed,
since /J. is here supposed to be expressed in seconds.
If we introduce logjp as one of the elements, from the equation
p = a cos* d -7 ~n - > and -7-= - by
d ) du -f- sin (I Q> ) di,
which determine the relations between the variations of the elements
of the orbit and those of the heliocentric longitude and latitude.
By differentiating the equations (88) w neglecting the latitude of
DIFFERENTIAL FORMULA. 165
the sun, and considering 1, /9, J, and O as variables, we derive, after
reduction,
r>
cos /? dX = _ cos (A O)dO,
R (94)
dp = -- sin sin A O
which determine the variation of the geocentric latitude and longitude
arising from an increment assigned to the longitude of the sun. It
appears, therefore, that an error in the longitude of the sun will
produce the greatest error in the computed geocentric longitude of a
heavenly body when the body is in opposition.
166 THEORETICAL ASTRONOMY.
CHAPTER III.
INVESTIGATION OP FORMULAE FOR COMPUTING THE ORBIT OF A COMET MOVING IN
A PARABOLA, AND FOR CORRECTING APPROXIMATE ELEMENTS BY THE VARIATION
OF THE GEOCENTRIC DISTANCE.
61. THE observed spherical co-ordinates of the place of a heavenly
body furnish each one equation of condition for the correction of the
elements of its orbit approximately known, and similarly for the
determination of the elements in the case of an orbit wholly unknown ;
and since there are six elements, neglecting the mass, which must
alw r ays be done in the first approximation, the perturbations not
being considered, three complete observations will furnish the six
equations necessary for finding these unknown quantities. Hence,
the data required for the determination of the orbit of a heavenly
body are three complete observations, namely, three observed longi-
tudes and the corresponding latitudes, or any other spherical co-
ordinates which completely determine three places of the body as
seen from the earth. Since these observations are given as made at
some point or at different points on the earth's surface, it becomes
necessary in the first place to apply the corrections for parallax. In
the case of a body whose orbit is wholly unknown, it i? impossible
to apply the correction for parallax directly to the place of the body ;
but an equivalent correction may be applied to the places of the
earth, according to the formula? which will be given in the next
chapter. However, in the first determination of approximate ele-
ments of the orbit of a comet, it will be sufficient to neglect entirely
the correction for parallax. The uncertainty of the observed places
of these bodies is so much greater than in the case of well-defined
objects like the planets, and the intervals between the observations
which will be generally employed in the first determination of the
orbit will be so small, that an attempt to represent the observed places
with extreme accuracy will be superfluous.
When approximate elements have been derived, we may find the
distances of the comet from the earth corresponding to the three
observed places, and hence determine the parallax in right ascension
DETERMINATION OF AN ORBIT. 167
and in declination for each observation by means of the usual formulae.
Thus, we have
?r/> cos and i, and these elements may refer to any fundamental
plane whatever. If we multiply the first of these equations by
sin (it" it'), the second by sin (u" it), and the third by
sin (u f it), and add the products, we find, after reduction,
- sin (u" u'} ~ sin (u" u) -f ^ sin (u' u) = 0,
which, by introducing the values of [rr'], [rr"], and [r'r"], becomes
[//'] x [rr"] x' + [rr'] x" = 0.
If we put
[r'r"] [rr'] , Q ,
n== D^J n FT
weget > " "_ r^
In precisely the same manner, we find
ny */' + n"f = 0, _
>' i ~"~"__n W
DETERMINATION OF AN ORBIT. 169
~ince the coefficients in these equations are independent of the posi-
tions of the co-ordinate planes, except that the origin is at the centre
of the sun, it is evident that the three equations are identical, ana
express simply the condition that the plane of the orbit passes through
the centre of the sun ; and the last two might have been derived
from the first by writing successively y and z in place of x.
Let ^, A', X" be the three observed longitudes, /9, /?', ft" the corre-
sponding latitudes, and A, J', A" the distances of the body from the
earth ; and let
A cos p = P , A' cos p = f>', A" cos p' = //',
which are called curtate distances. Then we shall have
x = p cos A R cos O , d = p' cos A' R' cos Q ',
y = p sin A R sin O , y' = p' sin X' R' sin O',
z p tan ft, z' = p' tan p,
of' = P " cos I" .fl"cosQ",
in which the latitude of the sun is neglected. The data may be so
transformed that the latitude of the sun becomes 0, as will be ex-
plained in the next chapter ; but in the computation of the orbit of
a comet, in which this preliminary reduction has not been made, it
will be unnecessary to consider this latitude which never exceeds 1",
while its introduction into the formulae would unnecessarily com-
plicate some of those which will be derived. If we substitute these
values of x, x> ', &c. in the equations (4) and (5), they become
= n (p cos A R cos O ) (p r cos / R r cos O ')
+ n" (p" cos A" R" cos O"),
= n (p sin A R sin O ) (p sin X' R' sin O') (6)
+ n" (p" sin A" R" sin "),
= np tan /9 p' tan p -f ri'p" tan ft".
These equations simply satisfy the condition that the plane of the
orbit passes through the centre of the sun, and they only become
distinct or independent of each other when n and n" are expressed
in functions of the time, so as to satisfy the conditions of undisturbed
motion in accordance with the law of gravitation. Further, they
involve five unknown quantities in the case of an orbit wholly
unknown, namely, w, n", p y p', and p" ; and if the values of n and
n" are first found, they will be sufficient to determine p, p', and p".
170 THEOEETICAL ASTRONOMY.
The determination, however, of n and n" to a sufficient degree of
accuracy, by means of the intervals of time between the obsei vations,
requires that p f should be approximately known, and hence, in
general, it will become necessary to derive first the values of n, n" y
and p r after which those of p and p tf may be found from equations
(6) by elimination. But since the number of equations will then
exceed the number of unknown quantities, we may combine them in
such a manner as will diminish, in the greatest degree possible, tho
effect of the errors of the observations. In special cases in which
the conditions of the problem are such that when the ratio of two
curtate distances is known, the distances themselves may be deter-
mined, the elimination must be so performed as to give this ratio
with the greatest accuracy practicable.
63. If, in the first and second of equations (6), we change the
direction of the axis of x from the vernal equinox to the place of the
sun at the time t' t and again in the second, from the equinox to the
second place of the body, we must diminish the longitudes in these
equations by the angle through which the axis of x has been moved,
and we shall have
= n(p cos(A 0') jRcos(O' O)) 0' cos (A' 0') #)
+ n" ( P " cos (A"_ 00 - R" cos (O" - 00),
= 7i 0> sin (A 0')+jRsm(0' O)) /sin(A' 0')
+ n" (," sin (A" _ 00 -R" Bin(0"- ')), (7)
= nO>sin(A' A) +jRsin(0 /)) R sin(O' A')
- n" (p" sin (A" A') R" sin (O " *')),
= np tan p p' tan p -f ri'p" tan p'.
If we multiply the second of these equations by tan/3', and the
fourth by sin (A' 0'), and add the products, we get
= n" P " (tan p' sin (A" O tan 0" sin (A' O 0)
n".R"sin(0" tan jfiH- n/ (tan ^ sin (A 0') tan/?sin(A' OO)
+ nJ2sin(' )tan/3 / . (8)
Let us now denote double the area of the triangle formed by the
sun and two places of the earth corresponding to R and R f fry
and we shall have
[#]= KR f sm(Q' O),
and similarly
[ JH2"] = RR" sin ( O " O ),
'] = R'R" sin(O" 00-
ORBIT OF A HEAVENLY BODY. 17l
Then, if we put
_ ,
[#']' ~ IEK"Y
we obtain
Substituting this in the equation (8), and dividing by the coefficient
of p", the result is
_ Q _
~ p n" ' tan ft" sin (A' Q ') tan ft' sin (A" Q ')
V tan ft' sin (A Q ') tan ft sin (A' Q ')
n (A' Q ') tan ft' sin (A" Q '
_!L _^M _ J?sin(0' Q)tanff' _
n" N" J tan /3" sin (A' - Q') tan /3' sin (A" 0')'
Let us also put
M' = tan ^ sin ( ; ~ Q ') ~ tan ^ sin (*' Q')
~~ tan ft" sin (A' 0') tan ft' sin (A" 0')'
j^ __ _ sin(0'-0)tanf __
tan/5" sin (A' 0') tan/3' sin (A" 0')'
and the preceding equation reduces to
-*, JTR (11)
"We may transform the values of 3/ r and M " so as to be better
adapted to logarithmic calculation with the ordinary tables. Thus,
if w' denotes the inclination to the ecliptic of a great circle passing
through the second place of the comet and the second place of the
sun, the longitude of its ascending node will be O', and we shall
have
sin (A' O') tan w' = tan ft'. (12)
Let /9 , ft Q " be the latitudes of the points of this circle corresponding
to the longitudes A and A", and we have, also,
tan ft = sin (A 0') tan w',
tan ft" = sin (A" Q') tan w'.
Substituting these values for tan/9', sin (A O 7 ) and sin (A" O ; )
in the expressions for M ' and M ", and reducing, they become
sin(/y- /?) cos p' cos ft "
sin (ft" /3 ") ' cos ft cos ft '
=^ S in( '-
172 THEORETICAL ASTRONOMY.
When the value of -77 has been found, equation (11) will give the
relation between p and p" in terms of known quantities. It is evi-
dent, however, from equations (14), that when the apparent path of
the comet is in a plane passing through the second place of the
sun, since, in this case, ft /? and ft" = /? ", we shall have M'= ^
and M fr = oo. In this case, theiefore, and also when /9 ft and
ft" ft Q " are very nearly 0, we must have recourse to some other
equation which may be derived from the equations (7), and which
does not involve this indetermination.
It will be observed, also, that if, at the time of the middle obser-
vation, the comet is in opposition or conjunction with the sun, the
values of M f and M " as given by equation (14) will be indeter-
minate in form, but that the original equations (10) will give the
values of these quantities provided that the apparent path of the
comet is not in a great circle passing through the second place of the
sun. These values are
,_ sin (A 0') _ s in(Q' Q)
" sin (A"')' " sm(A"_ 0') '
Hence it appears that whenever the apparent path of the body is
nearly in a plane passing through the place of the sun at the time of
the middle observation, the errors of observation will have great
influence in vitiating the resulting values of M f and M" and to
obviate the difficulties thus encountered, we obtain from the third of
equations (7) the following value of p" :
_ n sin (A' A)
p ==p f
^E sin (G - A') - ^ # sin (' - AO + 12" sin (" - X)
(15)
We may also eliminate p between the first and fourth of equa-
tions (7). If we multiply the first by tan/3', and the second by
cos (A' '), and add the products, we obtain
= n"p" (tan f cos (A" 0') tan /?" cos (A' '))
ri'R" tan p cos(O" 0') + w/> (tan /3' cos (A Q') tan/?cos(A' 00)
nR tan fi cos (O' O) + R' tan /3',
from which we derive
ORBIT OF A HEAVENLY BODY. 173
"_ JL tan ^ cos (A Q ') tan ff cos (A' Q')
p ~~ P n" ' tan /?" cos (A' 0') tan /5' cos (A" Q ') (16)
R" tan/?' cos(O" 00 + 7 12 tan ^ COB (0' 0) R'tan/?
tan /S" cos (A' Q ') tan $ cos (A" Q')
Let us now denote by I' the inclination to the ecliptic of a great
circle passing through the second place of the comet and that point
of the ecliptic whose longitude is 0' 90, which will therefore be
the longitude of its ascending node, and we shall have
cos (A' 0') tan I' = tan p ; (17)
and, if we designate by /9, and /? the latitudes of the points of this
circle corresponding to the longitudes ^ and X", we shall also have
tan /?, = cos (A 0')tan.T',
tan /? = cos (A" 0') tan I'.
Introducing these values into equation (16), it reduces to
_ n sin (/9, /9) cos /3" cos /?
~ p '
sin (/5" /?) cos /5 cos /?, (19)
tan r cos /5" cos/?,, / f ,. . n f . R' \
rinC9"- A( ) (^"^(e'-GO + ^^cosCG'-Q)-^).
from which it appears that this equation becomes indeterminate when
the apparent path of the body is in a plane passing through that
point of the ecliptic whose longitude is equal to the longitude of the
second place of the sun diminished by 90. In this case we may use
equation (11) provided that the path of the comet is not nearly in
the ecliptic. When the comet, at the time of the second observation,
is in quadrature with the sun, equation (19) becomes indeterminate
in form, and we must have recourse to the original equation (16),
which does not necessarily fail in this case.
When both equations (11) and (16) are simultaneously nearly in-
determinate, so as to be greatly affected by errors of observation, the
relation between p and p" must be determined by means of equation
(15), which fails only when the motion of the comet in longitude is
very small. It will rarely happen that all three equations, (14),
(15), and (16), are inapplicable, and when such a case does occur it
will indicate that the data are not sufficient for the determination of
the elements of the orbit. In general, equation (16) or (19) is to be
used when the motion of the comet in latitude is considerable, and
equation (15) when the motion in longitude is greater than in latitude.
174 THEORETICAL ASTRONOMY.
64. The formulae already derived are sufficient to determine the
relation between //' and p when the values of n and n" are known,
and it remains, therefore, to derive the expressions for these quan-
tities.
If we put
t(f -f) = r",
k (" - O = r, (20)
and express the values of x, y, z, x lf y y", z" in terms of re', y f , z r by
expansion into series, we have
f _dx^ T" J_ dV r^_ __ 1_ dV r" 8
x ~ , and the sign to be adopted must be determined
from the physical conditions of the problem.
If we suppose the chords g and x to be proportional to the linear
velocities of the earth and comet at the middle observation, we have,
the eccentricity of the earth's orbit being neglected,
= S r V
which shows that H is greater than g } and that d is positive, so long
as r f is less than 2. The comets are rarely visible at a distance from
the earth which much exceeds the distance of the earth from the sun,
and a comet whose radius-vector is 2 must be nearly in opposition in
order to satisfy this condition of visibility. Hence cases will rarely
occur in which d can be negative, and for those which do occur it
will generally be easy to determine which sign is to be used. How-
ever, if d is very small, it may be impossible to decide which of the
two solutions is correct without comparing the resulting elements
with other and more distant observations.
192 THEOKETICAL ASTKONOMY.
70. When the values of r and r" have been finally determined, as
just explained, the exact value of d may be computed, and then we
have
_ d + g cos / jy i \ \ sec 20'
(I (t*" -f- ) ,) = - - - -j =,
cos I (t* w) K rr
we get
0' = 22' 47".4, w = 115 40' 6' ; .3, log q = 9.887378.
Hence we have
* = 01 -f = 60 23' 17". S,
204 THEORETICAL ASTRONOMY.
and
i = u > = 27 12' 6".l, v" = ti" = 37 38' 43 '.1.
Then we obtain
log m = 9.9601277 | log 5 = 0.129061,
and, corresponding to the values of v and v", Table VI. gives
log M = 1.267163, log M" = 1.424152.
Therefore, for the time of perihelion passage, we have
and
T=t = t 13.74364,
m
T=t' =f 19.72836.
m
The first value gives T= 1863 Dec. 27.56473, and the second gives
T= Dec. 27.56463. The agreement between these results is the final
proof of the calculation of the elements from the adopted value of
M= f -.
P
If we find T by means of Table VIII., we have
log N = 0.021616, log N" = 0.018210,
and the equation
2 2
T = t 3 Nr% sin v = if' 3^ N"i"* sin v",
in which log ^ = 1.5883273, gives for T the values Dec. 27.56473
and Dec. 27.56469.
Collecting together the several results obtained, we have the fol-
lowing elements:
T 1863 Dec. 27.56471 Washington mean time.
= 6023'17".8 _._ . ,
ic and Mean
log q = 9.887378.
Motion Direct.
73. The elements thus derived will, in all cases, exactly represent
the extreme places of the comet, since these only have been used in
finding the elements after p and p" have been found. If, by means
NUMERICAL EXAMPLES. 205
of these elements, we compute n and n n ', and correct the \alue of M,
the elements which will then be obtained will approximate nearer
the true values; and each successive correction will furnish more
accurate results. When the adopted value of M is exact, the result-
ing elements must by calculation reproduce this value, and since the
computed values of A, A", /9, and ft" will be the same as the observed
values, the computed values of X and ft' must be such that when
substituted in the equation for M, the same result will be obtained
as when the observed values of X' and ft' are used. But, according
to the equations (13) and (14), the value of M depends only on the
inclination to the ecliptic of a great circle passing through the places
of the sun and comet for the time ', and is independent of the angle
at the earth between the sun and comet. Hence, the spherical co-
ordinates of any point of the great circle joining these places of the
sun and comet wilj, in connection with those of the extreme places,
give the same value of M, and when the exact value of M has been
used in deriving the elements, the computed values of A' and ft' must
give the same value for w r as that which is obtained from observa-
tion. But if we represent by if/ the angle at the earth between the
sun and comet at the time ', the values of if/ derived by observation
and by computation from the elements will differ, unless the middle
place is exactly represented. In general, this difference will be small,
and since w' is constant, the equations
cos 4/ = cos p cos (A' O')>
sin 4-' cos w' = cos ft' sin (A' '), (93)
sin 4/ sin w' = sin ft,
give, by differentiation,
cos p dl' = cos w' sec p d$ r ,
dp = smufcoa(X O'HV-
From these we get
cosft'dA' tan(*' Q')
dp sin p
which expresses the ratio of the residual errors in longitude and
latitude, for the middle place, when the correct value of M has been
used.
Whenever these conditions are satisfied, the elements will be
correct on the hypothesis of parabolic motion, and the magnitude
of the final residuals in the middle place will depend on the deviation
of the actual orbit of the comet from the parabolic form. Further,
206 THEORETICAL ASTRONOMY.
when elements have been derived from a value of M which has not
been finally corrected, if we compute X r and $' by means of these
elements, and then
the comparison of this value of tan w f with that given by observa-
tion will show whether any further correction of M is necessary, and
if the difference is riot greater than what may be due to unavoidable
errors of calculation, we may regard M as exact.
To compare the elements obtained in the case of the example
given with the middle place, we find
v f = 32 31' 13".5, u' = 148 IV 19".8, log / == 9.922836.
Then from the equations
tan (f Q> ) = cos i tan u',
tan b f = tan i sin (l f & ),
we derive
I' = 109 46' 48".3, V = 28 24' 56".0.
By means of these and the values of O' and R f , we obtain
A' = 302 57' 41".l, p = 57 39' 37".0 ;
and, comparing these results with the observed values of X' and /?',
the residuals for the middle place are found to be
Comp. Obs.
cos p AA' = -f 3".6, A/9 = + I'M.
The ratio of these remaining errors, after making due allowance for
unavoidable errors of calculation, shows that the adopted value of
M. is not exact, since the error of the longitude should be less than
that of the latitude.
The value of w 1 given by observation is
log tan w' = 0.966314,
and that given by the computed values of A' and ft' is
log tan w' = 0.966247.
The difference being greater than what can be attributed to errors of
calculation, it appears that the value of M requires further cor-
NUMERICAL EXAMPLES. 207
rection. Since the difference is small, we may derive the correct
value of M by using the same assumed value of ,, and, instead of
Yb
the value of tan&/ derived from observation, a value differing as
much from this in a contrary direction as the computed value differs.
Thus, in the present example, the computed value of log tan w f is
0.000067 less than the observed value, and, in finding the new value
of Mj we must use
log tan w' = 0.966381
in computing /? and /?/' involved in the first of equations (14). If
the first of equations (10) is employed, we must use, instead of tan/9'
as derived from observation,
tan p = tan vf sin (A' 0'),
or
log tan f = 0.966381 + log sin (A' 0') = 0.198559,
the observed value of X' being retained. Thus we derive
log M= 9.829586,
and if the elements of the orbit are computed by means of this
value, they will represent the middle place in accordance with the
condition that the difference between the computed and the observed
value of tan w r shall be zero.
A system of elements computed with the same data from
log M= 9.822906 gives for the error of the middle place,
C.-O.
cos f A;/ = 1' 26".2, &p = 40".l.
If we interpolate by means of the residuals thus found for two values
of M, it appears that a system of elements computed from
log M= 9.829586
will almost exactly represent the middle place, so that the data are
completely satisfied by the hypothesis of parabolic motion.
The equations (34) and (32) give
log -^- = 0.006955, log -^ = 0.006831,
/ f\/ JLV
and from (10) we get
log M ' = 9.822906, log M" = 9.663729...
208 THEORETICAL ASTRONOMY.
Then by means of the equation (33) we derive, for the corrected
value of My
log M= 9.829582,
which differs only in the sixth decimal place from the result obtained
by varying tanw' and retaining the approximate values t = -r f = j'
74. When the approximate elements of the orbit of a comet are
known, they may be corrected by using observations which include
a longer interval of time. The most convenient method of effecting
this correction is by the variation of the geocentric distance for the
time of one of the extreme observations, and the formulae which
may be derived for this purpose are applicable, without modification,
to any case in which it is possible to determine the elements of the
orbit of a comet on the supposition of motion in a parabola. Since
there are only five elements to be determined in the case of parabolic
motion, if the distance of the comet from the earth corresponding to
the time of one complete observation is known, one additional com-
plete observation will enable us to find the elements of the orbit.
Therefore, if the elements are computed which result from two or
more assumed values of A differing but little from the correct value,
by comparison of intermediate observations with these different sys-
tems of elements, we may derive that value of the geocentric distance
of the comet for which the resulting elements will best represent the
observations.
In order that the formulae may be applicable to the case of any
fundamental plane, let us consider the equator as this plane, and,
supposing the data to be three complete observations, let A, A f , A"
be the right ascensions, and D, D', D" the declinations of the sun
for the times t, t', t". The co-ordinates of the first place of the earth
referred to the third are
x R" cos D" cos A" R cos D cos A,
y = R" cos D" sin A" EcosD sin A,
z=R"smD" EsmD.
If we represent by g the chord of the earth's orbit between the places
for the first and third observations, and by G and K, respectively,
the right ascension and declination of the first place of the earth as
seen from the third, we shall have
x = g cos K cos G,
y = g cos K sin G,
z = g sin K.
VAKIATION OF THE GEOCENTKIC DISTANCE. 209
and, consequently,
g cos K cos ( G A) = R" cos D" cos (A" A) R cos D,
g cos K sin ( G A) = jR" cos D" sin (A" A), (96)
g sin iT = .#" sin Z>" jR sin D,
from which #, JT, and G may be found.
If we designate by x n y n z, the co-ordinates of the first place of
the comet referred to the third place of the earth, we shall havf
x, = A cos 5 cos a -f g cos K cos G,
y, A cos d sin a -f- g cos K sin (r,
2, = A sin 5 -f- # sin K.
Let us now put
x, = h' cos C' cos IT,
^ == A' cos :' sin #',
z, = h' sin C',
and we get
h f cos % cos (.ff ' G) = A cos fl cos (a G) + cos JT,
A' cos C' sin GET G) = J cos sin (a (f). C97)
A.' sin C' = J sin = 116 56' 33".7, o" = u" a> = 88 47' 55".l,
from which we get
T= 1858 Sept. 29.4274.
From these elements we find
log r' == 0.212844, v f = 107 7' 34".0, u' = 21 59' 12".3,
and from
tan (f &) cos i tan u',
tan b f = tan i sin (J! & ),
we get
I' = 154 56' 33".4, 6' = + 19 30' 22".l.
NUMERICAL EXAMPLE. 217
By means of these and the values of O' and R f , we obtain
A' = 137 39' 13".3, p = + 12 54' 45".3,
and comparing these results with observation, we have, for the error
of the middle place,
C. O.
/ = 27".2, A/5' = 23".7.
From the relative positions of the sun, earth, and cctnefc at the
time t" it is easily seen that, in order to diminish these residuals, the
geocentric distance must be increased, and therefore we assume, for
a second value of J,
log J = 0.398500,
from which we derive
H' = 153 44' 57".6, C' = -f 7 24' 26".l, log h' = 0.488026,
log C= 9.912587, logc = 0.472115, logr = 0.324207,
log A" = 0.311054, log r" = 0.054824, log x = 0.089922,
Then we find the heliocentric places
I = 159 40' 33".8, b = + 10 50' 8".6, logr = 0.324207,
I" = 144 17 12 .1, b" = + 35 8 37 .8, log r" = 0.054825,
and from these,
= 165 15' 41".l, i = 63 2' 49".2,
u = 12 10 30.8, w"=:40 13 26.0,
may be found.
Since
o , r r" - o f
cos 2/ = ; n, sm 2y = ; j
r + r ' r + r"
we have
g , 1> o^ _ 2 ^ o
r T F V / r/ / sin2/
p p _ 2p cot 2/
" '
and from equations (70),
. , sin i (it" w) tan 6r f -of cos /
cot 2/ == - - / - , sm 2/ = - T7-77 - r.
cos / cos ^ (u u)
Therefore the formulae (87) reduce to
* sin , - " t* = - == tan G'
(88)
cos (w (?/' + w)) = - x _ sec J (w v w;,
cos/F rr"
from which also 6 and co may be derived. Then
sin ) 8
~~
COS 2 , V f = U to, v" = u" ct> '
and if we compute r, r' 9 r" from these by means of the polar equa-
tion of the conic section, the results should agree with the values of
the same quantities previously obtained. According to the equation
(45)j, we have
tan \E = tan (45 ?) tan Jv,
tan \E' = tan (45 J?) tan jv', (90)
tan E" == tan (45 j?) tan W',
from which to find E, E' y and E" . The difference E" E should
agree with that derived from equation (84) within the limits of
accuracy afforded by the logarithmic tables. Then, to find the mean
anomalies, we have
M =E eswE,
M'=E'esmE', (91;
M" = E"esmE";
and, if M Q denotes the mean anomaly corresponding to any epoch T t
we have, also,
? n the application of which the values of t, t> ', and t" must be those
which have been corrected for the time of aberration. The agree-
262 THEORETICAL ASTRONOMY.
ment of the three values of M will be a final test of the accuracy of
the entire calculation. If the final values of P and Q are exact,
this proof will be complete within the limits of accuracy admitted
by the logarithmic tables.
When the eccentricity is such that the equations (91) cannot be
solved with the requisite degree of accuracy, we must proceed accord-
ing to the methods already given for finding the time from the peri-
helion in the case of orbits differing but little from the parabola.
For this purpose, Tables IX. and X. will be employed. As soon as
v, v f , and v" have been determined, we may find the auxiliary angle
V for each observation by means of Table IX. ; and, with V as the
argument, the quantities M 9 M f , M" (which are not the mean anoma-
lies) must be obtained from Table VI. Then, the perihelion distance
having been computed from
P
we shall have
in which log (7 = 9.96012771, for the determination of the time of
perihelion passage. The times t, t', t u must be those which have
been corrected for the time of aberration, and the agreement of the
three values of T is a final proof of the numerical calculation.
If Table X. is used, as soon as the true anomalies have been found,
the corresponding values of log B and log C must be derived from
the table. Then w is computed from
and similarly for w f and w" ; and, with these as arguments, we derive
M, M r , M" from Table VI. Finally, we have
MBq* _
(93)
for the time of perihelion passage, the value of O being the same as
in (92).
When the orbit is a parabola, e = I and p 2g, and the elements
Q and = 190 15' 39".57, log e = log sin ? = 9.2751434,
?> == 10 51 39 .62, TT = , -f- = 37 15' 40".29.
This value of ^ gives log cos , v" = u" cw,
according to which we have
v = 327 56' 39".97, v' = 330 27' 6".25, i" = 333 3' 27".57.
If we compute r, r f , and r" from these values by means of the polar
equation of the ellipse, we get
log r = 0.3048367, log r' = 0.3032586, log r" = 0.3017481,
and the agreement of these results with those derived directly from
p, p f , and p rr is a further proof of the calculation.
The equations
tan { E = tan (45 ?) tan v,
tan ^E' = tan (45 1$0 tan ^,
tan JJE" = tan (45 ?) tan W
give
E = 333 17' 28' .18, E' = 335 24' 38".00, E" = 337 36' 19".78.
NUMERICAL EXAMPLE. 275
The \alue of J (E" E) thus obtained differs only 0".003 from that
computed directly from x' '.
Finally, for the mean anomalies we have
M = E e sin E, M' = E' e sin E' y M" = E" e sin E'\
from which we get
M = 338 8' 36".71, M' = 339 54' 10".61, M" = 341 43' 6".97 ;
and if Jf denotes the mean anomaly for the date T=1863 Sept. 21.5
Washington mean time, from the formulae
M=M n(t r>
we obtain the three values 339 55' 25".97, 339 55' 25".96, and
339 55' 25".96, the mean of which gives
M 9 = 339 55' 25".96.
The agreement of the three results for Jf is a final proof of the
accuracy of the entire calculation of the elements.
Collecting together the separate results obtained, we have the fol-
lowing elements :
Epoch = 1863 Sept. 21.5 Washington mean time.
M = 339 55' 25".96
*= 37 15 40 .29)
ft - 207 72 V Ech P tlc and Mean
06 * J " ' v v I *i / . . - (99)
in which log JA = 8.8596330. We have, also, to the same degree of
approximation,
!< '=71' to^ = -5> (100)
For the values
log r = 9.0782331, log r f = 9.3724848, log r" ** 9.0645692,
log/ = 0.3032587,
these formulae give
log s = 0.0001277, log s' = 0.0004953, log " = 0.0001199,
which differ but little from the correct values 0.0001284, 0.0004954,
and 0.0001193 previously obtained.
Since
sec 8 / = 1 + 6 sin 2 tf + &c.,
the second of equations (65) gives
sin ' &c -
Substituting this value in the first of equations (66), we get
If we neglect terms of the fourth order with respect to the time, it
will be sufficient in this equation to put y r = f, according to (71), and
hence we have
4
3
and, since s f 1 is of the second order with respect to r', we have,
to terms of the fourth order,
280 THEORETICAL ASTRONOMY.
Therefore,
which, when the intervals are small, may be used to find s f from r
and r". In the same manner, we obtain
l^r,- (102)
For logarithmic calculation, when addition and subtraction loga-
rithms are not used, it is more convenient to introduce the auxiliary
angles , r , and ", by means of which these formulae become
K tW (103)
in which log |^ = 9.7627230. For the first approximation these
equations will be sufficient, even when the intervals are considerable,
to determine the values of s and s" required in correcting P and Q.
The values of r, r', r", and r" above given, in connection with
log r = 0.3048368, log r" = 0.3017481,
give
log 5 = 0.0001284, log s f = 0.0004951, log s" = 0.0001193.
These results for logs and logs" are correct, and that for logs' differs
only 3 in the seventh decimal place from the correct value.
ORBIT FEOM FOUK OBSERVATIONS.
CHAPTER V.
DETERMINATION OF THE ORBIT OF A HEAVENLY BODY FROM FOUR OBSERVATIONS,
OF WHICH THE SECOND AND THIRD MUST BE COMPLETE.
95. THE formulae given in the preceding chapter are not sufficient
to determine the elements of the orbit of a heavenly body when its
apparent path is in the plane of the ecliptic. In this case, however,
the position of the plane of the orbit being known, only four ele-
ments remain to be determined, and four observed longitudes will
furnish the necessary equations. There is no instance of an orbit
whose inclination is zero ; but, although no such case may occur, it may
happen that the inclination is very small, and that the elements
derived from three observations will on this account be uncertain,
and especially so, if the observations are not very exact. The diffi-
culty thus encountered may be remedied by using for the data in the
determination of the elements one or more additional observations,
and neglecting those latitudes which are regarded as most uncertain.
The formulae, however, are most convenient, and lead most expe-
ditiously to a knowledge of the elements of an orbit wholly unknown,
when they are made to depend on four observations, the second and
third of which must be complete ; but of the extreme observations
only the longitudes are absolutely required.
The preliminary reductions to be applied to the data are derived
precisely as explained in the preceding chapter, preparatory to a de-
termination of the elements of the orbit from three observations.
Let t, t', t", t'" be the times of observation, r, r', r", r'" the radii-
vectores of the body, u, u f , u n ', u'" the corresponding arguments of
the latitude, R, R', R", R" r the distances of the earth from the sun,
and O, O', O", O'" the longitudes of the sun corresponding to
these times. Let us also put
= r'r'"sm(u'" u'\
[rV"] = rV" sin (u" r u"\
and
282 THEORETICAL ASTRONOMY.
Then, according to the equations (5) 3 , we shall have
nx x' -f- ri'x" = 0,
ny _ + n y =0 ,
.Let ;, A', ;/', A'" be the observed longitudes, /9, /9', 0", /?'" the ob-
served latitudes corresponding to the times , t f , t' r , t" r , respectively,
and J, A', A", A'" the distances of the body from the earth. Further,
let
A'" cos f" = p'" 9
and for the last place we have
x f " = p'" cos A'" K" cos O'",
y m = p'" sin *'" "' sin O'".
Introducing these values of x'" and y f ", and the corresponding values
of x, x f , x", y, y f , y" into the equations (2), they become
= n (p cos I R cos Q ) (>' cos A' R' cos O')
+ n"G
= n sin A J2 sin O ) (f? sin A r R sin O'
+ w"
= n ' G/ cos / - R cos O') (p" cos A" #' cos O") (3)
-f n'" (p'" cos X'" R" cos 0'"),
= w' (p r sin A' ^ sin O') (p" sin X" R" sin 0")
4- n'" (/>"' sin A'" R" sin 0'").
If we multiply the first of these equations by sin ^, and the second
by cos ^, and add the products, we get
= nR sin (A Q) (p* sin (X X) + R f sin (A 0'))
+ n" (P" sin (A" A) -f- #' sin (A ")) ; (4)
and in a similar manner, from the third and fourth equations, we
find
= n' (p f sin (X" X ) R' sin (X" ')) (5j)
- (p" sin (/" A") K' sin (/" ")) - n'"R" sin (A'" 0'").
Whenever the values of n, n' 9 n", and n //r are known, or may be
determined in functions of the time so as to satisfy the conditions of
motion in a conic section, these equations become distinct or inde-
Dendent of each other ; and, since only two unknown quantities p 1
OEBIT FROM FOUR OBSERVATIONS. 283
and p n are involved in them, they will enable us to determine these
curtate distances.
Let us now put
cos p sin (/ /I) = A, cos ft" sin (A" X)=B,
cos ft" sin (/'" A") = C, cos p sin (*'" X ) = D,
and the preceding equations give
Ap' sec ^ Bn" P " sec /5" = wJ5 sin (A ) R r sin (A Q')
+ n"J2"Biny 0"),
!>>' sec /?' Q/' sec /9"= n'tf sin (/" ') R" sin (A'" Q") (7)
+ n'"R'" sin (/" Q'").
Tf we assume for n and n" their values in the case of the orbit of
the earth, which is equivalent to neglecting terms of the second order
in the equations (26) 3 , the second member of the first of these equa-
tions reduces rigorously to zero ; and in the same manner it can be
shown that when similar terms of the second order in the corre-
sponding expressions for n f and n" are neglected, the second member
of the last equation reduces to zero. Hence the second member of
each of these equations will generally differ from zero by a quantity
which is of at least the second order with respect to the intervals of
time between the observations. The coefficients of p' and p" are of
the first order, and it is easily seen that if we eliminate p" from
these equations, the resulting equation for // is such that an error of
the second order in the values of n and n" may produce an error of
the order zero in the result for p f , so that it will not be even an
approximation to the correct value ; and the same is true in the case
of p rr . It is necessary, therefore, to retain terms of the second order
in the first assumed values for n, n f , n", and ii' 1 ' \ and, since the
terms of the second order involve r f and r f/ , we thus introduce two
additional unknown quantities. Hence two additional equations in-
volving r f , r", p'j p" and quantities derived from observation, must
be obtained, so that by elimination the values of the quantities sought
may be found.
From equation (34) 4 we have
P' sec p = R' cos 4/ V r' 2 R r * sin 2 4/, (8)
which is one of the equations required; and similarly we find, for
the other eauation,
P " sec ft" = R" cos V ' db y r'" 2 R" 2 sin 2 4". (9)
284 THEORETICAL ASTRONOMY.
Introducing these values into the equations (7), and putting
x' = l/V 2 -.fl' 2 sin 2 47,
x''=V / r"* # 2 sin 2 4/',
we get
Ax' Bri'x" = nR sin (A O) R' sin (A 0')
-f n"R" sin (A O") AR' cos V + n"BR" cos *",
ZtoV Cx" = n'R' sin (X" 0') jR" sin (A"' 0")
-j- n'"R"' sin (A'" 0'") n'DR' cos 4' + CR" cos 4,".
Let us now put
B -h' D V>
A ht C =h >
or
, _ cos /3" sin (A" A) , _ cos ff sin (r r AQ
" cos f sin (A' A) ' ~ cos jt' sin (A'" A")'
and we have
x' = Kn"x" -f ndT of + n'V,
a;" = A" n y + n '"d" a" -f nV.
These equations will serve to determine x r and x", and hence r f and
r r/ , as soon as the values of n f n f , n ff , and n" 1 are known.
96. In order to include terms of the second order in the values of
n and n", we have, from the equations (26) 3 ,
and, putting
these give
ORBIT FROM FOUR OBSERVATIONS. 285
Lot us now put
T '" = k cr o, v = (*"' o
and, making the necessary changes in the notation in equations (26) s ,
we obtain
,r-"(r ' + r) , T"' (T"*+ T"'T - 1) dr" \
~* fn. -- - i^j --- ar*"f
r ** at i ,.
I . r ( r ' + ^'") , .r(r' + rr'"-r'"0 dr"
fs ^^ ~& 7r ~ "df
From these we get, including terms of the second order,
and hence, if we put
P" = ~, " = (' + '"-!)/', (17)
we shall have, since r ' = r + r' /r ,
= n.
When the intervals are equal, we have
P' P"
= 777 -pm
and these expressions may be used, in the case of an unknown orbit,
for the first approximation to the values of these quantities.
The equations (13) and (17) give
and, introducing these values, the equations (12) become
286 THEORETICAL ASTEOKOMY.
* = W-pr ( 1 + % } (*V + P'd' + c') -
1 ff'\
x " = TTpr ( 1 + ;* ) (A ' V + p " d " + c " }
Let us now put
P'cf+c'
An
c" -f
C 1 I TV/ J
(21;
1 + P" 1 +
and we shall have
(22)
+e ")-a".
We have, further, from equations (10),
r''=(*/'+P' 2 sinV)',
r f; =(a/'+JB"sinV / ) f ,
If we substitute these values of r' 3 and r" z in equations (22), the two
resulting equations will contain only two unknown quantities x' and
x". when P', P", ', and Q" are known, and hence they will be
sufficient to solve the problem. But if we effect the elimination of
either of the unknown quantities directly, the resulting equation
becomes of a high order. It is necessary, therefore, in the numerical
application, to solve the equations (22) by successive trials, which
may be readily effected.
If z' represents the angle at the planet between the sun and the
earth at the time of the second observation, and z" the same angle at
the time of the third observation, we shall have
, _ E' sin 4/
~'
(24)
,,_'
sin z"
Substituting these values of r f and r" in equations (10), we get
*'=r>z> (25)
x" = r" cos z ",
and hence
ORBIT FKOM FOUR OBSERVATIONS. 287
, R' sin 4,'
tan z== -- -, - ,
(26)
by means of which we may find z' and z" as ,soon as x f and x" shall
have been determined ; and then r' and r" are obtained from (24) or
(25). The last equations show that when x r is negative, z r must be
greater than 90, and hence that in this case r' is less than R f .
In the numerical application of equations (22), for a first approxi-
mation to the values of x r and x", since Q' and Q" are quantities of
the second order with respect to r or r r// , we may generally put
and we have
*'
*"
or, by elimination,
,
1 -/'/"
1 /'/"
With the approximate values of x' and x" derived from these equa-
tions, we compute first r r and r" from the equations (26) and (24),
and then new values of x' and x" from (22), the operation being
repeated until the true values are obtained. To facilitate these ap-
proximations, the equations (22) give
~/'(;+f ). >;
*/ = . " "^
Let an approximate value of x f be designated by # ', and let the
value of x" derived from this by means of the first of equations (27)
be designated by x ". With the value of x " for x" we derive a
new value of x f from the second of these equations, which we denote
by a?/. Then, recomputing x" and x', we obtain a third approximate
value of the latter quantity, which may be designated by a?,'; and,
if we put
a?/ x ' = a Q , #,' xj ',
288 THEORETICAL ASTRONOMY.
we shall have, according to the equation (67) 3 , the necessary changes
being made in the notation,
z' = < ?^- = x; PL-. (28)
< o o o
The value of x' thus obtained will give, by means of the first of
equations (27), a new value of x", and the substitution of this in the
last of these equations will show whether the correct result has been
found. If a repetition of the calculation be found necessary, the
three values of x' which approximate nearest to the true value will,
by means of (28), give the correct result. In the same manner, if
we assume for x" the value derived by putting Q f = and Q" = 0,
and compute x f , three successive approximate results for x" will
enable us to interpolate the correct value.
When the elements of the orbit are already approximately known,
the first assumed value of x' should be derived from
instead of by putting Q f and Q" equal to zero.
97. It should be observed that when A' = A or A'" = A", the equa-
tions (22) are inapplicable, but that the original equations (7) give,
in this case, either p" or p' directly in terms of n and n ff or of n'
and n'" and the data furnished by observation. If we divide the
first of equations (22) by h f , we have
The equations (21) give
f 1 /' P 17 ~r~ 17
h' ~~ 1 -I- P" A' ~ 1 + P'
and from (11) we get
of _ R' cos V B sin (I Q')
A' ~ A' B
1 = tf'cos*" + ^" sin ^- Q "), (29;
d' J?sin(A O)
A'~
Then, if we put
ORBIT FROM FOUR OBSERVATIONS.
its value may be found from the results for and ^ derived by
means of these equations, and we shall have
When A' = I, we have h' == oo, and this formula becomes
=( 1 + -J ) (*" + C ') - d + P'),
the value of j-, being given by the first of equations (29) This
equation and the second of equations (22) are sufficient to determine
x f and x" in the special case under consideration.
The second of equations (22) may be treated in precisely the same
manner, so that when X ffr = X", it becomes
o
=(
and this must be solved in connection with the first of these equations
in order to find x f and x ff .
98. As soon as the numerical values of x r and x n have been
derived, those of r f and r" may be found by means of the equations
(26) and (24). Then, according to (41) 4 , we have
X*& MO
The heliocentric places are then found from p f and p" by means of
the equations (71) 3 , and the values of r' and r" thus obtained should
agree with those already derived. From these places we compute
the position of the plane of the orbit, and thence the arguments of
the latitude for the times t f and t".
The values of r', r", u', u", n, n", n f , and n"' enable us to deter-
mine r, r" r , u t and u" f . Thus, we have
and, from the equations (1) and (3) 3 ,
19
MJK) THEORETICAL ASTRONOMY.
[rr"] =
Therefore,
n
r sin (u f u) = r" sin (u" w')
n
r sin (u" u) = - r' sin (it" it'),
(32)
i>" sin ("' ") = 4? r' sin (" '),
71
r'" sin (t*'" *') = ^7 r" sin (u" ')
From the first and second of these equations, by addition and sub-
traction, we get
r sin ((u f u) -J- 3 (w" w')) = sm 2 V w ')
/- re 'V" (33)
r cos ((i/ u) -f- 2 (w" i*')) = cos J (w" it'),
71
from which we may find r, u r it, and u = u r (u f it).
In a similar manner, from the third and fourth of equations (32),
we obtain
r m sin ((u m t*") + J (it" i/)) = V sin J (i*" - w'),
,/ w (34)
r'" cos ((i* w M ") + A (w" u')) = - m^- cos J (w" w'),
71
from which to find r'" and it'".
When the approximate values of r, r f , r", r f ", and w, it', it", u fn
have been found, by means of the preceding equations, from the
assumed values of P', P", Q f , and Q", the second approximation to
the elements may be commenced. But, in the case of an unknown
orbit, it will be expedient to derive, first, approximate values of
and r", using
p' _!_, p" = _L,
and then recompute P f and P" by means of the equations (14) and
ORBIT FROM FOUR OBSERVATIONS. 291
(18), before finding u f and u". The terms of the second order will
thus be completely taken into account in the first approximation.
99. If the times of observation have not been corrected for the
time of aberration, as in the case of an orbit wholly unknown, this
correction may be applied before the second approximation to the
elements is effected, or at least before the final approximation is com-
menced. For this purpose, the distances of the body from the earth
for the four observations must be determined ; and, since the curtate
distances p f and p ff are already given, there remain only p and p' tf to
be found. If we eliminate p f from the first two of equations (3), the
result is
nR sin (A' Q) R sin (A' Q') -f n" R" sin (X 0").
and, by eliminating p" from the last two of these equations, we also
obtain
, , n' sin (A" -A')
' : = / V'sm(A'"-/')
n 1 R' sin (A" - Q') R" sin (A" Q ") -f ri" R" f sin (A" "')
n'" sin (A'" /I")
by means of which p and p" f may be found. The combination of
the first and second of equations (3) gives
COS (X
nR cos a Q) R' cos (A 0') + n" R" cos (A 0")
/o = ^ cos (A' A) 2_- cos (A" A) (37)
n n
+
and from the third and fourth we get
/>'' cos( r'_A'o-$'
n' R 1 cos (A'" 0') R' cos (A'" 0") -f n'" R'" cos (A'" Q '")
.'" = ^7 cos (A'" - A") ^ cos (A'" A') (38)
Further, instead of these, any of the various formulae which have
been given for finding the ratio of two curtate distances, may be
employed; but, if the latitudes /?, /9', &c. are very small, the values
of p and p"' which depend on the differences of the observed longi-
tudes of the body must be preferred.
292 THEORETICAL ASTRONOMY.
T1ie values of p f and p rrf may also be derived by computing the
heliocentric places of the body for the times t and t' n by means of
the equations (82) 1? and then finding the geocentric places, or those
which belong to the points to which the observations have been
reduced, by means of (90) b writing p in place of Jcos/9. This
process affords a verification of the numerical calculation, namely,
the values of X and X' fr thus found should agree with those furnished
by observation, and the agreement of the computed latitudes ft and
ft'" with those observed, in case the latter are given, will show how
nearly the position of the plane of the orbit as derived from the
second and third observations represents the extreme latitudes. If
it were not desirable to compute X and X" in order to check the
calculation, even when ft and ft flf are given by observation, we mighi
derive p and p" r from the equations
p = r sin u sin i cot y5,
p'" = r'" sin u" f sin i cot p",
when the latitudes are not very small.
In the final approximation to the elements, and especially when
the position of the plane of the orbit cannot be obtained with the
required precision from the second and third observations, it will be
advantageous, provided that the data furnish the extreme latitudes
ft and ft" f , to compute p and p f " as soon as p f and p ff have been
found, and then find I, l f ", 6, and b'" directly from these by means
of the formulae (71) 3 . The values of & and i may thus be obtained
from the extreme places, or, the heliocentric places for the times t r
and t'" being also computed directly from p' and p", from those
which are best suited to this purpose. But, since the data will be
more than sufficient for the solution of the problem, when the extreme
latitudes are used, if we compute the heliocentric latitudes b f and b'"
from the equations
tan b' = tan i sin (I' & ),
tan b" = tan i sin (I" &),
they will not agree exactly with the results obtained directly from p 1
and /?", unless the four observations are completely satisfied by the
elements obtained. The values of r' and r n ', however, computed
directly from p r and p" by means of (71) 3 , must agree with those
derived from x r and x n .
The corrections to be applied to the times of observation on account
ORBIT FROM FOUR OBSERVATIONS. 293
of aberration may now be found. Thus, if t w t f , t Q ", and t Q '" are
the uncorrected times of observation, the corrected values will be
t = t Cp sec/5,
wherein log C== 7.760523, and from these we derive the corrected
values of r, r', r", r'", and r '.
100. To find the values of P', P", Q', and Q", which will be
exact when r, r', r", r f// , and w, u' 9 u", u'" are accurately known, we
have, according to the equations (47) 4 and (51) 4 , since Q f = |,
_ _ _
~" * ss" ' rr" cos J (t*" w') cos i (t*" M) cos j (' )'
In a similar manner, if we designate by s fff the ratio of the sector
formed by the radii-vectores r rr and r'" to the triangle formed by
the same radii-vectores and the chord joining their extremities, we
find
_ i rr'" r^
V : ~ 2 ^777 ' r y/, C o S l (>'" M ") COS i (U'" t*') COS J (u" 1*')'
The formulae for finding the value of s' rr are obtained from those for
s by writing / /r/ , f ff , G" f , &c. in place of /, /-, (r, &c., and using
r", r'", w //; u" instead of r', r", and it" M', respectively.
By means of the results obtained from the first approximation to
the values of P', P", Q', and ", we may, from equations (41) and
(42), derive new and more nearly accurate values of these quantities,
and, by repeating the calculation, the approximations to the exact
values may be carried to any extent which may be desirable. When
three approximate values of P' and ', and of P" and ", have
been derived, the next approximation will be facilitated by the use
of the formulae (82) 4 , as already explained.
When the values of P', P", Q f , and Q" have been derived with
sufficient accuracy, we proceed from these to find the elements of the
orbit. After &, i, r, r f , r" , r'", u, u', u" ', and u f/f have been found,
the remaining elements may be derived from any two radii-vectores
294 THEORETICAL ASTRONOMY.
and the corresponding arguments of the latitude. It will be most
accurate, however, to derive the elements from r, r /r/ , u, and u" f .
If the values of P', P", Q', and Q" have been obtained with great
accuracy, the results derived from any two places will agree with
those obtained from the extreme places.
In the first place, from
cos G = sin (u" f u), (43)
sin YQ sin O = cos (V" w) cos 2/ ,
cos r = cos | (u'" u) sin 2/ ,
we find p and 6r . Then we have
r '")S
"to (44 )
from which, by means of Tables XIII. and XIV., to find s and a? .
We have, further,
and the agreement of the value of p thus found with the separate
results for the same quantity obtained from the combination of any
two of the four places, will show the extent to which the approxima-
tion to P ; , P") Q f , and Q" has been carried. The elements are now
to be computed from the extreme places precisely as explained in the
preceding chapter, using r'" in the place of r" in the formula there
given and introducing the necessary modifications in the notation,
which have been already suggested and which will be indicated at
once.
101. EXAMPLE. For the purpose of illustrating the application
of the formulae for the calculation of an orbit from four observations,
let us take the following normal places of Eurynome derived by
comparing a series of observations with an ephemeris computed from
approximate elements.
Greenwich M. T. a {,
1863 Sept. 20.0 14 30' 35".6 + 9 23' 49".7,
Dec. 9.0 9 54 17 .0 2 53 41 .8,
1864 Feb. 2.0 28 41 34 .1 962 .8,
April 30.0 74 29 58 .9 -f- 19 35 41 .5.
NUMERICAL EXAMPLE. 295
These normals give the geocentric places of the planet referred to the
mean equinox and equator of 1864.0, and free from aberration. For
the mean obliquity of the ecliptic of 1864.0, the American Nautical
Almanac gives
e = 23 27' 24".49,
and, by means of this, converting the observed right ascensions and
declinations, as given by the normal places, into longitudes and lati-
tudes, we get
Greenwich M. T.
1863 Sept. 20.0
Dec. 9.0
1864 Feb. 2.0
April 30.0
A
16 59' 9".42
10 14 17 .57
29 53 21 .99
75 23 46 .90
ft
4- 2 56' 44".58,
1 15 48 .82,
2 29 57 .38,
3 4 44 .49.
These places are referred to the ecliptic and mean equinox of 1864.0,
and, for the same dates, the geocentric latitudes of the sun referred
also to the ecliptic of 1864.0 are
+ 0".60, -f-0".53, -f-0".36, + 0".19.
For the reduction of the geocentric latitudes of the planet to the
point in which a perpendicular let fall from the centre of the earth
to the plane of the ecliptic cuts that plane, the equation (6) 4 gives the
corrections 0".57, 0".38, 0".18, and 0".07 to be applied tc
these latitudes respectively, the logarithms of the approximate dis-
tances of the planet from the earth being
0.02618, 0.13355, 0.29033, 0.44990.
Thus we obtain
t = 0.0, A = 16 59' 9".42, ft = + 2 56' 44".01,
if = 80.0, X =10 14 17 .57, f = 1 15 49 .20,
" = 135.0, A" = 29 53 21 .99, f = 2 29 57 .56,
f" = 223.0, r = 75 23 46 .90, jf n = 3 4 44 .56 ;
and, for the same times, the true places of the sun referred to the
mean equinox of 1864.0 are
=177 0'58".6, logE =0.0015899.
0' =256 58 35.9, log# =9.9932638,
0" =312 57 49 .8, log-R" =9.9937748,
0'"= 40 21 26.8, log #" = 0.0035149,
296
THEORETICAL ASTRONOMY.
From the equations
tan/5'
tan w f = .
tan to" =
. ,., -j^>
sm(A' O)
tan"
. , ' ,
sm(A" O )
tan (A' Q')
tan * = -- - - 7=^1
cosw/
tan(A" O")
tan V = - - - \
cosw"
we obtain
4/ = 113 15' 20".10,
4,"= 76 5617.75,
cos V) = 9.5896777 n ,
log (R sin 4,') =9.9564624,
log (#' cos 4") = 9.3478848,
log (R" sin *") = 9.9823904.
The quadrant in which ty must be taken, is indicated by the condi-
tion that cosij/ and cos (A' O') must have the same sign. The
same condition exists in the case of i//' Then, the formulaB
A = cos ff sin (^ A),
C = cos /9" sin (A'" A"),
^-A'
A - h,
B = cos p' sin (A" jl),
Z> = cos f sin (A'" A'),
,,
O)
, ft
" Q'")
give the following results :
log A = 9.0699254 n ,
log B = 9.3484939,
log h' = 0.2785685 n ,
log a! = 0.8834880 n ,
log d = 0.9012910 n ,
log d' = 0.4650841,
log C = 9.8528803,
log D = 9.9577271,
log h" = 0.1048468,
log a" = 9.9752915 n ,
log c" = 9.7267348 w ,
log d" = 9.9096469 n .
We are now prepared to make the first hypothesis in regard to the
values of P', ', P f/ , and Q". If the elements were entirely un-
known, it would be necessary, in the first instance, to assume for these
quantities the values given by the expressions
NUMERICAL EXAMPLE. 297
then approximate values of r' and r" are readily obtained by means
of the equations (27), (26), and (24) or (25). The first assumed
value of x' to be used in the second member of the first of equations
(27), is obtained from the expression which results from (22) by
putting Q f = and Q" = 0, namely,
x' = *
after which the values of x f and x" will be obtained by trial from
(27). It should be remarked, further, that in the first determination
of an orbit entirely unknown, the intervals of time between the ob-
servations will generally be small, and hence the value of x f derived
from the assumption of Q f = and Q" will be sufficiently ap-
proximate to facilitate the solution of equations (27).
As soon as the approximate values of r f and r" have thus been
found, those of P' and P" must be recomputed from the expressions
With the results thus derived for P' and P", and with the values of
Q f and Q" already obtained, the first approximation to the elements
must be completed.
When the elements are already approximately known, the first
assumed values of P', P", Q', and Q f/ should be computed by means
of these elements. Thus, from
_rV'smQ/' t/) r/sinQ/ 1>)
~ rr" sin (" v) ' ~ rr" sin (v" v) '
,rV"sin(y" f/Q w _ rV ; sin(^ i;Q
~
r r sin (v v) TT sin (v v )
we find n, n' 3 n rf , and n"'. The approximate elements of Eurynom*
givo
v = 322 55' 9".3, logr =0.308327,
it =353 19 26 .3, log/ =0.294225,
v"= 14 45 8.5, log/' =0.296088,
v'"= 47 2332.8, log /" = 0.317278,
298 THEORETICAL ASTRONOMY.
and hence we obtain
log n = 9.653052, log n" = 9.806836,
log n 1 = 9.825408, logn'" = 9.633171.
Then, from
P" = ,
we get
log P' = 9.846216, log Q f = 9.840771,
log P"= 9.807763, log Qf' = 9.882480.
The values of these quantities may also be computed by means of the
equations (41) and (42).
Next, from
, _ P'd' + c' , _ h f
C ~~ = 1 + P" * ~ 1 -f- P''
,_Pd"+c" h"
c o = = x + P n > J - fqrp>
we find
log c ' = 0.541344 n , log/ = 0.047658.,
log C Q " = 9.807665 n , log/" = 9.889385.
Then we have
'
*= *"+"" <
, ,,
tan / = --, tan z" =
,_R' sin V _ a/ ,,_jrsinV_ a"
sin/ cos/' sin/ 7 cos/"
from which to find r' and r rr . In the first place, from
x' = v 1 jR"sinV,
we obtain the approximate value
log x f = 0.242737.
Then the first of the preceding equations gives
log *" = 0.237687.
NUMERICAL EXAMPLE. 299
From this we get
z" = 29 3' 11". 7, log r" = 0.296092 ;
and then the equation for x' gives
logo;' = 0.242768.
Hence we have
z' = 27 20' 59".6, log r' = 0.294249 ;
and, repeating the operation, using these results for x r and r', we get
log x" = 0.237678, log of = 0.242757.
The correct value of log a?' may now be found by means of equation
(28). Thus, in units of the sixth decimal place, we have
= 242768 242737 = + 31, a Q ' = 242757 242768 = 11,
and for the correction to be applied to the last value of log x', in
units of the sixth decimal place,
Therefore, the corrected value is
log af= 0.242760,
and from this we derive
log a" = 0.237681.
These results satisfy the equations for x' and x", and give
z' = 27 21' 1".2, log/ = 0.294242,
z" = 29 3 12 .9, log r" = 0.296087.
To find the curtate distances for the first and second observations,
the formulae are
oo o ,, = -
sin z' sm z"
which give
log p' = 0.133474, log p" = 0.289918.
Then, by means of the equations
300 THEORETICAL ASTRONOMY.
r' cos V cos (f 00 = p' cos (X K 9
r' cos V sin (? ') = p' sin (A' QO,
/ sin V = p' tan p,
r" cos ft" cos (r Q'O = P " cos (A" 0") J2",
r" cos b" sin (/" Q") = p" sin (A" "),
r" sin ft" = //'tan/S",
we find the following heliocentric places :
r = 37 35' 26".4, log tan ft' = 8.182861 n , log / = 0.294243,
r = 58 5815.3, logtanft" = 8.634209 n , log r" = 0.296087.
The agreement of these values of log r' and log r" with those obtained
directly from x' and x" is a partial proof of the numerical calcula-
tion.
From the equations
tan i sin ( (I" -f /') ) = J (tan ft" + tan ftO sec j (r 7),
tan i cos Q (J 1 -f I') & ) = J (tan ft" tan ftO cosec i (/" 0,
tan (^ &0 tan(r ft)
tan w = - - - r-^ , tan u =
cos i cos t
we obtain
& = 206 42' 24".0, i = 4 36' 47".2,
u r =190 55 6 .6 u" = 212 20 53 .5.
Then, from
we get
log n" = 9.806832, logw =9.653048,
log n' = 9.825408, log n'" = 9.633171,
and the equations
r sin ((u r - u) + i (u" - w')) = /+ J' V ' sin J (i*" - 0,
r cos ((i*' - u) + i (w" - tO) = ^= cos i (u" - tif ),
r"' sin ((*"' - w'O + J (u" - wO) - , sin j (t" - wO,
r"' cos ((u'" - i/O -f i (u" - 1*0) = cos i (u" -
NUMERICAL EXAMPLE. 301
give
logr = 0.308379, u = 160 30 ; 57".6,
log r"' = 0.317273, t*'"=244 5932.5.
Next, by means of the formulae
tan (I & ) = cos i tan u, tan 6 = tan i sin (/ ft ),
tan (Y" ) = cos i tan w'", tan 6'" = tan i sin ('" ft ),
^ cos (A o ) = f cos 6 cos (7 0) + -K,
jo sin (A Q ) = r cos 6 sin (7 O),
/o tan /? = r sin & ;
p ' C os (A"' O'") = r'" cos 6'" cos (f" "0 + 12"',
p' s in (A'" O'") = r'" cos 6"' sin (f" 0'"),
//"tan/5"' = r"'sin&'",
we obtain
J = 7 16' 51".8, r = 91 37' 40".0,
b = + 1 32 14 .4, 6"' = 4 10 47 .4,
A = 16 59 9 .0, /" 75 23 46 .9,
/? = + 2 56 40 .1, /?'" = 3 4 43 .4,
log P = 0.025707, log p'" = 0.449258.
The value of X" f thus obtained agrees exactly with that given by
observation, but / differs O r/ .4 from the observed value. This differ-
ence does not exceed what may be attributed to the unavoidable
errors of calculation with logarithms of six decimal places. The
differences between the computed and the observed values of /9 and
ft" show that the position of the plane of the orbit, as determined
by means of the second and third places, will not completely satisfy
the extreme places.
The four curtate distances which are thus obtained enable us, in
the case of an orbit entirely unknown, to complete the correction for
aberration according to the equations (40).
The calculation of the quantities which are independent of P',
P ff j Q f j and Q /f , and which are therefore the same in the successive
hypotheses, should be performed as accurately as possible. The
s*
value of -^-> required in finding x" from x', may be computed
directly from
jf .
the values of 77 and 77 being found by means of the equations (29);
302 THEORETICAL ASTRONOMY.
and a similar method may be adopted in the case of j, r > Further,
in the computation of x f and x", it may in some cases be advisable
to employ one or both of the equations (22) for the final trial. Thus,
in the present case, x" is found from the first of equations (27) by
means of the difference of two larger numbers, and an error in the
last decimal place of the logarithm of either of these numbers affects
in a greater degree the result obtained. But as soon as r" is known
Q"
so nearly that the logarithm of the factor 1 -f -^ remains unchanged,
the second of equations (22) gives the value of x" by means of the
sum of two smaller numbers. In general, when two or more for-
mulae for finding the same quantity are given, of those which are
otherwise equally accurate and convenient for logarithmic calculation,
that in which the number sought is obtained from the sum of smaller
numbers should be preferred instead of that in which it is obtained
by taking the difference of larger numbers.
The values of r, r f , r /f , r'", and u, u f y u ff , u fff , which result from
the first hypothesis, suffice to correct the assumed values of P f , P",
', and Q". Thus, from
"r 77 " '/ , /r 777 "
sin Y cos G = sin (u" u'), sin /' cos G" = sin A (u r u),
sin f sin G = cos J (u" u') cos 2/, sin /' sin G" cos (u r u) cos 2/',
cos f = cos | (u" u') sin 2/, cos f" = cos \ (V u) sin 2,%",
sin f cos G'" = sin I (u" r w"),
sin f sin G"' = cos (u" f u"} cos 2/"
cos /" = cos J (u"' u") sin 2/" ;
T 2 COS 6 / __ T" 2 COS 6 /'
m =
COS?'
m
r' /s cos 3 /'
t+j + ? *+/'+*
in connection with Tables XIII. and XIV. we find s, s", and '".
The results are
NUMEKICAL EXAMPLE.
303
log T = 9.9759441,
/ = 45 3'39".l,
r=-10 42 55 .9,
log m = 8.186217,
log; = 7.948097,
log s = 0.0085248,
log r"= 0.1386714,
7"=4432' 1".4,
/'IS 13 45 .0,
log m"= 8.516727,
log/'= 8.260013,
log "== 0.0174621,
log T"'= 0.1800641,
/"== 45 41' 55".2,
y"'=16 22 48 .5.
log m'"= 8.590596,
log/"= 8.325365,
log s'"= 0.020406?.
Then, by means of the formulae
a -*
V - 2 ,,
y, = , nl' _
2 ss "f r y // cog ^ (yftr u n^ cog ^ (jjin ^ CQg ^ ^/ % /y
we obtain
log P' = 9.8462100,
log P" = 9.8077615,
log ' = 9.8407536,
log " = 9.8824728,
with which the next approximation may be completed.
We now recompute c ', c/', /', /", x f , x n ', &c. precisely as already-
illustrated; and the results are
log c ' = 0.5413485 n ,
log/ = 0.0476614 n ,
log x' = 0.2427528,
z f = 27 21' 2".71,
log/ =0.2942369,
log />' = 0.1334635,
log n = 9.6530445,
log n' = 9.8254092,
log c " = 9.8076649 n ,
log/" = 9.8893851,
log x" = 0.2376752,
z" = 29 3' 14".09,
logr" =0.2960826,
log p" = 0.2899124,
log n" = 9.8068345,
log ri" = 9.6331707.
Then we obtain
/' = 37 35' 27".88,
I" 58 58 16 .48,
log tan V = 8.1828572 n ,
log tan b" = 8.6342073 n ,
log/ = 0.2942369,
log /' = 0.2960827.
These results for log r f and log r" agree with those obtained directly
from z f and z", thus checking the calculation of ty and ty f and of
the heliocentric places.
Next, we derive
206 42' 25".89,
55 6.27,
t = 4 36' 47".20,
2052.96,
304 THEORETICAL ASTRONOMY.
and from u ff u', r', r", n, n", n f , and n" f , we obtain
logr =0.3083734, u = 160 30' 55".45,
log/" =0.3172674, w'"=244 5931.98.
For the purpose of proving the accuracy of the numerical results,
we compute also, as in the first approximation,
/= 716'51".54, /'" = 91 37' 41".20,
b = + 1 32 14 .07, b'"= 4 10 47 .36,
A= 16 59 9.38, A'" = 75 2346.99,
= + 2 56 39 .54, p"= 3 4 43 .33,
log P = 0.0256960, log p" r = 0.4492539.
The values of ^ and \ tn thus found differ, respectively, only 0".04
and O r/ .09 from those given by the normal places, and hence the
accuracy of the entire calculation, both of the quantities which are
independent of P', P" ', ', and ", and of those which depend on
the successive hypotheses, is completely proved. This condition,
however, must always be satisfied whatever may be the assumed
values of P', P", Q f , and ".
From TJ r'j u, it', &c., we derive
log a = 0.0085254, log s" = 0.0174637, log s m = 0.0204076,
and hence the corrected values of P', P", ', and Q" become
logP' = 9.8462110, log q = 9.8407524,
log P" = 9.8077622, log " = 9.8824726.
These values differ so little from those for the second approximation,
the intervals of time between the observations being very large, that
a further repetition of the calculation is unnecessary, since the results
which would thus be obtained can differ but slightly from those
which have been derived. We shall, therefore, complete the deter-
mination of the elements of the orbit, using the extreme places.
Thus, from
sin YQ cos G = sin J (u tn u),
sin Y O sin G = cos J (u" f u) cos 2/ ,
cos YO cos (u"' u) sin 2/ ,
5 i i
~r Jo
? - _
"
m _ m a
~~ 7 ~
NUMERICAL EXAMPLE. 305
we get
log r = 0.5838863, log tan G Q = 8.0521953 W>
n = 42 14' 30".17, log m = 9.7179026,
log s 2 = 0.2917731, log x = 8.9608397.
The formula
s rr m sm(u m t*)\
gves
log p = 0.371 2401;
and if we compute the same quantity by means of
grV'sm(X' iQ \ 2 / g'Vr'sinK u) V I s'VV'sin (u'"-
u"} \
I*
the separate results are, respectively, 0.3712397, 0.3712418, and
0.3712414. The differences between these results are very small, and
arise both from the unavoidable errors of calculation and from the
deviation of the adopted values of P', P x/ , Q', and Q" from the
limit of accuracy attainable with logarithms of seven decimal places.
A variation of only O r/ .2 in the values of u' u and u fff u ff wil?
produce an entire accordance of the particular results.
From the equations
smi(V" u) /T,,
C 8 p = sin K"' -"
P
cos ) tan (45 ?),
t&n%E' = tani(X oO tan (45 ?>),
tan IE" = tan J O" - o) tan (45 ' p),
tan E'" = tan (u" r o) tan (45 ?>)
from which the results are
E = 329 11' 46".01, .E" 12 5' 33".63,
' = 354 29 11 .84, '" = 39 34 34 .65.
The value of J (E" f E) thus derived differs only 0".03 from that
obtained directly from x .
For the mean anomalies, we have
which give
, M" = E" esmE",
= E r e sin E', M'" = E'" e sin E'",
M = 334 55' 39".32, M" = 9 44' 52".82,
M' = 355 33 42 .97, JJf'" = 32 26 44 .74.
Finally, if M Q denotes the mean anomaly for the epoch T= 1864
Jan. 1.0 mean time at Greenwich, from
M Q = M fi (t T) =M t p(fT)
= M" fJL (f T} = M'" p. (f" T),
we obtain the four values
Jf = l29'39".40
39 .49
39 .40
39 .40,
the agreement of which completely proves the entire calculation of
the elements from the data. Collecting together the several results,
we have the following elements :
NUMERICAL EXAMPLE. 307
Epoch = 1864 Jan. 1.0 Greenwich mean time.
M = 1 29' 39".42
=
?= 11 15 52 .22
log a = 0.3881359
log ; = 2.9678027
A* = 928".54447.
102. The elements thus derived completely represent the four ob-
served longitudes and the latitudes for the second and third places,
which are the actual data of the problem ; but for the extreme lati-
tudes the residuals are, computation minus observation,
A = 4".47, A,?'" = + 1".23.
These remaining errors arise chiefly from the circumstance that the
position of the plane of the orbit cannot be determined from the
second and third places with the same degree of precision as from
the extreme places. It would be advisable, therefore, in the final
approximation, as soon as p f , p", n, n' f , n f , and n rrr are obtained, to
compute from these and the data furnished directly by observation
the curtate distances for the extreme places. The corresponding
heliocentric places may then be found, and hence the position of the
plane of the orbit as determined by the first and fourth observations.
Thus, by means of the equations (37) and (38), we obtain
log P = 0.0256953, log p'" = 0.4492542.
With these values of p and //", the following heliocentric places are
obtained :
I = 716'51".54, log tan b =8.4289064, logr =0.3083732,
l f " = 91 37 40 .96, log tan b'" = 8.8638549 n , log r"' = 0.3172678.
Then from
tan i sin (J (f" + 1) ft) = J (tan V" -f tan 6) sec J (f" Z),
tan i cos Q (f" -f & ) = i (tan b"' tan 6) cosec J (f I),
we get
^ = 206 42' 45".23, i = 4 36' 49".76.
For the arguments of the latitude the results are
u = 160 30' 35".99, u'" = 244 59' 12".53.
398 THEORETICAL ASTRONOMY.
The equations
tan V = tan i sin (I 1 Q ),
tan V = tan i sin (I" ),
give
log tan &' = 8.1827129^ log tan b" = 8,6342104 n ,
and the comparison of these results with those derived directly from
p' and p" exhibits a difference of -f 1".04 in &', and of 0".06 in
6". Hence, the position of the plane of the orbit as determined from
the extreme places very nearly satisfies the intermediate latitudes.
If we compute the remaining elements by means of these values
of r, r f ", and u, u" r , the separate results are :
log tan G = 8.0522282 n , log m = 9.7179026,
log sj = 0.2917731, log x = 8.9608397,
log j9 = 0.3712405, I (E" E) = 17 35' 42".12,
log (a cos ?) = 0.3796884, log cos 9 i, u, and u".
Then, if the value of a computed from the last result for u" u
differs from the last assumed value, a further repetition of the calcu-
CIKCULAB CEBIT. 313
latiori becomes necessary. But when three successive approximate
values of a have been found, the correct value may be readily inter-
polated according to the process already illustrated for similar cases.
As soon as the value of a has been obtained which completely
satisfies equation (5), this result and the corresponding values of &,
i, and the argument of the latitude for a fixed epoch, complete the
system of circular elements which will exactly satisfy the two observed
places. If we denote by u the argument of the latitude for the epoch
T, we shall have, for any instant t,
u being the mean or actual daily motion computed from
Jc
The value of u thus found, and r = a, substituted in the formulse for
computing the places of a heavenly body, will furnish the approxi-
mate ephemeris required.
The corrections for parallax and aberration are neglected in the
first determination of circular elements ; but as soon as these approxi-
mate elements have been derived, the geocentric distances may be
computed to a degree of accuracy sufficient for applying these cor-
rections directly to the observed places, preparatory to the determi-
nation of elliptic elements. The assumption of r f = a will also be
sufficient to take into account the term of the second order in the first
assumed value of P, according to the first of equations (98) 4 .
104. When approximate elements of the orbit of a heavenly body
have been determined, and it is desired to correct them so as to satisfy
as nearly as possible a series of observations including a much longer
interval of time than in the case of the observations used in finding
these approximate elements, a variety of -methods may be applied.
For a very long series of observations, the approximate elements
being such that the squares of the corrections which must be applied
to them may be neglected, the most complete method is to form the
equations for the variations of any two spherical co-ordinates which
fix the place of the body in terms of the variations of the six ele-
ments of the orbit; and the differences between the computed places
for different dates and the corresponding observed places thus furnish
equations of condition, the solution of which gives the corrections to
be applied to the elements But when the observations do not in-
314 THEORETICAL ASTRONOMY.
elude a very long interval of time, instead of forming the equations
for the variations of the geocentric places in terms of the variations
of the elements of the orbit, it will be more convenient to form the
equations for these variations in terms of quantities, less in number,
from which the elements themselves are readily obtained. If no as-
sumption is made in regard to the form of the orbit, the quantities
which present the least difficulties in the numerical calculation are
the geocentric distances of the body for the dates of the extreme
observations, or at least for the dates of those which are best adapted
to the determination of the elements. As soon as these distances are
accurately known, the two corresponding complete observations are
sufficient to determine all the elements of the orbit.
The approximate elements enable us to assume, for the dates t and
t n , the values of A and J"; and the elements computed from these
by means of the data furnished by observation, will exactly represent
the two observed places employed. Further, the elements may be
supposed to be already known to such a degree of approximation that
the squares and products of the corrections to be applied to the
assumed values of A and A" may be neglected, so that we shall have,
for any date,
.. da, . da
COS S Aa = COS d - A J -j- COS d -=--- rr A A ,
d d fR\
d* d3
** = dA^ + ^F AJ '
If, therefore, we compare the elements computed from A and A" with
any number of additional or intermediate observed places, each ob-
served spherical co-ordinate will furnish an equation of condition for
the correction of the assumed distances. But in order that the equa-
tions (6) may be applied, the numerical values of the partial differen-
tial coefficients of a and d with respect to A and d" must be found.
Ordinarily, the best method of effecting the determination of these is
to compute three systems of elements, the first from A and A", the
second from A + D and J", and the third from A and A" -f D", D
and D" being small increments assigned to A and A" respectively.
] f now, for any date t' y we compute a/ and d f from each system of
elements thus obtained, we may find the values of the differential
coefficients sought. Thus, let the spherical co-ordinates for the time
t f computed from the first system be denoted by a' and 3 f ; those
computed from the second system of elements, by a/ -f a sec 3 f and
6' + d; and those from the third system, by a'+ a!' sec d' and d'-\- d".
Then we shall have
VARIATION OF TWO GEOCENTRIC DISTANCES. 315
.. da! a dd d
,,a__ -..
dA"~ D'" dA"~ D"'
and the equations (6) give
COS d r Aa' = ~ A A -f -^77 A J",
d d" (8)
In the same manner, computing the places for various dates, for
which observed places are given, by means of each of the three systems
of elements, the equations for the correction of A and A" , as deter-
mined by each of the additional observations employed, may be
formed.
105. For the purpose of illustrating the application of this method,
let us suppose that three observed places are given, referred to the
ecliptic as the fundamental plane, and that the corrections for parallax,
aberration, precession, and nutation have all been duly applied. By
means of the approximate elements already known, we compute the
values of A and A" for the extreme places, and from these the helio-
centric places are obtained by means of the equations (71) 3 and (72) 3 ,
writing A cos /9 and A" cos ft" in place of p and p". The values of
&, i, u, and u" will be obtained by means of the formulae (76) 3 and
(77) 3 ; and from r, r" and u" u the remaining elements of the
orbit are determined as already illustrated. The first system of ele-
ments is thus obtained. Then we assign an increment to J, which
we denote by D, and with the geocentric distances A -f- D and A"
we compute in precisely the same manner a second system of ele-
ments. Next, we assign to A" an increment D", and from A and
A" + D" a third system of elements is derived. Let the geocentric
longitude and latitude for the date of the middle observation com-
puted from the first system of elements be designated, respectively,
by V and /9/ ; from, the second system of elements, by V and /? 2 ' ;
and from the third system, by ^/ and /9 3 '. Then from
It ( -\ I i l\ of fjff O f O f \** s
we compute a, a", d, and d" ', and by means of these and the values
of D and D" we form the equations
316 THEORETICAL ASTRONOMY.
= CM
for the determination of the corrections to be applied to the f.rst
assumed values of A and A n ', by means of the differences between
observation and computation. The observed longitude and latitude
being denoted by X and /?', respectively, we shall have
COS f A/ = (/ A/) COS /?',
ff = ff B' (11)
for finding the values of the second members of the equations (10),
and then by elimination we obtain the values of the corrections A J
and AJ" to be applied to the assumed values of the distances.
Finally, we compute a fourth system of elements corresponding to
the geocentric distances A -f A A and A" -f A A" either directly from
these values, or by interpolation from the three systems of elements
already obtained ; and, if the first assumption is not considerably in
error, these elements will exactly represent the middle place. It
should be observed, however, that if the second system of elements
represents the middle place better than the first system, A/ and /? 2 r
should be used instead of ^/ and /9/ in the equations (11), and, in
this case, the final system of elements must be computed with the
distances A -f- D + A A and A" + A A". Similarly, if the middles
place is best represented by the third system of elements, the cor-
rections will be obtained for the distances used in the third hy-
pothesis.
If the computation of the middle place by means of the final ele-
ments still exhibits residuals, on account of the neglected terms of
the second order, a repetition of the calculation of the corrections
A/f and AJ", using these residuals for the values of the second
members of the equations (10), will furnish the values of the dis-
tances for the extreme places with all the precision desired. The
increments D and D" to be assigned successively to the first assumed
values of A and A" may, without difficulty, be so taken that the
true elements shall differ but little from one of the three systems
computed ; and in all the formulae it will be convenient to use, in-
stead of the geocentric distances themselves, the logarithms of these
distances, and to express the variations of these quantities in units
of the last decimal place of the logarithms.
These formulae will generally be applied for the correction of
VARIATION OF T\VO GEOCENTRIC DISTANCES. 317
approximate elements by means of several observed places, which
may be either single observations or normal places, each derived from
several observations, and the two places selected for the computation
of the elements from A and A" should not only be the most accurate
possible, but they should also be such that the resulting elements are
not too much affected by small errors in these geocentric places.
They should moreover be as distant from each other as possible, the
other considerations not being overlooked. When the three systems
of elements have been computed, each of the remaining observed
places will furnish two equations of condition, according to equations
(10), for the determination of the corrections to be applied to the
assumed values of the geocentric distances ; and, since the number
of equations will thus exceed the number of unknown quantities,
the entire group must be combined according to the method of least
squares. Thus, we multiply each equation by the coefficient of AJ
in that equation, taken with its proper algebraic sign, and the sum
of all the equations thus formed gives one of the final equations
required. Then we multiply each equation by the coefficient of AJ"
in that equation, taken also with its proper algebraic sign, and the
sum of all these gives the second equation required. From these
two final equations, by elimination, the most probable values of AJ
and AZ/ ;/ will be obtained ; and a system of elements computed with
the distances thus corrected will exactly represent the two funda-
mental places selected, while the sum of the squares of the residuals
for the other places will be a minimum. The observations are thus
supposed to be equally good; but if certain observed places are
entitled to greater influence than the others, the relative precision
of these places must be taken into account in the combination of the
equations of condition, the process for which will be fully explained
in the next chapter.
When a number of observed places are to be used for the correction
of the approximate elements of the orbit of a planet or comet, it w* ! ll
be most convenient to adopt the equator as the fundamental plane.
In this case the heliocentric places will be computed from the assumed
values of A and J", and the corresponding geocentric right ascensions
and declinations by means of the formulae (106) 3 and (107) 3 ; and the
position of the plane of the orbit as determined from these by means
of the equations (76) 3 will be referred to the equator as the funda-
mental plane. The formation of the equations of condition for the
corrections AJ and A A" to be applied to the assumed values of the
distances will then be effected precisely as in the case of I and /?, the
318 THEORETICAL ASTRONOMY.
necessary changes being made in the notation. In a similar manner,
the calculation may be effected for any other fundamental plane which
may be adopted.
It should be observed, further, that when the ecliptic is taken as
the fundamental plane, the geocentric latitudes should be corrected
by means of the equation (6) 4 , in order that the latitudes of the sun
phall vanish, otherwise, for strict accuracy, the heliocentric places
must be determined from A and A" in accordance with the equations
(39):-
106. The partial differential coefficients of the two spherical co-
ordinates with respect to A and A" may be computed directly by
means of differential formula; but, except for special cases, the
numerical calculation is less expeditious than in the case of the indi-
rect method, while the liability of error is much greater. If we
adopt the plane of the orbit as determined by the approximate values
of A and A" as the fundamental plane, and introduce / as one of the
elements of the orbit, as in the equations (72) 2 , the variation of the
geocentric longitude 6 measured in this plane, neglecting terms of the
second order, depends on only four elements; and in this case the
differential formulae may be applied with facility. Thus, if we ex-
press r and v in terms of the elements ^, M Q) and /*, we shall have
and
or
In like
dr
dA
dv
dA
_dr
dtp
dv
dtp
dr
dM
dr
h d,u
dv
dp.
dp.
' dA
'dA
->
dv
dA
' dA'
dp.
dx
dM Q
d
'dA
dM
dA
1 dp.
__ _
d? ' dA dM ' dA " dfi ' dA'
dA ~ d TT> an( i - n -- are
dA dA dA dA
known, the equations necessary for finding the differential coefficients
of the elements , = o.
dA aA
To find -7-7- and - 7 . , from the equations
dA dA
A cos ^ cos = x -f- X,
A cos >? sin = y -{- F",
in which ^ is the geocentric latitude in reference to the plane of the
orbit computed from A and A" as the fundamental plane, and X, Y
the geocentric co-ordinates of the sun referred to the same plane, we
get
dx = cos 7] cos 6 dA,
dy = cos T? sin 6 dA,
or, substituting for dx and dy their values given by (73) 2 ,
cos t] cos 6 dA = cos u dr r sin u d (v -f- /),
cos r) sin O d A = sin udr -\- r cos it c? (v -f~ /).
Eliminating, successively, c? (v + #) and c?r, we get
dr , n x
- = cos ^ cos (6 u),
d
I . , .
= - cos ^ sm (6 u).
7 -
dA r
Therefore, we shall have
dy , dv dv dv
dip dA ' dM Q dA dp. dA ,_
d X _ _dif_ d?_ d^_ dM^,dv^_ ^_
dA H " ' dA + ^jf o jj ' ' ^ ' ^j
_dr^_ ^__
' ~
^ dA dif dA dfji dA
and if we compute the numerical values of the differential coefficients
of r, r", v, and v tr with respect to the elements Sm(t U } '
dr" d dr" dM , dr"
dA" r dft dA"
Since the geocentric latitude ^ is affected chiefly by a change of the
position of the plane of the orbit, while the variation of the longitude
6 is independent of & and i when the squares and products of the
variations of the elements are neglected, if we determine the elements
which exactly represent the places to which A and A" belong, as well
as the longitudes for two additional places, or, if we determine those
which satisfy the two fundamental places and the longitudes for any
number of additional observed places, so that the sum of the squares
of their residuals shall be a minimum, the results thus obtained will
very nearly satisfy the several latitudes.
Let 0' denote the geocentric longitude of the body, referred to the
plane of the orbit computed from A and A" as the fundamental plane,
for the date t r of any one of the observed places to be used for cor-
recting these assumed distances. Then, to find the partial differential
coefficients of 6 f with respect to A and A", we have
, do' .do' dy .do' dy . dtf dM Q
COS rj -7-7- = COS TJ -j ~- -f COS y 3 -j-j- + COS ff -7-=- j-f
dA d/ dA d , and , respectively, we compute the values of 0, 37, and f for the
dates of the several places to be employed. Then the residuals for
each of the observed places are found from the formulae
cos f) A0 = sin Y A<5 -f- cos f cos 8 Aa,
Aiy = cos Y A# sin f cos d Aa,
the values of Aa and A for each place being found by subtracting
from the observed right ascension and declination, respectively, the
right ascension and declination computed by means of the elements
derived from A and A". The values of 9 f], and f being required
only for finding cos y A#, A/7, and the differential coefficients of d and
ry, with respect to the elements of the orbit, need not be determined
with great accuracy.
Next, we compute -^ and - k , . from equations (12), and from
,+ n ^ .1 i dr dr" dv dv" dr , / i i
(16), the values of , , , -= , -r^, &c,, by means of which,
d d , A!/ O , and
may be obt
we compute
may be obtained. If, from the values of -- - and -p
VARIATION OF TWO GEOCENTRIC DISTANCES. 323
and apply these corrections to the values of v and v" found from J
and J", we obtain the true anomalies corresponding to the distances
A -}- A J and J" + A J". The corrections to be applied to the values
of r and r" derived from A and A" are given by
If A J and AJ" are expressed in seconds of arc, the corresponding
values of Ar and Ar" must be divided by 206264.8. The corrected
results thus obtained should agree with the values of r and r" com-
puted directly from the corrected values of v, v", p, and e by means
of the polar equation of the conic section. Finally, we have
dz = sin rj dd,
and similarly for dz ff and the last of equations (73) 2 gives
r sin u Ai' r cos u sin i' A & ' == sin 7 A J,
r" sin u" Ai' r" cos u" sin i' A ft ' sin ^ A J",
from which to find AI V and A & ', u and u" being the arguments of
the latitude in reference to the equator. We have also, according to
(72)*
A a/ = A/ COS i' A & ',
4^ = A/ + 2 Sin 1 t* A &',
from which to find the corrections to be applied to a) 1 and n f . The
elements which refer to the equator may then be converted into those
for the ecliptic by means of the formula which may be derived from
(109)! by interchanging & and &' and 180 i' and i.
The final residuals of the longitudes may be obtained by substi-
tuting the adopted values of A J and A A" in the several equations of
condition, or, which affords a complete proof of the accuracy of the
entire calculation, by direct calculation from the corrected elements ;
and the determination of the remaining errors in the values of rj will
show how nearly the position of the plane of the orbit corresponding
to the corrected distances satisfies the intermediate latitudes.
Instead of and i, the equations may be formed by
means of which the corrections to be applied to the assumed values
of the two geocentric distances, or to those of & and i, will be
obtained.
110. The formulae which have thus far been given for the correc-
tion of an approximate orbit by varying the geocentric distances,
depend on two of these distances when no assumption is made in
regard to the form of the orbit, and these formulae apply with equal
facility whether three or more than three observed places are used.
But when a series of places can be made available, the problem may
be successfully treated in a manner such that it will only be necessary
to vary one geocentric distance. Thus, let x, y, z be the rectangular
heliocentric co-ordinates, and r the radius-vector of the body at the
time t, and let X, F, Z be the geocentric co-ordinates of the sun at
the same instant. Let the geocentric co-ordinates of the body be
designated by x 09 y Q , Z Q , and let the plane of the equator be taken as
the fundamental plane, the positive axis of x being directed to the
vernal equinox. Further, let p denote the projection of the radius-
vector of the body on the plane of the equator, or the curtate dis-
tance with respect to the equator; then we-shall have
x = p cos a, y = p sin a, z p tan d. (32)
If we represent the right ascension of the sun by A, and its declina-
tion by _D, we also have
VARIATION OF ONE GEOCENTRIC DISTANCE. 329
Z=RsmD. (33)
The fundamental equations for the undisturbed motion of the planet
or comet, neglecting its mass in comparison with that of the sun, are
but since
and, neglecting also the mass of the earth,
these become
Substituting for a? , y , and 2; their values in terms of a and d, and
putting
we get
+ ^co Sa + , = 0,
+ sina + , = 0, ' (36)
Differentiating the equations (32) with respect to t, we find
dx n dp . da
-d?= c
f
31)0 THEORETICAL ASTRONOMY.
Differentiating again with respect to t, and substituting in the equa
tions (36) the values thus found, the results are
If we multiply the first of these equations by sin a, and the second
by cos a, and add the products, we obtain
d' 2 a
j gBin.-?COB.-/-
dt 2 rff
57
Now, from (35) we get
sin a ^ cos a = & 2 1 -= -- 3 I -K cos D sin (a A),
and the preceding equation becomes
The value of -vr thus found is independent of the differential co-
dp
efficients of d with respect to t. To find another value of -J-, using
dt
all three of equations (38), we multiply the first of these equations
by sin A tan 8, the second by cos A tan , and the third by
sin (a A). Then, adding the products, since sin A = TJ cos A,
the result is
from which we get
% - cot (a - A) % + **( 2 ' + cot % ) + i cot*
a/> _ , . In this case, and H will be referred to the
equator as the fundamental plane. The angles ^ and i//' will be
obtained from the equations (102) 3 , or from equations of the form
350 THEORETICAL ASTRONOMY.
of (26), and finally the auxiliary quantities A, B, B", &c. will be
obtained from (51) 3 , writing d and d /f in place of ft and ft", respect-
ively.
As soon as these auxiliary quantities have been determined, by
means of (52) 3 the value of % must be found which will exactly
satisfy equation (65). To effect this, we first compute e from
sin ^e =
and, if it be required, we also find 3 from
using approximate values of r + r" and x. Then we find Q from
(66), and Ar ' from (76) or from (78), the logarithms of the auxiliary
quantities B Q and N being found by means of Table XV. with the
argument e. The value of r/ having been found from (77), the
equations (73) and (74), in connection with Table XI., enable us to
obtain a closer approximation to the correct value of x. With this
we compute new values of r and r ff , and repeat the determination
of x. A few trials will generally give the correct result, and these
trials may be facilitated by the use of the formula (67) 3 . It will be
observed, also, that Q and Ar ' are very slightly changed by a small
change in the values of r -j- r fl and x, so that a repetition of the
calculation of these quantities pnly becomes necessary for the final
trial in finding the value of x which completely satisfies the equa-
tions (52) 3 and (65). When the value of a is such that the values
of Q and ^V exceed the limits of Table XV., the equation (61) may
be employed, and, in the case of hyperbolic motion, when Q and '
exceed the limits of Table XVI., we may employ the complete ex-
pression for the time r' in terms of m and n as given by (79).
The values of r, r /r , and x having thus been found, the equations
will determine the curtate distances p and />". When the equator is
the fundamental plane, we have
p = A cos 9, p" = J" cos d".
From p, p", and the corresponding geocentric spherical co-ordinates,
the radii- vectores and the heliocentric spherical co-ordinates I, I", 6,
and b" will be obtained, and thence &, i, u, u ff , and the remaining
VARIATION OF THE SEMI-TRANSVERSE AXIS. 351
elements of the orbit, as already illustrated. In the case of elliptic
motion, if we compute the auxiliary quantities e and d by means of
the equations (60), we shall have
e cos J (E" + J) = cos J (e + d),
from which e and \(E fr -\~ E} may be found, and hence, since
i(J" J) = i(e ), we derive .# and J0". The values of q and
i? may then be found directly from these and quantities already
obtained. Thus, the last of equations (43)! gives
cos ^v cos \E cos ^v" cos ^E"
V~q ' 1/r Vq 1/7'
Multiplying the first of these expressions by sin Ju", and the second
oy sin J0, adding the products, and reducing, we obtain
smi(v" v)sin%v _ cos j (v" v) cos %E cos \E"
l/q Vr V7' '
Therefore, we shall have
_1_ g . n _ cosjJg cos E"
V~q Sm 3V " ~ 1/r tan J (M" *) T/7' sin J (t*" uj
1 cos IE
-= cos Iv = -,
VQ V
from which q and v may be found as soon as cos \E and cos \E" are
known. In the case of parabolic motion the eccentric anomaly is
equal to zero, and these equations become identical with (92) 3 . The
angular distance of the perihelion from the ascending node will be
obtained from
(a = u V.
Since r = a ae co$E, and q = a(l e), we have
and hence
1-1
a
(89)
352 THEORETICAL ASTRONOMY.
When the eccentricity is nearly equal to unity, the value of q given
by approximate elements will be sufficient to compute cos^ and
cos^E" by means of these equations, and the results thus derived
will be substituted in the equations (88), from which a new value of
q results. If this should differ considerably from that used in com-
puting cos \E and cos \E", a repetition of the calculation will give
the correct result.
In the case of hyperbolic motion, although E and E f> are imagi-
nary, we may compute the numerical values of cos^ and cos^E''
from the equations (89), regarding a as negative, and the results will
be used for the corresponding quantities in (88) in the computation
of q and v for the hyperbolic orbit.
Next, we compute a second system of elements from M and/ -f Sf,
and a third system from M -f- dM and/, df and dM denoting the
arbitrary increments assigned to / and M respectively. The com-
parison of these three systems of elements with additional observed
places of the comet, will enable us to form the equations of condition
for the determination of the most probable values of the corrections
&M and A/ to be applied to M and /respectively. The formation of
these equations is effected in precisely the same manner as in the case
of the variation of the geocentric distances or of & and i, and it does
not require any further illustration. The final elements will be ob-
tained from M-\- &M, and/-f- A/, either directly or by interpolation.
We may remark, further, that it will be convenient to use log M as
the quantity to be corrected, and to express the variations of log M
in units of the last decimal place of the logarithms.
When the orbit differs very little from the parabolic form, it will
be most expeditious to make two hypotheses in regard to M, putting
in each case = 0, and only compute elliptic or hyperbolic elements
in the third hypothesis, for which we use J/and f=df. The first
and second systems of elements will thus be parabolic.
120. Instead of M and - we may use A and - as the quantities to
a a
be corrected. In this case we assume an approximate value of J by
means of elements already known, and by means of (96) 3 , (98) 3 , (102),,
and (103) 3 , we compute the auxiliary quantities C, B, B", &c., re-
quired in the solution of the equations (104) 3 . We assume, also, an
approximate value of J r/ and compute the corresponding value of r",
the value of r having been already found from the assumed value of
J. Then, by trial, we find the value of x which, in connection witn
EQUATIONS OF CONDITION. 353
the assumed value of -, will satisfy the equations (104) 3 and (65) or
(61). The corresponding value of A" is given by
J" = e i/x 2 C\
When A" has thus been determined, the heliocentric places will be
obtained by means of the equations (106) 3 and (107) 3 , and, finally,
the corresponding elements of the orbit will be computed. If the
ecliptic is taken as the fundamental plane, we put D = Q, A = O,
and write ^ and /9 in place of a and d respectively.
If we now compute a second system of elements from A + d J and
f = -, and a third system from J and /+/", the comparison of the
three systems of elements with additional observed places will furnish
the equations of condition for the determination of the corrections
A J and A/ to be applied to J and - respectively.
When the eccentricity is very nearly equal to unity, we may as-
sume /=0 for the first and second hypotheses, and only compute
elliptic or hyperbolic elements for the third hypothesis.
121. The comparison of the several observed places of a heavenly
body with one of the three systems of elements obtained by varying
the two quantities selected for correction, or, when the required dif-
ferential coefficients are known, with any other system of elements
such that the squares and products of the corrections may be neg-
lected, gives a series of equations of the form
mx -f ny p,
m'x -f- n'y =p', &c.,
in which x and y denote the final corrections to be applied to the two
assumed quantities respectively. The combination of these equations
which gives the most probable values of the unknown quantities, is
effected according to the method of least squares. Thus, we multiply
each equation by the coefficient of x in that equation, and the sum
of all the equations thus formed gives the first normal equation.
Then we multiply each equation of condition by the coefficient of y
in that equation, and the sum of all the products gives the second
normal equation. Let these equations be expressed thus:-
[mm] x + [mn\ y = [mp],
\mri\ x + \nri\ y = \np],
23
354 THEORETICAL ASTRONOMY.
in which [mra]=m 2 -{-m' 2 -|-m'' 2 -}-&c., [mn]=mw+mV-f ra'W'-f&c.,
and similarly for the other terms. These two final equations give,
by elimination, the most probable values of x and y, namely, those
for which the sum of the squares of the residuals will be a minimum.
It is, however, often convenient to determine x in terms of y, or y
in terms of x, so that we may find the influence of a variation of one
of the unknown quantities on the differences between computation
and observation when the most probable value of the other unknown
quantity is used. Thus, if it be desired to find x in terms of y, the
most probable value of x will be
[mp~\ \mn\
[mm] [mm]
If we substitute this value of x in the original equations of condition,
the remaining differences between computation and observation will
be expressed in terms of the unknown quantity y, or in the form
A0 = m + n y. (90)
Then, by assigning different values to y, we may find the correspond-
ing residuals, and thus determine to what extent the correction y may
be varied without causing these residuals to surpass the limits of the
probable errors of observation.
In the determination of the orbit of a comet there must be more
or less uncertainty in the value of a, and if y denotes the correction
to be applied to the assumed value of -, we may thus determine the
CL
probable limits within which the true value of the periodic time
must be found. In the case of a comet which is identified, by the
similarity of elements, with one which has previously appeared, if
we compute the system of elements which will best satisfy the series
of observations, the supposition being made that the comet has per-
formed but one revolution around the sun during the intervening
interval, it will be easy to determine whether the observations are
better satisfied by assuming that two or more revolutions have been
completed during this interval. Thus, let T denote the periodic
time assumed, and the relation between T and a is expressed by
k '
in which n denotes the semi-circumference of a circle whose radius
ORBIT OF A COMET. 355
is unity. Let the periodic time corresponding to - -f- y be denoted
m * Q>
by ; then we shall have
1 , 1
y = i*~7
and the equations for the residuals are transformed into the form
A f f\ I f \ f(\-\ \
A0 = (m nj ) -f- n j z*. (91)
If we now assign to 2, successively, the values 1, 2, 3, &c., the re-
siduals thus obtained will indicate the value of z which best satisfies
the series of observations, and hence how many revolutions of the
comet have taken place during the interval denoted by T.
122. In the determination of the orbit of a comet from three ob-
served places, a hypothesis in regard to the semi-transverse axis may
w r ith facility be introduced simultaneously with the computation of
the parabolic elements. The numerical calculation as far as the form-
ation of the equations (52) s will be precisely the same for both the
parabolic and the elliptic or hyperbolic elements. Then in the one
case we find the values of r, r f/ , and x which will satisfy equation
(56) 3 , and in the other case we find those which will satisfy the equa-
tion (65), as already explained. From the results thus obtained, the
two systems of elements will be computed. Let /= -> then in the
case of the system of parabolic elements Ave have/=0, and the com-
parison of the middle place with these and also with the elliptic or
hyperbolic elements will give the value of
in which 6 l denotes the geocentric spherical co-ordinate computed
from the parabolic elements, and 2 that computed from the other
system of elements. Further, let A# denote the difference between
computation and observation for the middle place, and the correction
to be applied to /, in order that the computed and the observed
values of 6 may agree, will be given by
Hence, the two observed spherical co-ordinates for the middle place
will give two equations of condition from which A/ may be found,
356 THEORETICAL ASTRONOMY.
and the corresponding elements will be those which best represent
the observations, assuming the adopted value of M to be correct.
123. The first determination of the approximate elements of the
orbit of a comet is most readily effected by adopting the ecliptic as
the fundamental plane. In the subsequent correction of these ele-
ments, by varying - and M or J, it will often be convenient to use
a
the equator as the fundamental plane, and the first assumption in
regard to M will be made by means of the values of the distances
given by the approximate elements already known. But if it be
desired to compute M directly from three observed places in reference
to the equator, without converting the right ascensions and declina-
tions into longitudes and latitudes, the requisite formula may be
derived by a process entirely analogous to that employed when the
curtate distances refer to the ecliptic. The case may occur in which
only the right ascension for the middle place is given, so that the
corresponding longitude cannot be found. It will then be necessary
to adopt the equator as the fundamental plane in determining a
system of parabolic elements by means of two complete observations
and this incomplete middle place. If we substitute the expressions
for the heliocentric co-ordinates in reference to the equator in the
equations (4) 3 and (5) s , we shall have
n (p cos a R cos D cos A} (p r cos a' R f cos U cos J/)
-f n" (P" sin a" R" cos D" cos A"\
= n (p sin a R cos D sin A) (p f sin a' E' cos D' sin A') (92)
-f n" (p" sin a" R" cos D" sin A"),
Q = n(pttmd R sin D} (p tan S' R' sin Z>')
-f- n" (p" tan a" R" sin IT),
in which p y //, p tf denote the curtate distances with respect to the
equator, A, A f , A" the right ascensions of the sun, and D, Z>', I/'
its declinations. These equations correspond to (6) 3 , and may be
treated in a similar manner.
From the first and second of equations (92) we get
= n (p sin (a' a) RcosD sin (ofA)) -f R cos D' sin (a' 4')
n" (P" sin (a" a') -f R" cosZ>" sin (a' A")),
and hence
=? = ^ >4 na-a
p n" Sin (a" a')
?i ft cos D sin (a' A) R r cosD' sin (a'- A'}-\-ri'R" cos D" sin (a' A")
pri'sm(a" a')
VARIATION OF TWO RADII- VECTORES. 357
This formula, being independent of the decimation d', may be used
to compute M when only the right ascension for the middle place ia
given. For the first assumption in the case of an unknown orbit,
we take
1f sin (a' a)
M=
if t sin(a' / a')'
and, by means of the results obtained from this hypothesis, the com-
plete expression (93) may be computed. By a process identical with-
that employed in deriving the equation (36) 3 , we derive, from (93),
the expression
p" = p--j$~$- ) (94)
j TT f^\tt\i^- \ \R r cos D' sin (a' A!}
-g^F^ ^ ; \r* #/ sin (a" a')
and, putting
,_. n sin (a' a)
n" ' sin (a" a') '
F I ' rr * (?' T ^ cos D' sin (a' A') ^/1__1_\
^ g n ' r" ^ T sin (a' a) p \r' 3 R' 3 f'
we have
(95)
The calculation of the auxiliary quantities in the equations (52) 3
will be effected by means of the formulae (96) 3 , (86), (87), (102) 3 , and
(51) 3 . The heliocentric places for the times t and t" will be given
by (106) 3 and (107) 3 , and from these the elements of the orbit will
be found according to the process already illustrated.
124. The methods already given for the correction of the approxi-
mate elements of the orbit of a heavenly body by means of additional
observations or normal places, are those which will generally be
applied. There are, however, modifications of these which may be
advantageous in rare and special cases, and which will readily suggest
themselves. Thus, if it be desired to correct approximate elements
by varying two radii-vectores r and r", we may assume an approxi-
mate value of each of these, and the three equations (88)j will con-
tain only the three unknown quantities d, b, and /. By elimination,
these unknown quantities may be found, and in like manner the
S58 THEORETICAL ASTRONOMY.
values of J", b"j and I". It will be most convenient to compute
the angles ^ an< ^ Vj anc ^ then find z and z" from
R sin * R" sin V
sm z = - , sm z = - ^ ,
or, putting # 2 = r 2 J? 2 sin 2 ^, and x" 2 = r" 2 -R" 2 sin 2 ^", from
sin 4 " sin V
tan z = - , tan z = - .
x x"
The curtate distances will be given by the equations (3), and the
heliocentric spherical co-ordinates by means of (4), writing r in place
of a. From these u" u may be found, and by means of the values
of r, r n ', and u" u the determination of the elements of the orbit
may be completed. Then, assigning to r an increment dr, we com-
pute a second system of elements, and from r and r" -f- dr ff a third
system. The comparison of these three systems of elements with an
additional or intermediate observed place will furnish the equations
for the determination of the corrections Ar and Ar /r to be applied to
r and r", respectively. The comparison of the middle place may be
made with the observed geocentric spherical co-ordinates directly, or
with the radius-vector and argument of the latitude computed directly
from the observed co-ordinates; and in the same manner any number
of additional observed places may be employed in forming the equa-
tions of condition for the determination of Ar and Ar r/ .
Instead of r and r", we may take the projections of these radii-
vectores on the plane of the ecliptic as the quantities to be corrected.
Let these projected distances of the body from the sun be denoted
by r and r ", respectively ; then, by means of the equations (88)^
we obtain
(96)
from which I may be found ; and in a similar manner we may find
I". If we put
we have
tan(;-,>) = * sin(A - 0) . (97,
X Q
Let 8 denote the angle at the sun between the earth and the place
of the planet or comet projected on the plane of the ecliptic ; then
we shall have
VARIATION OF TWO RADII-VECTORES. 359
=180 + O I,
==: R8in(Z- L 0) (98)
and
p tan ft
tan 6 = , (99)
r o
by means of which the heliocentric latitudes b and b" may be found.
The calculation of the elements and the correction of r and r " are
then effected as in the case of the variation of r and r" .
In the case of parabolic motion, the eccentricity being known, we
may take q and T as the quantities to be corrected. If we assume
approximate values of these elements, r, r', r ff , and v, v f , v" will be
given immediately. Then from r, r r , r' 1 and the observed spherical
co-ordinates of the body we may compute the values of u" u' and
u' u. In the same manner, by means of the observed places, wo
compute the angles u" u' and u f u corresponding to q-\-dq and T y
and to q and T-+- d T } 3q and dT denoting the arbitrary increments
assigned to q and T, respectively. The comparison of the helio-
centric motion, during the intervals t" t 1 and t' i, thus obtained,
in the case of each of the three systems of elements, from the ob-
served geocentric places with the corresponding results given by
w" u' = v" v',
enables us to form the equations by which we may find the cor-
rections A^ and &T to be applied to the assumed values of q and T,
respectively, in order that the values of u" u f and u' u computed
by means of the observed places shall agree with those given by the
true anomalies computed directly from q and T.
360 THEORETICAL ASTRONOMY.
CHAPTER VII.
METHOD OP LEAST SQUARES, THEORY OF THE COMBINATION OP OBSERVATIONS, AND
DETERMINATION OF THE MOST PROBABLE SYSTEM OF ELEMENTS FROM A SERIES
OF OBSERVATIONS.
125. WHEN the elements of the orbit of a heavenly body are known
to such a degree of approximation that the squares and products of
the corrections which should be applied to them may be neglected,
by computing the partial differential coefficients of these elements
with respect to each of the observed spherical co-ordinates, we may
form, by means of the differences between computation and observa-
tion, the equations for the determination of these corrections. Three
complete observations will furnish the six equations required for the
determination of the corrections to be applied to the six elements of
the orbit ; but, if more than three complete places are given, the
number of equations will exceed the number of unknown quantities,
and the problem will be more than determinate. If the observed
places were absolutely exact, the combination of the equations of
condition in any manner whatever would furnish the values of these
corrections, such that each of these equations would be completely
satisfied. The conditions, however, which present themselves in the
actual correction of the elements of the orbit of a heavenly body by
means of given observed places, are entirely different. When the
observations have been corrected for all known instrumental errors,
and when all other known corrections have been duly applied, there
still remain those accidental errors which arise from various causes,
such as the abnormal condition of the atmosphere, the imperfections
of vision, and the imperfections in the performance of the instrument
employed. These accidental and irregular errors of observation cannot
be eliminated from the observed data, and the equations of condition
for the determination of the corrections to be applied to the elements
of an approximate orbit cannot be completely satisfied by any system
of values assigned to the unknown quantities unless the number of
equations is the same as the number of these unknown quantities.
It becomes an important problem, therefore, to determine the par-
ticular combination of these equations of condition, by means of which
METHOD OF LEAST SQUARES. 361
the resulting values of the unknown quantities will be those which,
while they do not completely satisfy the several equations, will afford
the highest degree of probability in favor of their accuracy. It will
be of interest also to determine, as far as it may be possible, the
degree of accuracy which may be attributed to the separate results.
But, in order to simplify the more general problem, in which the
quantities sought are determined indirectly by observation, it will be
expedient to consider first the simpler case, in which a single quantity
is obtained directly by observation.
126. If the accidental errors of observation could be obviated, the
different determinations of a magnitude directly by observation would
be identical ; but since this is impossible when an extreme limit of
precision is sought, we adopt a mean or average value to be derived
from the separate results obtained. The adopted value may or may
not agree with any individual result, since it is only necessary that
the residuals obtained by comparing the adopted value with the
observed values shall be such as to make this adopted value the most
probable value. It is evident, from the very nature of the case, that
we approach here the confines of the unknown, and, before we pro-
ceed further, something additional must be assumed.
However irregular and uncertain the law of the accidental errors
of observation may be, we may at least assume that small errors are
more probable than large errors, and that errors surpassing a certain
limit will not occur. We may also assume that in the case of a large
number of observations, errors in excess will occur as frequently as
errors in defect, so that, in general, positive and negative residuals
of equal absolute value are equally probable. It appears, therefore,
that the relative frequency of the occurrence of an accidental error A
in the observed value will depend on the magnitude of this error,
and may be expressed by (n' x)
Q = (n x) 7 - ^rjr, -- ^ + (ri x) 7-? & \\, , , 4- &c.,
(n x) d (n x) (ri x) d (ri x) '
and the comparison of this with (6) shows that
d log (ri x)
(n x)d(n x) (ri x)d(ri x) "
k being a constant quantity. Hence we derive
d loge ( J) = kA d J,
the integration of which gives
log e c being the constant of integration. From this equation tiiere
results
, ,, ^*A ,
^ ( J) = ce , (8)
in which e is the base of Naperian logarithms. Since 1 O ^ 1 2 S ' ^ 1 2 ^J 4 *'
which converges rapidly when T is small. To find the value of T
which corresponds to the value 0.44311 assigned to the integral, we
compute the value of the series (17) for the values 0.45, 0.47, and
0.49 assigned to T, successively, and from the results thus obtained
it is easily seen that when the sum of the terms of the series is
0.44311, we have
T=Jir = 0.47694,
or
0.47694
r = -j , (18)
which determines the relation between the probable error and the
measure of precision.
The probability that the error of an observation, without regard to
sign, does not exceed nr, is expressed by
, > I <-/ wi/j ( J.t7 )
VTT^O
and this integral, therefore, indicates the ratio of the number of obser-
vations affected with an error which does not exceed nr to the whole
number of observations. Hence, if we assign different values to n,
the integral (19) computed for the several assumed values of
nhr = 0.47694n
will give the relative number of errors of a given magnitude. Thus,
if we put n = \) we obtain
0.2385
-4: JV" eft = 0.264.
from which it appears that in a series of 1000 observations there
ought to be 264 observations in which the error does not exceed \r.
It has been found, in this manner, that in the case of an extended
series of observations the number of errors of a given magnitude
assigned by theory agrees very closely with that actually given by
the series of observations ; and hence we conclude that the error com-
mitted in extending the limits of the summation in the expression (1)
to oo and -f oo, instead of the finite limits which it is presumed
that the actual errors cannot exceed, is very slight, so that the form
368 THEORETICAL ASTRONOMY.
of the function (p (J) which has been derived may be regarded as that
which best satisfies all the conditions of the problem.
129. The relative accuracy of different series of observations may
also be indicated by means of what are called the mean error and the
mean of the errors for each series, the former being the error whose
square is equal to the mean of the squares of all the errors of the
series, and the latter the mean of these errors without reference to
their algebraic sign.
Let denote the mean error ; then, since the number of observa-
tions having the error J is m
V x m
and we have
In the case of a single observation, if P denotes the probability of
the error zero, and P f the probability of the error d, we have
Hence it appears that if h denotes the measure of precision of the
arithmetical mean of m observations, the relation between h and A,
the measure of precision of an observation, is given by
h * = mh*; (27)
and if r is the probable error of the arithmetical mean, and its
mean error, we have, according to the equations (18) and (20),
(28)
These expressions determine the probable and the mean error of the
arithmetical mean of a number of observations when these errors in
the case of a single observation are known.
131. The expressions for the relation between the mean and pro-
bable errors have been derived for the case of a very large number
of observations, a number so great that the error of the arithmetical
mean becomes equal to zero. In the case of a limited number of
observed values of #, the residuals given by comparing the arith-
metical mean with the several observations will not, in general, give
the true errors of the observations ; but the greater the number of
observations, the nearer will these residuals approach the absolute
errors. If J, J r , A n ', &c. are the actual errors of the observations,
and v, v f , v ff , &c. those which result from the most probable value of
x, we shall have, denoting the arithmetical mean by x w and the true
value by X Q -f 3,
A=v S, A' = v' 3, A" = v " d, &c.;
METHOD OF LEAST SQUABES. 371
hence
me 2 = [ A J] = [ W ] + m<* 2 . (29)
This equation will enable us to determine the mean error of an ob-
servation when 3 is given ; but, since this is necessarily unknown,
some assumption in regard to its value must be made. If we assume
it to be equal to the mean error of the arithmetical mean, the re-
maining error will be wholly insensible, and hence the equation (29)
becomes
we 2 = [iw] + me 2 = [vv~\ -f- 3 .
Therefore, we shall have
<
and, according to (21),
r = 0.6745
These equations give the values of the mean and probable errors of
a single .observation in terms of the actual residuals found by com-
paring the arithmetical mean with the several observed values.
The probable and the mean error of the arithmetical mean will be
given by
' ^ (32;
r = 0.6745^ |^j_ .
When the number of observations is very large, the probable error
of an observation and also that of the arithmetical mean may be de-
termined by means of the mean of the errors. If we suppose the
number of positive errors to be the same as the number of negative
errors, the mean of the errors without reference to the algebraic sign
gives
and hence we have, according to (23),
r = 0.8453 &i (33)
m
For the mean error of an observation we have
e = l/ii = 1.2533 4 (34)
372 THEORETICAL ASTRONOMY.
If the number of observations is very great, the results given oy
these equations will agree with those given by (30) and (31); but for
any limited series of observed values, the results obtained by means
of the mean error will afford the greatest accuracy.
132. The relative accuracy of two or more observed values of a
quantity may be expressed by means of what are called their weights.
If the observations are made under precisely similar circumstances,
so that there is no reason for preferring one to the other, they are said
to have the same weight. The weight must therefore depend on the
measure of precision of the observations, and hence on their probable
errors. The unit of the weight is entirely arbitrary, since only the
relative weights are required, and if we denote the weight by p, the
value of p indicates the number of observations of equal accuracy
which must be combined in order that their arithmetical mean may
have the same degree of precision as the observation whose weight is
p. Hence, if the weight of a single observation is 1, the arithmetical
mean of m such observations will have the weight m. Let the pro-
bable error of an observation of the weight unity be denoted by r,
and the probable error of that whose weight is p f by r' then, ac-
cording to the first of equations (28), we shall have
or
For the case of an observation whose weight is p rr and whose pro-
bable error is r", we have
from which it appears that the weights of two observations are to each
other inversely as the squares of their probable or mean errors, and f
according to (18), directly as the squares of their measures of precision.
Let us now consider two values of #, which may be designated by
x r and a?", the mean errors of these values being, respectively, e' and
e": then, if we put
X=x'x"
and suppose that both a;' and x tr have been derived from a large num-
ber m of observations (and the same number in each case), so that the
residuals v 9 v/, v n ', &c. in the case of x f and the residuals v,, v/, v/' t
&c. in the case of x n may be regarded as the actual errors of obser-
METHOD OF LEAST SQUARES.
vatiou, the errors of the value of X t as determined from the several
observations, will be
V V,, V' V,', V" V,", &G.
Let the mean error of X be denoted by E; then we have
mE* =S(v v,y = [w] zh 2 [>,] + [v,v t ] ;
and since the number of observed values is supposed to be so great
that the frequency of negative products w, is the same as that of the
similar positive products, so that [vvj] = 0, this equation gives
or
E 2 = e" -I- e" 2 .
Combining X with a third value x" r whose mean error is e'", the
mean error of x f x" x 1 " will be found in the same manner to be
equal to e /2 + e //2 -j- s ///2 ; and hence we have, for the algebraic sum
of any number of separate values,
E = * + e' 2 + e" 2 + Ac., (35)
and, according to the last of equations (21),
R = ir 2 -hr' 2 + r" 2 + &c., (36)
R being the probable error of the algebraic sum. If the probable
errors of the several values are the same, we have
and the probable error of the sum of m values will be given by
E = rl/m.
Hence tne probable error of the arithmetical mean of m observed
values will be
R r
*= =77
m Vm
which agrees with the first of equations (28).
Let P denote the weight of the sum X, p f the weight of a?', and p /f
that of x" ; then we shall have
p _
"
374 THEORETICAL ASTRONOMY.
from which we get
Since the unit of weight is arbitrary, we may take
and hence we have, for the weight of the algebraic sum of any
number of values,
= ^ == t /1 + r" 2 -f-r'" 2 -{-&c.'
or, whatever may be the unit of weight adopted,
. I * I __ I
p' "" p" ' p'" '
In the case of a series of observed values of a quantity, if we
designate by r f the probable error of a residual found by comparing
the arithmetical mean with an observed value, by r the probable
error of the observation, by X Q the arithmetical mean, and by n any
observed value, the probable error of
n=x +v,
according to (36), will be
.+**=; + ",
r being the probable error of the arithmetical mean. Hence we derive
m
m 1
and if we adopt the value
/ = 0.8453 ^3 ;
m'
the expression for the probable error of an observation becomes
r = 0.8453 M (40)
l/m(m 1)
in which [v] denotes the sum of the residuals regarded as positive,
and m the number of observations.
133. Let n, n f , n", &c. denote the observed values of x, and let p,
p', p ff } &c. be their respective weights ; then, according to the defi-
METHOD OF LEAST SQUARES. 375
nition of the weight, the value n may be regarded as the arithmetical
mean of p observations whose weight is unity, and the same is true
in the case of n f , n f/ , &c. We thus resolve the given values into
p + p f + p" + observations of the weight unity, and the arith-
metical mean of all these gives, for the most probable value of x,
_pn + p'ri + p"n" + Ac. _ [pn\ ,
' ' &c. ' '
The unit of weight being entirely arbitrary, it is evident that the
relation given by this equation is correct as well when the quantities
P) p'y P"> & c - are fractional as when they are whole numbers. The
weight of X Q as determined by (41) is expressed by the sum
and the probable error of x is given by
r ' (42)
when r, denotes the probable error of an observation whose weight
is unity. The value of r, must be found by means of the observa-
tions themselves. Thus, there will be p residuals expressed by
n x Q , p' residuals expressed by n' o? , and similarly in the case of
n", n" r , &c. Hence, according to equation (31), we shall have
r, = 0.6745
in which m denotes the number of values to be combined, or the
number of quantities n, n f , n", &c. For the mean error of # , we
have the equations
(44;
If different determinations of the quantity x are given, for which
the probable errors are r, r', r", &c., the reciprocals of the squares
of these probable errors may be taken as the weights of the respective
values n, n', n rf , &c., and we shall have
!L L ?L- I ^ L
.,2 ~r ^/ 2 ~r v //2 T *
376 THEORETICAL ASTRONOMY.
with the probable error
= ' (46)
The mean errors may be used in these equations instead of the pro-
bable errors.
134. The results thus obtained for the case of the direct observa-
tion of the quantity sought, are applicable to the determination of
the conditions for finding the most probable values of several un-
known quantities when only a certain function of these quantities is
directly observed. In the actual application of the formulae it will
always be possible to reduce the problem to the case in which the
quantity observed is a linear function of the quantities sought. Thus,
let V be the quantity observed, and , y, , &c. the unknown quan-
tities to be determined, so that we have
Let , y w f , &c. be approximate values of these quantities supposed
to be already known by means of previous calculation, and let x, y,
z, &c. denote, respectively, the corrections which must be applied to
these approximate values in order to obtain their true values. Then,
if we suppose that the previous approximation is so close that the
squares and products of the several corrections may be neglected, we
have
_ T , dV , dV . dV
v - v =d! x +^y+~dt z +-->
and thus the equation is reduced to a linear form. Hence, in general,
if we denote by n the difference between the computed and the ob-
served value of the function, and similarly in the case of each obser-
vation employed, the equations to be solved are of the following
form :
ax -f- by "h cz ~f~ du -\- ew -f- / -j~ n 0,
a'x + b'y -f c'z + d'u + e'w +ft + ri = 0, (47)
a " x 4. y'y 4- c " z 4. d"u -f e"w +ft + n" = 0,
&c. &c.
which may be extended so as to include any number of unknown
quantities. If the number of equations is the same as the number
of unknown quantities, the resulting values of these will exactly
satisfy the several equations; but if the number of equations exceeds
the number of unknown quantities, there will not be any system of
METHOD OF LEAST SQUARES. 377
values for these which will reduce the second members absolutely to
zero, and we can only determine the values for which the errors for
the several equations, which may be denoted by v, v', v", &c., will be
those which we may regard as belonging to the most probable values
of the unknown quantities.
Let J, J r , J", &c. be the actual errors of the observed quantities;
then the probability that these occur in the case of the observations
used in forming the equations of condition, will be expressed by
and the most probable values of the unknown quantities will be those
which make P a maximum. The form of the function
shall be a minimum. Hence it appears that when the observations are
equally precise, the most probable values of the unknown quantities
are those which render the sum of the squares of the residuals a
minimum, and that, in general, if each error is multiplied by its
measure of precision, the sum of the squares of the products thus
formed must be a minimum.
If we denote the actual residuals by v, v f , v", &c., and regard the
observations as having the same measure of precision, the condition
that the sum of their squares shall be a minimum gives
dM =0 ^M = 0> ^W^o.&c.,
dx * dy dz
378 THEOEETICAL ASTRONOMY.
or
dv . ,dv' ,,dv" ,
V ~T- + v -J- + *> -T- + 0,
eta ' eta n da? ~
dv . ,dv r dv"
I p **/ . ff UU . -
V ^ + V ^ + V dF +---- = >
&c. &c.
If we differentiate the equations
ax -{-by + cz -f- du -j- ew -(-/ + n = v,
a'z + 6'y -f- c'z -f d'w + e'w -\-ft + ^' =v', (49)
a "a; 4. 5" y _j_ c " z 4_ rf" M + e "w; +f"t + n" = i;' r ,
&c. &c.
with respect to x, y, z, &c., successively, we obtain
dv _ dv' f dv" o
(50)
dv , dv f dv lf
dy dy dy
&c. &c. &c.
Introducing these values into the equations (48), and substituting for
0, v r , v ff , &c. their values given by (49), we get
[aa] x -j- [a&] y -f [ac] z -j- [ad] w + [ae] w + [a/] < + [a^] = 0,
[a6] x + [65] y + [6c] 2 + [ W] w + \be] w + [6/] t + [6/1] = 0,
[ac] a; + [6c] y + [cc] z -f [cd] u + [ce] ti; + [c/] * + [c/i] = 0,
\_ad~] x + [6d] y + [cd] 3 + [dd] u + [de] w + [d/] * + [dw] = 0,
[ae] a; + {be] y + [ce] 2 + \de\ u + [ee] w -f [e/] < + [en] = 0,
[/]* -h P/]y + ['/]* + [d/] + [>/]* + [//] * + [>] - 0,
in which
[aa] = aa -f- a'a' -j- a" a" -f- . . . .
[ac] =ac-f-aV+a'V' + ....
[66] = 66 -f W -f V'b" + ____
&c. &c.
The equations of condition are thus reduced to the same number as
the number of the unknown quantities, and the solution of these
will give the values for which the sum of the squares of the residuals
will be a minimum. These final equations are called normal equations.
When the observations are not equally precise, in accordance with
the condition that AV + h /2 v' 2 -\- h tt2 v m -\- &c. shall be a minimum,
METHOD OF LEAST SQUARES. 379
each equation of condition must be multiplied by the measure of
precision of the observation; or, since the weight is proportional to
the square of the measure of precision, each equation of condition
must be multiplied by the square root of the weight of the observa-
tion, and the several equations of condition, being thus reduced to
the same unit of weight, must be combined as indicated by the equa-
tions (51).
135. It will be observed that the formation of the first normal
equation is eifected by multiplying each equation of condition by
the coefficient of x in that equation and then taking the sum of all
the equations thus formed. The second normal equation is obtained
in the same manner by multiplying by the coefficient of y; and thus
by multiplying by the coefficient of each of the unknown quantities
the several normal equations are formed. These equations will gene-
rally give, by elimination, a system of determinate values of the
unknown quantities #, y, z, &c. But if one of the normal equations
may be derived from one of the others by multiplying it by a con-
stant, or if one of the equations may be derived by a combination of
two or more of the remaining equations, the number of distinct rela-
tions will be less than the number of unknown quantities, and the
problem will thus become indeterminate. In this case an unknown
quantity may be expressed in the form of a linear function of one or
more of the other unknown quantities. Thus, if the number of
independent equations is one less than the number of unknown
quantities, the final expressions for all of these quantities except one,
will be of the form
X = a + Pt, y = a' -f fit, S = a" -f fi't, &C. (53^)
The coefficients a, /9, a/, /9', &c. depend on the known terms and co-
efficients in the normal equations, and if by any means t can be de-
termined independently, the values of x, y, z, &c. become determinate.
It is evident, further, that when two of the normal equations may be
rendered nearly identical by the introduction of a constant factor, the
problem becomes so nearly indeterminate that in the numerical appli-
cation the resulting values of the unknown quantities will be very
uncertain, so that it will be necessary to express them as in the equa-
tions (53).
The indeterrnination in the case of the normal equations results
necessarily from a similarity in the original equations of condition,
and when the problem becomes nearly indeterminate, the identity of
380 THEORETICAL ASTRONOMY.
the equations will be closer in the normal equations than in the equa-
tions of condition from which they are derived. It should be observed,
also, that when we express x } y, z, &c. in terms of t, as in (53), the
normal equation in t, which is the one formed by multiplying by the
coefficient of t in each of the equations of condition, is not required.
136. The elimination in the solution of the equations (51) is most
conveniently effected by the method of substitution. Thus, the first
of these equations gives
lab'] [ae] [ad] [ae] [a/] [an] .
it/ - 7/ ~ "p ^r Z ' ~p - ^ Zv ~" "* r- -. IV ~^~* 1^ ^^ .
[aa] [aa] [aa] \_aa\ \_aa] [aa]
and if we substitute this for x in each of the remaining normal equa-
tions, and put
[66] - |g| [at] = [M.1], [fc] - [gj [ao] = [Jc.1],
(54)
r i r i
[] - gj [ae] = [ce.ll, [ffl - jgj- [/] = [c/1] ;
[dd] - [afl = [ may be found directly, terms of the third order
being neglected.
In the case of the orbits of comets for which e differs but little
from unity, instead of dM we compute by means of the formula
(142) the value of dT, and since we have
ddT _ _!_ dSM
dt n dt
the equation for p becomes
(jxm \2
l-^f) (! + ")'; (186)
and for a we have
Then e, t/, and q will be found by means of the equations
514 THEORETICAL ASTRONOMY.
e sin (Vf v) = a sin v, /5 cos v,,
e cos (v, v) = e Q -f a (cos v, + e ) + /5 sin v,, (188)
P
=IT7
and the time of perihelion passage will be derived from e and v by
means of Table IX. or Table X.
There remain yet to be found the elements and ?'.
Let us now differentiate the equation
regarding the elements as variable, and we get
2rdrl_ _1_ da 2V l
r L dt J " a 2 ' dt + k 2 (1 + m) ' dt
T
dt tf (1 + m) dt
The differential coefficient is here the increment of the accele-
dt
rating force, in the direction of the tangent to the orbit at the given
point, due to the action of the disturbing force; and if we designate
the angle which the tangent makes with the prolongation of the
radius-vector by ^ , we shall have
_ = E cos 4' + S sin 4' v
dt
-Substituting this value in the preceding equation, we obtain
VARIATION OF CONSTANTS. 519
But we have, according to the equations (50) 6 ,
dv
in wnich v denotes the true anomaly in the instantaneous orbit; and
hence there results
e sn
by means of which the variation of a may be found.
If we introduce the mean daily motion ft, we shall have
^r = ^'-^p ( 199 )
and hence
? sin vR + 2- ), (200)
for the determination of dp.
The first of the equations (97) gives
and hence we obtain
d (i/p) =
) It
f (2 cos 2 v cos v cos E} cot y S] (t O-> />
sin v ^ dt
or
dM 1
- . = (( p cot cos vR
d/
-f cos Sfi.
ftjf
1863 Dec.
12.0,
+ 0'
'.01
0"
.00
+ 8".43
+ 0"
.12 -
f 0".0007
5".48
1864 Jan.
21.0,
.04
.01
8
.49
.38
.0040
+ 5 .72
March
1.0,
.24
.03
26
.78
1
.80
.0216
19 .15
April
10.0,
.66
.06
48
.01
3
.88
.0502
37 .11
May
20.0,
1
.35
.08
72
.82
6
.27
.0875
60 .91
June
29.0,
2
.28
.10
100
.83
8
.61
.1299
90 .73
Aug.
8.0,
3
.40
.09
130
.93
10
.65
.1751
125 .79
Sept.
17.0,
4
.63
.07
161
.66
12
.26
.2200
164 .79
Oct.
27.0,
5
.84
.03
191
.48
13
.48
.2624
206 .19
Dec.
6.0,
6
.93
4-0
.03
219
.27
14
.44
.3004
248 .72
1865 Jan.
15.0,
y
.81
-f
.10
244
.24
15
.37
+
.3323
+ 291 .33
Applying the variations of the elements thus obtained to the oscu-
lating elements for 1864 Jan. 1.0, as given in Art. 166, the osculating
elements for the instant 1865 Jan. 15.0 are found to be the following:
Epoch = 1865 Jan. 15.0 Berlin mean time.
M= 99 34' 48".81
TT = 44 13 7 .93
i = 4 36 52. 21
? = 11 15 35 .65
log a = 0.3880283
ft = 928".8897.
1860.0.
NUMERICAL, EXAMPLE. 531
In order to compare the results thus derived with the nerturbations
computed by the other methods which have been given, let us com-
pute the heliocentric longitude and latitude, in the case of the dis-
turbed orbit, for the date 1865 Jan. 15.0, Berlin mean time. Thus,
by means of the new elements, we find
M= 99 34' 48".81, E= 110 5' 14".15,
logr= 0.4162182, v = 120 19 18 .01,
I = 164 37' 59".04, I = 3 5 32 .54,
agreeing completely with the results already obtained by the other
methods. The heliocentric place thus found is referred to the ecliptic
and mean equinox of 1860.0, to which the elements ;r, 2, and i are
referred ; and it may be reduced to any other ecliptic and equinox by
means of the usual formulae. Throughout the calculation of the per-
turbations it will be convenient to adopt a fixed equinox and ecliptic,
the results being subsequently reduced by the application of the cor-
rections for precession and nutation.
In the determination of dM, if we denote by AM the value which
7 Ttr
is obtained when we neglect the last term of the equation for jr> we
shall have
which form is equally convenient in the numerical calculation. Thus,
for 1865 Jan. 15.0, we find
AM = + 234".74,
and from the several values of 1600^- we obtain, for the same date v
by means of the formula for double integration,
Hence we derive
dM = + 234".74 + 56".59 = + 291".33,
agreeing with the result already obtained.
If we compute the variation of the mean anomaly at the epoch, by
means of equation (209), we find, in the case under consideration,
dM = + 165".29,
532 THEORETICAL ASTKONOMY.
and since the place of the body in the case of the instantaneous orbit
is to be computed precisely as if the planet had been moving con-
stantly in that orbit, we have, for 1865 Jan. 15.0,
and hence
dM = 3M + (* 4,) 8?= + 9.91" M.
The error of this result is 0".23, and arises chiefly from the in-
crease of the accidental and unavoidable errors of the numerical cal-
culation by the factor t 1 09 which appears in the last term of the
equation (209). Hence it is evident that it will always be preferable
to compute the variation of the mean anomaly directly; and if the
variation of the mean anomaly at a given epoch be required, it may
easily be found from dM by means of the equation
If the osculating elements of one of the asteroid planets are thus
determined for the date of the opposition of the planet, they will
suffice, without further change, to compute an ephemeris for the brief
period included by the observations in the vicinity of the opposition,
unless the disturbed planet shall be very near to Jupiter, in which
case the perturbations during the period included by the ephemeris
may become sensible. The variation of the geocentric place of the
disturbed body arising from the action of the disturbing forces, may
be obtained by substituting the corresponding variations of the ele-
ments in the differential formulae as derived from the equation (1) 2 ,
whenever the terms of the second order may be neglected. It should
be observed, however, that if we substitute the value of dM directly
in the equations for the variations of the geocentric co-ordinates, the
coefficient of d/j. must be that which depends solely on the variation
of the semi-transverse axis. But when the coefficient of dp. has been
computed so as to involve the effect of this quantity during the in-
terval t t w the value of dM Q must be found from dM and substi-
tuted in the equations.
200. It will be observed that, on account of the divisor e in the
expressions for -> -TT-, and -^p these elements will be subject to large
perturbations whenever e is very small, although the absolute effect
on the heliocentric place of the disturbed body may be small ; and on
VARIATION OF CONSTANTS. 533
account of the divisor sin i in the expression for 5 the variation
of & will be large whenever i is very small. To avoid the difficul-
ties thus encountered, new elements must be introduced. Thus, in
the case of &, let us put
a" = sin i sin & , 0" = sin i cos & ; (224)
then we shall have
da" . di . d
7 = sm & cos i-j- -f- sin t cos & 7 >
at ac at
- cos a cos t- - sm * sin
Introducing the values of -^- and given by the equations (212),
and introducing further the auxiliary constants a, b, A, and B com-
puted by means of the formulae (94) x with respect to the fundamental
plane to which & and i are referred, we obtain
7 tl -j
rZsin a cos (.4 + w).
(225)
6 cos (B + w),
(1 -f m)
by means of which the variations of a." and /9" may be found. If
the integrals are put equal to zero at the beginning of the integration,
the values of da" and dp" will be obtained, so that we shall have
sin i sin & = sin i sin & -f- a",
sin i cos & = sin i cos & -f &&' *
or
sin i sin ( & ) = cos #" sin & <5y?",
sin i cos (& $^ ) = sin i + sin ^ ^a" -J- cos ^ dp", (226)
by means of which i and & & may be found.
In the case of #, let us put
y" = e sin /, C" = e cos /, (227 )
and we have
de . d x
d?' de d x
534 THEORETICAL ASTRONOMY.
Substituting for -=- and -j- the values given by the equations (203)
and (205), and reducing, we obtain
-{- m) *
-f-ersin/j/S),
d" 1 I (228)
p Sin + ^ ^ + ^ + ') CQS
ercos/j/Sl,
by means of which the values of dy" and ^f" may be found. Then
we shall have
e sin x = e sin TT O -f &/',
e cos / = e cos TT O -f- <5C",
or
e sin (/ TT O ) = cos TT O &/' sin TT O (*>
and, also,
dy m' dtf
~dt = 1-fm' ' "3T
m' fr
1 -f- m' ' eft '
If, therefore, from the elements of the orbit of the disturbing planet
we compute the auxiliary constants for the adopted fundamental
plane by means of the equations (94) x or (99) 1? and also V and U r
from
V P '
(e r sin a/ + sin w') =7' sin J7',
l m ' (e r cos / + cos ti') = F cos U',
'
Vp'
the equations (100)! and (49), in connection with (232) and (233),
give
9x = ' r , r' sin a' sin (A* + i*'), (234^
PERTURBATIONS OF COMETS. 539
dy = f r ' sin b f sin ( f -f u') t
m f
Sz = -T r r sm c sm ( (J -f- w') ;
X -J- Wi
*" m ' i' + U'), (234)
"* j,., . , ^ ~, TT f \
If we add the values of dx 9 %, z, 3-j-> J-TT and (5 -yr- to the cor-
w* wt dt
responding co-ordinates and velocities of the comet in reference to
the centre of gravity of the sun, the results will give the co-ordinates
and velocities of the comet in reference to the common centre of
gravity of the sun and disturbing planet, and from these the new
elements of the orbit may be determined as explained in Art. 168.
The time at which the elements of the orbit of the comet may be
referred to the common centre of gravity of the sun and planet, can
be readily estimated in the actual application of the formulae, by
means of the magnitude of the disturbing force. In the case of Mer-
cury as the disturbing planet, this transformation may generally be
effected when the radius-vector of the comet has attained the value
1.5, and in the case of Venus when it has the value 2.5. It should
be remarked, however, that the distance here assigned may be in-
creased or diminished by the relative position of the bodies in their
orbits. The motion relative to the common centre of gravity of
the sun and planet disregarding the perturbations produced by the
other planets, which should be considered separately may then be re-
garded as undisturbed until the comet has again arrived at the point
at which the motion must be referred to the centre of the sun, and at
which the perturbations of this motion by the planet under consider-
ation must be determined. The reduction to the centre of the sun
will be effected by means of the values obtained from (234), when the
second member of each of these equations is taken with a contrary
sign.
204. In the cases in which the motion of the comet will be referred
to the common centre of gravity of the sun and disturbing planet,
the resulting variations of the co-ordinates and velocities will be so
smalJ that their squares and products may be neglected, and, there-
540 THEORETICAL ASTRONOMY.
fore, instead of using the complete formulae in finding the new ele-
ments, it will suffice to employ differential formulae. The formulas
(100), give
dx . . f . dr . . , A dv
-Tr = 8iiia sin ( A -f u) -%- -f r sm a cos (A + u) -JT*
Hj- = sin b sin (B + u) ^ + r sin b cos ( + u) ^-, (235)
(it ut (it
dz . . , dr . f ~ . dv
-jr = sm c sm ( C + w ) -57 + r sm c cos ( (7 -f- w) -TI-
at at at
If we multiply the first of these equations by 3x, the second by dy }
and the third by dz; then multiply the first by d> the second by
7 7 Ctf
^-^Y* and the third by S-j-> and put
at dt
P= sin a sin (A -{-u) 3x -{- sin 6 sin (JB -}- w
-f- sin c sin ( C -}- w) ^2,
Q = sin a cos ( J. + w) 5 -f sin 6 cos (JB -f w
4- sin c cos ((7 -f- u) dz;
P' = sin a sin (4 + i*) *- + sin 6 sin (B + ti) d-- (236)
g' = sin a cos (J. + w) d-^ + sin 6 cos (B + w) J-~
-f- sin c cos ( C + u) djr*
we shall have, observing that - = -^= e sin t? and that -^ = -^r-
-dt -dT--dr -df- r -dt v ^t
From the equations
dr dx
_
~ df ^ dt' ^ de'
PERTUKBATIONS OF COMETS. 541
we get
which by means of (237) become
'+ ! +v (238)
=-= esinvP'H- ^ O^
Xf>
From the equation
V-W-TF *
we get
Substituting the values given by (238), observing also that P=dr,
this becomes
dk dp __ V*r re* sin 2 v esi
T' h ^ ^ J "7T-
and, since
F 8 = -(l + 2ecosv+e),
we obtain
flr-., (239,
by means of which the variation of \/p may be found.
The equation
^. = 2P_ F
a t
gives
from which we derive
542 THEORETICAL ASTRONOMY.
from which the new value of the semi-transverse axis a may be
found. To find dp we have
1 3k
dfj. = Ifiad J- /A -, (241 ")
a k
or
Next, to find 8e, we have, from p = a (1 e 2 ),
**! '<** (243 >
or
jocose smv smv~l/p' f , ^P ,
to = -j P H H r-^ 1 -P H 77- (cos v + cos E}
The equation (12) 2 gives
(2 + e cos i>) to, (245)
v
2
a 2 cos? a'cos'p
and from = 1 + e cos v we get
e sin i; re sm v re sin -y
Substituting this value of dv in (245), and reducing, we find
(24 =
Vp
Jc Jct/p ds
V cos (f> = -7- e sin v, Fsin ^ = - -, V = -rr>
r
\ve
have
EESISTING MEDIUM IN SPACE. 553
Substituting these values of R and 8 in the equation (205), it reduces to
= 2Ky ( - 1 sin v ds.
Now, since
we have
V= -- (I -f- 2e cos v + e,
Vp
ds = Vdt = (l + 2e cos v
and hence
' X \ f* ' V^OOJ
If we suppose the function
2e cos v -f- e*) 3 ,
( i j r 3 (1
the value of which is always positive, to be developed in a series
arranged in reference to the cosines of v and of its multiples, so that
we have
^ ( 7 ) ** C 1 + 2e cos v + e rf = A + ^ cos v + c cos 2v + &c -> ( 267 >
in which A, B, &c. are positive and functions of e, the equation (266)
becomes
2
ed% = -- ( A -j- -B cos v -f- . . . .) sin v dv.
Hence, by integrating, we derive
+ ---- ), (268)
from which it appears that / is subject only to periodic perturbations
on account of the resisting medium.
In a similar manner it may be shown that the second term of the
second member of equation (210) produces only periodic terms in the
value of dM, so that if we seek only the secular perturbations due to
the action of the ethereal fluid, the first and second terms of the
second member of (210) will not be considered, and only the secular
perturbations arising from the variation of // will be required.
Let us next consider the elements a and e. Substituting in the
554 THEORETICAL ASTRONOMY.
equations (198) and (202) the values of E and S given by (265), and
reducing, we get
2o / 1 \ 2 .. | ,
da = r- K
n\
\ r I
(270)
It remains now to make an assumption in regard to the law of the
density of the resisting medium. In the case of Encke's comet it
has been assumed that
and this hypothesis gives results which suffice to represent the obser-
vations at its successive returns to the perihelion. Substituting for F
its value in terms of r and a, the equations (270) thus become
dt r'r a
- OM TT
- LKr \J
a COS COS E
r
dt v*
l 2 *
1 ---
\ r a
by means of which djy. and d -j- = a cos cos v I sin v
in which s = 206264". 8, /* being expressed in seconds of arc. Com-
bining the results thus obtained with the differential coefficients of
the geocentric spherical co-ordinates with respect to r and v, as indi-
cated by the equations (42) 2 , we obtain the required coefficients of x
and y to be introduced into the equations of condition. The solution
of all the equations of condition by the method of least squares will
then furnish the most probable values of y and #, or of the secular
variations of the eccentricity and mean motion, without any assump-
tion being made in rei jrence either t ) the density of the ethereal fluid
or to the modifications of the resistance on account of the changes in
the form and dimensions of the comet, and the results thus derived
may be employed in determining the values of J!f, //, and f
9476
9
99
50
2.15
'5
J 3
4807
40
12
13
14
15
16
17
4 40.06
5 1-85
5 23.28
5 44-33
6 4.95
6 25.14
21.79
21.43
21.05
20.62
20.19
19.72
9-999 9377
9271
9'57
935
8905
8768
106
114
122
I 3
137
144
37
10
20
30
40
50
II 3.28
4-39
5-47
6-54
7-58
8-59
.11
.08
.07
.04
I.OI
1. 00
9-999 4767
4726
4686
4645
4604
45 6 3
41
40
4 1
4 1
4 1
4'
18
19
20
21
22
23
6 44.86
7 4.09
7 22.80
7 40.99
7 58.61
8 15.66
19.23
18.71
18.19
17.62
17.05
16.44
9.999 8624
8472
8314
8149
7977
7799
165
172
I 7 8
185
38
10
20
30
40
50
ii 9.59
10.56
11.51
12.44
13-34
14.22
0.97
0.95
0.93
0.90
0.88
0.86
9-999 45"
4481
4440
4399
4358
43'7
4 1
4 1
4 1
4 1
24
25
26
27
28
8 32.10
8 47.93
9 3- 12
9 17.65
9 3 x -5
I5-83
15.19
14-53
I3-85
9-999 7614
7424
7228
7027
6820
190
196
201
20 7
39
10
20
30
40
ii 15.08
15.92
16.73
17.52
18.29
0.84
0.8 1
0.79
0.77
9-999 4276
4234
4*93
4152
4110
42
42
29
9 44.66
13.16
12.46
6608
212
216
50
19.04
0.75
0.72
4069
42
30
10
20
30
40
50
9 57.12
9 59- 12
10 I. II
3-7
5.02
6.94
2.00
1-99
1.96
1.95
1.92
1.91
9-999 6392
6 355
6319
6282
6245
6208
P
37
37
37
37
40
10
20
30
40
50
ii 19.76
20.46
21.13
21.79
22.42
23.02
0.70
0.67
0.66
0.63
0.60
9.999 4027
3985
3944
3902
3860
42
4 1
42
42
4 1
42
31
10
20
30
40
50
10 8.85
10.73
12.59
14.44
16.26
1 8.06
1.88
1.86
I.g 5
1.82
1.80
1.78
9.999 6171
6096
6059
6021
5984
I
38
3 8
41
10
20
30
40
50
ii 23.61
24.17
24.70
25.22
25.71
26.18
0.56
o-53
0.52
0.49
0.47
0.44
9-999 3777
3735
3 6 93
3651
3609
3567
42
42
42
42
42
42
32
10
20
30
40
50
10 19.84
21. 60
23-34
25.05
26.75
28.43
1.76
1.74
1.71
1.70
1.68
1.65
9-999 5946
5908
5870
5832
5794
5755
38
38
38
38
39
3 8
42
10
20
30
40
50
ii 26.62
27.04
27.44
27.82
28.17
28.50
0.42
0.40
0.38
o-35
o-33
0.30
9-999 3525
3483
3399
3357
33'5
42
42
42
42
42
33
10
20
cO
40
10 30.08
33-32
34-9J
36.48
1.63
1.61
i-59
i-57
9-999 5717
5678
5640
5601
39
38
39
39
43
10
20
30
40
ii 28.80
29.08
29-34
29.58
29.79
0.28
0.26
0.24
0.21
9-999 3273
3 270
3188
3H6
3*4
43
42
42
42
42
50
38.03
"55
1.52
5523
39
39
50
29.98
o.i 6
3062
T-*
43
34
10
20
30
40
50
10 39-55
41.06
42.54
44.00
JIU
1.51
1.48
1.46
i-44
1.42
I. -JO
9-999 5484
5445
5406
53 6 7
5327
5288
39
39
39
40
39
4
44
10
10
SO
40
50
ii 30.14
30.29
30.41
30.50
3-57
30.62
0.15
0.12
0.09
0.07
0.05
0.03
9-999 3 OI 9
2977
2935
2892
2850
2808
42
42
43
42
42
42
35
10 48.25
J 7
9-999 5248
45
ii 30.65
9.999 2766
661
TABLE I, Angle of the Vertical and Logarithm of the Earth's Eadius,
= Geocentric Latitude.
P = Earth's Eadius.
0-0'
Diff.
logp
Diff.
0-0'
Diff.
logp
Diff.
O /
45
10
20
30
40
50
II 30.65
30.65
30.63,
30.58
30.42
0.00
0.02
O.O5
0.07
0.0 9
Oil
9.999 2766
2723
2681
2639
2596
43
42
42
43
42
O /
55
10
20
30
40
50
10 49.74
48.36
46.97
45-55
2%
II
1.38
1.39
1.42
1.44
1.46
9-999 *75
0235
0195
0153
OIK
0076
40
40
40
39
40
42
1.49
39
46
10
20
30
40
50
II 30.31
30.17
30.01
29.82
29.61
29.38
0.14
0.16
0.19
0.21
0.23
O.26
9.999 2512
2470
2427
2385
2343
2300
42
43
42
42
43
42
56
10
20
30
40
50
10 41.16
39-65
38.13
36.5!
35-oJ
33-4 1
1.51
1.52
'I 7
1. 60
1.61
9-999 37
9.998 9998
9958
9919
9880
9841
39
40
39
39
39
39
4T
10
20
30
II 29.12
28.85
28.54
28.22
0.27
0.31
0.32
9.999 225*
2216
2174
2132
42
42
42
57
10
20
30
10 31.80
30.16
28.50
26.83
1.64
1.66
1.67
9.998 9802
9764
9725
9686
38
39
39
40
50
27.87
27.50
-35
0-37
0.40
2089
2047
43
42
42
40
50
23.40
1.70
1.73
1.74
9648
9610
39
48
10
20
II 27.10
26.69
26.24
0.41
0.45
9.999 2005
1963
1921
42
42
58
10
20
10 21.66
19.90
18.11
1.76
9.998 9571
9533
9495
1
30
40
50
25.78
25.29
24.78
0.40
0.49
0.51
0.54
1879
1837
'795
42
42
42
42
30
40
50
16.31
14.48
12.63
I'll
1.85
1.86
9457
9419
9382
|
49
10
20
30
40
50
ii 24.24
23.69
23.11
22.50
21.87
21.22
0.55
0.58
0.61
0.63
0.65
0.67
9-999 753
1711
1669
1627
1586
1544
42
42
42
42
42
59
10
20
30
40
50
10 10.77
8.88
6.97
5.04
3-o8
IO I. II
1.89
1.91
1.93
1.96
1.97
1.99
9.998 9344
9307
9269
9232
9158
11
37
37
37
37
50
10
20
30
40
50
II 20.55
19.85
18.39
17.63
16.84
0.70
0.72
0.74
0.76
0.79
0.82
9.999 1502
1460
1419
'377
^335
1294
42
42
42
4'
42
60
61
62
63
64
65
9 59-i*
9 46.74
9 33-65
9 19.85
9 5-36
8 50.21
12.38
13.09
13.80
14.49
15.15
15.81
9.998 9121
8902
8688
8479
8275
8077
219
214
209
204
198
51
10
20
30
40
50
II l6.02
15.19
14-33
13-45
12.55
11.62
0-.83
0.86
0.90
-93
0.95
9.999 1252
I2II
II 7
1128
1087
1046
42
4 1
4 1
66
67
68
69
70
71
8 34-40
8 17-97
8 0.92
7 43- 2 9
7 25.08
7 6.33
16.43
17.05
17.63
18.21
18.75
19.27
9.998 7884
7697
75'7
734*
7174
7013
187
1 80
III
161
J 54
52
10
20
30
II 10.67
9.70
8.71
7.69
0.97
0.99
.02
9-999 1005
0963
0922
O88l
42
72
73
74
75
6 47.06
6 27.28
6 7.03
5 46-33
19.78
20.25
20.70
9.998 6859
6713
6573
6441
146
140
132
40
50
6.66
5.60
3
.09
0840
0800
40
4 1
76
77
5 3-67
21.13
2I -53
21.90
6317
6201
Til
1 08
53
10
20
30
40
50
ii 4.51
3.40
2.27
II 1. 12
10 59.94
58.74
.11
13
.1C
.18
.20
.22
9-999 759
0718
0677
0637
0596
0556
40
4i
40
4 1
78
79
80
81
82
83
4 41-77
4 19-53
3 56.96
3 34-io
3 10.98
2 47.63
22.24
22.57
22.86
23.12
2 3-35
23.56
9.998 6093
5993
5901
5818
5676
100
ii
57
54
10
20
30
40
50
10 57.52
56.28
55.02
53-73
52.42
51.09
.29
33
35
9.999 0515
0475
435
395
355
03*5
40
40
4
40
40
4
84
85
86
87
88
89
2 24.07
2 0.33
I 36.44
I 12.43
o 4.8.3!
o 24.18
23:89
24.01
24.09
24.16
24.18
9.998 5619
5570
553
5498
5476
5463
49
40
3*
22
5
55
10 49.74
J J
9.999 0275
90
O O.OO
9-998 5458
662
TABLE II.
For converting intervals of Mean Solar Time into equivalent intervals of Sidereal
Hours.
Minutes.
Seconds.
Decimals.
Mean T.
Sidereal Time.
Mean T.
Sidereal Time.
lean T.
Sidereal Time.
Mean T.
Sidereal Time.
h
A m s
m
m s
g
1
s
s
I
I o 9.856
I
I 0.164
I
1.003
0.02
0.02O
2
2 19.713
2
2 0.329
2
2.005
0.04
0.04.0
3
3 o 29.569
3
3 0-493
3
3.008
0.06
O.OOO
4
4 o 39.426
4
4 0.657
4
4.011
0.08
0.080
c
5 o 49.282
5 0.821
5
5.014.
0.10
O. IOO
i
6 o 59-139
6
6 0.986
6
6.016
0.12
O.I 20
1
7 8.995
8 18.852
I
7 -15
8 -3H
i
7.019
8.022
O.IJ
o.i 6
0.140
O.I 60
9
9 28.708
9
9 .478
9
9.025
0.18
0.180
10
10 38.565
10
10 643
10
10.027
0.2O
0.201
ii
n 48.421
ii
ii .807
ii
I 1.030
0.22
0.221
12
12 58.278
12
12 .971
12
12.033
0.24
0.241
'3
4
13 2 8.134
14 2 17.991
14
13 2.136
I 4 2.300
13
4
I3-036
14.038
0.26
0.28
0.26l
0.28l
i5
15 2 27.847
'5
'5 2.464
15.041
0.3.0
0.301
16
16 2 37.704
if
16 2.628
1 6
16.044
O.J2
0.321
18
17 2 47.560
1 8 2 5,7.416
18
17 2.793
18 2.957
ii
17.047
18.049
0-34
0.36
0.341
0.361
19
19 3 7.273
19
19 3.121
19
19.052
0.38
0.381
20
20 3 17.129
20
20 3.285
20
20.055
0.40
0.401
21
21 3 26.986
21
21 3.450
21
21.057
0.42
0.421
22
22 3 36.842
22
22 3.614
22
22.000
o-44
0.441
23
24
23 3 46.699
*4 3 56.555
2 4
24 3-943
23
24
24.066
. f O
0.46
o. 4 S
0.461
0.481
11
25 4.107
26 4-271
11
25 ooo
26.071
0.50
0.52
0.501
0.521
oJ
27
27 4-435
27
27.074
o-54
0.541
Q C
28
28 4.600
28
28.077
0.56
0.562
EH 8
2 9
30
29 4-764
30 4.928
2 9
3
29.079
30.082
0.58
0.60
0.582
0.602
g i
3 1
31 5.092
3 1
31.085
0.62
O.622
rQ fl
32
S 2 5-257
32.088
0.64
0.642
OQ bfi
33
33 5-421
33
33.090
0.66
O.662
2 .5 o5
34
34 5-585
34
34-093
0.68
0.682
8 15
35
35 5-750
35
35.096
0.70
0.702
S fj
S 6
3 6 5-9H
3 6
36.099
0.72
0.722
ill
37
37 6.078
37
37.101
0-74
0.742
UU 4) 2
PJ rfi H
38 '
38 6.242
38
38.104
0.76
0.762
* 1 s
39
39 6.407
39
39.107
0.78
0.782
40
40 6.571
40
40.110
0.80
0.802
H a "So
4 1 6-735
41.112
0.82
0.822
*o .5 a)
42
42 6.899
42
42.115
0.84
0.842
Ii!
43
44
43 7-064
44 7.228
43
44
43.118
44.120
086
0.88
0.862
0.882
ill
O i, c3
46
45 7-392
46 7-557
Jl
45.123
46.126
0.90
0.92
0.902
0.923
.
47
47 7-721
47
47.129
0.94
0.943
Ill
48
48 7.885
48
48.131
0.96
0.963
5-1 O
O *J3 ^
49
49 8-049
49
49.134
0.98
0.983
5
50 8.214
5
50.137
1. 00
1.003
3 * i
51 8.378
5 1
51.140
<3
5*
52 8.542
52
52.142
j S
53
53 8.707
53
53-145
.s
54
54 8.871
54
54.148
<2 ^
II
55,
57
55 9-035
5 6 9 '99
57 9-364
I?
57
55-I5I
56.153
57.156
e
58
58 9.528
58
58.159
M
I 9
60
59 9.692
60 9.856
59
60
59.162
60.164
563
TABLE III.
For converting intervals of Sidereal Time into equivalent intervals of Mean Solar Time,
Hours.
Minutes.
Seconds.
Decimals.
Sid. T.
Mean Time.
Sid. T.
Mean Time.
Sid. T.
Mean Time.
Sid. T.
Mean Time.
h
Am
m
m 5
s
8
s
s
I
o 59 50.170
I
o 59.836
I
0.997
O.02
0.020
a
i 59 40.341
a
I 59.672
2
1.995
0.04
0.040
3
a 59 30.511
3
a 59.509
3
2.992
O.O6
0.060
4
3 59 20.682
4
3 59-345
4
3-989
0.08
0.080
5
4 59 IO - 8 5 2
5
4 59.181
5
4.986
O.IO
O.I 00
6
5 59 i- oa 3
6
5 59017
6
5.984
0.12
0.120
7
6 58 51.193
7
6 58.853
7
6.981
0.14
0.140
8
7 58 4i-3 6 3
8
7 58.689
8
7-97 8
0.16
0.160
9
8 58 31.534
9
8 58.526
9
8.975
0.18
0.180
10
9 58 21.704
10
9 58-362
10
9-973
0.20
0.199
ii
10 58 11.875
ii
10 58.198
ii
30.970
0.22
0.219
la
ji 5$ 2.045
la
ii 58.034
12
11.967
0.24
0.239
J 3
la 57 52.216
'3
12 57.870
13
12.964
0.26
0.259
H
13 57 42.386
4
13 57.706
H
13.962
0.28
0.279
II
H 57 3 2 -557
15 57 22.727
11
*4 57-543
15 57-379
!i
14-959
I5-956
0.30
0.32
0.299
0.319
17
16 57 12.897
17
16 57.215
17
16.954
0-34
o-339
18
17 57 3.068
18
17 57.051
18
17.953
0.36
o-359
19
18 56 53.238
'9
18 56.887
19
18.948
0.38
0:379
ao
19 5 6 43-49
20
19 56.723
20
19.945
0.40
o-399
ai
20 56 33.579
21
20 56.560
21
20.943
0.42
0.419
aa
21 56 23.750
22
21 56.396
22
21.940
o-44
0.439
a 3
22 56 13.920
a 3
22 56.232
2 3
22.937
0.46
0,459
24
43 5 6 4-09 1
24
23 56.068
2 4
a 3-934
0.48
0-479
a 5
24 55.904
2 5
24.932
0.50
0.499
26
25 55.740
26
25.929
0.52
0.519
C)
a 7
26 55.577
27
26.926
0.54
0.539
a
a8
27 55-4I3
28
27.924
0.56
0.558
EH
a 9
28 55.249
29
28.921
0.58
0-578
83 73
30
2 9 55.085
3
29.918
0.60
0-598
'o g
3 1
30 54.921
3 s
30.915
0.62
0.618
rj
3 a
3 1 54-758
32
3 I -9 I 3
0.64
0.638
o3 "OD
33
32 54.594
3.3
32.910
0.66
0.658
S C? S
34
33 54-43
34
33.907
0.68
0.678
Vn '&
3S
34 54.266
35
34.904
0.70
0.698
C 8 73
36
35 54.102
36
35-9 oa
0.72
0.718
i II
H
36 53.938
37 53-775
P
36.899
37.896
0.74
0.76
0-738
0.758
coirections to be applied to the corresponding elements are , y^,
and T7y^7- In the same manner, we may adopt as the unknown
quantity, instead of the actual variation of any one of the elements
of the orbit, n times that variation, in which case its coefficient in
the equations must be divided by n.
The value of Aa, derived by taking the difference between the
computed and the observed place, is affected by the uncertainty
necessarily incident to the determination of a by observation. The
unavoidable error of observation being supposed the same in the case
of a as in the case of d, when expressed in parts of the same unit,
it is evident that an error of a given magnitude will produce a
greater apparent error in a than in o, since in the case of a it is
measured on a small circle, of which the radius is cos d and hence,
in order that the difference between computation and observation in
a and d may have the same influence in the determination of the
corrections to be applied to the elements, we introduce COS^AOC
instead of Aa. The same principle, is applied in the case of the
longitude and of all corresponding spherical co-ordinates.
DIFFERENTIAL FORMULAE. 143
52. The formulae already given will determine also the variations
of the geocentric longitude and latitude corresponding to small in-
crements assigned to the elements of the orbit of a heavenly body.
In this case we put e = 0, and compute the values of A, I>, sin a,
and sin b by means of the equations (94) r We have also (7=0,
sin c = sin i, and, in place of a and d, respectively, we write ^ and /?.
But when the elements are referred to the same fundamental plane
as the geocentric places of the body, the formulae which depend on
the position of the plane of the orbit may be put in a form which is
more convenient for numerical application.
If we differentiate the equations
x f = r cos u cos ft r sin u sin ft cos i t
y' = r cos u sin ft -f- r sin u cos ft cos *,
z' = r sin u sin i,
we obtain
dx r = dr r (sin u cos ft -f- cos u sin ft cos i) du
r (cos u sin ft -f sin u cos ft cos i) dft -f- r sin u sin ft sin i cfo',
dtf = dr r (sin u sin & cos u cos & cos i) du
-f r (cos tt cos ft sin w sin ft cos i) c?ft r sin w cos ft sin i eft, (46)
dz 1 = - dr -{ r cos u sin idu -\- r sin it cos i eft,
r
in which #', /', z' are the heliocentric co-ordinates of the body in
reference to the ecliptic, the positive axis of x being directed to the
vernal equinox. Let us now suppose the place of the body to be
referred to a system of co-ordinates in which the ecliptic remains as
the plane of xy, but in which the positive axis of x is directed to the
point whose longitude is ft ; then we shall have
dx = dx' cos & + dtf sin & ,
dy = dx' sin & + dtf cos ft ,
dz = dz',
and the preceding equations give
dx = - dr r sin u du r sin u cos i d& ,
T
dy = - dr + r cos u cos i du -f- r cos u d& r sin it sin i di, (47)
dz = _ dr -f r cos u sin idu ~{-r sin ?f cos i eft.
144 THEORETICAL ASTRONOMY.
This transformation, it will be observed, is equivalent to diminishing
the longitudes in the equations (46) by the angle & through which
the axis of x has been moved.
Let X n Y,, Z, denote the heliocentric co-ordinates of the earth
referred to the same system of co-ordinates, and we have
x -f- X, = A cos/3 cos 0* Q),
V + F,= Jcos08in(J ),
z -f Z, = A sin /?,
in which ^ is the geocentric longitude and /9 the geocentric latitude.
In differentiating these equations so as to find the relation between
the variations of the heliocentric co-ordinates and the geocentric lon-
gitude and latitude, we must regard Q, as constant, since it indicates
here the position of the axis of x in reference to the vernal equinox,
and this position is supposed to be fixed. Therefore, we shall have
cos sn
dy= cos sin (A &)d/f A sin /5 sin (A )d-|- A cos/? cos (A
cfe =sin /? dJ H- J cos y9 d/9,
from which, by elimination, we find
sin/? cos (A a) , sin /9 sin (A q) cos /?
-- _
These equations give
COS p = ==
r
= -- -. - . j
dx A dx
Q ft
cosj? -j- = 0, -,-
C?3 ^2
If we introduce the distance o between the ascending node and the
place of the perihelion as one of the elements of the orbit, we have
du = dv -{- d<0,
and the equations (47) give
dx x dy y dz z . . .
~-=- = cosu, /- = -== smu cos i, = =- = smw smt;
dr r dr r dr r
dx dx dy dy . dz dz
= = -_ = rsmu, . y = -. y - =r cos u cost, = = -=
av aw dv dw dv dot
DIFFERENTIAL FORMULA. 145
^-=rcos, -*l=0j (49)
dx A dy . . c?z
7 . = 0. =r- = ?* sin w sin ^. ^r- = r sin u cos i.
at at di
If we introduce TT, the longitude of the perihelion, we have
du = dv -{- dit a* & ,
and hence the expressions for the partial differential coefficients of
the heliocentric co-ordinates with respect to it and & become
dx dy dz
r= r sm w, 7 r cos u cos ^, , = r cos w sin i ;
f f t ()
dx _ . , , . aw . . t . az
- = 2r sm w sm 2 Jt, ~- = 2r cos w sm 2 J*, -7 = r cos w sin i.
When the direct inclination exceeds 90 and the motion is regarded
as being retrograde, we find, by making the necessary distinctions in
regard to the algebraic signs in the general equations,
dx dy dz
-rr = 0, -T- = f sin u sin i, j-r = r sin u cos i ; (51)
di di di
dx dx dx dy ,
and the expressions for -7-, -j-, -TQ-> -f-, &c. are derived directly
from (49) by writing 180 i in place of i. If we introduce the
longitude of the perihelion, we have, in this case,
du = dv oV + d Q ,
and hence
dx dy dz . .
/- = r cos i* cos i, -j = rcosttsmi;
(52)
-r~- = 2r sin w sin 2 i, -~- = 2r cos tt sin 2 Ji, -^ - = r cos w sin i.
But, to prevent confusion and the necessity of using so many for-
mulae, it is best to regard i as admitting any value from to 180,
and to transform the elements which are given with the distinction
of retrograde motion into those of the general case by taking
180 i instead of i, and 2& n instead of TT, the other elements
remaining the same in both cases.
53. The equations already derived enable us to form those for the
differential coefficients of ^ and ft with respect to r, v, & , i, and co or
T, by writing successively X and ft in place of 0, and &, i, &c. in
10
146 THEORETICAL, ASTRONOMY.
place of it in equation (2). The expressions for the differential coeffi-
cients of r and v, with respect to the elements which determine the
form of the orbit and the position of the body in its orbit, being
independent of the position of the plane of the orbit, are the same as
those already given; and hence, according to (42) and (43), we may
derive the values of the partial differential coefficients of A and /9
with respect to these elements. The numerical application, however,
is facilitated by the introduction of certain auxiliary quantities.
Thus, if we substitute the values given by (48) and (49) in the
equations
. c?A ctt d x . Q ctt dy
cos /9 -j- = cos /? -=- -j f- cos /5- -/-,
dv dx dv l dy dv
&L&L ^L\~ ty-+dP-
dv dx dv dy dv dz dv'
and put
cos i cos (A & ) = A Q sin A,
sin (A & ) = A Q cos A,
sin i = n sin N,
sin (A 2 ) cos i = n cos N,
in which A Q and n are always positive, they become
dv dw A
-~ = -== (sin /9 cos (A & ) sin u -f- n cos u sin (N + ft ).
Let us also put
n sin ( N + ft = -B g i n -^>
sin /5 cos (A & ) = J? cos jB,
and we have
cos /? -j = cos /? -j = ^4 sin (A -\- u -,
* ,! CM)
7 TO
The expressions for cos ft -j- and -j- give, by means of the same
auxiliary quantities,
cos -= = ^ cos (J. + u),
dr (56)
f = -^eos(* + ).
In the same manner, if we put
DIFFERENTIAL FORMULA. 147
cos (/I )= C sin C,
cos i sin (A & ) = (7, cos (7;
1 67)
cos i = D sin Z>
sin f A &) sin i = D cos D;
we obtain
-^- = = ^ (7 sin
A r
cos /9 JT- = -j sin i sin w cos (A & ),
-2(L = -D fl sin w sin (D + ).
Cd T
If we substitute the expressions (55) and (56) in the equations
fl . . fl
cos /9 3 = cos /? -=- 3 -- ^- cos /? -j- 3
d? dr dy dv d'= cos C' cos H' cos <5" cos a"-\- cos C' sin JET' cos 5" sin a"-}- sin C' sin <5",
or
cos ?/ = cos ' cos M gi>H
fe
Ret
Direct.
Retrog
O
vo
Tj- <* ON VO d
O co co TJ- vo ONVO O
r^OooOd d vo M o t-^
vo oo d rj- M ^t-ddt^vo
d vo d d ON TJ- cooo vo vo
oo covo co d vo O f-- w vo
O t^- t*^ vo d M co o w ^J*
o" c?\ o\ o\ o" o oV o" o* o\
i O VO oo ON r-N
oo ^~ cooo oo T^ vo ^ vo
M t^% M ^- M vo t-^ co O
vo M Ooo co oo ^-<4-dvo
ONTJ-VOCOM M COTJ-OOO
t^vo i-i d co ONOO O O d
OvONONOjNVO ONOJNO\O\O
ON ON ON ON ON ON ON ON ON O
vo M T- - TJ- f^ ON OOvOMwd VO ONVO OO OO OO VO M ONOO
TJ-rj-Tj-coTi- cp co cp cp T}- r p rr| '*"'^"' : J" rj- co cp cp co c P T ^" T }" T i" r P
66666 66666 66666 66666 66666
ON f^ r^- co
ON CO vo -rt- O OOOMVOOO
covo OMd oo *!- d vo ON
d rj-oo d ON vo oo M >H ON
rJ-Ocor^co co d r~- ON r>
d t^ *}- O t- ON M vo d t-~
^J" cp TJ- TJ- cp r J" T *" T t" r t" T l"
66666 66666
VO
oo vo vo O
M VD OO OO d
T}- ONVO CO vo
d O O coo
q oo q vo vo
M ONTJ- "
vooooooo vodMTl-
d vo ON ONOO vo t^-vo d O
d O O oo M co T$- ONOO vo
CO OO CO CO II VO CO r*.* O VO
Ovococoro M oo oo d ^h
co oo oot~ vor~>oooooo rj-oo t^
covo oo t^. r^ cooo r>--i-d vorj-d^j-co OOONONCOM
O VO ONOO OO rl- COOO COCO ^-OdTj-Tj- OVOVOVOTJ-
vo vo ** t^ vo vo M o vo d oo ON d co t*^ vo t^*co vooo
r^* OMd vo vo vo M f^oo vo ON M co d oo vovo oo t^
M l-^. vo
ON N vo d -4-
COVO OO VO CO
OsONO ONVO
T$- T}-
c> ON
CO VO
vo vo d
ON ONOO OO OO
ONd O O oo
O w d COrJ-
t^ ONO ON ON
CO rj- COVO CO
CO 11 CO CO VO
ON f~-oo vo ON
ONVO VOONVO codMVoO
MCOVOVOt^ COVOCOQOO
OO ON t^OO O> VO OO t^-OO VO
OO VOOO COM OOCNjCOrl-Tl- MONM
d t^ vovo 6 4- O vo d vo
M d CO CO
CO IO VO
OO
O d d d co 4-oo oo d 06
M d COM ^4- M VOVO
VO CO CO CO I
COOO Tj- M
Q
co rj- vo O M co ON t-^ vo o d Tj- ON O
i O O ON O O co d
d d ONOO vo
vo T}-VO M cp
ON l^-vo d vo
t--OO CO CO O\
^cT^d 5^"
q q > vo TJ- vo vo vo co
t^ co O I s * "*J"
co d co ^
<1-00 vo ONM
CO vo O M O
VOVOM
d vo oo M co cooo vo d oo
t^ vovo d OO t^ vooo vo w
vocodM Mvorhco
I -4- t~~ <* ON ONVO t~- ON t^ d
cp vo Cfr d d
ON M ON O ON
ON vo ON vooo
. vo O d vo r^oo M O ON M oo vooo M
O vo O w vo
vo d vo M O
O vo r^vo <*
00 vo VO CO w
d
VO CO rj-VO OO
rl- ON rt- vo O
M CO CO
M O ONTj-VO
co d oo vo
d M co
tJ-OO ONVO w
00 VO M 00 00
M cod -*
i- q q oo vo
M VO VOVO 6
co d M
vocoti-coO T}-ONVOI
d w co vo co
vo rf- vo t~~- ON
rj- M VO VO Tj-
00 vo M vd CO
OVOVOOOON vOMt^dO
to vo ONOO ON
OO d rj- ON vo
VO M N M VO
W TJ- CO t^ CO
H
vo O vo O oo
M VO n CO ON
d
t^OO coO d
d vod rj-0
tO M IH CO
vo t^ d vo M
coco d vo co
d M
CO r}- d O M
ON ON T*- voO
M CO M d
t^ O ONOO vo
VOO*NOO vod
co TJ- d co
_ t^ M ONVO Tj-
O m M d M co
vo q vo t^ q
6 6* co 4- d
H VO 11
MONCpclvo vOMQcpd MO_NCO t^oo vo cooo cp O O O
vo ON M ON M vo 4" l"^ co 4" vovo vo O M vo M r^ d vo d 4"
llCOrl- CO T}-dCO rj-
rh d T}- vo vo O t r}- vo co
= 1 e-*
il -sf.pl HI'L
1 3 m'j 0^ "S'djajJ
ifq PL|HS^S H^OHHPn
lMl.2 Ii
II Is
3 S* tn
P- S ^> 3
3 c o .-
Ol^0004O r*
646
J. xi JJ J_LU WiYV i
Elements of the Orbits of the Minor Planets.
i 1: 1 1 I aa a a d . a a a I *i . d a
-g^SS-g -goc-g-g a ri S $ i! j 4 I IB i 1 Ss" c -3 "3 g.'o
TO O C 5 OT O O X W O M Jg TO 0> O 3 TO . * 0. O ( - g
2-52 2 Softs 2 So:222? 25213 2:2:2 S ** Seg-gg SS)5|2
00000 0000
nJOi-^Oi
CO T-l
,d
t- CO iO OS OS T*
iii-31 SSJrS* 9fttt
O O i-s fM S S ^ <; S >-s - d vo
f^oo t^ O
OQ M t-^ N
U~,
O
VO CO T^-
vo O oo r- d
<* o vo vo vo
COTJ-O M t-.
d r-- co rj- T|-
rt-00 rt-00 CO
TJ- CO CO CO 5J-
66666
vo vovO ON O
SON N OO CO
vo ONOO d
r}-vo d ONOO d M vo d
ON M M vo vo M ON ONOO ON
dd^J--^-ON ONTj-ddO
<4- ON M co t~- Tt- vo ONOO Tj-
M ON co cooo t~- ON t^ d CO
; . , , T *" T *"i" r *' t r T *"i" r r' t *" i o
ddddd ddddd ddddd
o o
M t^OO ^ ^J-
vo d O O oo
T}" M M O M
co rj- rj- co r>
r4 oo T*- t M
<* "*
d c d d d
VOMCOVOOOVO ^ONON
vo oo O oo d oo f-- vovo d oo t^ vooo co
CO <$- O vo CO vo vovo ON d COOO ON O **
ft ON 4 vo vo vo ONO co vo co vo ON d d **
VOVOdOO-^- ONONONCOTj- VOTJ-MVOON
ONVO dOON ONO4-MO covot^-voci
t^ d oo t^- co vo cooo xj- ON oo d TJ- vo d
r-^oo t^-r^O t^ONOONt^. oo t-~vo vo oo
TJ- t^ M VO VO Tj-VO CO O M
cooo O d M cooo co f>- t^
M o oo co d covo vo ON T!-
O O d O O
VOMTj-QVO MON OO
MONQVOt^ MOCOQVO
oovodOr^ dONdcoo
O 00 co cooo ON Tj- covo ON
i d ON
.VOCOONS. TJ-COONONVO oo^J-voOvb d -"t-vo ON co
ONVO OOt>t^ OOVOt^t^ON vovo ONOO vo OO ON f-vo OO
ON cooo
ON ON vo
O r^oo O MM
-)vo d rj-'
M r-~ d vo r--vo ooONd t-^t^-ONcovo
VO M c< COOO
. O oo vo d vo co ONOO tj- ON f^ ^- covo
ONVO OO co vo
to co cp ON rt-
ON 4 VOOO VO VO 00 vo 00 I
VO Tj- VO CO Tj- VO <$- CO
M Tj-M d VO
ONO OOO
00 CO CO M 00 MO vpvo O
4-dvovo4 - dcodt~-vo
MMd VO MCO^t-VOd
M vo d 'T rj- co TJ- TJ-
ONVO dVOVO VOVpCpOOO MONONdd
d oo co cooo vo r-'-vo d ON vo vo 4 ON 6
rf- VO ^4 M d M M VO M CO CO M CO
^ d
Tj-
00 COVO O
VOOO CO COVO d VOVO COd ONl*VOMt^ OOVO VOOO CO OOdvOMCO COVO t^OO M
847
1
jj -^ ** a *>
9
IUJI
O c- 3 o
a s i S is 9
I/ate of
)iscovery.
CO OS i OS C0 01
ec
being expressed in parts of the equatorial radius as the unit. These quantities are required in the determination of the parallax of a heavenly body. The formula for the parallax in right ascension and in declination are given in Art. 61. TABLE II. gives the intervals of sidereal time corresponding to given intervals of mean time. It is required for the conversion of mean solar into sidereal time. TABLE III. gives the intervals of mean time corresponding to given intervals of sidereal time. It is required for the conversion t)f sidereal into mean solar time. TABLE IV. furnishes the numbers required in converting hours, minutes, and seconds into decimals of a day. Thus, to convert I3h 19m 43.5s into the decimal of a day, we find from the Table 13A =0.5416667 19m =0.0131944 43s = 0.0004977 0.5s = 0.0000058 Therefore ISA 19m 43.5s = 0.5553646 861 652 THEORETICAL ASTRONOMY. The decimal corresponding to 0.5s is found from that for 5s by changing the place of the decimal point. TABLE V. serves to find, for any instant, the number of days from the beginning of the year. Thus, for 1863 Sept. 14, 15h 53m 37.2s, we have Sept. 0.0 = 243.00000 days from the beginning of the year. Ud 15/i 53m 37.2s = 14.66224 Required number of days = 257.66224 TABLE VI. contains the values of M = 75 tan %v -f 25 tan 3 \v for values of v at intervals of one minute from to 180. For an ex- planation of its construction and use, see Articles 22, 27, 29, 41, and 72. In the case of parabolic motion the formulae are m== f M mO T\ wherein log C = 9.9601277. From these, by means of the Table, v may be found when t T is given, or t T when v is known. From v = 30 to v = !SO the Table contains the values of log M. TABLE VII., the construction of which is explained in Art. 23, serves to determine, in the case of parabolic motion, the true anomaly or the time from the perihelion when v approaches near to 180. The formulae are 8 /200 =Y- w being taken in the second quadrant. The Table gives the values of A O with w as the argument. As an example, let it be required to find the true anomaly corresponding to the values t T= 22.5 days and log g = 7.902720. From these we derive log M = 4.4582302. Table VI. gives for this value of log M, taking into account the second differences, i> = 16859'32".49; but, using Table VII., we have w = 168 59' 29".ll, A O = 3".37, EXPLANATION OF THE TABLES. 653 and hence v = w + A O = 168 59' 32".48, the two results agreeing completely. TABLE VIII. serves to find the time from the perihelion in the case of parabolic motion. For an explanation of its construction and use, see Articles 24, 69, and 72. TABLE IX. is used in the determination of the true anomaly 01 the time from the perihelion in the case of orbits of great eccen- tricity. Its construction is fully explained in Art. 28, and its use in Art. 41. TABLE X. serves to find the value of v or of t T in the case of elliptic or hyperbolic orbits. The construction of this Table is ex- plained in Art. 29. The first part gives the values of log B and log C 9 with A as the argument, for the ellipse and the hyperbola. In the case of log C there are given also log I. Diff. and log half II. Diif., expressed in units of the seventh decimal place, by means of which the interpolation is facilitated. Thus, if we denote by log (C) the value which the Table gives directly for the argument next less than the given value of A, and by &A the difference between this argument and the given value of A, expressed in units of the second decimal place, we have, for the required value, log 0= log (0) + Aj. X I. Diff. + AJ. 2 X half II. Diff. For example, let it be required to find the value of log C correspond- ing to A = 0.02497 944, and the process will be: (1) (2) Arg. 0.02, log (C) = 0.0034986 log I. Diff. =4.24585 log half II.Diff. = 1.778 ^ _ 8770.6 log* A =9.69718 21ogA4 9.394 A4 = 0.497944, (2)= 14.8 3.94303 1.172 log = 0.0043771 The second part of the Table gives the values of A corresponding to given values of r. TABLE XI. serves to determine the chord of the orbit when the extreme radii-vectores and the time of describing the parabolic arc are given. For an explanation of the construction and use of this Table, see Articles 68, 72, and 117. 654 THEORETICAL ASTRONOMY. TABLE XII. exhibits the limits of the real roots of the equation sin (/ C) = m sin* /. The construction and use of this table are fully explained in Articles 84 and 93. TABLES XIII. and XIV. are used in finding the ratio of the sector included by two radii-vectores to the triangle included by the same radii-vectores and the chord joining their extremities. For an explanation of the construction and use of these Tables, see Articles 88, 89, 93, and 101. TABLE XV. is used in the determination of the chord of the part of the orbit described in a given time in the case of very eccentric elliptic motion, and in the determination of the interval of time whenever the chord is known. For an explanation of its construc- tion and use, see Articles 116, 117, and 119. TABLE XVI. is used in finding the chord or the interval of time in the case of hyperbolic motion. See Articles 118 and 119 for an explanation of the use of the Table, and also the explanation of Table X. for an illustration of the use of the columns headed log I. Diff. and log half II. Diff. TABLE XVII. is used in the computation of special perturbations when the terms depending on the squares and higher powers of the masses are taken into account. For an explanation of its construc- tion and use, see Articles 157, 165, 166, 170, and 171. TABLE XVIII. contains the elements of the orbits of the comets which have been observed. These elements are: T } the time of peri- helion passage (mean time at Greenwich) ; TT, the longitude of the perihelion; &, the longitude of the ascending node; i, the inclina- tion of the orbit to the plane of the ecliptic; e, the eccentricity of the orbit; and q, the perihelion distance. The longitudes for Nos. 1, 2, 12, 16, 91, 92, 115, 127, 138, 155, 156, 159, 160, 162, 171, 173-175 180, 181, 185, 191, 192, 195-199, 201, 203, 204, 207, 208, 212-215, 217-219, 221-228, 230, 233, 234, 237-248, 251-258, 261-267, 269-275, 277-279, are in each case measured from the mean equinox of the beginning of the year. In the case of Nos. 134, 146, 172, 182, 189, 190, 205, 231, 232, 236, 259, and 268, the longitudes are EXPLANATION OF THE TABLES. 655 pleasured from the mean equinox of the beginning of the next year. The longitudes for Nos. 19 and 27 are measured from the me*\n equinox of 1850.0; for No. 186, from the mean equinox of July 3; for No. 187, from the mean equinox of Nov. 9; for No. 200, from the mean equinox of July 1 ; for No. 202, from the mean equinox of Oct. 1 ; for No. 206, from the mean equinox of Oct. 7 ; for No. 211, from the mean equinox of 1848.0; for No. 216, from the mean equi- nox of Feb. 20 ; for No. 229, from the mean equinox of April 1 ; foi No. 250, from the mean equinox of Oct. 1 ; and for No. 276, from the mean equinox of 1865 Oct. 4.0. Nos. 1, 2, 11, 12, 20, 23, 29, 41, 53, 80, and 177 give the elements for the successive appearances of Halley's comet; Nos. 104, 116, 126, 143, 149, 157, 167, 170, 176, 178, 183, 194, 210, 220, 235, 249, and 260, those for Encke's comet, the longitudes being measured from the mean equinox for the instant of the perihelion passage. Nos. 92, 127, 159, 172, 196, and 222 give the elements for the successive ap- pearances of Biela's comet; Nos. 187, 216, 250, and 276, those for Faye's comet; Nos. 197 and 238, those for Brorsen's comet; Nos. 217 and 243, those for D' Arrest's comet; and Nos. 145 and 245, those for Winnecke's comet. For epochs previous to 1583 the dates are given according to the old style. This Table is useful for identifying a comet which may appear with one previously observed, by means of a similarity of the ele- ments, its periodic character being otherwise unknown or at least un- certain. The elements given are those which appear to represent the observations most completely. For a collection of elements by vari- ous computers, and also for information in regard to the observations made and in regard to the place and manner of their publication, consult Carl's Repertorium der Cometen-Astronomie (Munich, 1864), or Galle's Cometen-Verzeichniss appended to the latest edition of Olbers's Meihode die Bahn eines Cometen zu berechnen. TABLE XIX. contains the elements of the orbits of the minor planets, derived chiefly from the Berliner Aslronomisches Jahrbiivh fur 1868. The epoch is given in Berlin mean time ; M denotes the mean anomaly,
, r , \