MathGroup Archive 2001

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: kepler equation

  • To: mathgroup at
  • Subject: [mg31241] Re: kepler equation
  • From: "Hermann Meier" <hermann.ernst.meier at>
  • Date: Tue, 23 Oct 2001 04:53:19 -0400 (EDT)
  • References: <9qrd3f$c17$>
  • Sender: owner-wri-mathgroup at

Kepler's equation: mean anomaly (input) and eccentric anomaly (output) are
both in radians:

keplerEquation[meanAnomaly_,eccentricity_] :=

The following solution may also be used (for moderate values of eccentricity
only; quite slow) (Taff, Celestial Mechanics, p. 34)

kep[M_,e_] :=
   1/m*(((e*m)^m*BesselJ[m, Sqrt[e^2*m^2]])/(e^2*m^2)^(-(-m/2)))*
 {m,1,Infinity}, Method -> Fit])

M = 199.7855*Degree; e = 0.047837;

keplerEquation[M,e] / Degree // InputForm or kep[M,e]/Degree // InputForm
both give 198.89778838060755

Hermann Meier
hermann.ernst.meier at

"ammar" <ammarsak at> wrote in message
news:9qrd3f$c17$1 at
> Dear sir
> would you send to me some source files of Kepler equation based on
> or any related  topics
> thanks
> Ammar Jamil Sakaji
> P.O.Box 45247
> Abu Dhabi - UAE

  • Prev by Date: Re: font problem
  • Next by Date: Re: Limit and Abs
  • Previous by thread: kepler equation
  • Next by thread: Mathematica question