  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).
