Direct Solution of Renormalization Group Equations of QCD in $x-$space: NLO Implementations at Leading Twist

We illustrate the implementation of a method based on the use of recursion relations in (Bjorken) $x-$space for the solution of the evolution equations of QCD for all the leading twist distributions. The algorithm has the advantage of being very fast. The implementation that we release is written in C and is performed to next-to-leading order in $\alpha_s$.