- Kepler's laws of planetary motion
In

astronomy ,**Kepler's Laws of Planetary Motion**are three mathematical laws that describe the motion ofplanet s in theSolar System . Germanmathematician andastronomer (Johannes Kepler 1571 –1630 ) discovered them.Kepler studied the

observation s (theRudolphine tables ) ofTycho Brahe . Around 1605, Kepler found that Brahe's observations of the planets' positions followed three relatively simple mathematical laws.Kepler's laws challenged Aristotelean and Ptolemaic astronomy and physics. His assertion that the Earth moved, his use of ellipses rather than

epicycle s, and his proof that the planets' speeds varied, changedastronomy andphysics . Almost a century laterIsaac Newton was able to deduce Kepler's laws from Newton's own laws of motion and his law of universal gravitation, using classical Euclidean geometry.In modern times, Kepler's laws are used to calculate approximate orbits for

artificial satellite s, and bodies orbiting theSun of which Kepler was unaware (such as the outer planets and smallerasteroid s). They apply where any relatively small body is orbiting a larger, relatively massive body, though the effects ofatmospheric drag (e.g. in a low orbit),relativity (e.g.Perihelion precession of Mercury ), and other nearby bodies can make the results insufficiently accurate for a specific purpose.Kepler's laws are formulated below using heliocentric

