A C++ code to solve the DGLAP equations applied to ultra high energy cosmic rays. R. Toldra.

PROGRAM SUMMARY
Title of program: evolve
Catalogue identifier: ADPQ
Ref. in CPC: 143(2002)287
Distribution format: tar gzip file
Operating system: SunOS 5.7, OSF1 4.0, Linux Mandrake 7.2
Number of lines in distributed program, including test data, etc: 2995
Keywords: UHECR, Fragmentation functions, DGLAP equations, Laguerre method, Supersymmetry, Elementary particle physics, Phenomenological model.
Programming language used: C++
Computer: Sun , Alpha , Dell-PC .

Nature of physical problem:
In order to predict the spectra of UHECR produced by the decay of a dark matter superheavy particle with mass MX one needs to calculate fragmentation functions at the energy scale MX. These can be calculated from measured low energy fragmentation functions using the DGLAP equations.

Method of solution:
The DGLAP equations are solved by expanding them in Laguerre polynomials which reduces their integration to the computation of a set of coefficients. These coefficients are given by algebraic recursive relations.

Typical running time:
A few seconds.

Restrictions:
Gluon coherence at x < 10^-3 is not included.