The simulation scenario allows for easily adding more objects with varied orbit patterns.
We can use a script like this:
import rebound # setup simulation sim = rebound.Simulation() sim.getWidget() # add star sim.add(m=1) # add ten 'planets' for i in range(10): sim.add(m=1e-3, a=0.4+0.1*i, inc=0.03*i, omega=5.*i) # center all particles sim.move_to_com() # advance the simulation sim.integrate(500) # get our final display on screen fig = rebound.OrbitPlot(sim) fig
This results in a display like: