eph = EphemerisProvider("de440.bsp")
eph = EphemerisProvider(["de440.bsp", "moon_pa_de440.bpc"])

# from/to are NAIF integer IDs, time is TDB seconds since J2000
r     = ephem_vector3(eph, from, to, time)    # SVector{3}  position
rv    = ephem_vector6(eph, from, to, time)    # SVector{6}  pos + vel
rva   = ephem_vector9(eph, from, to, time)    # SVector{9}  pos + vel + acc
rvaj  = ephem_vector12(eph, from, to, time)   # SVector{12} pos + vel + acc + jerk

ang   = ephem_rotation3(eph, from, to, time)   # SVector{3}  angles
angd  = ephem_rotation6(eph, from, to, time)   # SVector{6}  angles + rates
angdd = ephem_rotation9(eph, from, to, time)   # SVector{9}

ephem_get_points(eph)        # Available NAIF IDs
ephem_get_axes(eph)          # Available frame IDs
ephem_timescale_id(eph)      # 1=TDB, 2=TCB
ephem_spk_records(eph)       # Vector{EphemRecordSPK}
ephem_spk_timespan(eph)      # (first_time, last_time, continuity)