A C++ Code to Solve the DGLAP Equations Applied to Ultra High Energy Cosmic Rays

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.