The directed bond percolation process is an important model in statistical physics. By now its universal properties are known only up to the second-order of the perturbation theory. Here, our aim is to put forward a numerical technique with anomalous dimensions of directed percolation to higher orders of perturbation theory and is focused on the most complicated Feynman diagrams with problems in calculation. The anomalous dimensions are computed up to three-loop order in ϵ = 4-d.