Love seeing stuff like this. The corrections for an oblate spheroid threw me for a loop at first, until I realized "yeah, of course". I've only ever played around with ideal bodies when simulating the n-body problem (sounds a bit raunchy...) so never even considered the fact that a rotating planet isn't perfectly spherical.