The three-body recombination reaction, or ternary association, is a termolecular reaction leading to a molecule after a three-body encounter that plays a vital role in many relevant scenarios in chemical physics. Here, we introduce the Python 3-Body Recombination program, which is dedicated to the computation of atomic three-body recombination rate coefficients. The software is based on a classical trajectory approach in hyperspherical coordinates after mapping the three-body problem as a single particle in a higher-dimensional space. This theoretical approach is fully general and applicable to any ion-atom-atom or atom-atom-atom three-body process. The predictive power of the methodology has been tested in several different experimental scenarios, reaching a good description of every system. The code structure is presented alongside examples and tests to ensure the software's capacity. In addition, the performance of the software after parallelization is shown.