def test_hyperbolic_orbits_1(self):
o = Orbit(
i = 0.14170287439640022,
Ω = 3.978394514307273,
ω = 0.22281479333802098,
e = 1.077961355413604,
a = -8856935.204227254,
M0 = -6.003136275330101,
t0 = 34056451.522458464,
body = CelestialBody(
name = 'kerbin',
equatorial_radius = 600000.0,
gravitational_parameter = 3531600035840.0,
rotational_speed = 0.0002908894093707204,
),
)
t = o.epoch
tpe = o.time_to_periapsis_at_epoch
npoints = 10
tmin = t + tpe - 30
tmax = t + tpe + 30
tt = np.linspace(tmin,tmax,npoints)
lat = o.latitude_at_time(tt) / deg
lon = o.longitude_at_time(tt) / deg
r = o.radius_at_time(tt)
