QCDNUM numerically solves the DGLAP evolution equations on a discrete grid in x and Q2. You can evolve unpolarised parton density functions in NNLO, and polarised pdfs or fragmentation functions in NLO. The jobs example.f and exampleCxx.cc show how to do this in a few lines of code.
QCDNUM supports evolution in the fixed and variable number schemes. To study the scale uncertainties, the renormalisation scale can be varied with respect to the factorisation scale.
QCDNUM also provides a large toolbox where you can enter your own splitting functions and then solve N-fold coupled DGLAP evolution equations (this is available in the pre-releases 17-01/xx). Also provided are routines to convolute pdfs with each other, or with user-defined kernels. With these tools you can write your own evolution routine (e.g. QCD-QED) or calculate parton luminosities and structure functions in both the massless and generalised mass schemes.
With the toolbox you can extend the functionality of QCDNUM as is done by the structure function add-on packages ZMSTF (unpolarised zero-mass structure functions) and HQSTF (heavy quark structure functions in the fixed-flavour number scheme). Both these packages are included in the distribution and described in the write-up.
Also included in the distribution is the MBUTIL collection of FORTRAN routines which should be linked to QCDNUM; apart from this the program is completely stand-alone.
If you use QCDNUM, please refer to
M. Botje, Comput. Phys. Commun. 182 (2011) 490, arXiv:1005.1481
The stable version is QCDNUM-17-00/08, latest release on May 1, 2017.
The latest QCDNUM-18 pre-release is QCDNUM-17-01/14 (21-Dec-2017).
For the evolution of fragmentation functions please do not use versions prior to 17-00/07 and 17-01/12 (see erratum).