Hello, I am not sure if this is the right place to ask this, but I don't readily see a "numerical methods" forum here so I assumed this would be the place to go. Sorry if I overlooked another place to post this! Anyway, I have some points interacting via potentials that are dependent only on distances (IE, not time or velocity), and I am wondering what would be the best ODE algorithm to use to simulate the motions of these points. The number of interacting points is fairly large, from 150 to 300. Let's say they're interacting via an additive pairwise potential similar to Newton's law of gravitation. What would be a decent ODE algorithm to use that is both fairly accurate (and not subject to huge error propagation) and not too complicated for solving the equations of motion? Any help would be greatly appreciated.