In this presentation we are going to discuss in detail the relevant aspects of the computation of the corrections to the lightest CP-even Higgs boson mass at three-loop accuracy in the SUSYQCD sector. We have adopted two different renormalization schemes: the first one is based on a variant of dimensional regularization where the so called dimensional reduction is performed in order to preserve supersymmetry to all perturbative orders, while the second scheme includes an on-shell renormalization of the squark masses. This calculation extends the region of validity of previous studies to the whole parameters space. The emerging loop integrals were computed by exploiting new proposed approaches for the numerical calculation of multi-scale multi-loop Feynman integrals. With the help of dispersion relations, three-loop Feynman integrals with arbitrary mass scales can be written in terms of one-dimensional and two-dimensional integrals of functions that are suitable for an efficient numerical evaluation.