added routine to compute sky positions and RVs for 2-body Keplerian system