polar coordinate s. However, they can alternatively be formulated and derived usingCartesian coordinates . [*Hyman, Andrew. [*]*http://www.andrewhyman.com/articles/newton.pdf "A Simple Cartesian Treatment of Planetary Motion"*] , "European Journal of Physics", Vol. 14, pp. 145-147 (1993).**First law**The first law says:"The

orbit of everyplanet is anellipse with the sun at one of the foci."The mathematics of the ellipse is as follows. The equation is :$r=frac\{p\}\{1+epsiloncdotcos\; u\}$ where ("r","ν") are heliocentricical polar coordinates for the planet, "p" is the "

semi-latus rectum ", and "ε" is the "eccentricity", which is greater than or equal to zero, and less than one. For "ν"=0 the planet is at theat minimum distance: :$r\_mathrm\{min\}=frac\{p\}\{1+epsilon\}$perihelion for "ν"=90°: "r"="p", and for "ν"=180° the planet is at the

at maximum distance: :$r\_mathrm\{max\}=frac\{p\}\{1-epsilon\}$aphelion The

semi-major axis is thearithmetic mean between "r"_{min}and "r"_{max}::$a=frac\{p\}\{1-epsilon^2\}$ Thesemi-minor axis is thegeometric mean between "r"_{min}and "r"_{max}::$b=frac\; p\{sqrt\{1-epsilon^2$and it is also thegeometric mean between the semimajor axis and the semi latus rectum::$frac\; a\; b=frac\; b\; p$**econd law**The second law: "A line joining a planet and the sun sweeps out equal areas during equal intervals of time." [

*" [*]*http://demonstrations.wolfram.com/KeplersSecondLaw/ Kepler's Second Law*] " by Jeff Bryant with Oleksandr Pavlyk,The Wolfram Demonstrations Project .This is also known as the law of equal areas. It is a direct consequence of the law of conservation of angular momentum; see the derivation below.

Suppose a planet takes one day to travel from point "A" to "B". The lines from the Sun to "A" and "B", together with the planet orbit, will define a (roughly triangular) area. This same amount of area will be formed every day regardless of where in its orbit the planet is. This means that the planet moves faster when it is closer to the sun.

This is because the sun's gravity accelerates the planet as it falls toward the sun, and decelerates it on the way back out, but Kepler did not know that reason.

The two laws permitted Kepler to calculate the position, ("r","ν"), of the planet, based on the time since

perihelion , "t", and the orbital period, "P". The calculation is done in four steps.:1. Compute the**mean anomaly**"M" from the formula::$M=frac\{2pi\; t\}\{P\}$:2. Compute the"E" by numerically solvingeccentric anomaly **Kepler's equation**:::$M=E-epsiloncdotsin\; E$:3. Compute the"ν" by the equation:::$anfrac\; u\; 2\; =\; sqrt\{frac\{1+epsilon\}\{1-epsiloncdot\; anfrac\; E\; 2$:4. Compute thetrue anomaly **heliocentric distance**"r" from the first law:::$r=frac\; p\; \{1+epsiloncdotcos\; u\}$The proof of this procedure is shown below.**Third law**The third law : "The squares of the

orbital period s of planets are directly proportional to the cubes of the axes of the orbits." Thus, not only does the length of the orbit increase with distance, theorbital speed decreases, so that the increase of theorbital period is more than proportional.:$P^2\; propto\; a^3$:$P$ = orbital period of planet:$a$ = semimajor axis of orbit

So the expression "P"

^{2}·"a"^{–3}has the same value for all planets in theSolar System as it has forEarth . When certain units are chosen, namely "P" is measured insidereal year s and "a" inastronomical unit s, "P"^{2}·"a"^{–3}has the value 1 for all planets in the Solar System.In

SI units : $frac\{P^\{2\{a^\{3\; =\; 3.00\; imes\; 10^\{-19\}\; frac\{s^\{2\{m^\{3\; pm\; 0.7\%,$.The law, when applied to

circular orbit s where theacceleration is proportional to "a"·"P"^{−2}, shows that the acceleration is proportional to "a"·"a"^{−3}= "a"^{−2}, in accordance withNewton's law of gravitation .The general equation, which was derived from Newton's law of gravity, is

:$left(\{frac\{P\}\{2pi\; ight)^2\; =\; \{a^3\; over\; G\; (M+m)\},$where $G$ is the

gravitational constant , $M$ is the mass of the sun, and $m$ is the mass of planet. The latter appears in the equation since the equation of motion involves thereduced mass .Note that "P" is time perorbit and "P"/2π is time perradian .See the actual figures: attributes of major planets.

This law is also known as the

**harmonic law**.**Position as a function of time**The Keplerian problem assumes an elliptical orbit and the four points:

*"s" the sun (at one focus of ellipse);

*"z" the perihelion

*"c" the center of the ellipse

*"p" the planetand :$a=|cz|,$ distance from center to perihelion, the**semimajor axis**,:$varepsilon=\{|cs|over\; a\},$ the**eccentricity**,:$b=asqrt\{1-varepsilon^2\},$ the**semiminor axis**,:$r=|sp|\; ,$ the distance from sun to planet.and the angle:$u=angle\; zsp,$ the planet as seen from the sun, the.true anomaly The problem is to compute the

polar coordinates ("r","ν") of the planet from the**time since perihelion**, "t".It is solved in steps. Kepler began by adding the orbit's auxiliary circle (that with the major axis as a diameter) and defined these points:

*"x" is the projection of the planet to the auxiliary circle; then the area $|zsx|=frac\; a\; b\; cdot|zsp|$

*"y" is a point on the auxiliary circle such that the area $|zcy|=|zsx|$and :$M=angle\; zcy$, "y" as seen from the centre, the.mean anomaly The area of the

circular sector $|zcy|\; =\; frac\{a^2\; M\}2$, and the area swept since perihelion, :$|zsp|=frac\; b\; a\; cdot|zsx|=frac\; b\; a\; cdot|zcy|=frac\; b\; acdotfrac\{a^2\; M\}2\; =\; frac\; \{a\; b\; M\}\{2\}$ , is by Kepler's second law proportional to time since perihelion. So the mean anomaly, "M", is proportional to time since perihelion, "t".:$M=\{2\; pi\; t\; over\; P\},$ where "P" is theorbital period .The mean anomaly "M" is first computed. The goal is to compute the true anomaly "ν". The function "ν"="f"("M") is, however, not elementary. Kepler's solution is to use :$E=angle\; zcx$, "x" as seen from the centre, the

as an intermediate variable, and first compute "E" as a function of "M" by solving Kepler's equation below, and then compute the true anomaly "ν" from the eccentric anomaly "E". Here are the details.:$|zcy|=|zsx|=|zcx|-|scx|$eccentric anomaly :$frac\{a^2\; M\}2=frac\{a^2\; E\}2-frac\; \{avarepsiloncdot\; asin\; E\}2$Division by "a"²/2 gives

**Kepler's equation**:$M=E-varepsiloncdotsin\; E$.The catch is that Kepler's equation cannot be rearranged to isolate "E". The function "E"="f"("M") is not an elementary formula, but Kepler's equation is solved either iteratively by aroot-finding algorithm or, as derived in the article oneccentric anomaly , by aninfinite series .Having computed the eccentric anomaly "E" from Kepler's equation, the next step is to calculate the true anomaly "ν" from the eccentric anomaly "E".

Note from the geometry of the problem that:$acdotcos\; E=acdotvarepsilon+rcdotcos\; u.$Dividing by "a" and inserting from Kepler's first law:$frac\; r\; a\; =frac\{1-varepsilon^2\}\{1+varepsiloncdotcos\; u\}$ to get:$cos\; E=varepsilon+frac\{1-varepsilon^2\}\{1+varepsiloncdotcos\; u\}cdotcos\; u=frac\{varepsiloncdot(1+varepsiloncdotcos\; u)+(1-varepsilon^2)cdotcos\; u\}\{1+varepsiloncdotcos\; u\}=frac\{varepsilon\; +cos\; u\}\{1+varepsiloncdotcos\; u\}.$The result is a usable relationship between the eccentric anomaly "E" and the true anomaly "ν".

A computationally more convenient form follows by substituting into the

trigonometric identity ::$an^2frac\{x\}\{2\}=frac\{1-cos\; x\}\{1+cos\; x\}.$Get:$an^2frac\{E\}\{2\}=frac\{1-cos\; E\}\{1+cos\; E\}=frac\{1-frac\{varepsilon+cos\; u\}\{1+varepsiloncdotcos\; u\{1+frac\{varepsilon+cos\; u\}\{1+varepsiloncdotcos\; u=frac\{(1+varepsiloncdotcos\; u)-(varepsilon+cos\; u)\}\{(1+varepsiloncdotcos\; u)+(varepsilon+cos\; u)\}=frac\{1-varepsilon\}\{1+varepsilon\}cdotfrac\{1-cos\; u\}\{1+cos\; u\}=frac\{1-varepsilon\}\{1+varepsilon\}cdot\; an^2frac\{\; u\}\{2\}.$Multiplying by (1+ε)/(1−ε) and taking the square root gives the result:$anfrac\; u2=sqrtfrac\{1+varepsilon\}\{1-varepsilon\}cdot\; anfrac\; E2.$We have now completed the third step in the connection between time and position in the orbit.

One could even develop a series computing "ν" directly from "M". [

*http://info.ifpan.edu.pl/firststep/aw-works/fsII/mul/mueller.html*]The fourth step is to compute the heliocentric distance "r" from the true anomaly "ν" by Kepler's first law::$r=acdotfrac\{1-varepsilon^2\}\{1+varepsiloncdotcos\; u\}.$

**Derivation from Newton's laws**[

thumb|right|300px|Two_bodies_with_similar_mass_orbiting_around_a_common_barycenter with elliptic orbits.] Kepler's laws are about the motion of the planets around the sun, whileNewton's laws more generally are about the motion of point particles attracting each other by the force ofgravitation .In the special case where there are only two particles, and one of them is much lighter than the other, and the distance between the particles remains limited, then the lighter particle moves around the heavy particle as a planet around the Sun according to Kepler's laws, as shown below. Newton's laws however also admit other solutions, where the trajectory of the lighter particle is a

parabola or ahyperbola . These solutions show that there is a limitation to the applicability of Kepler's first law, which states that the trajectory will always be an ellipse.In any two-body system (the

two body problem ), the objects always orbit around their common center of mass (barycenter ). When one object is very heavy compared to the other, the common center of mass is simply assumed to be the same as the center of mass of the heavy object (which is nearly correct). However, Kepler's laws can be used even in the case of two objects of similar mass, by decomposing the motion of each body as an ellipse around the barycenter.While Kepler's laws are expressed either in geometrical language or as equations connecting the coordinates of the planet and the time variable with the

orbital elements , Newton's second law is adifferential equation . So the derivations below involve the art of solving differential equations. Kepler's second law is derived first, as the derivation of the first law depends on the derivation of the second law.**Deriving Kepler's second law**Newton's law of gravitation says that "every object in the universe attracts every other object along a line of the centers of the objects, proportional to each object's mass, and inversely proportional to the square of the distance between the objects," and his second law of motion says that "the mass times the acceleration is equal to the force." So the mass of the planet times the acceleration vector of the planet equals the mass of the sun times the mass of the planet, divided by the square of the distance, times minus the radial

unit vector , times a constant of proportionality. This is written::$mcdotddotmathbf\{r\}\; =\; frac\{Mcdot\; m\}\{r^2\}cdot(-hat\{mathbf\{r)cdot\; G$

where a dot on top of the variable signifies differentiation with respect to time, and the second dot indicates the second derivative.

Assume that the planet is so much lighter than the sun that the acceleration of the sun can be neglected.

:$dothat\{mathbf\{r\; =\; dot\; u\; hat\{oldsymbol\; u\}$where $hat\{oldsymbol\; u\}$ is the tangential unit vector, and:$dothat\{oldsymbol\; u\}\; =\; -dot\; u\; hat\{mathbf\{r.$

So the position vector:$mathbf\{r\}\; =\; r\; hat\{mathbf\{r$is differentiated twice to give the velocity vector and the acceleration vector:$dotmathbf\{r\}\; =dot\; r\; hatmathbf\{r\}\; +\; r\; dothatmathbf\{r\}=dot\; r\; hat\{mathbf\{r\; +\; r\; dot\; u\; hat\{oldsymbol\; u\},$:$ddotmathbf\{r\}\; =\; (ddot\; r\; hat\{mathbf\{r\; +dot\; r\; dothat\{mathbf\{r\; )+\; (dot\; rdot\; u\; hat\{oldsymbol\; u\}\; +\; rddot\; u\; hat\{oldsymbol\; u\}+\; rdot\; u\; dothat\{oldsymbol\; u\})=\; (ddot\; r\; -\; rdot\; u^2)\; hat\{mathbf\{r\; +\; (rddot\; u\; +\; 2dot\; r\; dot\; u)\; hat\{oldsymbol\; u\}.$

Note that for constant distance, $r$, the planet is subject to the

centripetal acceleration , $rdot\; u^2$, and for constant angular speed, $dot\; u$, the planet is subject to thecoriolis acceleration , $2dot\; r\; dot\; u$.Inserting the acceleration vector into Newton's laws, and dividing by "m", gives the vector

equation of motion :$(ddot\; r\; -\; rdot\; u^2)\; hat\{mathbf\{r\; +\; (rddot\; u\; +\; 2dot\; r\; dot\; u)\; hat\{oldsymbol\; u\}=\; -GMr^\{-2\}hat\{mathbf\{r$Equating component, we get the two

ordinary differential equation s of motion, one for the radial acceleration and one for the tangential acceleration::$ddot\; r\; -\; rdot\; u^2\; =\; -GMr^\{-2\},$:$rddot\; u\; +\; 2dot\; rdot\; u\; =\; 0.$In order to derive Kepler's second law only the tangential acceleration equation is needed. Divide it by $r\; dot\; u:$:$frac\{ddot\; u\}\{dot\; u\}\; +2frac\{dot\; r\}\{r\}=0$and integrate::$logdot\; u\; +2log\; r\; =\; logell,$where $logell$ is a constant of integration, and exponentiate::$r^2dot\; u\; =ell\; .$This says that the

specific angular momentum $r^2\; dot\; u$ is aconstant of motion , even if both the distance $r$ and theangular speed $dot\; u$ vary.The area swept out from time "t"

_{1}to time "t"_{2},:$int\_\{t\_1\}^\{t\_2\}frac\; 1\; 2\; cdot\; basecdot\; d(height)\; =\; int\_\{t\_1\}^\{t\_2\}frac\; 1\; 2\; cdot\; rcdot\; rdot\; u\; dt=frac\; 1\; 2\; cdotell\; cdot(t\_2-t\_1)$depends only on the duration "t"_{2}−"t"_{1}. This is Kepler's second law.**Deriving Kepler's first law**The expression :$p=ell\; ^2\; G^\{-1\}M^\{-1\}$has the dimension of length and is used to make the equations of motion dimensionless. We define:$u\; =pr^\{-1\}$and get:$-GMr^\{-2\}=-ell^2\; p^\{-3\}u^\{2\}$and:$dot\; u\; =ell\; r^\{-2\}=ell\; p^\{-2\}u^2.$

Differentiation with respect to time is transformed into differentiation with respect to angle::$dot\; X=frac\; \{dX\}\{dt\}=frac\; \{dX\}\{d\; u\}cdotfrac\; \{d\; u\}\{dt\}=frac\; \{dX\}\{d\; u\}cdot\; dot\; u=frac\; \{dX\}\{d\; u\}cdotell\; p^\{-2\}u^2.$

Differentiate:$r\; =pu^\{-1\}$twice::$dot\; r\; =\; frac\{d(pu^\{-1\})\}\{d\; u\}cdotell\; p^\{-2\}u^\{2\}\; =\; -pu^\{-2\}frac\{du\}\{d\; u\}cdotell\; p^\{-2\}u^\{2\}=\; -ell\; p^\{-1\}frac\{du\}\{d\; u\}$:$ddot\; r\; =\; frac\{ddot\; r\}\{d\; u\}cdotell\; p^\{-2\}u^\{2\}=\; frac\{d\}\{d\; u\}(-ell\; p^\{-1\}\; frac\{du\}\{d\; u\})cdotell\; p^\{-2\}u^\{2\}=\; -ell^2\; p^\{-3\}u^\{2\}frac\{d^2\; u\}\{d\; u^2\}$Substitute into the radial equation of motion:$ddot\; r\; -\; rdot\; u^2\; =\; -GMr^\{-2\}$and get:$(-ell^2\; p^\{-3\}u^2frac\{d^2u\}\{d\; u^2\})\; -\; (pu^\{-1\})(ell\; p^\{-2\}u^2)^2\; =\; -ell\; ^2\; p^\{-3\}\; u^2$Divide by $-ell^2\; p^\{-3\}u^2$ to get a simple non-homogeneous linear differential equation for the orbit of the planet:

:$frac\{d^2u\}\{d\; u^2\}\; +\; u\; =\; 1\; .$

An obvious solution to this equation is the circular orbit:$u\; =\; 1.$Other solutions are obtained by adding solutions to the homogeneous linear differential equation with constant coefficients:$frac\{d^2u\}\{d\; u^2\}\; +\; u\; =\; 0$These solutions are:$u\; =\; epsiloncdotcos(\; u-\; u\_0)$where $epsilon$ and $u\_0,$ are arbitrary constants of integration. So the result is:$u\; =\; 1+\; epsiloncdotcos(\; u-\; u\_0)$

Choosing the axis of the

coordinate system such that $u\_0=0$, and inserting $u=pr^\{-1\}$, gives::$pr^\{-1\; \}\; =\; 1+\; epsiloncdotcos\; u\; .$

If $epsilon<1\; ,$ this is Kepler's first law.

**Deriving Kepler's third law**The area of the planetary orbit ellipse is:$A=pi\; a\; b=pi\; a(asqrt\{1-epsilon^2\})=pi\; a^2sqrt\{1-epsilon^2\},$The area speed of the radius vector sweeping the orbit area is:$dot\; A=ell/\; 2\; ,$where:$ell^2=pGM=a(1-epsilon^2)GM\; ,.$The period of the orbit is:$P=frac\; A\; dot\; A\; =frac\; \{pi\; a^2sqrt\{1-epsilon^2\{ell/2\}=2pifrac\; \{a^2sqrt\{1-epsilon^2\{ell\},$satisfying:$left(frac\; P\{2\; pi\}\; ight)^2=left(\; frac\; \{a^2sqrt\{1-epsilon^2\{ell\}\; ight)^2=frac\; \{a^4(1-epsilon^2)\}\{ell^2\}=frac\; \{a^4(1-epsilon^2)\}\{a(1-epsilon^2)GM\; \}=frac\{a^3\}\{GM\},$which is Kepler's third law.

**Notes****References***Kepler's life is summarized on pages 627-623 and Book Five of his "magnum opus", "

Harmonice Mundi " ("harmonies of the world"), is reprinted on pages 635-732 of "On the Shoulders of Giants": The Great Works of Physics and Astronomy (works byCopernicus , Kepler,Galileo , Newton, and Einstein).Stephen Hawking , ed. 2002 ISBN 0-7624-1348-4*A derivation of Kepler's third law of planetary motion is a standard topic in engineering mechanics classes. See, for example, pages 161-164 of Citation|first=J. L. |last=Meriam|year=1971|date=1966, 1971|title=Dynamics, 2nd ed. |location=New York|publisher=John Wiley|ISBN=0-471-59601-9.

*Murray and Dermott, Solar System Dynamics, Cambridge University Press 1999, ISBN-10 0-521-57597-4

**ee also***

Kepler orbit

*Kepler problem

*Circular motion

*Gravity

*Two-body problem

*Free-fall time

*Laplace-Runge-Lenz vector **External links*** B.Surendranath Reddy; animation of Kepler's laws: [

*http://www.surendranath.org/Applets/Dynamics/Kepler/Kepler1Applet.html applet*]

* Crowell, Benjamin, "Conservation Laws", [*http://www.lightandmatter.com/area1book2.html http://www.lightandmatter.com/area1book2.html*] , an online book that gives a proof of the first law without the use of calculus. (see section 5.2, p.112)

* David McNamara and Gianfranco Vidali, "Kepler's Second Law -JAVA Interactive Tutorial", [*http://www.phy.syr.edu/courses/java/mc_html/kepler.html http://www.phy.syr.edu/courses/java/mc_html/kepler.html*] , an interactive JAVA applet that aids in the understanding of Kepler's Second Law.

* University of Tennessee's Dept. Physics & Astronomy: Astronomy 161 page on Johannes Kepler: The Laws of Planetary Motion [*http://csep10.phys.utk.edu/astr161/lect/history/kepler.html*]

* Equant compared to Kepler: interactive model [*http://people.scs.fsu.edu/~dduke/kepler.html*]

* Kepler's Third Law:interactive model [*http://people.scs.fsu.edu/~dduke/kepler3.html*]

*Wikimedia Foundation.
2010.*

### Look at other dictionaries:

**Kepler's laws of planetary motion**— ▪ astronomy in astronomy and classical physics, laws describing the motions of the planets in the solar system. They were derived by the German astronomer Johannes Kepler, whose analysis of the observations of the 16th century Danish… … Universalium**Kepler's law of planetary motion**— noun (astronomy) one of three empirical laws of planetary motion stated by Johannes Kepler • Syn: ↑Kepler s law • Topics: ↑astronomy, ↑uranology • Hypernyms: ↑law, ↑law of nature … Useful english dictionary**laws of planetary motion**— kepler s laws … Useful english dictionary**Kepler’s laws**— of planetary motion, named after their discoverer Johannes Kepler (1571–1630), state: Time to travel from A to X = C to Y Time to travel (E Q E)2 = (EO)3 i. The path of a planet around the sun is an ellipse, and the sun is at one of the foci of… … Aviation dictionary**Kepler's laws**— Astron. three laws governing planetary motion: each planet revolves in an ellipse, with the sun at one focus; the line connecting a planet to the sun sweeps out equal areas in equal periods of time (law of areas); the square of the period of… … Universalium**Kepler's laws**— noun the three laws of planetary motion discovered by Kepler in the early 17th century See Also: Keplers equation … Wiktionary**Kepler's laws**— /ˈkɛpləz lɔz/ (say kepluhz lawz) plural noun three laws of planetary motion stating: 1. that the planets move in elliptical orbits about the sun, which is situated at one focus of the ellipses. 2. that the radius vectors joining each planet to… … Australian English dictionary**Kepler's laws**— Law Law (l[add]), n. [OE. lawe, laghe, AS. lagu, from the root of E. lie: akin to OS. lag, Icel. l[ o]g, Sw. lag, Dan. lov; cf. L. lex, E. legal. A law is that which is laid, set, or fixed; like statute, fr. L. statuere to make to stand. See… … The Collaborative International Dictionary of English**Motion (physics)**— Motion involves change in position, such as in this perspective of rapidly leaving Yongsan Station In physics, motion is a change in position of an object with respect to time. Change in action is the result of an unbalanced force. Motion is… … Wikipedia**Kepler (disambiguation)**— Kepler may refer to: *Johannes Kepler, a key figure in the scientific revolution. Named after him are:Ideas: *Kepler s laws of planetary motion, a set of calculations for the orbits of planets *The Kepler Poinsot polyhedra, a set of geometrical… … Wikipedia