def setup_orbit():
b = CelestialBody(
equatorial_radius = 6378*km,
gravitational_parameter = 398600*km**3)
o = Orbit(body=b)
o.epoch = 0
o.periapsis = 6600*km
o.speed_at_periapsis = 1.2 * b.escape_speed(6600*km)
return o
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