- #1
DivergentSpectrum
- 149
- 15
is there a numerical method to determine whether two bodies will stay bounded forever in an n body simulation? i know if the energy of a particle orbitting the origin is negative, then it is bounded, where -∫(force)dr+(dr/dt)2/2=energy.
but I am curious about a genereral case, where there are more than 2 bodys, and even forces that are not inverse square proportional where i don't necessarily have a potential to subtract kinetic energy from, just the ability to numerically do the line integral of force
∫ [Fx,fy,fz]*[dx/dt,dy/dt,dz/dt] dt
thanks!
but I am curious about a genereral case, where there are more than 2 bodys, and even forces that are not inverse square proportional where i don't necessarily have a potential to subtract kinetic energy from, just the ability to numerically do the line integral of force
∫ [Fx,fy,fz]*[dx/dt,dy/dt,dz/dt] dt
thanks!