The SU(3) beta function from numerical stochastic perturbation theory

The SU(3) beta function is computed from Wilson loops to 20th order numerical stochastic perturbation theory. An attempt is made to include massless fermions, whose contribution is known analytically to 4th order. The question whether the theory admits an infrared stable fixed point is addressed.