A numerical solution of the coupled Dyson-Schwinger equations for the ghost and gluon propagators in Yang-Mills theory is presented in Landau gauge. Aimed at investigating the infrared behavior of the propagators, the equations are simplified by neglecting the gluon loops, according to the ghost dominance hypothesis motivated by the Gribov-Zwanziger scenario. The equations are solved with an iterative method, eliminating the ultraviolet divergence through a continuous regulator function depending on the cut off scale. The solutions, derived for different values of the Euclidean space-time dimension, present scaling (the infrared exponents are obtained) or decoupling behavior, depending on whether the horizon condition is or not implemented. Moreover, it is shown that the running coupling constant approaches a constant value in the IR, corresponding to an attractive fixed point. © Published under licence by IOP Publishing Ltd.