NBodySimulation

NBodySimulation[law, {state1, …, staten}, t]

NBodySimulation[law, < | body1 -> state1, …, bodyn -> staten | >, t]