A program package for calculating regularized classical trajectories for Coulomb n−body problem is developed. The Coulomb singularities from the equations of motion are removed by transformations of variables including the time. This effectively conserves the energy of the time-independent systems to a high accuracy for long time propagation. Sample calculations are shown for the cases of 2,3,4, and 5 particle systems giving comparisons with the un-regularized trajectories. The program can be used for general purposes including the classical-trajectory Monte-Carlo simulations for charged-particle collisions in free or laser environments.