Incorporate some basic oblateness

The oblateness force is for now radial.  Only Jupiter has oblateness.

Also restructure the test to be table-driven and adjust the numbers for the Galilean satellites (hint: they look better).
let's rename this to `expected_offset_from_parent`, similarly for `vector_error`.