It is widely known that in QCD the infrared region is not available to perturbative calculations due to the presence of a Landau pole in the running of the coupling constant. We present non perturbative solutions for the propagators in Yang-Mills theory by numerically solving a truncated system of Dyson-Schwinger equations and compare them with analytical results found by applying Callan-Symanzik renormalization group equations in an epsilon expansion near the infrared fixed point, when a mass term for the gluon is added to the action. Copyright © 2012 by INFN.