If we use this script: import reboundsim = rebound.Simulation()sim.add(m=1)sim.add(a=1)fig = rebound.OrbitPlot(sim)fig We have the resulting figure: