In completely generic four-dimensional gauge-Yukawa theories, the renormalization group $$ \beta $$
β
-functions are known to the 3–2–2 loop order in gauge, Yukawa, and quartic couplings, respectively. It does, however, remain difficult to apply these results to realistic models without the use of dedicated computer tools. We describe a procedure for extracting $$ \beta $$
β
-functions using the general results and introduce , a dedicated Mathematica package for extracting the $$ \overline{\text {MS}} $$
MS
¯
$$ \beta $$
β
-functions in broad classes of models. The package and example notebooks are available from the GitHub repository.