It is possible to obtain the true anomaly of a planet this way:

In[1]:= PlanetData["Earth", "TrueAnomaly"]
Out[1]= Quantity[MixedMagnitude[{7, 18.8561}], MixedUnit[{"AngularDegrees", "ArcMinutes"}]]

However, the true anomaly is actually a function of time. So I suppose it gives the anomaly for the current moment. Is it possible to obtain this with a DateObject (or similar) parameter? I already have an implementation of the Kepler equations, but I am curious if there is a simpler (and perhaps more precise) way to obtain the anomalies.


1 Answer 1


So I found this solution:

                          {"Date" -> DateObject[{2020, 1, 1}]}]]

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.