from multistar import multi as M
from multistar.direct import Direct as D
from physconst import DAY, YEAR
m = M('^arrokoth')
m.rund(1300*YR)
m.rund(1., cont=(m.t[-1] - 10 * DAY) * 1j)
d = D(m, orbit=0)
e = D(d)
e.set_flags(sd=False, td=False, tf=False, pqm=False, pqq=True)
e.rund()
e.plot_jctc()
e.flags
e.plot_et()