# A C++ code to solve the DGLAP equations applied to Ultra High Energy Cosmic Rays

## Description

Abstract We solve numerically the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations for the evolution of fragmentation functions using the Laguerre method. We extend this method to include supersymmetric evolution. The solution to the DGLAP equations is particularly interesting to calculate the expected spectra of Ultra High Energy Cosmic Rays in models where they are produced by the decay of a massive particle X, M_X> 10^(12)GeV. Title of program: evolve Catalogue Id: ADPQ_v1_0 Nature of 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. Versions of this program held in the CPC repository in Mendeley Data ADPQ_v1_0; evolve; 10.1016/S0010-4655(01)00458-1 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)