def METHOD_NAME(self):
    # discrete epochs - astropy.time.Time objects
    epochs = Time(['2018-01-02', '2018-01-05'], format='iso', scale='tdb')
    data = Orbit.from_horizons('Ceres', epochs=epochs)
    assert len(data.table) == 2