Feature request: Add mean motion (n) approximant

[duetosymmetry]

The Newtonian definition of mean motion is n = 2\pi/T where T is the radial period. I think this would be a nice smooth function, even to extrapolate past the ends of the earliest/latest detectable peri/apo-center passages.