Optimal smoothing for N-body codes

In any collisionless N-body code, there is an optimal choice for the smoothing parameter that minimizes the average error in the force evaluations. We show how to compute the optimal softening length in a direct-summation code and demonstrate that it varies roughly as 1/N^(1/3).