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)