We discuss precision Monte Carlo (MC) calculations for solving the QCD evolution equations up to the next-to-leading-order (NLO) level. They employ forward Markovian Monte Carlo algorithms, which provide rigorous solutions of the above equations. These algorithms are implemented in the form of the Monte Carlo program  EvolFMC. This program has been cross-checked with independent, non-MC, programs ( QCDNum16 and  APCheb33) and the numerical agreement at the level of 0.1\% has been found.

PACS numbers: 12.38.--t, 12.38.Bx, 12.38.Cy