The periodic time varying inversely as n, this equation expresses Kepler’s third law. This reasoning tacitly supposes the orbit to be a circle of radius a, and the mass of the planet to be negligible. The rigorous relation is expressed by a slight modification of the law. Putting M and m for the respective masses of the sun and planet, a for the semi-major axis of the orbit, and n for the mean angular motion in unit of time, the relation then is

a³n² = M + m.

What is noteworthy in this theorem is that this relation depends only on the sum of the masses. It follows, therefore, that were any portion of the mass of the sun taken from it, and added to the planet, the relation would be unchanged. Kepler’s third law therefore expresses the fact that the mass of the sun is the same for all the planets, and deviates from the truth only to the extent that the masses of the latter differ from each other by quantities which are only a small fraction of the mass of the sun.

Problem of Three Bodies.—As soon as the general law of gravitation was fully apprehended, it became evident that, owing to the attraction of each planet upon all the others, the actual motion of the planets must deviate from their motion in an ellipse according to Kepler’s laws. In the Principia Newton made several investigations to determine the effects of these actions; but the geometrical method which he employed could lead only to rude approximations. When the subject was taken up by the continental mathematicians, using the analytical method, the question naturally arose whether the motions of three bodies under their mutual attraction could not be determined with a degree of rigour approximating to that with which Newton had solved the problem of two bodies. Thus arose the celebrated “problem of three bodies.” Investigation soon showed that certain integrals expressing relations between the motions not only of three but of any number of bodies could be found. These were:—

First, the law of the conservation of the centre of gravity. This expresses the general fact that whatever be the number of the bodies which act upon each other, their motions are so related that the centre of gravity of the entire system moves in a straight line with a constant velocity. This is expressed in three equations, one for each of the three rectangular co-ordinates.

Secondly, the law of conservation of areas. This is an extension of Kepler’s second law. Taking as the radius vector of each body the line from the body to the common centre of gravity of all, the sum of the products formed by multiplying each area described, by the mass of the body, remains a constant. In the language of theoretical mechanics, the moment of momentum of the entire system is a constant quantity. This law is also expressed in three equations, one for each of the three planes on which the areas are projected.

Thirdly, the entire vis viva of the system or, as it is now called, the energy, which is obtained by multiplying the mass of each body into half the square of its velocity, is equal to the sum of the quotients formed by dividing the product of every pair of the masses, taken two and two, by their distance apart, with the addition of a constant depending on the original conditions of the system. In the language of algebra putting m1, m2, m3, &c. for the masses of the bodies, r1.2, r1.3, r2.3, &c. for their mutual distances apart; v1, v2, v3, &c., for the velocities with which they are moving at any moment; these quantities will continually satisfy the equation

½(m1v1² +m2v2² + ...) =m1m2 +m1m3 +m2m3 + ... + a constant.
r1.2r1.3r2.3

The theorems of motion just cited are expressed by seven integrals, or equations expressing a law that certain functions of the variables and of the time remain constant. It is remarkable that although the seven integrals were found almost from the beginning of the investigation, no others have since been added; and indeed it has recently been shown that no others exist that can be expressed in an algebraic form. In the case of three bodies these do not suffice completely to define the motion. In this case, the problem can be attacked only by methods of approximation, devised so as to meet the special conditions of each case. The special conditions which obtain in the solar system are such as to make the necessary approximation theoretically possible however complex the process may be. These conditions are:—(1) The smallness of the masses of the planets in comparison with that of the sun, in consequence of which the orbit of each planet deviates but slightly from an ellipse during any one revolution; (2) the fact that the orbits of the planets are nearly circular, and the planes of their orbits but slightly inclined to each other. The result of these conditions is that all the quantities required admit of development in series proceeding according to the powers of the eccentricities and inclinations of the orbits, and the ratio of the masses of the several planets to the mass of the sun.

Perturbations of the Planets.—Kepler’s laws do not completely express the motion of a planet around a central body, except when no force but the mutual attraction of the two bodies comes into play. When one or more other bodies form a part of the system, their action produces deviations from the elliptic motion, which are called perturbations. The problem of determining the perturbations of the heavenly bodies is perhaps the most complicated with which the mathematical astronomer has to grapple; and the forms under which it has to be studied are so numerous that they cannot be easily arranged under any one head. But there is one conception of perturbations of such generality and elegance that it forms the common base of all those methods of determining these deviations which have high scientific interest. This conception is embodied in the method of “variation of elements,” originally due to J.L. Lagrange. The simplest method of presenting it starts with the second view of the elliptic motion already set forth.