On The Integral-Conservative Numerical Solution Of Few-Body Gravitational Problems