  We solve CCFM evolution equation numerically using the CohRad program based
on Monte Carlo methods. We discuss the effects of removing soft emissions and
non-Sudakov form factor by comparing the obtained distributions as functions of
accumulated transverse momenta or fractions of proton's longitudinal momenta.
We also compare the solution of the CCFM with the DGLAP equation in the gluonic
channel. Finally, we analyze the infra-red behaviour of solutions using the
so-called diffusion plots.
