April 2018
Beginner
238 pages
7h 13m
English
We can use a script like this:
import rebound# setup simulationsim = rebound.Simulation()sim.getWidget()# add starsim.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 particlessim.move_to_com()# advance the simulationsim.integrate(500)# get our final display on screenfig = rebound.OrbitPlot(sim)fig
This results in a display like:
