[ ] Add an initial condition with the Cartesian coordinates and velocities of N-1 bodies with respect to the first body and/or with respect to the center of mass of the system or the Cartesian coordinates of each of the nested Keplerians (e.g. planets surrounding a star).
[ ] Add an initial condition which is two-dimensional (plane-parallel orbits).