Monte Carlo solutions of the QCD evolution equations

We present various methods of solving the QCD evolution equations by means of Monte Carlo techniques developed in Cracow. In particular we present two algorithms that simulate the Constrained Markovian evolution. The second algorithm works for both GLAP- and HERWIG-type evolution.