def setup_orbit(): b = Bunch( equatorial_radius = 6378*km, gravitational_parameter = 398600*km**3, ) o = Orbit(body=b) o.eccentricity = 1 o.periapsis = 7000*km return o